Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Interpolation on pressure surfaces

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.


New member
Dear WRF Users,

I would like to interpolate the geopotential height on 500 hPa pressure surface.
To do that I read the PB and P variables from NETCDF output file and PH and PHB for geopotential. I divide the geopotential (both PH and PHB) by gravity to have the geopotential height.
To interpolate the geopotential height to the 500 hPa I assume that for a grid point i,j,k the pressure is PB(i,j,k)+P(i,j,k) and the geopotential is (ph(i,j,k)+phb(i,j,k))/gravity. Is that correct?

Hi Stefano,

You can use the funtion wrf_user_vert_interp to do this automatically!

Let me know if you have further questions.