<!doctype html>

<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">

</head>

<body>

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">

<p>Hi,</p>

<p><br>

</p>

<p>Attached there are two scripts that I used to compute Geopotential for CMIP6 model MPI-ESM1.2-LR and to interpolate after to pressure levels.</p>

<p>Also attached there are 3 figures that compare the geopotential computed and interpolated to pressure levels with NCL with that provided by MPI. You may see that they compare fairly well.</p>

<p><br>

</p>

<p>You should check that every variable is what NCL's expects (units, coordinates ordering, etc. ). I think that for example you P0 is in hPa, but NCL expects that to be in Pa.</p>

<p><br>

</p>

<p>I don't think that you can compute and interpolate to pressure levels without <span>

the surface geopotential</span>.</p>

<p><br>

</p>

<p>Hope that helps.<br>

</p>

<p><br>

</p>

<div id="Signature">

<div style="font-family:Tahoma; font-size:13px">

<div style="font-family:Tahoma; font-size:13px"><br>

C. A. F. Marques<br>

CESAM & Physics Department<br>

<span lang="pt">

<div style="margin:0pt"><font size="2" face="Calibri,sans-serif"><span style="font-size:11pt">University of Aveiro</span></font></div>

<div style="margin:0pt"><font size="2" face="Calibri,sans-serif"><span style="font-size:11pt">3810-193 Aveiro. Portugal<br>

</span></font></div>

</span><a href="http://www.cesam.ua.pt/cafmarques" target="_blank" id="LPNoLP">http://www.cesam.ua.pt/cafmarques</a>

</div>

</div>

</div>

</div>

<hr style="display:inline-block;width:98%" tabindex="-1">

<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> ncl-talk <ncl-talk-bounces@mailman.ucar.edu> on behalf of markus mingel via ncl-talk <ncl-talk@mailman.ucar.edu><br>

<b>Sent:</b> Monday, November 14, 2022 12:42:27 PM<br>

<b>To:</b> ncl-talk@ucar.edu<br>

<b>Subject:</b> [ncl-talk] Computing geopotential height</font>

<div> </div>

</div>

<div>

<p class="default-style">Hi All,<br>

</p>

<p class="default-style">I am trying to derive the geopotential height of a given CMIP6 model from its air temperature and humidity.

<br>

</p>

<p class="default-style">Firstly, I compute the geopotential in the native hybrid coordinates using the NCL function cz2ccm, then I interpolate the just computed geopotential to some pressure levels using vinth2p_ecmwf.</p>

<p class="default-style">The comparison of the geopotential computed with my script (attached) with the zg downloaded from ESGF shows large differences, so I was wondering if anyone could suggest me what's wrong in my script? </p>

<p class="default-style">Besides, when the variable phis (surface geopotential (m^2/sec^2)) is not available as input, which values I should use? Is it possible to compute it in anyway with NCL?<br>

</p>

<p class="default-style">Thanks into advance for help,<br>

</p>

<p class="default-style">Best,<br>

</p>

<p class="default-style">Markus<br>

</p>

</div>

</body>

</html>