Pressure level interpolation problems

This post was from a previous version of the WRF&MPAS-A Support Forum. Please do not add new replies here and if you would like the thread moved out of the Historical / Archive section then contact us, making sure to include the link of the thread to be moved.

rclare2

New member
Hello,

I'm using NCL to create a netCDF file with height data at pressure levels as opposed to eta levels. This has involved the use of the wrf_user_intrp3d() function to interpolate height fields at 19 different pressure levels. The code I've used can be found here: https://github.com/rclare2/ncl-Scripts/blob/master/PressureLevelHeights_netCDF.ncl

The attachment shows my colleague's efforts to plot the data which show expected missing values at lower levels where terrain is in the way. We're particularly interested, however, in the 950 mb level where large areas of missing values appear over eastern Ontario. On this day (March 6, 2005), there is an ordinary low pressure system in that region but we were surprised to find the missing values at 950 mb in an area with such low topography.

Is this a problem stemming from vertical interpolation too close to the lower boundary or something else entirely? Can it be fixed? Is there another, more reliable way to convert WRF data to pressure levels?

Thank you,
Ryan Clare
University of Wisconsin-Madison
 

Attachments

  • case12_hght_950_900hPa.jpg
    case12_hght_950_900hPa.jpg
    75.1 KB · Views: 962

abjaye

New member
Hmmm, I'm not exactly sure why that is happening. My suggestion is to try using a newer ncl interpolator called wrf_user_vert_interp. It has an option to "extrapolate below the ground". Here is the documentation:

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

I have also attached the Fortran source code for wrf_user_interp3d for your information.
 

Attachments

  • wrf_user.f
    21 KB · Views: 87

D_Fora

New member
abjaye said:
Hmmm, I'm not exactly sure why that is happening. My suggestion is to try using a newer ncl interpolator called wrf_user_vert_interp. It has an option to "extrapolate below the ground". Here is the documentation:

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

I have also attached the Fortran source code for wrf_user_interp3d for your information.

I tried to compare the difference between the interpolated 10m wind speed and the 10m wind speed in the WRF model output files (wrfout*).

However I got missing values using wrf_user_vert_interp to get 10m wind speed (agl) with opts@extrapolate set to True already.

Is there any explaination? Thanks a lot.
 
Top