
<!doctype html>

<html>

 <head> 

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

 </head>

 <body> 

  <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> 

 </body>

</html>