winds from WRF not matching RAOB

From: steve kirby <thin_runner_at_nyahnyahspammersnyahnyah>
Date: Tue Sep 20 2011 - 15:59:45 MDT

Sorry  - forgot to note, in the attached figure, WRF winds are in orange; RAOB winds are in black   --Steve ----- Forwarded Message ----- From: steve kirby <thin_runner_at_yahoo.com> To: "ncl-talk_at_ucar.edu" <ncl-talk_at_ucar.edu> Sent: Tuesday, September 20, 2011 3:54 PM Subject: [ncl-talk] winds from WRF not matching RAOB Hi, 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" ... a=addfile("./wrfout_d02_2010-07-16_20_00_00.nc","r");         ;get the variables  wrftc=wrf_user_getvar(a,"tc",0)  wrftd=wrf_user_getvar(a,"td",0)  wrfp=wrf_user_getvar(a,"pressure",0)  wrfz=wrf_user_getvar(a,"z",0)  uvm=wrf_user_getvar(a,"uvmet",0)   ;extract u and v from uvm array and convert to knots  u=uvm(0,:,:,:)*1.94386  v=uvm(1,:,:,:)*1.94386   ;;;;  ; 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?   locX=loc(0)-1  locY=loc(1)-1         dataOpts_at_colTemperature = "red"  dataOpts_at_colDewPt = "magenta"  ;new location for winds (default is 45)  dataOpts_at_xpWind = 40.  dataOpts_at_colWindP = "orange"   skewT_data = skewT_PlotData (wks,skewt_bkgd,wrfp(:,locY,locX),\                                             wrftc(:,locY,locX),\                                             wrftd(:,locY,locX),\                                             wrfz(:,locY,locX),\                                             u(:,locY,locX),\                                             v(:,locY,locX),\                                             dataOpts)   draw(skewT_data) ===================================================================================================== 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. Thanks, Steve _______________________________________________ ncl-talk mailing list List instructions, subscriber options, unsubscribe: http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Received on Tue Sep 20 15:59:56 2011

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