<div dir="ltr">I have storm I am analyzing but the vectors are not displaying appropriately. I am trying to show the appropriate updrafts and downdrafts of the storm.<div><br></div><div><br></div><div>My code:</div><div><br>
</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">&#39;reinit&#39;<br><blockquote>
&#39;open 20110416_d03.ctl&#39;<br>&#39;set display color white&#39;<br>&#39;clear&#39;<br>&#39;set map 2 1 5&#39;<br>&#39;set annot 1 12&#39;<br>&#39;set xlopts 1 6 .17&#39;<br>&#39;set ylopts 1 6 .16&#39;<br>&#39;set mpdset hires&#39;<br>
&#39;set grid off&#39;<br>&#39;set timelab off&#39;<br>&#39;set csmooth on&#39;<br>*&#39;set gxout shaded&#39;<br>&#39;set t 84&#39;<br>*light green to dark green<br>&#39;set rgb 31 230 255 225&#39;<br>&#39;set rgb 32 200 255 190&#39;<br>
&#39;set rgb 33 180 250 170&#39;<br>&#39;set rgb 34 150 245 140&#39;<br>&#39;set rgb 35 120 245 115&#39;<br>&#39;set rgb 36  80 240  80&#39;<br>&#39;set rgb 37  55 210  60&#39;<br>&#39;set rgb 38  30 180  30&#39;<br>&#39;set rgb 39  15 160  15&#39;<br>
*light beige to dark brown<br>&#39;set rgb 71 250 240 230&#39;<br>&#39;set rgb 72 240 220 210&#39;<br>&#39;set rgb 73 225 190 180&#39;<br>&#39;set rgb 74 200 160 150&#39;<br>&#39;set rgb 75 180 140 130&#39;<br>&#39;set rgb 76 160 120 110&#39;<br>
&#39;set rgb 77 140 100  90&#39;<br>&#39;set rgb 78 120  80  70&#39;<br>&#39;set rgb 79 100  60  50&#39;<br>&#39;set xlint 1.5&#39;<br>********************************************************<br>*<br>*          Color Shades<br>
*******************************************************<br>&#39;set rgb 59 0 0 0&#39;<br>&#39;set rgb 60 15 15 15&#39;<br>&#39;set rgb 61 30 30 30&#39;<br>&#39;set rgb 62 45 45 45&#39;<br>&#39;set rgb 63 60 60 60&#39;<br>
&#39;set rgb 64 75 75 75&#39;<br>&#39;set rgb 65 90 90 90&#39;<br>&#39;set rgb 66 105 105 105&#39;<br>&#39;set rgb 67 120 120 120&#39;<br>&#39;set rgb 68 135 135 135&#39;<br>&#39;set rgb 69 150 150 150&#39;<br>&#39;set rgb 70 165 165 165&#39;<br>
&#39;set rgb 71 180 180 180&#39;<br>&#39;set rgb 72 195 195 195&#39;<br>&#39;set rgb 73 210 210 210&#39;<br>&#39;set rgb 74 225 225 225&#39;<br>&#39;set rgb 75 245 245 245&#39;<br>&#39;set rgb 76 255 255 255&#39;<br>*light yellow to dark red<br>
&#39;set rgb 20 255 250 200&#39;<br>&#39;set rgb 21 255 250 170&#39;<br>&#39;set rgb 22 255 232 120&#39;<br>&#39;set rgb 23 255 192  60&#39;<br>&#39;set rgb 24 255 160   0&#39;<br>&#39;set rgb 25 255  96   0&#39;<br>&#39;set rgb 26 255  50   0&#39;<br>
&#39;set rgb 27 225  20   0&#39;<br>&#39;set rgb 28 192   0   0&#39;<br>&#39;set rgb 29 165   0   0&#39;<br>*light blue to dark blue<br>&#39;set rgb 81 225 255 255&#39;<br>&#39;set rgb 82 180 240 250&#39;<br>&#39;set rgb 83 150 210 250&#39;<br>
&#39;set rgb 84 120 185 250&#39;<br>&#39;set rgb 85  80 165 245&#39;<br>&#39;set rgb 86  60 150 245&#39;<br>&#39;set rgb 87  40 130 240&#39;<br>&#39;set rgb 88  30 110 235&#39;<br>&#39;set rgb 89  20 100 210&#39;<br>&#39;set rgb 40 0 0 0&#39;<br>
&#39;set rgb 41 0 236 236&#39;<br>&#39;set rgb 42 1 160 246&#39;<br>&#39;set rgb 43 0 0 246&#39;<br>&#39;set rgb 44 0 255 0&#39;<br>&#39;set rgb 45 0 200 0&#39;<br>&#39;set rgb 46 0 144 0&#39;<br>&#39;set rgb 47 255 255 0&#39;<br>
&#39;set rgb 48 231 192 0&#39;<br>&#39;set rgb 49 255 144 0&#39;<br>&#39;set rgb 50 255 0 0&#39;<br>&#39;set rgb 51 214 0 0&#39;<br>&#39;set rgb 52 192 0 0&#39;<br>&#39;set rgb 53 255 0 255&#39;<br>&#39;set rgb 54 153 85 201&#39;<br>
&#39;set rgb 55 255 255 255&#39;<br>* These are the BLUE shades<br>&#39;set rgb 16   0   0 255&#39;<br>&#39;set rgb 17  55  55 255&#39;<br>&#39;set rgb 18 110 110 255&#39;<br>&#39;set rgb 19 165 165 255&#39;<br>&#39;set rgb 20 220 220 255&#39;<br>
&#39;c&#39;<br>&#39;set background 0&#39;<br>&#39;set line 0&#39;<br>&#39;draw rec 0 0 11 8.5&#39;<br>&#39;set lev 1000 100&#39;<br>*&#39;define prs = 100&#39;<br>*&#39;qv = qvapor&#39;<br>*&#39;th = theta&#39;<br>*&quot;define e=0.01*prs*qv/(0.6219718+qv)+1e-20&quot;<br>
*&quot;define t=th*pow((prs/100000.0),(287.04/1005.7))&quot;<br>*&quot;define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805)&quot;<br>*&quot;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))&quot;<br>
&quot;define vvort= hcurl(u,v)&quot;<br>*******************************************************************************************<br>&#39;set x 1&#39;<br>&#39;set y 10&#39;<br>*&#39;set lev 1000 100&#39;<br>*&#39;set lat 35.65 35.75&#39;<br>
*&#39;set lon 280.95 281.3&#39;<br>*&#39;set gxout shaded&#39;<br>&#39;set parea 1 9.25 1 7.75&#39;<br>* ARBITRARY CROSS-SECTION defined by lon1-lat1 lon2-lat2<br><b>lon1=281.2 (78.8)<br></b><b>lon2=281.36 (78.64)<br></b><b>lat1=35.78<br>
</b><b>lat2=35.86<br></b>*Vector winds projected in the plane<br>m=(lat2-lat1)/(lon2-lon1)<br>angle=math_atan(m)<br>xproj = math_cos(angle)<br>yproj = math_sin(angle)<br>lon=lon1<br>lat=lat1<br>&#39;collect 1 free&#39;<br>
&#39;collect 2 free&#39;<br>&#39;collect 3 free&#39;<br>&#39;collect 4 free&#39;<br>*&#39;collect 5 free&#39;<br>&#39;collect 6 free&#39;<br>&#39;collect 7 free&#39;<br>while (lon&lt;=lon2)<br>     &#39;collect 1 gr2stn(w,&#39;lon&#39;,&#39;lat&#39;)&#39;<br>
     &#39;collect 2 gr2stn(u,&#39;lon&#39;,&#39;lat&#39;)&#39;<br>     &#39;collect 3 gr2stn(v,&#39;lon&#39;,&#39;lat&#39;)&#39;<br>     &#39;collect 4 gr2stn(dbz,&#39;lon&#39;,&#39;lat&#39;)&#39;<br>*    &#39;collect 5 gr2stn(the,&#39;lon&#39;,&#39;lat&#39;)&#39;<br>
    &#39;collect 6 gr2stn(mag(u,v),&#39;lon&#39;,&#39;lat&#39;)&#39;<br>    &#39;collect 7 gr2stn(vvort,&#39;lon&#39;,&#39;lat&#39;)&#39;<br>     say lon<br>     say lat<br>     lon=lon+0.013<br>     lat=lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)<br>
endwhile<br>&#39;set arrlab off&#39;<br>&#39;set grads off&#39;<br>&#39;set x 1 50&#39;<br>&quot;* Display vector winds projected into plane<br>&#39;set gxout vector&#39;<br>&#39;set ccolor 1&#39;<br>&#39;set arrscl 0.2 7&#39; <br>
&#39;set cthick 6&#39;<br>&#39;d skip(coll2gr(2,-u)*&#39;%xproj%&#39;+coll2gr(3,-u)*&#39;%yproj%&#39;,2,2);coll2gr(1,-u)*10&#39;&quot;<br>So I dont see what I did wrong. Here&#39;s the image. </blockquote></blockquote></div>
<div><div><br></div><div>Thanks for the help. </div><div><br></div><div><br></div><div>Also the range that is being displayed is incorrect</div><div><br></div><div>It should be from 78.64 W to 78.8 W NOT 84.05W to 83.55 W.</div>
<div><br></div><div>I dont understand why it&#39;s not displaying that.   </div><div><br></div><div><br></div><div>Thanks! </div><div><br></div><div><br></div>-- <br><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">
Nimrod Micael
</blockquote></div></div>