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.

F

federico

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?

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.
 
You must log in or register to reply here.
Back
Top