[ncl-talk] Computing geopotential height

[markus.mingel at libero.it]

Hi All,


I am trying to derive the geopotential height of a given CMIP6 model from its air temperature and humidity. 


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.

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? 

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?


Thanks into advance for help,


Best,


Markus

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20221114/c4b325b3/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Compute_zg.ncl
Type: text/x-ncl
Size: 5603 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20221114/c4b325b3/attachment-0001.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Surface_Geopotential.png
Type: image/png
Size: 248298 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20221114/c4b325b3/attachment-0001.png>