[gradsusr] Problem with arbitrary cross section

[Jeff Duda]

I wrote a script that basically alters the given script in the section for
"arbitrary cross sections" in the Grads user documentation:

'clear'
'set grads off' ; 'set timelab off'
'set zlog off'
'set x 1'
'set y 1'
'set z 1 30'
lon1 = -101.2
lon2 = -99.0
lat1 = 36.75
lat2 = 35.6
lon = lon1
'collect 1 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(w,'lon','lat')'
  lon = lon + 0.01
endwhile

'set lon 'lon1' 'lon2
'set xaxis 'lon1' 'lon2
'set ylab off'
'set gxout shaded'
'run dewrgb 1 25 1'
'd coll2gr(1,-u)'

Here, w is a field from a WRF history dump file and displays just fine in
any 2D cross section. The data certainly exists within the bounds of the
lat, lon, and height levels. However, the output from this script just
gives me "cannot plot data - all undefined values."

Why is this happening?

Jeff Duda

-- 
Jeff Duda
Graduate research assistant
University of Oklahoma School of Meteorology
Center for Analysis and Prediction of Storms
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20151111/0c094722/attachment-0001.html