Interpolation on pressure surfaces

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?

Thanks.
Stefano

Hi Stefano,

You can use the funtion wrf_user_vert_interp to do this automatically!

http://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_vert_interp.shtml

Let me know if you have further questions.