<div dir="ltr"><div><div><div>I wrote a script that basically alters the given script in the section for &quot;arbitrary cross sections&quot; in the Grads user documentation:<br><br>&#39;clear&#39;<br>&#39;set grads off&#39; ; &#39;set timelab off&#39;<br>&#39;set zlog off&#39;<br>&#39;set x 1&#39;<br>&#39;set y 1&#39;<br>&#39;set z 1 30&#39;<br>lon1 = -101.2<br>lon2 = -99.0<br>lat1 = 36.75<br>lat2 = 35.6<br>lon = lon1<br>&#39;collect 1 free&#39;<br>while (lon &lt;= lon2)<br>  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)<br>  &#39;collect 1 gr2stn(w,&#39;lon&#39;,&#39;lat&#39;)&#39;<br>  lon = lon + 0.01<br>endwhile<br><br>&#39;set lon &#39;lon1&#39; &#39;lon2<br>&#39;set xaxis &#39;lon1&#39; &#39;lon2<br>&#39;set ylab off&#39;<br>&#39;set gxout shaded&#39;<br>&#39;run dewrgb 1 25 1&#39;<br>&#39;d coll2gr(1,-u)&#39;<br><br></div>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 &quot;cannot plot data - all undefined values.&quot;<br><br></div>Why is this happening?<br clear="all"><br></div>Jeff Duda<br><div><div><div><div><div><br>-- <br><div class="gmail_signature">Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br></div>
</div></div></div></div></div></div>