Re: winds from WRF not matching RAOB

From: Tardif, Elliot M <elliot.tardif_at_nyahnyahspammersnyahnyah>
Date: Wed Sep 21 2011 - 07:48:00 MDT

Hi Steve,

I recently had to produce Skew-T plots of WRF output and RAOB data overlaid together in NCL, similar to what you are trying to do. I've attached a template of the NCL script I used as well as a RAOB file that the script read in. Please let me know if this helps you.


Elliot Tardif.

Elliot Tardif, Meteorologist I
NC DENR, Division of Air Quality
Planning Section, Attainment Planning Branch
1641 Mail Service Center
Raleigh, NC 27699-1641
Phone: 919-715-9773
Fax : 919-715-7476
Web :<>

Email correspondence to and from this address is subject to the North Carolina Public Records Law and may be disclosed to third parties unless the content is exempt by statute or other regulation.

From: [] On Behalf Of steve kirby
Sent: Tuesday, September 20, 2011 6:00 PM
Subject: winds from WRF not matching RAOB

Sorry - forgot to note, in the attached figure, WRF winds are in orange; RAOB winds are in black


----- Forwarded Message -----
From: steve kirby <>
To: "" <>
Sent: Tuesday, September 20, 2011 3:54 PM
Subject: winds from WRF not matching RAOB

I'm doing some plots of RAOB's vs WRF and the WRF-extracted winds aren't matching the RAOB. Here's the WRF portion --

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"



;get the variables

;extract u and v from uvm array and convert to knots

 ; get ij point in model domain for location of interest
 ; loc(1) is touth-north (y) and loc(0) is west-east (x)
 ; ensure this point is in the domain, as we are not checking
 ; we are interested in NCL array pointers, so subtract 1
  loc=wrf_user_ll_to_ij(a,-106.93,32.28,True) ; it should read (a, lon,lat,True) is that right?

dataOpts@colTemperature<mailto:dataOpts@colTemperature> = "red"
 dataOpts@colDewPt<mailto:dataOpts@colDewPt> = "magenta"
 ;new location for winds (default is 45)
 dataOpts@xpWind<mailto:dataOpts@xpWind> = 40.
 dataOpts@colWindP<mailto:dataOpts@colWindP> = "orange"

skewT_data = skewT_PlotData (wks,skewt_bkgd,wrfp(:,locY,locX),\


QUESTION - Do I have it right, that in the wrf_user_ll_to_ij( ) call that I use longitude, latitude (in that order) to locate the point of interest?

I have done a couple of other plots of WRF winds vs. RAOB and the winds are not matching up with those either, so I was wondering if I had something wrong in code.


ncl-talk mailing list
List instructions, subscriber options, unsubscribe:

ncl-talk mailing list
List instructions, subscriber options, unsubscribe:

Received on Wed Sep 21 07:48:15 2011

This archive was generated by hypermail 2.1.8 : Mon Oct 10 2011 - 11:43:35 MDT