[gradsusr] (no subject)

[Jonathan Wynn Smith]

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.


On Wed, May 21, 2014 at 10:50 PM, James T. Potemra <jimp at hawaii.edu> wrote:

> Jonathan,
>
> The syntax for plotting vectors is such that the length/direction are
> based on the first
> two arguments while the optional third argument is for color.  Thus,
> this line in your script:
>
> 'd skip(a,3,1);skip(b,3,1);skip(c,3,1)'
>
> will plot vectors (length/direction) based on (u,v) that are colored
> based on the values of w.
> If you want vectors of vertical motion, and given that you have a
> north-south section,
> I believe you want to plot (v,w), or
>
> 'd skip(b,3,1);skip(c,3,1)'
>
> Note however that you are likely to have serious issues with scale since
> the horizontal
> extent (lat) is much larger than the vertical.
>
> Jim
>
> On 5/21/14 3:25 PM, Jonathan Wynn Smith wrote:
> > Hello
> >
> > I would like overlay vectors of vertical motion (W) on top of ozone
> > contours in a N-S cross section from WRF model output.
> >
> > Here is my script:
> > * B Vs -VTL
> > 'reinit'
> > *'open /Volumes/JWS/GRADS_OUTPUT/NEW/lnox_bb_p.ctl'
> > 'open /Volumes/JWS/GRADS_OUTPUT/NEW/lnox_only_p.ctl'
> > *'open /Volumes/JWS/GRADS_OUTPUT/NEW/bb_p.ctl'
> > 'open /Volumes/JWS/GRADS_OUTPUT/NEW/ctl_p.ctl'
> >
> > 'set mpdset hires'
> > 'set mpt 1 1 1 4'
> > "set mapdraw on"
> > 'set display color white'
> > 'run rgbset.gs'
> >
> >
> > 'c'
> > "set parea 1 10 1 7.5"
> > "set LON -20"
> > "set LAT -10 10"
> > "set xlopts 1 10 0.25"
> > "set ylopts 1 10 0.25"
> > "set LEV 1000 100"
> > 'set grads off'
> > * Use 14 for blue/purple color
> > "set rgb 23 0 0 255"
> > "set rgb 24 20 20 255"
> > "set rgb 25 40 40 255"
> > "set rgb 26 60 60 255"
> > "set rgb 27 80 80 255"
> > "set rgb 28 100 100 255"
> > "set rgb 29 120 120 255"
> > "set rgb 30 140 140 255"
> > "set rgb 31 160 160 255"
> > "set rgb 32 180 180 255"
> > "set rgb 33 200 200 255"
> > "set rgb 34 220 220 255"
> > "set rgb 35 240 240 255"
> >
> > * Red light to dark for anamolies
> > "set rgb 36 255 240 240"
> > "set rgb 37 255 230 230"
> > "set rgb 38 255 220 220"
> > "set rgb 39 255 200 200"
> > "set rgb 40 255 180 180"
> > "set rgb 41 255 160 160"
> > "set rgb 42 255 140 140"
> > "set rgb 43 255 120 120"
> > "set rgb 44 255 100 100"
> > "set rgb 45 255 80 80"
> > "set rgb 46 255 60 60"
> > "set rgb 47 255 40 40"
> > "set rgb 48 255 20 20"
> > "set rgb 49 255 0 0"
> >
> > "set gxout shaded"
> > *Ozone color
> > "set clevs -6 -2 0 5 7 10 11 13 15 18 24 30 33"
> > "set ccols 23 28 31 36 38 40 42 44 47 6 14 8 12 7"
> > 'd ave((o3.1-o3.2)*1000,t=289, t=1008)'
> >
> > "set gxout contour"
> > "set line 1 1 2"
> > "set ccols 1"
> > *cl = "-0.01 0.02 0.03"
> > cl = "0 7 10 11 13"
> > 'set clevs 'cl
> > 'set clopts -1 10 0.25'
> > 'set cthick 10'
> > 'set clab masked'
> >
> >
> > 'a = ave(U.1,t=289,t=1008,1)'
> > 'b = ave(V.1,t=289,t=1008,1)'
> >
> > 'c = ave(W.1,t=289, t=1008)'
> > "xcbar 1 10 0.3 0.5 -fw 0.175 -fh 0.175 -line on"
> >
> >
> > "set gxout vector"
> > 'd skip(a,3,1);skip(b,3,1);skip(c,3,1)'
> >
> > What am I doing wrong?
> > Jonathan
> >
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>



-- 
*Jonathan Wynn Smith, Ph.D.*
jonathan.smith at noaa.gov
301-683-0038
Postdoctoral Research Associate
National Academy of Science/National Research Council Associateship
NOAA/NESDIS/STAR/SMCD, #2844
5830 University Research Ct.
College Park, MD 20740-3818


*Lord, you establish peace for us; all that we have accomplished you have
done for us. - Isaiah 26:12*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140523/57fd8f9f/attachment-0001.html