Push notifications in your browser are not yet configured.

You are not logged in, you may not see all content and functionalities. If you have an account, please login .

Push notifications in your browser are not yet configured.

## Retreive the geopotential height at 500 hPa from lbfd files

Hello,

I would like to compute the geopotential height at 500 hPa from lbfd files that I have generated using the INT2LM software. In principle, it should be possible using PP, but I do not know how to compute the pressure profile of the reference state. irefatm=2 with delta_t=75, h_scal=10000, posl=100000 and tosl=288.15.

Then I can approximatively estimate the height for which P is 500 hPa. As the height will likely fall between 2 vertical levels, which kind of interpolation do you suggest if I want to refine this first estimate? Using CDO for doing that?

Thank you for your input,

Sylvain

I think the easiest way to get the reference atmosphere is just to run CCLM for a few timesteps and write out both full pressure (P) and perturbation pressure (PP) on model levels. The reference atmosphere would then be obtained by subtracting PP from P. Alternatively, you can copy the code from vgrid_refatm_utils.f90.