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

Pressure level interpolation problems

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.

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: 969
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: 97
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