<div dir="ltr">Is there any other way. I just want to show the simple up/down pointing arrow for upward/downward motion alon the north-south transect.<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, May 21, 2014 at 10:50 PM, James T. Potemra <span dir="ltr"><<a href="mailto:jimp@hawaii.edu" target="_blank">jimp@hawaii.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Jonathan,<br>
<br>
The syntax for plotting vectors is such that the length/direction are<br>
based on the first<br>
two arguments while the optional third argument is for color. Thus,<br>
this line in your script:<br>
<br>
'd skip(a,3,1);skip(b,3,1);skip(c,3,1)'<br>
<br>
will plot vectors (length/direction) based on (u,v) that are colored<br>
based on the values of w.<br>
If you want vectors of vertical motion, and given that you have a<br>
north-south section,<br>
I believe you want to plot (v,w), or<br>
<br>
'd skip(b,3,1);skip(c,3,1)'<br>
<br>
Note however that you are likely to have serious issues with scale since<br>
the horizontal<br>
extent (lat) is much larger than the vertical.<br>
<br>
Jim<br>
<div class="HOEnZb"><div class="h5"><br>
On 5/21/14 3:25 PM, Jonathan Wynn Smith wrote:<br>
> Hello<br>
><br>
> I would like overlay vectors of vertical motion (W) on top of ozone<br>
> contours in a N-S cross section from WRF model output.<br>
><br>
> Here is my script:<br>
> * B Vs -VTL<br>
> 'reinit'<br>
> *'open /Volumes/JWS/GRADS_OUTPUT/NEW/lnox_bb_p.ctl'<br>
> 'open /Volumes/JWS/GRADS_OUTPUT/NEW/lnox_only_p.ctl'<br>
> *'open /Volumes/JWS/GRADS_OUTPUT/NEW/bb_p.ctl'<br>
> 'open /Volumes/JWS/GRADS_OUTPUT/NEW/ctl_p.ctl'<br>
><br>
> 'set mpdset hires'<br>
> 'set mpt 1 1 1 4'<br>
> "set mapdraw on"<br>
> 'set display color white'<br>
> 'run <a href="http://rgbset.gs" target="_blank">rgbset.gs</a>'<br>
><br>
><br>
> 'c'<br>
> "set parea 1 10 1 7.5"<br>
> "set LON -20"<br>
> "set LAT -10 10"<br>
> "set xlopts 1 10 0.25"<br>
> "set ylopts 1 10 0.25"<br>
> "set LEV 1000 100"<br>
> 'set grads off'<br>
> * Use 14 for blue/purple color<br>
> "set rgb 23 0 0 255"<br>
> "set rgb 24 20 20 255"<br>
> "set rgb 25 40 40 255"<br>
> "set rgb 26 60 60 255"<br>
> "set rgb 27 80 80 255"<br>
> "set rgb 28 100 100 255"<br>
> "set rgb 29 120 120 255"<br>
> "set rgb 30 140 140 255"<br>
> "set rgb 31 160 160 255"<br>
> "set rgb 32 180 180 255"<br>
> "set rgb 33 200 200 255"<br>
> "set rgb 34 220 220 255"<br>
> "set rgb 35 240 240 255"<br>
><br>
> * Red light to dark for anamolies<br>
> "set rgb 36 255 240 240"<br>
> "set rgb 37 255 230 230"<br>
> "set rgb 38 255 220 220"<br>
> "set rgb 39 255 200 200"<br>
> "set rgb 40 255 180 180"<br>
> "set rgb 41 255 160 160"<br>
> "set rgb 42 255 140 140"<br>
> "set rgb 43 255 120 120"<br>
> "set rgb 44 255 100 100"<br>
> "set rgb 45 255 80 80"<br>
> "set rgb 46 255 60 60"<br>
> "set rgb 47 255 40 40"<br>
> "set rgb 48 255 20 20"<br>
> "set rgb 49 255 0 0"<br>
><br>
> "set gxout shaded"<br>
> *Ozone color<br>
> "set clevs -6 -2 0 5 7 10 11 13 15 18 24 30 33"<br>
> "set ccols 23 28 31 36 38 40 42 44 47 6 14 8 12 7"<br>
> 'd ave((o3.1-o3.2)*1000,t=289, t=1008)'<br>
><br>
> "set gxout contour"<br>
> "set line 1 1 2"<br>
> "set ccols 1"<br>
> *cl = "-0.01 0.02 0.03"<br>
> cl = "0 7 10 11 13"<br>
> 'set clevs 'cl<br>
> 'set clopts -1 10 0.25'<br>
> 'set cthick 10'<br>
> 'set clab masked'<br>
><br>
><br>
> 'a = ave(U.1,t=289,t=1008,1)'<br>
> 'b = ave(V.1,t=289,t=1008,1)'<br>
><br>
> 'c = ave(W.1,t=289, t=1008)'<br>
> "xcbar 1 10 0.3 0.5 -fw 0.175 -fh 0.175 -line on"<br>
><br>
><br>
> "set gxout vector"<br>
> 'd skip(a,3,1);skip(b,3,1);skip(c,3,1)'<br>
><br>
> What am I doing wrong?<br>
> Jonathan<br>
><br>
<br>
</div></div><div class="HOEnZb"><div class="h5">_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr"><span style="color:rgb(0,0,255)"><font><b>Jonathan Wynn Smith, Ph.D.</b></font></span><br><font size="1"><a href="mailto:jonathan.smith@noaa.gov" target="_blank">jonathan.smith@noaa.gov</a><br>
301-683-0038<br>Postdoctoral Research Associate <br>National Academy of Science/National Research Council Associateship<br>NOAA/NESDIS/STAR/SMCD, #2844<br>5830 University Research Ct.<br>College Park, MD 20740-3818</font><br>
<font size="1"><i><br>Lord, you establish peace for us; all that we have accomplished you have <br>done for us. - Isaiah 26:12</i></font><br><br><br><br><br><br><br></div>
</div>