[gradsusr] Make a Skew-T diagram from WRF output

[Ryan Meiszberg]

Hello,

I'm trying to create a Skew-T diagrams using....

"reinit"^M
"open WRFPRS_2011031712_00.ctl"^M
"set t 2"^M
"set lat 18.275"^M
"set lon -67.189"^M
"set lev 1000 200"^M
"define tt=TMPprs-273"^M
"define qq=SPFHprs*1000"^M
"define vel=mag(UGRDprs,VGRDprs)"^M
"define dir=atan2(VGRDprs,UGRDprs)"^M
*rc=plotskew(tt,qq,vel,dir)^M
rc=plotskew(tt,qq,-1,-1)

this is followed by http://moe.met.fsu.edu/~rhart/plotskew.gs

The script seems to run properly for the first half and does plot temperature, but somewhere it is changing the Z dimensions and causes undefined regions for dewpoint and wind. (see output below)

Data Request Warning:  Request beyond file limits
Data Request Warning:  Request beyond file limits
Data Request Warning:  Request beyond file limits
Data Request Warning:  Request beyond file limits
Data Request Warning:  Request beyond file limits
Data Request Warning:  Request beyond file limits
Data Request Warning:  Request beyond file limits
Data Request Warning:  Request beyond file limits
Data Request Warning:  Request beyond file limits
Data Request Warning:  Request beyond file limits
Cannot contour grid - all undefined values 
Drawing temperature sounding.
Drawing dewpoint sounding.
Non-numeric args to numeric operation
  Error occurred on line 2237
  In file plotskew2_new.gs
  Error occurred on line 2252
  In file plotskew2_new.gs
  Error occurred on line 1058
  In file plotskew2_new.gs
  Error occurred on line 12
  In file plotskew2_new.gs


ga-> q dims
Default file number is: 1 
X is fixed     Lon = -67.1913  X = 95.0002
Y is fixed     Lat = 18.2755  Y = 138
Z is fixed     Lev = -9.99e+33  Z = 3.33e+33
T is fixed     Time = 13Z17MAR2011  T = 2

Thanks, 

Ryan Meiszberg
Senoir/Atmospheric Science
Purdue University