<div dir="ltr"><div><b>"Collect Error: Z or T must be the only varying dimension</b></div><div><b>Collect Error: Z or T must be the only varying dimension</b></div><div><b>Collect Error: Z or T must be the only varying dimension</b></div>
<div><b>Collect Error: Z or T must be the only varying dimension</b></div><div><b>Collect Error: Z or T must be the only varying dimension</b></div><div><b>Collect Error: Z or T must be the only varying dimension</b></div>
<div><b>281.275</b></div><div><b>35.7428571428571</b></div><div><b>Collect Error: Z or T must be the only varying dimension</b></div><div><b>Collect Error: Z or T must be the only varying dimension</b></div><div><b>Collect Error: Z or T must be the only varying dimension</b></div>
<div><b>Collect Error: Z or T must be the only varying dimension</b></div><div><b>Collect Error: Z or T must be the only varying dimension</b></div><div><b>Collect Error: Z or T must be the only varying dimension</b></div>
<div><b>281.288</b></div><div><b>35.7465714285714</b></div><div><b>Operation error: Invalid dimension environment</b></div><div><b> Looping dimension does not vary</b></div><div><b>Operation error: Invalid dimension environment</b></div>
<div><b> Looping dimension does not vary</b></div><div><b>Cannot plot color bar: No shading information"</b></div><div><b><br></b></div><div>Part of the error.</div><div><br></div><div><br></div><div><b><u>Script used: </u></b></div>
<div><br></div><div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">'reinit'</blockquote><div>'open 20110416_d03.ctl'</div>
<div>'set display color white'</div><div>'clear'</div><div>'set map 2 1 5'</div><div>'set annot 1 12'</div><div>'set xlopts 1 6 .17'</div><div>'set ylopts 1 6 .16'</div><div>
'set mpdset hires'</div><div>'set grid off'</div><div>'set timelab off'</div><div>'set csmooth on'</div><div>*'set gxout shaded'</div><div>'set t 83'</div><div><br></div><div>
*light green to dark green</div><div>'set rgb 31 230 255 225'</div><div>'set rgb 32 200 255 190'</div><div>'set rgb 33 180 250 170'</div><div>'set rgb 34 150 245 140'</div><div>'set rgb 35 120 245 115'</div>
<div>'set rgb 36 80 240 80'</div><div>'set rgb 37 55 210 60'</div><div>'set rgb 38 30 180 30'</div><div>'set rgb 39 15 160 15'</div><div><br></div><div>*light beige to dark brown</div>
<div>'set rgb 71 250 240 230'</div><div>'set rgb 72 240 220 210'</div><div>'set rgb 73 225 190 180'</div><div>'set rgb 74 200 160 150'</div><div>'set rgb 75 180 140 130'</div><div>'set rgb 76 160 120 110'</div>
<div>'set rgb 77 140 100 90'</div><div>'set rgb 78 120 80 70'</div><div>'set rgb 79 100 60 50'</div><div><br></div><div>'set xlint 1.5'</div><div><br></div><div>********************************************************</div>
<div>*</div><div>* Color Shades</div><div>*******************************************************</div><div>'set rgb 59 0 0 0'</div><div>'set rgb 60 15 15 15'</div><div>'set rgb 61 30 30 30'</div>
<div>'set rgb 62 45 45 45'</div><div>'set rgb 63 60 60 60'</div><div>'set rgb 64 75 75 75'</div><div>'set rgb 65 90 90 90'</div><div>'set rgb 66 105 105 105'</div><div>'set rgb 67 120 120 120'</div>
<div>'set rgb 68 135 135 135'</div><div>'set rgb 69 150 150 150'</div><div>'set rgb 70 165 165 165'</div><div>'set rgb 71 180 180 180'</div><div>'set rgb 72 195 195 195'</div><div>'set rgb 73 210 210 210'</div>
<div>'set rgb 74 225 225 225'</div><div>'set rgb 75 245 245 245'</div><div>'set rgb 76 255 255 255'</div><div><br></div><div><br></div><div>*light yellow to dark red</div><div>'set rgb 20 255 250 200'</div>
<div>'set rgb 21 255 250 170'</div><div>'set rgb 22 255 232 120'</div><div>'set rgb 23 255 192 60'</div><div>'set rgb 24 255 160 0'</div><div>'set rgb 25 255 96 0'</div><div>'set rgb 26 255 50 0'</div>
<div>'set rgb 27 225 20 0'</div><div>'set rgb 28 192 0 0'</div><div>'set rgb 29 165 0 0'</div><div><br></div><div>*light blue to dark blue</div><div>'set rgb 81 225 255 255'</div>
<div>'set rgb 82 180 240 250'</div><div>'set rgb 83 150 210 250'</div><div>'set rgb 84 120 185 250'</div><div>'set rgb 85 80 165 245'</div><div>'set rgb 86 60 150 245'</div><div>'set rgb 87 40 130 240'</div>
<div>'set rgb 88 30 110 235'</div><div>'set rgb 89 20 100 210'</div><div><br></div><div>'set rgb 40 0 0 0'</div><div>'set rgb 41 0 236 236'</div><div>'set rgb 42 1 160 246'</div><div>
'set rgb 43 0 0 246'</div><div>'set rgb 44 0 255 0'</div><div>'set rgb 45 0 200 0'</div><div>'set rgb 46 0 144 0'</div><div>'set rgb 47 255 255 0'</div><div>'set rgb 48 231 192 0'</div>
<div>'set rgb 49 255 144 0'</div><div>'set rgb 50 255 0 0'</div><div>'set rgb 51 214 0 0'</div><div>'set rgb 52 192 0 0'</div><div>'set rgb 53 255 0 255'</div><div>'set rgb 54 153 85 201'</div>
<div>'set rgb 55 255 255 255'</div><div><br></div><div>* These are the BLUE shades</div><div>'set rgb 16 0 0 255'</div><div>'set rgb 17 55 55 255'</div><div>'set rgb 18 110 110 255'</div>
<div>'set rgb 19 165 165 255'</div><div>'set rgb 20 220 220 255'</div><div><br></div><div><br></div><div>'c'</div><div>'set background 0'</div><div>'set line 0'</div><div>'draw rec 0 0 11 8.5'</div>
<div><br></div><div>'set lev 1000 100'</div><div>*'define prs = 100'</div><div>*'qv = qvapor'</div><div>*'th = theta'</div><div>*"define e=0.01*prs*qv/(0.6219718+qv)+1e-20"</div><div>
*"define t=th*pow((prs/100000.0),(287.04/1005.7))"</div><div>*"define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805)"</div><div>*"define the=t*pow(100000.0/prs,0.2854*(1.0-0.28*qv))*exp(((3376.0/tlcl)-2.54)*qv*(1.0+0.81*qv))"</div>
<div>"define vvort= hcurl(u,v)"</div><div>*******************************************************************************************</div><div><b>'set x 1 50'</b></div><div><b>'set y 10'</b></div>
<div>*'set lev 1000 100'</div><div>'set lat 35.65 35.75'</div><div>'set lon 280.95 281.3'</div><div>*'set gxout shaded'</div><div>'set parea 1 9.25 1 7.75'</div><div><br></div><div>
* ARBITRARY CROSS-SECTION defined by lon1-lat1 lon2-lat2</div><div>lon1=280.95</div><div>lon2=281.3</div><div>lat1=35.65</div><div>lat2=35.75</div><div><br></div><div>*Vector winds projected in the plane</div><div>m=(lat2-lat1)/(lon2-lon1)</div>
<div>angle=math_atan(m)</div><div>xproj = math_cos(angle)</div><div>yproj = math_sin(angle)</div><div><br></div><div>lon=lon1</div><div>lat=lat1</div><div>'collect 1 free'</div><div>'collect 2 free'</div><div>
'collect 3 free'</div><div>'collect 4 free'</div><div>*'collect 5 free'</div><div>'collect 6 free'</div><div>'collect 7 free'</div><div><br></div><div>while (lon<=lon2)</div><div>
'collect 1 gr2stn(w,'lon','lat')'</div><div> 'collect 2 gr2stn(u,'lon','lat')'</div><div> 'collect 3 gr2stn(v,'lon','lat')'</div><div> 'collect 4 gr2stn(dbz,'lon','lat')'</div>
<div>* 'collect 5 gr2stn(the,'lon','lat')'</div><div> 'collect 6 gr2stn(mag(u,v),'lon','lat')'</div><div> 'collect 7 gr2stn(vvort,'lon','lat')'</div>
<div> say lon</div><div> say lat</div><div> lon=lon+0.013</div><div> lat=lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)</div><div>endwhile</div><div><br></div><div>'set arrlab off'</div><div>'set grads off'</div>
<div>'set x 0 48'</div><div><br></div><div>* Horizotal windspeed and potential temparature</div><div>*'set lon 'lon1' 'lon2</div><div>*'set xaxis 'lon1-360' 'lon2-360' 0.2'</div>
<div>*'set parea 1 9.25 1 7.75'</div><div>*'set vrange 10 'height</div><div>*'set yaxis 0 'hgt' 0.5'</div><div>*'set clab off'</div><div>*'set xlopts 1 7 0.18'</div><div>*'set ylopts 1 7 0.18'</div>
<div>*'set clevs -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5'</div><div>*'set ccols 49 47 46 45 44 43 42 41 21 22 24 25 27 29'</div><div><br></div><div>*THETA-E Calculation</div><div>*'set clevs 350 355 360 365 367 370 375 380 385 390 395'</div>
<div>*'set ccols 18 19 20 21 22 23 24 25 26 27 28 29'</div><div>*'d coll2gr(5)'</div><div>*'draw title 8/29 1200UTC Theta-E(K) (29.8,-89.8 to 29.0,-88.0)'</div><div><br></div><div>* Display reflectivity</div>
<div>'set gxout shaded'</div><div>'set ccols 18 19 20 42 41 44 45 46 47 48 49 50 51 52 53 54'</div><div>'set clevs -10 -5 0 5 10 15 20 25 30 35 40 45 50 55 60'</div><div>'d coll2gr(4)'</div>
<div>*'draw title 8/29 1200UTC Ref.(dBz) (29.8,-89.8 to 29.0,-88.0)'</div><div><br></div><div>*W-wind</div><div>*'set gxout contour'</div><div>*'set cint 5'</div><div>*'set ccolor 1'</div><div>
*'set cmin 0'</div><div>*'d coll2gr(1)'</div><div><br></div><div>* Winds</div><div>*'set ccols 21 22 23 24 25 26 27 28 29 30'</div><div>*'set clevs 20 25 30 35 40 45 50 55 60 65'</div><div>
*'d coll2gr(6)'</div><div><br></div><div><br></div><div>*Vertical Vorticity</div><div>'set gxout contour'</div><div>'d coll2gr(7)'</div><div><br></div><div><br></div><div>* Display vector winds projected into plane</div>
<div>'set ccolor 1'</div><div>*'set arrscl 0.2 4' </div><div>*'set cthick 6'</div><div>*'d skip(umotion-coll2gr(2,-u)*'%xproj%'+coll2gr(3,-u)*'%yproj%',2,2);vmotion-coll2gr(1,-u)*10'</div>
<div><br></div><div><br></div><div>'run cbarn'</div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
'set grads off'</blockquote><div><br></div></div><div><br></div><div><br></div><div>I speculate what's highlighted maybe the error in the x range may not be the appropriate range. Everything else seems ok. I am trying to display vertical vorticity in contours.</div>
<div><br></div><div>I know it states in the coll2gr function that the "<span style="background-color:rgb(224,240,255);color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium">The X axis of the output grid will contain the equally spaced station profiles and will span the range of the current X dimension environment.</span><span style="background-color:rgb(224,240,255);color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium"> " </span></div>
<div><br></div><div><br></div><div><br></div><div>At any rate...I would really like your assistance on how to resolve this problem because I am totally stuck. </div><div><br></div><div>Thanks!</div><div><br></div><div><br>
</div>-- <br>Nimrod Micael
</div>