<div dir="ltr">Hello GrADS community,<div><br></div><div>    I am trying to plot an arbitrary cross section of RH from lev = 900 to lev = 800. My .ctl file has the following levels in between: 900mb,850mb,825mb,820mb,810mb,800mb (Im working in a high elevation and trying to get as close to just above the ground as possible without going into the ground). If I plot RH individually say, for example, at 820mb by doing 'd rh(z = 11)', where z = 11 corresponds to 820mb, it works fine but if I use the arbitrary cross section method and have z vary from (in my case) z = 8 to z = 13 I get the error "World Coordinates convert to non-integer grid coordinates". Does anyone have a solution to this or a workaround?I am still new to GrADS so perhaps the answer is obvious here.<br clear="all"><div><br></div><div><br></div><div><br></div><div>Here is my script if it is of any use:</div><div><br></div><div>'open (file).ctl'<br>'set display white'<br>*sets the time<br>'set t 1'</div><div>*sets the cross section (note the exact values were modified for the problem purposes)</div><div>lon1 = -112<br>lon2 = -110<br>lat1 = 35.5<br>lat2 = 34.5<br>lon = lon1</div><div>*frees collection 1 from memory<br>'collect 1 free'<br>while (lon <= lon2)</div><div>*sets the vertical range<br>'set z 8 13'</div><div>*iterates the latitude <br>lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)<br>'set lon 'lon<br>'set lat 'lat</div><div>*Creates collection<br>'collect 1 gr2stn(rh,'lon','lat')'</div><div>*iterates the longitude by an increment based on the file resolution<br>lon = lon + 0.008<br>endwhile<br>'set lon 'lon1' 'lon2   </div><div>*displays "station" data as shaded gridded data<br>'set gxout shaded'<br>'d coll2gr(1,-u)'<br>'cbarn'<br>'close 1'<br></div><div><br></div><div><br></div><div>Any help/suggestions would be greatly appreciated.</div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><br></div><div dir="ltr">-Thank you,<div>Jon</div></div></div></div></div>