Meredith,<br>I also use standard grads background maps for most of my work, but have increasingly taken advantage of the new shapefiles capabilities. Attached script file includes a simple example of using a shapefile to create a basic background map. There are detailed instructions for using shapefiles in the GrADS Documentation Index.<br>
<br>Unfortunately, I don't have any experience with google maps. Hopefully someone else who does will furnish a more helpful response. You might also check out the "basemap" script which should be included in your library.<br>
<br>Best wishes...<br>Stephen Mc<br><br><div class="gmail_quote">On Wed, Nov 10, 2010 at 4:26 PM, Meredith Croke <span dir="ltr"><<a href="mailto:mcroke@airdat.com">mcroke@airdat.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">Stephen,</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">I typically just use the standard grads background, however I
am very interested in how you create other base/background maps and use
them. I followed the previous email chain and was able to plot my data on
nwcoast.png, but now I'm curious how to create my own background images from
google maps. </font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">Thanks,</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">Meredith</font></span></div>
<div dir="ltr" align="left">
<hr>
</div>
<div dir="ltr" align="left"><font face="Tahoma" size="2"><b>From:</b>
<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>] <b>On
Behalf Of </b>Stephen McMillan<br><b>Sent:</b> Tuesday, November 09, 2010 2:25
PM<br><b>To:</b> GrADS Users Forum<br><b>Subject:</b> Re: [gradsusr] GrADS help,
plotting shaded contour precip/ptype products over a
basemap/background<br></font><br></div>
<div></div>Jeff,
<div>I will provide a response later on this evening when I have more time
outside my work hours. It might also help if you send me a copy of your
google basemap, if you didn't already attach earlier.</div>
<div>Stephen Mc<br><br>
<div class="gmail_quote">On Tue, Nov 9, 2010 at 1:55 PM, Jeff Chabot <span dir="ltr"><<a href="mailto:jsc219@gmail.com" target="_blank">jsc219@gmail.com</a>></span>
wrote:<br>
<blockquote style="border-left: 1px solid rgb(204, 204, 204); margin: 0px 0px 0px 0.8ex; padding-left: 1ex;" class="gmail_quote">Hello Stephen,<br><br>Thanks for your code and sample
images. That is definitely leading me down the right path. I have
attached my attempt to reproduce your work, called test_4.png. Please
note, that is only a work in progress, the image isn't perfect, but it proves
that this works. Using the same code, it didn't work on my google map
basemap, same problem white map covering the basemap, but it worked on your
GrADS produced basemap. Can I ask you how you produced the basemap and
could you share your code to produce that basemap called
USNE.png?<br><br>Thanks again,<br><br><br>Jeff<br><br><br>
<div class="gmail_quote">
<div>On Tue, Nov 9, 2010 at 7:06 AM, <span dir="ltr"><<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>></span> wrote:<br></div>
<blockquote style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;" class="gmail_quote">
<div>Send gradsusr mailing list submissions to<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br><br>To subscribe or unsubscribe
via the World Wide Web, visit<br> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>or, via
email, send a message with subject or body 'help' to<br>
<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br><br>You can reach the
person managing the list at<br> <a href="mailto:gradsusr-owner@gradsusr.org" target="_blank">gradsusr-owner@gradsusr.org</a><br><br>When replying, please
edit your Subject line so it is more specific<br>than "Re: Contents of
gradsusr digest..."<br><br><br>Today's Topics:<br><br> 1. Re: GrADS
help, plotting shaded contour precip/ptype products<br></div>
over a basemap/background (Stephen
McMillan)<br><br><br>----------------------------------------------------------------------<br><br>Message:
1<br>Date: Tue, 9 Nov 2010 11:11:33 -0500<br>From: Stephen McMillan <<a href="mailto:smcmillan@planalytics.com" target="_blank">smcmillan@planalytics.com</a>>
<div><br>Subject: Re: [gradsusr] GrADS help, plotting shaded
contour<br> precip/ptype products over a
basemap/background<br>To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>Message-ID:<br></div>
<AANLkTinozhUHTUAXzsO78rxVciLwtO3M=<a href="mailto:p9xrGEgU-Zq@mail.gmail.com" target="_blank">p9xrGEgU-Zq@mail.gmail.com</a>>
<div>
<div></div>
<div><br>Content-Type: text/plain;
charset="iso-8859-1"<br><br>Jeff,<br>Thanks for sending the image. I
tried something very similar using the same<br>colors, but only with precip.
Resulting image attached, along with the<br>short script to produce
it. I also attached the "USNE.png" background map<br>used.<br><br>One
of the differences I see between your and my image is that your<br>displayed
(foreground) image borders appear to be transparent while the<br>plotted
white is not. My image sample appears to be true in reverse, but
I<br>don't know why. Perhaps my sample script will give you a
clue.<br><br>Stephen Mc<br><br>On Mon, Nov 8, 2010 at 8:42 PM, Jeff Chabot
<<a href="mailto:jsc219@gmail.com" target="_blank">jsc219@gmail.com</a>>
wrote:<br><br>> Stephen,<br>><br>> Awesome idea, but GrADS still
produced a map with a white background. I<br>> did confirm that the
white is now a pure white: rgb 255 255 255.<br>><br>>
Sample:<br>><br>> <a href="http://jeffsweatherservice.com/grads/nam/test_2.png" target="_blank">http://jeffsweatherservice.com/grads/nam/test_2.png</a><br>><br>>
Printim line used:<br>><br>> printim test_2.png png -t 200 x570 y435
-b NE.png<br>><br>> Colors:<br>> set rgb 200 255 255
255<br>><br>> set rgb 30 210 180 140<br>> set rgb 16 107 159
91<br>> set rgb 17 90 142 75<br>> set rgb 18 74 125 60<br>> set rgb
19 58 109 45<br>> set rgb 20 41 92 19<br>> set rgb 21 25 75 14<br>>
set rgb 22 253 248 2<br>> set rgb 23 229 188 0<br>> set rgb 24 253 149
0<br>> set rgb 25 253 0 0<br>> set rgb 26 212 0 0<br>> set rgb 27
255 0 153<br>><br>> set clevs 0.01 0.05 0.10 0.25 0.50 1.00 2.00 3.00
4.00 5.00<br>> set ccols 200 16 18 19 20 22 23 24 25 26
27<br>><br>> I used 200 because I already had 16 and several other
numbers assigned to<br>> different colors. I tried both png and gif
format to see if it would make a<br>> difference, but it did
not.<br>><br>> Any other ideas? I did use the full/local path to
my images above.<br>><br>> Sincerely,<br>><br>><br>>
Jeff<br>><br>> On Mon, Nov 8, 2010 at 12:16 PM, <<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>>
wrote:<br>><br>>> Send gradsusr mailing list submissions
to<br>>> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>>><br>>> To subscribe
or unsubscribe via the World Wide Web, visit<br>>>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>>>
or, via email, send a message with subject or body 'help' to<br>>>
<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>>><br>>> You
can reach the person managing the list at<br>>>
<a href="mailto:gradsusr-owner@gradsusr.org" target="_blank">gradsusr-owner@gradsusr.org</a><br>>><br>>> When
replying, please edit your Subject line so it is more specific<br>>>
than "Re: Contents of gradsusr
digest..."<br>>><br>>><br>>> Today's
Topics:<br>>><br>>> 1. Re: gradsusr Digest, Vol 9, Issue
18 (Jeff Chabot)<br>>><br>>><br>>>
----------------------------------------------------------------------<br>>><br>>>
Message: 1<br>>> Date: Mon, 8 Nov 2010 16:21:29
-0500<br>>><br>>> From: Jeff Chabot <<a href="mailto:jsc219@gmail.com" target="_blank">jsc219@gmail.com</a>><br>>> Subject: Re: [gradsusr]
gradsusr Digest, Vol 9, Issue 18<br>>><br>>> To: <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>>> Message-ID:<br>>>
<AANLkTi=<a href="mailto:hSyCi59RifNAgKttG0ZFeyMmz79rTGAQ8GUdp@mail.gmail.com" target="_blank">hSyCi59RifNAgKttG0ZFeyMmz79rTGAQ8GUdp@mail.gmail.com</a>><br>>><br>>>
Content-Type: text/plain; charset="iso-8859-1"<br>>><br>>>
Jennifer, Jeff,<br>>><br>>> I was the one that discovered the
NE.png looks to the file ne.png<br>>> (basically<br>>> looks to
lowercase files only). I resolved that problem by renaming
my<br>>> file<br>>> NE.png to ne.png. I can use a
background image now, but when I draw<br>>> shaded<br>>>
contours over that background, I either get a black or white map area
for<br>>> precip =0. I can still see my background outside of
the map area, but I<br>>> would like to show that background inside
the map area for areas that<br>>> precip<br>>> = 0. I
thought -t 0 or -t 00 would allow me to do just that.
Instead,<br>>> when<br>>> I used -t 0 or -t 00, I the map
area is black for precip equals zero.<br>>><br>>> #1: Yes I am
using the file ne.png.<br>>> #2: My map area background is black using
the color 0 and -t 0 or -t 00.<br>>> But, in the printim -t 00, that
should equal transparent, not black. If I<br>>> change -t 0 to
white, then it the black area is now white. I am trying to<br>>>
use a background map behind the shaded contours for precip = 0 instead
of<br>>> white or black.<br>>><br>>> Ultimately, my main
goal here is to use a colorful basemap to my precip<br>>> contour
maps, instead of using just black or white, or any other solid<br>>>
color<br>>> for that fact. And, if anyone knows how to do that,
I would love to know.<br>>><br>>> Once again, thanks for all
your input.<br>>><br>>>
Sincerely,<br>>><br>>><br>>> Jeff
Chabot<br>>><br>>><br>>> On Mon, Nov 8, 2010 at 11:50 AM,
<<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>>
wrote:<br>>><br>>> > Send gradsusr mailing list submissions
to<br>>> > <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>>> ><br>>> > To
subscribe or unsubscribe via the World Wide Web, visit<br>>> >
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>>>
> or, via email, send a message with subject or body 'help'
to<br>>> > <a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>>> ><br>>>
> You can reach the person managing the list at<br>>> >
<a href="mailto:gradsusr-owner@gradsusr.org" target="_blank">gradsusr-owner@gradsusr.org</a><br>>> ><br>>>
> When replying, please edit your Subject line so it is more
specific<br>>> > than "Re: Contents of gradsusr
digest..."<br>>> ><br>>> ><br>>> > Today's
Topics:<br>>> ><br>>> > 1. Re: GrADS help, plotting
shaded contour precip/ptype products<br>>> >
over a basemap/background (Jeffrey Duda)<br>>> ><br>>>
><br>>> >
----------------------------------------------------------------------<br>>>
><br>>> > Message: 1<br>>> > Date: Mon, 8 Nov 2010
14:55:31 -0600<br>>><br>>> > From: Jeffrey Duda <<a href="mailto:jdduda@iastate.edu" target="_blank">jdduda@iastate.edu</a>><br>>> > Subject: Re:
[gradsusr] GrADS help, plotting shaded contour<br>>> >
precip/ptype products over a
basemap/background<br>>> > To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>>> >
Message-ID:<br>>> > <AANLkTi=<a href="mailto:UuhCprPO4JG9zSvZdHJ6fo7e2gLVtP4oLUyFu@mail.gmail.com" target="_blank">UuhCprPO4JG9zSvZdHJ6fo7e2gLVtP4oLUyFu@mail.gmail.com</a>><br>>><br>>>
> Content-Type: text/plain; charset="iso-8859-1"<br>>>
><br>>> > Jeff,<br>>> > 1) Did you try Jennifer's
suggestion? If that doesn't work,<br>>> > 2) what is your
background color?<br>>> ><br>>> > Jeff<br>>>
><br>>><br>>> > On Mon, Nov 8, 2010 at 2:31 PM, Jeff
Chabot <<a href="mailto:jsc219@gmail.com" target="_blank">jsc219@gmail.com</a>> wrote:<br>>> ><br>>>
> > Jeff,<br>>> > ><br>>> > > I liked the
idea, but that didn't help. I tried the following, both<br>>>
> > resulting in a black background for the map area:<br>>> >
><br>>> > > printim precip_NE.png png x570 y435
-b<br>>><br>>> >
/home/jeff/public_html/images/NE.png<br>>> > > -t 00<br>>>
> ><br>>> > > printim precip_NE.png png x570 y435 -t 00
-b<br>>> > > /home/jeff/public_html/images/NE.png<br>>>
> ><br>>> > > Example: <a href="http://jeffsweatherservice.com/grads/nam/precip_NE_2.png" target="_blank">http://jeffsweatherservice.com/grads/nam/precip_NE_2.png</a><br>>>
> ><br>>> > > I even tried using "00" in the color table,
but that had the same<br>>> impact.<br>>> > ><br>>>
> > Now, if I use the following command, the black background
becomes<br>>> white:<br>>> > ><br>>> > >
printim precip_NE.png png x570 y435 -b<br>>><br>>> >
/home/jeff/public_html/images/NE.png<br>>> > > white<br>>>
> ><br>>> > > Thanks for the idea though.<br>>> >
><br>>> > ><br>>> > >
Jeff<br>>><br>>> > ><br>>> > ><br>>>
> ><br>>> > ><br>>> > ><br>>> >
><br>>> > ><br>>> > ><br>>> > > On
Mon, Nov 8, 2010 at 10:59 AM, <<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>><br>>>
wrote:<br>>> > ><br>>> > >> Send gradsusr mailing
list submissions to<br>>> > >> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>>> >
>><br>>> > >> To subscribe or unsubscribe via the World
Wide Web, visit<br>>> > >> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>>>
> >> or, via email, send a message with subject or body 'help'
to<br>>> > >> <a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>>> >
>><br>>> > >> You can reach the person managing the
list at<br>>> > >> <a href="mailto:gradsusr-owner@gradsusr.org" target="_blank">gradsusr-owner@gradsusr.org</a><br>>> >
>><br>>> > >> When replying, please edit your Subject
line so it is more specific<br>>> > >> than "Re: Contents of
gradsusr digest..."<br>>> > >><br>>> >
>><br>>> > >> Today's Topics:<br>>> >
>><br>>> > >> 1. Re: GrADS help, plotting shaded
contour precip/ptype products<br>>> > >>
over a basemap/background (Jeffrey Duda)<br>>> >
>><br>>> > >><br>>> > >><br>>>
----------------------------------------------------------------------<br>>>
> >><br>>> > >> Message: 1<br>>> > >>
Date: Mon, 8 Nov 2010 14:04:13 -0600<br>>><br>>> >
>><br>>> > >> From: Jeffrey Duda <<a href="mailto:jdduda@iastate.edu" target="_blank">jdduda@iastate.edu</a>><br>>> > >> Subject:
Re: [gradsusr] GrADS help, plotting shaded contour<br>>> > >>
precip/ptype products over a
basemap/background<br>>> > >> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>>> > >>
Message-ID:<br>>> > >>
<AANLkTinO=D7vgkJzwo=+60vLO5VgOZjp=<a href="mailto:4Hhr7P93qo%2B@mail.gmail.com" target="_blank">4Hhr7P93qo+@mail.gmail.com</a><<a href="mailto:4Hhr7P93qo%252B@mail.gmail.com" target="_blank">4Hhr7P93qo%2B@mail.gmail.com</a>><br>
>> <<a href="mailto:4Hhr7P93qo%252B@mail.gmail.com" target="_blank">4Hhr7P93qo%2B@mail.gmail.com</a> <<a href="mailto:4Hhr7P93qo%25252B@mail.gmail.com" target="_blank">4Hhr7P93qo%252B@mail.gmail.com</a>>><br>
>> >
<<a href="mailto:4Hhr7P93qo%252B@mail.gmail.com" target="_blank">4Hhr7P93qo%2B@mail.gmail.com</a> <<a href="mailto:4Hhr7P93qo%25252B@mail.gmail.com" target="_blank">4Hhr7P93qo%252B@mail.gmail.com</a>>
<<br></div></div>>> <a href="mailto:4Hhr7P93qo%25252B@mail.gmail.com" target="_blank">4Hhr7P93qo%252B@mail.gmail.com</a> <<a href="mailto:4Hhr7P93qo%2525252B@mail.gmail.com" target="_blank">4Hhr7P93qo%25252B@mail.gmail.com</a>>>>
<div>
<div></div>
<div><br>>><br>>> > >> ><br>>> >
>><br>>> > >> Content-Type: text/plain;
charset="iso-8859-1"<br>>> > >><br>>> > >>
Jeff,<br>>> > >> Try using a two-digit color number (i.e.,
00) in the printim command<br>>> as<br>>> > >>
per<br>>> > >> the recommendation on the index page for the
printim command:<br>>> > >> <a href="http://www.iges.org/grads/gadoc/gradcomdprintim.html" target="_blank">http://www.iges.org/grads/gadoc/gradcomdprintim.html</a><br>>>
> >><br>>> > >> Jeff<br>>> >
>><br>>> > >> On Mon, Nov 8, 2010 at 1:41 PM, Jeff
Chabot <<a href="mailto:jsc219@gmail.com" target="_blank">jsc219@gmail.com</a>><br>>> wrote:<br>>> >
>><br>>> > >> > I am using two separate color
tables. For precip:<br>>> > >> ><br>>> >
>> > *Categorical Rain<br>>> > >> > set rgb 30
210 180 140<br>>> > >> > set rgb 16 107 159 91<br>>>
> >> > set rgb 17 90 142 75<br>>> > >> > set
rgb 18 74 125 60<br>>> > >> > set rgb 19 58 109
45<br>>> > >> > set rgb 20 41 92 19<br>>> >
>> > set rgb 21 25 75 14<br>>> > >> > set rgb 22
253 248 2<br>>> > >> > set rgb 23 229 188 0<br>>>
> >> > set rgb 24 253 149 0<br>>> > >> > set
rgb 25 253 0 0<br>>> > >> > set rgb 26 212 0 0<br>>>
> >> > set rgb 27 255 0 153<br>>> > >>
><br>>> > >> > set clevs 0.01 0.05 0.10 0.25 0.50 1.00
2.00 3.00 4.00 5.00<br>>> > >> > set ccols 0 16 18 19 20
22 23 24 25 26 27<br>>> > >> ><br>>> > >>
> For precip type:<br>>> > >> ><br>>> >
>> > *Categorical Rain<br>>> > >> > set rgb 16
107 159 91<br>>> > >> > set rgb 17 90 142 75<br>>>
> >> > set rgb 18 74 125 60<br>>> > >> > set
rgb 19 58 109 45<br>>> > >> > set rgb 20 41 92
19<br>>> > >> > set rgb 21 25 75 14<br>>> >
>> > set rgb 22 253 248 2<br>>> > >> > set rgb 23
229 188 0<br>>> > >> > set rgb 24 253 149 0<br>>>
> >> > set rgb 25 253 0 0<br>>> > >> > set rgb
26 212 0 0<br>>> > >> > set rgb 27 255 0 153<br>>>
> >> ><br>>> > >> > *Categorical Freezing
Rain<br>>> > >> > set rgb 28 255 188 188<br>>> >
>> > set rgb 29 255 172 172<br>>> > >> > set rgb
30 255 156 156<br>>> > >> > set rgb 31 253 140
140<br>>> > >> > set rgb 32 255 124 124<br>>> >
>> > set rgb 33 255 108 108<br>>> > >> > set rgb
34 255 96 96<br>>> > >> > set rgb 35 255 80 80<br>>>
> >> > set rgb 36 255 64 56<br>>> > >> > set
rgb 37 240 32 128<br>>> > >> > set rgb 38 240 16
255<br>>> > >> ><br>>> > >> >
*Categorical Ice Pellets / Mix<br>>> > >> > set rgb 39 255
200 0<br>>> > >> > set rgb 40 255 180 0<br>>> >
>> > set rgb 41 255 160 0<br>>> > >> > set rgb 42
255 140 0<br>>> > >> > set rgb 43 255 124 0<br>>>
> >> > set rgb 44 255 108 0<br>>> > >> > set
rgb 45 255 96 0<br>>> > >> > set rgb 46 255 80
0<br>>> > >> > set rgb 47 255 64 0<br>>> >
>> > set rgb 48 255 32 0<br>>> > >> > set rgb 49
255 16 0<br>>> > >> ><br>>> > >> >
*Categorical Snow<br>>> > >> > set rgb 50 4 233
231<br>>> > >> > set rgb 51 0 173 255<br>>> >
>> > set rgb 52 0 148 255<br>>> > >> > set rgb 53
0 123 255<br>>> > >> > set rgb 54 0 104 255<br>>>
> >> > set rgb 55 0 85 255<br>>> > >> > set
rgb 56 4 67 245<br>>> > >> > set rgb 57 0 38
255<br>>> > >> > set rgb 58 0 14 255<br>>> >
>> > set rgb 59 0 0 255<br>>> > >> > set rgb 60 0
0 223<br>>> > >> > set rgb 61 255 255 255<br>>> >
>> ><br>>> > >> > set clevs 0.01 0.05 0.1 0.25
0.5 1 2 3 4 5 10 10.01 10.05 10.1<br>>> 10.25<br>>> >
>> 10.5<br>>> > >> > 11 12<br>>> > >>
> 13 14 15 20 20.01 20.05 20.1 20.25 20.5 21 22 23 24 25 30
30.01<br>>> 30.05<br>>> > >> 30.1<br>>> >
>> > 30.25<br>>> > >> > 30.5 31 32 33 34
35<br>>> > >> > set ccols 0 16 18 19 20 22 23 24 25 26 27
0 39 40 41 42 43 44 45 46<br>>> 47<br>>> > 0<br>>> >
>> 0<br>>> > >> > 28 29<br>>> > >>
> 30 31 32 33 34 35 36 0 0 50 51 52 53 54 55 56 57 58
59<br>>> > >> ><br>>> > >> > So, Jeff, I
believe that I am using #0 correctly. But, when I do,<br>>>
it<br>>> > >> just<br>>> > >> > uses white
or black as the background as in this example:<br>>> > >>
><br>>> > >> > <a href="http://jeffsweatherservice.com/grads/nam/ptype_NE_2.gif" target="_blank">http://jeffsweatherservice.com/grads/nam/ptype_NE_2.gif</a><br>>>
> >> ><br>>> > >> > Thanks,<br>>> >
>> ><br>>> > >> ><br>>> > >> >
Jeff<br>>> > >> ><br>>> > >>
><br>>> > >> ><br>>> > >> > On Mon,
Nov 8, 2010 at 10:18 AM, <<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>><br>>> >
wrote:<br>>> > >> ><br>>> > >> >>
Send gradsusr mailing list submissions to<br>>> > >> >>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>>> > >>
>><br>>> > >> >> To subscribe or unsubscribe via
the World Wide Web, visit<br>>> > >> >><br>>>
> >> >> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>>>
> >> >> or, via email, send a message with subject or body
'help' to<br>>> > >> >> <a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>>> > >>
>><br>>> > >> >> You can reach the person
managing the list at<br>>> > >> >>
<a href="mailto:gradsusr-owner@gradsusr.org" target="_blank">gradsusr-owner@gradsusr.org</a><br>>> > >>
>><br>>> > >> >> When replying, please edit your
Subject line so it is more<br>>> specific<br>>> > >>
>> than "Re: Contents of gradsusr digest..."<br>>> > >>
>><br>>> > >> >><br>>> > >>
>> Today's Topics:<br>>> > >> >><br>>> >
>> >> 1. GrADS help, plotting shaded
contour precip/ptype<br>>> products<br>>> > >> >>
over a basemap/background (Jeff
Chabot)<br>>> > >> >> 2. Re: GrADS help,
plotting shaded contour precip/ptype products<br>>> > >>
>> over a basemap/background (Jeffrey
Duda)<br>>> > >> >> 3. Problem with skewt plot
using WRF output (Cristian ...)<br>>> > >>
>><br>>> > >> >><br>>> > >>
>><br>>> >
----------------------------------------------------------------------<br>>>
> >> >><br>>> > >> >> Message:
1<br>>> > >> >> Date: Mon, 8 Nov 2010 13:51:33
-0500<br>>> > >> >> From: Jeff Chabot <<a href="mailto:jsc219@gmail.com" target="_blank">jsc219@gmail.com</a>><br>>> > >> >>
Subject: [gradsusr] GrADS help, plotting shaded contour<br>>>
precip/ptype<br>>> > >> >>
products over a basemap/background<br>>> > >> >>
To: <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>>> > >> >>
Message-ID:<br>>> > >> >>
<AANLkTinKn=<br>>> <a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%2BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf+FgLPz@mail.gmail.com</a><<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%2BFgLPz@mail.gmail.com</a>><br>
>>
<<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%2BFgLPz@mail.gmail.com</a><<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%25252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%252BFgLPz@mail.gmail.com</a>><br>
>>
><br>>> > <<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%2BFgLPz@mail.gmail.com</a><<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%25252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%252BFgLPz@mail.gmail.com</a>><br>
>>
<<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%25252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%252BFgLPz@mail.gmail.com</a><<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%2525252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%25252BFgLPz@mail.gmail.com</a>><br>
>>
><br>>> > ><br>>> > >> <<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%2BFgLPz@mail.gmail.com</a><<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%25252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%252BFgLPz@mail.gmail.com</a>><br>
>>
<<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%25252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%252BFgLPz@mail.gmail.com</a><<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%2525252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%25252BFgLPz@mail.gmail.com</a>><br>
>>
><br>>> > <<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%25252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%252BFgLPz@mail.gmail.com</a><<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%2525252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%25252BFgLPz@mail.gmail.com</a>><br>
</div></div>>>
<<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%2525252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%25252BFgLPz@mail.gmail.com</a><<a href="mailto:eTr6myE00LdycRWRxdyy5e8hpBtf%252525252BFgLPz@mail.gmail.com" target="_blank">eTr6myE00LdycRWRxdyy5e8hpBtf%2525252BFgLPz@mail.gmail.com</a>>
<div>
<div></div>
<div><br>>> ><br>>><br>>> >
><br>>> > >> ><br>>> > >><br>>>
> >> >> ><br>>> > >> >> Content-Type:
text/plain; charset="iso-8859-1"<br>>> > >>
>><br>>> > >> >><br>>> > >>
>> Dear GrADS Users,<br>>> > >> >> I have been
trying to apply a basemap to my precip/ptype shaded<br>>> >
contour<br>>> > >> >> products on my personal website
for some time now. I have<br>>> > brainstormed<br>>>
> >> >> and<br>>> > >> >> searched the
GrADS online documents and have come up with ideas,<br>>>
but<br>>> > >> they<br>>> > >> >> have
never worked. My latest idea was to use Google Maps
terrain<br>>> > data<br>>> > >> as<br>>> >
>> >> a<br>>> > >> >> background to the
map. I thought it was a good idea, but I<br>>> >
discovered<br>>> > >> >> that<br>>> > >>
>> I still cannot resolve the issue in that precip = 0 shows
either<br>>> all<br>>> > >> white<br>>> >
>> >> or all black depending on what I choose for the
background,<br>>> covering<br>>> > up<br>>> >
>> my<br>>> > >> >> basemap. I really
thought the following command would have<br>>> resolved<br>>>
> >> that<br>>> > >> >> issue, but it did
not:<br>>> > >> >><br>>> > >> >>
printim precip_NE_1.png png x570 y435 -b NE.png -t 0<br>>>
Versions:GrADS:<br>>> > >> 2.0a9<br>>> > >>
>> OS: Fedora 13<br>>> > >> >> Sample: <a href="http://jeffsweatherservice.com/grads/nam/ptype_NE_2.png" target="_blank">http://jeffsweatherservice.com/grads/nam/ptype_NE_2.png</a><br>>>
> >> >><br>>> > >> >> I know this can be
done because of the example:<br>>> > >> >><br>>>
> >> >> <a href="http://wxmaps.org/pix/ez.east.html" target="_blank">http://wxmaps.org/pix/ez.east.html</a><br>>> >
>> >><br>>> > >> >> Any assistance here
would be much appreciated.Sincerely,<br>>> > >>
>><br>>> > >> >><br>>> > >>
>> Jeff Chabot<br>>> > >> >> Email: <a href="mailto:jsc219@gmail.com" target="_blank">jsc219@gmail.com</a><br>>> > >> >> Web:
<a href="http://jeffsweatherservice.com" target="_blank">http://jeffsweatherservice.com</a><br>>> > >>
>> -------------- next part --------------<br>>> > >>
>> An HTML attachment was scrubbed...<br>>> > >>
>> URL:<br>>> > >> >><br>>> >
>><br>>> ><br>>> <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20101108/93f97642/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20101108/93f97642/attachment-0001.html</a><br>
>>
> >> >><br>>> > >> >>
------------------------------<br>>> > >>
>><br>>> > >> >> Message: 2<br>>> >
>> >> Date: Mon, 8 Nov 2010 12:57:30 -0600<br>>> >
>> >> From: Jeffrey Duda <<a href="mailto:jdduda@iastate.edu" target="_blank">jdduda@iastate.edu</a>><br>>> > >> >>
Subject: Re: [gradsusr] GrADS help, plotting shaded contour<br>>> >
>> >> precip/ptype products over a
basemap/background<br>>> > >> >> To: GrADS Users Forum
<<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>>> > >>
>> Message-ID:<br>>> > >> >>
<<br>>> <a href="mailto:AANLkTikPmH-Og0-6uNAi4tpUMCoZbAZ3YTqm_HXkndpV@mail.gmail.com" target="_blank">AANLkTikPmH-Og0-6uNAi4tpUMCoZbAZ3YTqm_HXkndpV@mail.gmail.com</a><br>>>
> ><br>>> > >> >> Content-Type: text/plain;
charset="iso-8859-1"<br>>> > >> >><br>>> >
>> >><br>>> > >> >> Jeff,<br>>> >
>> >> What color table are you using? Are you sure you
have a #0 color?<br>>> I<br>>> > >>
know<br>>> > >> >> that's supposed to be the
background, but perhaps you are<br>>> overwriting<br>>> >
>> it<br>>> > >> >> with some other color that
has the same color as your background<br>>> but<br>>> >
a<br>>> > >> >> different color number.<br>>>
> >> >><br>>> > >> >> Jeff
Duda<br>>> > >> >><br>>> > >>
><br>>> > >> ><br>>> > >>
><br>>> > >> >><br>>> > >> >>
On Mon, Nov 8, 2010 at 12:51 PM, Jeff Chabot <<a href="mailto:jsc219@gmail.com" target="_blank">jsc219@gmail.com</a>><br>>> > wrote:<br>>>
> >> >><br>>> > >> >> > Dear GrADS
Users,<br>>> > >> >> > I have been trying to apply a
basemap to my precip/ptype shaded<br>>> > >>
contour<br>>> > >> >> > products on my personal
website for some time now. I have<br>>> > >>
brainstormed<br>>> > >> >> and<br>>> >
>> >> > searched the GrADS online documents and have come up
with ideas,<br>>> > but<br>>> > >> >>
they<br>>> > >> >> > have never worked. My
latest idea was to use Google Maps<br>>> terrain<br>>> >
>> data<br>>> > >> >> as a<br>>> >
>> >> > background to the map. I thought it was a good
idea, but I<br>>> > >> discovered<br>>> > >>
>> that<br>>> > >> >> > I still cannot resolve
the issue in that precip = 0 shows either<br>>> > all<br>>>
> >> >> white<br>>> > >> >> > or all
black depending on what I choose for the background,<br>>> >
covering<br>>> > >> up<br>>> > >> >>
my<br>>> > >> >> > basemap. I really thought
the following command would have<br>>> > resolved<br>>> >
>> >> that<br>>> > >> >> > issue, but it
did not:<br>>> > >> >> ><br>>> > >>
>> > printim precip_NE_1.png png x570 y435 -b NE.png -t
0<br>>> > Versions:GrADS:<br>>> > >> >>
2.0a9<br>>> > >> >> > OS: Fedora 13<br>>> >
>> >> > Sample: <a href="http://jeffsweatherservice.com/grads/nam/ptype_NE_2.png" target="_blank">http://jeffsweatherservice.com/grads/nam/ptype_NE_2.png</a><br>>>
> >> >> ><br>>> > >> >> > I know
this can be done because of the example:<br>>> > >> >>
><br>>> > >> >> > <a href="http://wxmaps.org/pix/ez.east.html" target="_blank">http://wxmaps.org/pix/ez.east.html</a><br>>> >
>> >> ><br>>> > >> >> > Any
assistance here would be much appreciated. Sincerely,<br>>> >
>> >> ><br>>> > >> >> ><br>>>
> >> >> > Jeff Chabot<br>>> > >> >>
> Email: <a href="mailto:jsc219@gmail.com" target="_blank">jsc219@gmail.com</a><br>>> > >> >> >
Web: <a href="http://jeffsweatherservice.com" target="_blank">http://jeffsweatherservice.com</a><br>>> > >>
>> ><br>>> > >> >> ><br>>> >
>> >> >
_______________________________________________<br>>> > >>
>> > gradsusr mailing list<br>>> > >> >> >
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>>> > >> >>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>>>
> >> >> ><br>>> > >> >>
><br>>> > >> >><br>>> > >>
>><br>>> > >> >> --<br>>> > >>
>> Jeff Duda<br>>> > >> >> Iowa State
University<br>>> > >> >> Meteorology Graduate
Student<br>>> > >> >> 3134 Agronomy Hall<br>>>
> >> >> <a href="http://www.meteor.iastate.edu/%7Ejdduda" target="_blank">www.meteor.iastate.edu/~jdduda</a><<a href="http://www.meteor.iastate.edu/%7Ejdduda" target="_blank">http://www.meteor.iastate.edu/%7Ejdduda</a>><br>
>>
<<a href="http://www.meteor.iastate.edu/%7Ejdduda" target="_blank">http://www.meteor.iastate.edu/%7Ejdduda</a>><br>>><br>>>
> <<a href="http://www.meteor.iastate.edu/%7Ejdduda" target="_blank">http://www.meteor.iastate.edu/%7Ejdduda</a>><<br>>>
> >> <a href="http://www.meteor.iastate.edu/%7Ejdduda" target="_blank">http://www.meteor.iastate.edu/%7Ejdduda</a>><br>>>
> >><br>>> > >> >> -------------- next part
--------------<br>>> > >> >> An HTML attachment was
scrubbed...<br>>> > >> >> URL:<br>>> >
>> >><br>>> > >><br>>> ><br>>> <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20101108/159827b2/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20101108/159827b2/attachment-0001.html</a><br>
>>
> >> >><br>>> > >> >>
------------------------------<br>>> > >>
>><br>>> > >> >> Message: 3<br>>> >
>> >> Date: Mon, 8 Nov 2010 16:23:09 -0300<br>>> >
>> >> From: "Cristian ..." <<a href="mailto:galenso85@hotmail.com" target="_blank">galenso85@hotmail.com</a>><br>>> > >>
>> Subject: [gradsusr] Problem with skewt plot using WRF
output<br>>> > >> >> To: <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>>> > >>
>> Message-ID:
<SNT112-W63DF5D09C452E1C6FB44FBF4F0@phx.gbl><br>>> > >>
>> Content-Type: text/plain; charset="iso-8859-1"<br>>> >
>> >><br>>> > >> >><br>>> >
>> >> Hello,<br>>> > >> >><br>>> >
>> >> My name is Cristian. I am having some troubles to perform
some<br>>> skewt<br>>> > >> >> graphics using WRF
output.<br>>> > >> >> I have created a binary and .ctl
archives using ARWpost.<br>>> > >> >> The file is OK,
and no troubles are expected.<br>>> > >> >><br>>>
> >> >> The problem starts when I want to create a skewt
graphic.<br>>> > >> >> The temperature and dewpoint
profiles are OK, but no information<br>>> > about<br>>> >
>> >> wind is shown. Neither the Hodograpf, nor the wind speed
and<br>>> > direction<br>>> > >> >> profiles
are shown.<br>>> > >> >><br>>> > >>
>> No message of an error apears, and i have no idea what
is<br>>> happening.<br>>> > >> >> Even the undef
value is th same as the .ctl archive.<br>>> > >>
>><br>>> > >> >> I am using Grads version 1.a8,
in UBUNTU 9.04.<br>>> > >> >><br>>> > >>
>> Here is the script. I only changed parameters at the top of
the<br>>> > >> archive.<br>>> > >>
>><br>>> > >> >> I really appreciate your
help.<br>>> > >> >><br>>> > >> >>
Thanks<br>>> > >> >><br>>> > >>
>><br>>> > >> >><br>>> >
>><br>>> ><br>>>
-----------------------------------------------------------------------------------------------------------------------------------------<br>>>
> >> >> * Script para construir radiosondeos<br>>> >
>> >><br>>> > >> >> 'reinit'<br>>>
> >> >> 'run
/home/cristian/Escritorio/Experimentos/lib/jaecol'<br>>> > >>
>> 'set display color white'<br>>> > >> >>
'reset'<br>>> > >> >><br>>> > >>
>> * Defino el archivo a utilizar<br>>> > >>
>><br>>> > >> >> 'open
20081004_YSU.ctl'<br>>> > >> >> *path=''<br>>>
> >> >><br>>> > >> >> * Latitud Longitud
y tiempo<br>>> > >> >><br>>> > >>
>> lat=-46<br>>> > >> >> lon=-68<br>>> >
>> >> time=12z04Oct2008<br>>> > >>
>><br>>> > >> >><br>>> > >>
>> 'set time 'time<br>>> > >> >> 'set lon
'lon<br>>> > >> >> 'set lat 'lat<br>>> >
>> >><br>>> > >> >> *Entre que niveles hace
el perfil hacerlo entre 1000 y 600<br>>> > >>
>><br>>> > >> >> levmax=925<br>>> >
>> >> levmin=500<br>>> > >> >><br>>>
> >> >> 'set lev 'levmax ' ' levmin<br>>> >
>> >><br>>> > >> >> *calculo Td a partir de
HR, P y T<br>>> > >> >> *'define
et=(-2937.4/(t))-(4.9283*log10(t))+22.5518'<br>>> > >>
>> *'define es=pow(10,et)*10'<br>>> > >> >>
*'define e=es*0.01*hr'<br>>> > >> >> *'define
q=0.622*e/(lev-e)'<br>>> > >> >> *'define
aux1=q*1000*(1.+0.81*q)'<br>>> > >> >> *'define
aux2=0.2854*(1.-0.28*q)'<br>>> > >> >> *'define
aux3=(1000/lev)'<br>>> > >> >> *'define
td=243.5/((17.67/(log(e/6.11)))-1)'<br>>> > >> >>
*'define tc=t-273'<br>>> > >> >><br>>> >
>> >> 'set ylopts 1 5 0.18'<br>>> > >>
>><br>>> > >> >> plotskew(tc,td,u,v)<br>>>
> >> >><br>>> > >> >> 'q
dims'<br>>> > >> >> line1=sublin(result,4)<br>>>
> >> >> line2=sublin(result,5)<br>>> > >>
>> itime1=subwrd(line1,6)<br>>> > >> >>
itit=substr(itime1,1,6)<br>>> > >> >>
itime1=subwrd(line2,6)<br>>> > >> >>
itime=substr(itime1,1,12)<br>>> > >> >>
d=substr(itime1,4,2)<br>>> > >> >>
m=substr(itime1,6,3)<br>>> > >> >>
y=substr(itime1,9,4)<br>>> > >> >>
h=substr(itime1,1,2)<br>>> > >> >>
lati=-1*lat<br>>> > >> >> loni=-1*lon<br>>> >
>> >><br>>> > >> >> 'draw title Perfil
vertical de T y Td en 'lati'S y 'loni'W\'itime<br>>> > >>
>> 'printim 'd'.'h'sondeo.png png x600 y800 white'<br>>> >
>> >><br>>> > >> >> function
plotskew(sndtemp,snddewp,u,v)<br>>> > >> >><br>>>
> >> >> *El valor Undef que se coloca ac? debe coincidir con
el valor<br>>> undef<br>>> > del<br>>> > >>
>> CTL, sin? hay problemas<br>>> > >> >> *en la
construcci?n de la hod?grafa y del perfil vertical de<br>>>
viento.<br>>> > >> >> undef=1.e30<br>>> >
>> >><br>>> > >><br>>> ><br>>>
*************************************************************************<br>>>
> >> >> *<br>>> > >> >> * GrADS Script
to Plot a SkewT/LogP Diagram<br>>> > >> >>
*<br>>> > >> >> * Bob Hart<br>>> > >>
>> * Penn State University / Dept of Meteorology<br>>> >
>> >> * Last Update: January 23, 2001<br>>> >
>> >> *<br>>> > >> >> * Recent
Changes:<br>>> > >> >> *<br>>> > >>
>> * 01/23/01 - Fixed a small bug in the theta-e
calculation.<br>>> > >> >> *
Errors averaged 0.5-3K. Thank you George
Bryan.<br>>> > >> >> *<br>>> > >>
>> * 11/10/99 - Change in calculation method for CAPE/CIN.
Trapezoid<br>>> > >> >> *
integration method is now used. Speeds up
execution<br>>> > >> >> *
by 25%, and increases accuracy by 5-10%.<br>>> >
>> >> *<br>>> > >> >> * 10/18/99 - Minor
glitch fixed that occasionally caused crash.<br>>> > >>
>> *<br>>> > >> >> * 8/26/99 - Datasets
with missing data can now be used.<br>>> > >> >>
*<br>>> > >> >> * Features:<br>>> > >>
>> * - All features of standard skewt/logp plot<br>>>
> >> >> * - RH sounding<br>>> > >>
>> * - LCL location<br>>> > >> >> *
- Parcel trajectory for both sfc based convection and elevated<br>>>
> from<br>>> > >> >> * most unstable
level (highest theta-e level reported)<br>>> > >> >> *
- Stability indices and precipitable water calculations<br>>>
> >> >> * - CAPE & CIN Calculations<br>>>
> >> >> * - Wind Profile<br>>> > >>
>> * - Hodograph / Hodograph scaling<br>>> > >>
>> * - Helicity and SR Helicity Calculations and
Display<br>>> > >> >> * - Color aspects of
output<br>>> > >> >> * - Line Thickness, style
aspects of output<br>>> > >> >> * - Can be run
in either PORTRAIT or LANDSCAPE mode.<br>>> > >> >>
*<br>>> > >> >> * There are numerous tunable parameters
below to change the<br>>> structure<br>>> > >> >>
* and output for the diagram.<br>>> > >> >>
*<br>>> > >> >> * Function Arguments:<br>>> >
>> >> * sndtemp - temperature data (Celsius) as a
function of<br>>> pressure<br>>> > >> >> *
snddewp - dewpoint data (Celsius) as a function of
pressure<br>>> > >> >> * sndspd -
wind speed data (knots) as a function of pressure<br>>> > >>
>> * snddir - wind direction data as a function of
pressure<br>>> > >> >> *<br>>> > >>
>> * Use '-1' for any of the above 4 arguments to indicate that
you<br>>> > >> >> * are not passing that variable.
The appropriate options will<br>>> > >> >> * be
ignored based on your specifying '-1' for that variable.<br>>> >
>> >> *<br>>> > >> >> * NOTE: Make
sure to set the vertical range of the plot before<br>>> > >>
running.<br>>> > >> >> *
I.e., "SET LEV 1050 150", for example. This does
not<br>>> have<br>>> > to<br>>> > >> >>
* be limited to the pressure range of your
data.<br>>> > >> >> *<br>>> > >>
>> * Labelling: Pressure/Height is labelled along left
side.<br>>> > Temperature<br>>> > >>
is<br>>> > >> >> *
labelled along bottom. Mixing ratio is labelled<br>>>
along<br>>> > >> right<br>>> > >> >> *
side/top.<br>>> >
>> >> *<br>>> > >> >> *<br>>> >
>> >> * PROBLEMS: First check out the web page for the
script (which<br>>> also<br>>> > >> >> *
has a link to a FAQ with answers to many
common<br>>> > questions<br>>> > >> >> *
about using the
script):<br>>> > >> >> *
<a href="http://www.ems.psu.edu/%7Ehart/skew.html" target="_blank">http://www.ems.psu.edu/~hart/skew.html</a><<a href="http://www.ems.psu.edu/%7Ehart/skew.html" target="_blank">http://www.ems.psu.edu/%7Ehart/skew.html</a>><br>
>>
<<a href="http://www.ems.psu.edu/%7Ehart/skew.html" target="_blank">http://www.ems.psu.edu/%7Ehart/skew.html</a>><br>>>
> <<a href="http://www.ems.psu.edu/%7Ehart/skew.html" target="_blank">http://www.ems.psu.edu/%7Ehart/skew.html</a>><br>>>
> >> <<a href="http://www.ems.psu.edu/%7Ehart/skew.html" target="_blank">http://www.ems.psu.edu/%7Ehart/skew.html</a>><br>>><br>>>
> >><br>>> > >> >> *<br>>> > >>
>> * Please send any further problems, comments, or suggestions
to<br>>> > >> >> * <<a href="mailto:hart@ems.psu.edu" target="_blank">hart@ems.psu.edu</a>><br>>> > >> >>
*<br>>> > >> >> * ACKNOWLEDGMENTS: Thanks go to
the innumerable users who have<br>>> > helped<br>>> >
>> >> * fine tune the script from the horrible mess from which
it began.<br>>> > >> >> * In particular, thanks go out
to Steve Lord (NCEP), Mike Fiorino<br>>> > >>
(ECMWF),<br>>> > >> >> * George Bryan (PSU), Davide
Sacchetti (CMIRL), and Enrico<br>>> Minguzzi<br>>> > >>
>> (CMIRL).<br>>> > >> >> *<br>>> >
>> >><br>>> > >><br>>> ><br>>>
**************************************************************************<br>>>
> >> >> * !!!!!
BEGINNING OF USER-SPECIFIED OPTIONS !!!!!!<br>>> > >>
>><br>>> > >><br>>> ><br>>>
**************************************************************************<br>>>
> >> >> *<br>>> > >> >> *
--------------------- Initialization options<br>>> >
----------------------<br>>> > >> >>
*<br>>> > >> >> * ClrScrn = Whether to clear the screen
before drawing diagram<br>>> > >> >> *
[1 = yes, 0 = no]<br>>> > >>
>><br>>> > >> >> ClrScrn = 1<br>>> >
>> >><br>>> > >> >> *<br>>> >
>> >> * ------------------- Define Skew-T Diagram<br>>>
> >> Shape/Slope-----------------<br>>> > >>
>> *<br>>> > >> >> * (P1,T1) = Pres, Temp of some
point on left-most side<br>>> > >> >> * (P2,T2) = Pres,
Temp of some point on right-most side<br>>> > >> >> *
(P3,T3) = Pres, Temp of some point in diagram which is mid-point<br>>>
> >> >> * in the
horizontal between 1 and 2.<br>>> > >> >> *<br>>>
> >> >> * P1, P2, P3 are in mb ; T1, T2, T3 are in
Celsius<br>>> > >> >> *<br>>> > >>
>> * These define the SLOPE and WIDTH of the diagram as you see
it<br>>> but<br>>> > DO<br>>> > >>
NOT<br>>> > >> >> * DEFINE THE HEIGHT of the diagram as
you see it. In other words,<br>>> > >> >> * 1 and
2 do NOT necessarily need to be at the bottom of the<br>>>
diagram<br>>> > >> and<br>>> > >> >> * 3
does NOT necessarily need to be at the top. THE VERTICAL<br>>>
> PRESSURE<br>>> > >> >> * RANGE OF THE SKEWT AS YOU
SEE IT IS DETERMINED BY YOUR 'SET Z<br>>> ...'<br>>> >
>> >> * COMMAND OR THE 'SET LEV ...' COMMAND BEFORE RUNNING THIS
SCRIPT.<br>>> > >> >> *<br>>> > >>
>> * _______________________<br>>> > >>
>> * |
|<br>>> > >> >> * |
|<br>>> > >> >> * |
3 |<br>>> > >>
>> * |
|<br>>> > >> >> * |
|<br>>> > >> >> * |
|<br>>> > >>
>> * |
|<br>>> > >> >> * |
|<br>>> > >> >> * |
|<br>>> > >>
>> * |
|<br>>> > >> >> * |
|<br>>> > >> >> * |1
2|<br>>> > >>
>> * |
|<br>>> > >> >> *
|_______________________|<br>>> > >> >> *<br>>>
> >> >> *<br>>> > >> >> * A good set of
defining points are given below. Feel free<br>>> > >>
>> * to experiment with variations.<br>>> > >>
>><br>>> > >> >><br>>> > >>
>> *P1 = 1000<br>>> > >> >> *T1 = -40<br>>>
> >> >><br>>> > >> >> *P2 =
1000<br>>> > >> >> *T2 = 40<br>>> > >>
>><br>>> > >> >> *P3 = 200<br>>> >
>> >> *T3 = 0<br>>> > >> >><br>>>
> >> >> * Another good set of defining points suggested by
Juan Ruiz<br>>> > (Emagrama)<br>>> > >> >> *
are:<br>>> > >> >> *<br>>> > >> >>
P1 = 1000<br>>> > >> >> T1 =
-80<br>>> > >> >> *<br>>> > >> >>
P2 = 1000<br>>> > >> >> T2 = 40<br>>>
> >> >> *<br>>> > >> >> P3 =
500<br>>> > >> >> T3 = -20<br>>> >
>> >><br>>> > >> >> * -------------------
Contour Intervals / Levels<br>>> > >> >>
--------------------------<br>>> > >> >> *<br>>>
> >> >> * All variables below are contour intervals/levels
for diagram<br>>> > >> >> *<br>>> > >>
>> * Thetaint = interval for potential temperature lines<br>>>
> >> >> * Thetwint = interval for moist pseudo
adiabats<br>>> > >> >> * tempint = interval for
temperature lines<br>>> > >> >> * wsclevs =
contour LEVELS for mixing ratio lines<br>>> > >> >>
*<br>>> > >> >> *<br>>> > >> >>
thetaint= 20<br>>> > >> >> thetwint= 10<br>>>
> >> >> tempint = 10<br>>> > >> >>
wsclevs = ".1 .5 1 3 6 10 15 20 25 30"<br>>> > >> >>
*<br>>> > >> >> *<br>>> > >> >> *
------------------------ Output Options<br>>> > >>
--------------------------------<br>>> > >> >>
*<br>>> > >> >> * All variables below are logical ..
1=yes, 0=no, unless otherwise<br>>> > >> >> *
specified.<br>>> > >> >> *<br>>> > >>
>> * DrawBarb = Draw wind barbs along right side of plot<br>>>
> >> >> * DrawThet = Draw dry adiabats<br>>> >
>> >> * DrawThtw = Draw moist pseudo-adiabats<br>>> >
>> >> * DrawTemp = Draw temperature lines<br>>> >
>> >> * DrawMix = Draw mixing ratio lines<br>>> >
>> >> * DrawTSnd = Draw temperature sounding<br>>> >
>> >> * DrawDSnd = Draw dewpoint sounding<br>>> >
>> >> * DrawRH = Draw relative humidity
sounding<br>>> > >> >> * DrawPrcl = Draw parcel path
from surface upward<br>>> > >> >> * DrawPMax = Draw
parcel path from most unstable level upward<br>>> > >>
>> * DrawIndx = Display stability indices & CAPE<br>>> >
>> >> * DrawHeli = Calculate and display absolute and
storm-relative<br>>> > helicity<br>>> > >> >>
* DrawHodo = Draw hodograph<br>>> > >> >> * DrawPLev =
Draw Pressure Levels<br>>> > >> >> * DrawZLev = Draw
height levels and lines<br>>> > >> >> *
0 = no lines<br>>> > >> >> *
1 = above ground level
(AGL)<br>>> > >> >> *
2 = above sea level (ASL)<br>>> > >> >> *
DrawZSTD = Draw Height levels using standard atm lapse rate<br>>> >
>> >> * LblAxes = Label the x,y axes (temperature,
pressure,mixing<br>>> ratio)<br>>> > >> >>
*<br>>> > >> >> * ThtwStop = Pressure level at which to
stop drawing Theta-w lines<br>>> > >> >> * MixStop
= Pressure level at which to stop drawing Mixratio<br>>>
lines<br>>> > >> >><br>>> > >> >>
DrawBarb= 1<br>>> > >> >> DrawThet= 1<br>>> >
>> >> DrawThtw= 1<br>>> > >> >> DrawTemp=
1<br>>> > >> >> DrawMix = 1<br>>> > >>
>> DrawTSnd= 1<br>>> > >> >> DrawDSnd=
1<br>>> > >> >> DrawRH = 0<br>>> >
>> >> DrawPrcl= 1<br>>> > >> >> DrawPMax=
1<br>>> > >> >> DrawIndx= 1<br>>> > >>
>> DrawHeli= 0<br>>> > >> >> DrawHodo=
1<br>>> > >> >> DrawPLev= 1<br>>> > >>
>> DrawZLev= 0<br>>> > >> >> DrawZSTD=
0<br>>> > >> >> LblAxes = 1<br>>> > >>
>><br>>> > >> >> ThtwStop = 200<br>>> >
>> >> MixStop = 400<br>>> > >>
>><br>>> > >> >><br>>> > >>
>> *<br>>> > >> >> * -----------------
Sounding Geography options<br>>> > >>
------------------------<br>>> > >> >> *<br>>>
> >> >> * SfcElev = Elevation above sea-level (meters) of
lowest level<br>>> > reported<br>>> > >> >> *
in sounding. Used only if DrawZLev
= 2<br>>> > >> >><br>>> > >> >>
SfcElev = 0<br>>> > >> >><br>>> > >>
>><br>>> > >> >> *<br>>> > >>
>> * ------------------ Thermodynamic Index Options<br>>>
--------------------<br>>> > >> >> *<br>>> >
>> >> * All variables here are in inches. Use -1 for the
default<br>>> values.<br>>> > >> >> *<br>>>
> >> >> * Text1XC = X-location of midpoint of K,TT,PW
output box<br>>> > >> >> * Text1YC = Y-location
of midpoint of K,TT,PW output box<br>>> > >> >> *
Text2XC = X-Location of midpoint of surface indices output
box<br>>> > >> >> * Text2YC = Y-location of
midpoint of surface indices output box<br>>> > >> >> *
Text3XC = X-Location of midpoint of most unstable
level-based<br>>> > >> indices<br>>> > >>
>> * output box<br>>>
> >> >> * Text3YC = Y-location of midpoint of most
unstable level-based<br>>> > >> indices<br>>> >
>> >> * output
box<br>>> > >> >><br>>> > >> >>
Text1XC = -1<br>>> > >> >> Text1YC = -1<br>>>
> >> >> Text2XC = -1<br>>> > >> >>
Text2YC = -1<br>>> > >> >> Text3XC = -1<br>>>
> >> >> Text3YC = -1<br>>> > >>
>><br>>> > >> >> *<br>>> > >>
>> * ----------------- Wind Barb Profile Options<br>>> >
>> ----------------------------<br>>> > >> >>
*<br>>> > >> >> * All variables here are in units of
inches, unless otherwise<br>>> > specified<br>>> >
>> >> *<br>>> > >> >> * barbint =
Interval for plotting barbs (in units of levels)<br>>> > >>
>> * poleloc = X-Location of profile. Choose -1 for the
default.<br>>> > >> >> * polelen = Length of
wind-barb pole<br>>> > >> >> * Len05 =
Length of each 5-knot barb<br>>> > >> >> * Len10
= Length of each 10-knot barb<br>>> > >> >> *
Len50 = Length of each 50-knot flag<br>>> > >>
>> * Wid50 = Width of base of 50-knot flag<br>>>
> >> >> * Spac50 = Spacing between 50-knot flag
and next ror occurred on<br>>> > >> >>
libarb/flag<br>>> > >> >> * Spac10 =
Spacing between 10-knot flag and next flag<br>>> > >>
>> * Spac05 = Spacing between 5-knot flag and next
flag<br>>> > >> >> * Flagbase= Draw flagbase
(filled circle) for each windbarb<br>>> [1=yes,<br>>> >
0<br>>> > >> >> =no]<br>>> > >> >>
* Fill50 = Solid-fill 50-knot flag [1=yes, 0=no]<br>>>
> >> >> * barbline= Draw a vertical line connecting all
the wind barbs<br>>> > [1=yes,<br>>> > >> >>
0=no]<br>>> > >> >> *<br>>> > >>
>> barbint = 1<br>>> > >> >> poleloc =
-1<br>>> > >> >> polelen = 0.35<br>>> >
>> >> len05 = 0.07<br>>> > >> >>
len10 = 0.15<br>>> > >> >> len50 =
0.15<br>>> > >> >> wid50 = 0.06<br>>> >
>> >> spac50 = 0.07<br>>> > >> >>
spac10 = 0.05<br>>> > >> >> spac05 =
0.05<br>>> > >> >> Fill50 = 1<br>>> >
>> >> flagbase= 1<br>>> > >> >> barbline=
1<br>>> > >> >><br>>> > >> >>
*<br>>> > >> >> *<br>>> > >> >>
*---------------- Hodograph Options<br>>> > >>
-------------------------------------<br>>> > >> >>
*<br>>> > >> >> * All variables here are in units of
inches, unless otherwise<br>>> > specified<br>>> >
>> >> *<br>>> > >> >> * HodXcent=
x-location of hodograph center. Use -1 for default<br>>> >
>> location.<br>>> > >> >> * HodYcent= y-location
of hodograph center. Use -1 for default<br>>> > >>
location.<br>>> > >> >> * HodSize = Size of hodograph
in inches<br>>> > >> >> * NumRing = Number of rings to
place in hodograph (must be at<br>>> least<br>>> >
1)<br>>> > >> >> * HodRing = Wind speed increment of
each hodograph ring<br>>> > >> >> * HodoDep = Depth
(above lowest level in mb) of end of hodograph<br>>> >
trace<br>>> > >> >> * TickInt = Interval (in kts) at
which tick marks are drawn along<br>>> the<br>>> > >>
axes<br>>> > >> >> *
Use 0 for no tick marks.<br>>> > >> >> * TickSize= Size
of tick mark in inches<br>>> > >> >> * Text4XC =
X-location of midpoint of hodograph text output. Use<br>>>
-1<br>>> > for<br>>> > >> >>
default.<br>>> > >> >> * Text4YC = Y-location of
midpoint of hodograph text output. Use<br>>> -1<br>>> >
for<br>>> > >> >> default.<br>>> > >>
>><br>>> > >> >> HodXcent= 6<br>>> >
>> >> HodYcent= 9<br>>> > >> >> HodSize =
2<br>>> > >> >> NumRing = 3<br>>> > >>
>> HodRing = 10<br>>> > >> >> HodoDep =
500<br>>> > >> >> TickInt = 5<br>>> > >>
>> TickSize= 0.1<br>>> > >> >> Text4XC =
-1<br>>> > >> >> Text4YC = -1<br>>> > >>
>><br>>> > >> >> *--------------- Helicity
Options<br>>> > >>
---------------------------------------<br>>> > >> >>
*<br>>> > >> >> * MeanVTop = Top pressure level (mb) of
mean-wind calculation<br>>> > >> >> * MeanVBot = Bottom
pressure level (mb) of mean-wind calculation<br>>> > >>
>> * HelicDep = Depth in mb (above ground) of helicity
integration<br>>> > >> >> * StormMot = Type of storm
motion estimation scheme. Use<br>>> following:<br>>> >
>> >> * 0 = No
departure from mean wind.<br>>> > >> >> *
1 = Davies-Jones (1990) approach<br>>> >
>> >> * FillArrw = Whether to fill the arrowhead of the storm
motion<br>>> vector<br>>> > >> >> *
[1 = yes, 0 = no]<br>>> > >>
>><br>>> > >> >> MeanVTop= 300<br>>> >
>> >> MeanVBot= 850<br>>> > >> >> HelicDep=
300<br>>> > >> >> StormMot= 0<br>>> > >>
>> FillArrw= 1<br>>> > >> >><br>>> >
>> >> *<br>>> > >> >> *----------------
Color Options<br>>> > >>
------------------------------------------<br>>> > >>
>> *<br>>> > >> >> * ThetCol = Color of dry
adiabats<br>>> > >> >> * TempCol = Color of temperature
lines<br>>> > >> >> * MixCol = Color of mixing
ratio lines<br>>> > >> >> * ThtwCol = Color of moist
adiabats<br>>> > >> >> * TSndCol = Color of Temperature
Sounding<br>>> > >> >> * DSndCol = Color of Dewpoint
Sounding<br>>> > >> >> * RHCol = Color of RH
Sounding<br>>> > >> >> * PrclCol = Color of parcel
trace<br>>> > >> >> * BarbCol = Color of wind barbs
(choose -1 for color according to<br>>> > >>
speed)<br>>> > >> >> * HodoCol = Color of hodograph
trace<br>>> > >> >><br>>> > >> >>
ThetCol = 23<br>>> > >> >> TempCol = 79<br>>>
> >> >> MixCol = 38<br>>> > >> >>
ThtwCol = 39<br>>> > >> >> TSndCol = 29<br>>>
> >> >> DSndCol = 49<br>>> > >> >> RHCol
= 3<br>>> > >> >> PrclCol = 1<br>>> >
>> >> BarbCol = 1<br>>> > >> >> HodoCol =
2<br>>> > >> >><br>>> > >> >>
*<br>>> > >> >> *-------------------- Line Style
Options<br>>> > >> >>
------------------------------------<br>>> > >> >>
*<br>>> > >> >> * GrADS Styles: 1=solid;2=long
dash;3=short dash;4=long,short<br>>> dashed;<br>>> > >>
>> * 5=dotted;6=dot
dash;7=dot dot dash<br>>> > >> >> *<br>>> >
>> >> * ThetLine = Line Style of dry adiabats<br>>> >
>> >> * TempLine = Line Style of temperature lines<br>>>
> >> >> * MixLine = Line Style of mixing ratio
lines<br>>> > >> >> * ThtwLine = Line Style of moist
adiabats<br>>> > >> >> * TSndLine = Line Style of
Temperature Sounding<br>>> > >> >> * DSndLine = Line
Style of Dewpoint Sounding<br>>> > >> >> * RHLine
= Line Style of RH sounding<br>>> > >> >> *
PrclLine = Line Style of parcel trace<br>>> > >> >> *
HodoLine = Line Style of hodograph trace<br>>> > >> >>
*<br>>> > >> >><br>>> > >> >>
ThetLine = 1<br>>> > >> >> TempLine = 1<br>>>
> >> >> MixLine = 5<br>>> > >> >>
ThtwLine = 2<br>>> > >> >> TSndLine = 1<br>>>
> >> >> DSndLine = 2<br>>> > >> >>
RHLine = 1<br>>> > >> >> PrclLine =
3<br>>> > >> >> HodoLine = 1<br>>> > >>
>><br>>> > >> >> *<br>>> > >>
>> *------------------- Line Thickness<br>>> > >>
>> Options---------------------------------<br>>> > >>
>> * GrADS Line Thickness: increases with increasing
number.<br>>> Influences<br>>> > >> >> *
hardcopy output more strongly than screen<br>>> > >>
output.<br>>> > >> >> *<br>>> > >>
>> *<br>>> > >> >> * ThetThk = Line Thickness of
dry adiabats<br>>> > >> >> * TempThk = Line Thickness
of temperature lines<br>>> > >> >> * MixThk =
Line Thickness of mixing ratio lines<br>>> > >> >> *
ThtwThk = Line Thickness of moist adiabats<br>>> > >>
>> * TSndThk = Line Thickness of temperature sounding<br>>> >
>> >> * DSndThk = Line thickness of dewpoint
sounding<br>>> > >> >> * RHThk = Line thickness
of RH sounding<br>>> > >> >> * PrclThk = Line thickness
of parcel trace<br>>> > >> >> * HodoThk = Line
thickness of hodograph trace<br>>> > >> >> * BarbThk =
Line thickness of wind barbs<br>>> > >> >><br>>>
> >> >> ThetThk = 3<br>>> > >> >>
TempThk = 1<br>>> > >> >> MixThk = 7<br>>>
> >> >> ThtwThk = 3<br>>> > >> >>
TSndThk = 12<br>>> > >> >> DSndThk = 12<br>>>
> >> >> RHThk = 8<br>>> > >> >>
PrclThk = 6<br>>> > >> >> HodoThk = 6<br>>> >
>> >> BarbThk = 2<br>>> > >> >><br>>>
> >> >> *<br>>> > >> >>
*------------------- Data Point Marker Options<br>>> > >>
>> -----------------------------<br>>> > >> >> *
GrADS Marker Types: 0 = none ; 1 = cross ; 2 = open circle ;<br>>>
> >> >> *
3 = closed circle ; 4 = open square ; 5 =<br>>>
> closed<br>>> > >> >> square<br>>> >
>> >> *
6 = X ; 7 = diamond ; 8 = triangle ; 9 =<br>>>
none<br>>> > >> >> *
10 = open circle with vertical line ; 11
=<br>>> open<br>>> > >> oval<br>>> > >>
>> *<br>>> > >> >> * TSndMrk = Mark type of data
point marker for temperature<br>>> sounding<br>>> > >>
>> * DSndMrk = Mark type of data point marker for dewpoint
sounding<br>>> > >> >> * RHMrk = Mark type of
data point marker for relative humidity<br>>> > >>
sounding<br>>> > >> >> * MrkSize = Mark size (inches)
of each data marker<br>>> > >> >><br>>> >
>> >> TSndMrk = 3<br>>> > >> >> DSndMrk =
3<br>>> > >> >> RHMrk = 0<br>>> >
>> >> MrkSize = 0.1<br>>> > >>
>><br>>> > >> >><br>>> > >>
>> * !!!!! YOU SHOULD NOT NEED TO CHANGE ANYTHING BELOW HERE
!!!!!<br>>> > >> >><br>>> > >>
>><br>>> > >><br>>> ><br>>>
****************************************************************************<br>>>
> >> >><br>>> > >> >>
*-------------------------------------------<br>>> > >>
>> * grab user-specified environment dimensions<br>>> >
>> >> *-------------------------------------------<br>>>
> >> >><br>>> > >> >> "q
dims"<br>>> > >> >> rec=sublin(result,2)<br>>>
> >> >> _xtype=subwrd(rec,3)<br>>> > >>
>> _xval=subwrd(rec,9)<br>>> > >> >>
rec=sublin(result,3)<br>>> > >> >>
_yval=subwrd(rec,9)<br>>> > >> >>
_ytype=subwrd(rec,3)<br>>> > >> >>
rec=sublin(result,4)<br>>> > >> >>
_ptype=subwrd(rec,3)<br>>> > >> >>
_pmax=subwrd(rec,6)<br>>> > >> >>
_pmin=subwrd(rec,8)<br>>> > >> >>
_zmin=subwrd(rec,11)<br>>> > >> >>
_zmax=subwrd(rec,13)<br>>> > >> >>
rec=sublin(result,5)<br>>> > >> >>
_ttype=subwrd(rec,3)<br>>> > >> >>
_tval=subwrd(rec,9)<br>>> > >> >><br>>> >
>> >> "q file"<br>>> > >> >>
rec=sublin(result,5)<br>>> > >> >>
_zmaxfile=subwrd(rec,9)<br>>> > >> >><br>>> >
>> >>
*-------------------------------------------------------------<br>>>
> >> >> * Check to ensure that dimensions are valid.
Warn & exit if not.<br>>> > >> >>
*--------------------------------------------------------------<br>>>
> >> >><br>>> > >> >>
dimrc=0<br>>> > >> >> If (_xtype !=
"fixed")<br>>> > >> >> say "X-Dims Error:
Not fixed. Use 'set lon' or 'set x' to<br>>>
specify<br>>> > a<br>>> > >> >>
value."<br>>> > >> >> dimrc=-1<br>>> >
>> >> Endif<br>>> > >> >><br>>> >
>> >> If (_ytype != "fixed")<br>>> > >> >>
say "Y-Dims Error: Not fixed. Use 'set lat' or 'set y'
to<br>>> specify<br>>> > a<br>>> > >> >>
value"<br>>> > >> >> dimrc=-1<br>>> >
>> >> Endif<br>>> > >> >><br>>> >
>> >> If (_ptype != "varying")<br>>> > >>
>> say "Z-Dims Error: Not varying. Use 'set lev' or
'set z' to<br>>> > specify<br>>> > >> a<br>>>
> >> >> range."<br>>> > >> >>
dimrc=-1<br>>> > >> >> Endif<br>>> > >>
>><br>>> > >> >> If (_ttype !=
"fixed")<br>>> > >> >> say "Time Error:
Not fixed. Use 'set time' or 'set t' to<br>>> >
specify<br>>> > >> a<br>>> > >> >>
value"<br>>> > >> >> dimrc=-1<br>>> >
>> >> Endif<br>>> > >> >><br>>> >
>> >><br>>> > >> >> If (dimrc <
0)<br>>> > >> >> Return(-1)<br>>> >
>> >> Endif<br>>> > >> >><br>>> >
>> >><br>>> > >> >> *<br>>> >
>> >> * A few global variables used in units
conversion<br>>> > >> >> *<br>>> > >>
>><br>>> > >> >> _pi=3.14159265<br>>> >
>> >> _dtr=_pi/180<br>>> > >> >>
_rtd=1/_dtr<br>>> > >> >> _ktm=0.514444<br>>>
> >> >> _mtk=1/_ktm<br>>> > >>
>><br>>> > >> >> * A few global constants used in
thermo calcs<br>>> > >> >><br>>> > >>
>> _C0=0.99999683<br>>> > >> >>
_C1=-0.90826951/100<br>>> > >> >> _C2=
0.78736169/10000<br>>> > >> >>
_C3=-0.61117958/1000000<br>>> > >> >> _C4=
0.43884187/pow(10,8)<br>>> > >> >>
_C5=-0.29883885/pow(10,10)<br>>> > >> >> _C6=
0.21874425/pow(10,12)<br>>> > >> >>
_C7=-0.17892321/pow(10,14)<br>>> > >> >> _C8=
0.11112018/pow(10,16)<br>>> > >> >>
_C9=-0.30994571/pow(10,19)<br>>> > >> >><br>>>
> >> >> *Calculo la direcci?n e intensidad del viento usando
las funciones<br>>> > >> GetWdir<br>>> > >>
>> y GetWsp<br>>> > >> >><br>>> >
>> >><br>>> > >> >><br>>> >
>> >><br>>> > >> >><br>>> >
>> >><br>>> > >> >> * A pressure array of
power calculations which should be performed<br>>> > >>
>> * only once to reduce execution time.<br>>> > >>
>><br>>> > >> >> zz=1100<br>>> >
>> >> while (zz > 10)<br>>> > >> >>
subscr=zz/10<br>>> > >> >>
_powpres.subscr=pow(zz,0.286)<br>>> > >> >>
zz=zz-10<br>>> > >> >>
endwhile<br>>> > >> >><br>>> > >>
>> *<br>>> > >> >> * Turn off options not
available due to user data limitations<br>>> > >> >>
*<br>>> > >> >><br>>> > >> >> If
(ClrScrn = 1)<br>>> > >> >> "clear"<br>>>
> >> >> Endif<br>>> > >> >><br>>>
> >> >> If (sndspd = -1 | snddir = -1)<br>>> >
>> >> DrawBarb = 0<br>>> > >> >>
DrawHodo = 0<br>>> > >> >> DrawHeli =
0<br>>> > >> >> Endif<br>>> > >>
>><br>>> > >> >> If (snddewp = -1)<br>>>
> >> >> DrawDSnd = 0<br>>> > >> >>
DrawRH = 0<br>>> > >> >> DrawPrcl =
0<br>>> > >> >> DrawPMax = 0<br>>> >
>> >> DrawIndx = 0<br>>> > >> >>
Endif<br>>> > >> >><br>>> > >> >>
If (sndtemp = -1)<br>>> > >> >> DrawTSnd =
0<br>>> > >> >> DrawRH = 0<br>>>
> >> >> DrawPrcl = 0<br>>> > >> >>
DrawPMax = 0<br>>> > >> >> DrawIndx =
0<br>>> > >> >> DrawZLev = 0<br>>> >
>> >> Endif<br>>> > >> >><br>>> >
>> >> If (NumRing < 1)<br>>> > >> >>
DrawHodo = 0<br>>> > >> >> Endif<br>>> >
>> >><br>>> > >> >> "q gxinfo"<br>>>
> >> >> rec=sublin(result,2)<br>>> > >>
>> xsize=subwrd(rec,4)<br>>> > >> >><br>>>
> >> >> If (xsize = 11)<br>>> > >> >>
PageType = "Landscape"<br>>> > >> >>
Else<br>>> > >> >> PageType =
"Portrait"<br>>> > >> >> Endif<br>>> >
>> >><br>>> > >> >>
*------------------------------------------------------<br>>> >
>> >> * calculate constants determining slope/shape of
diagram<br>>> > >> >> * based on temp/pressure values
given by user<br>>> > >> >>
*-------------------------------------------------------<br>>> >
>> >><br>>> > >> >> "set x 1"<br>>>
> >> >> "set y 1"<br>>> > >> >> "set z
1"<br>>> > >> >> "set t 1"<br>>> > >>
>> _m1=(T1+T2-2*T3)/(2*log10(P2/P3))<br>>> > >>
>> _m2=(T2-T3-_m1*log10(P2/P3))/50<br>>> > >> >>
_m3=(T1-_m1*log10(P1))<br>>> > >> >><br>>> >
>> >> "set z "_zmin" "_zmax<br>>> > >> >>
"set zlog on"<br>>> > >> >> "set xlab off"<br>>>
> >> >><br>>> > >> >>
*-------------------------------------------------<br>>> > >>
>> * perform coordinate transformation to Skew-T/LogP<br>>> >
>> >>
*-------------------------------------------------<br>>> > >>
>><br>>> > >> >> "set gxout stat"<br>>>
> >> >> "set x "_xval<br>>> > >> >> "set
y "_yval<br>>> > >> >> "set t "_tval<br>>> >
>> >> "define
tempx=("sndtemp"-"_m1"*log10(lev)-"_m3")/"_m2<br>>> > >>
>> "define dewpx=("snddewp"-"_m1"*log10(lev)-"_m3")/"_m2<br>>>
> >> >><br>>> > >> >> If (PageType =
"Portrait")<br>>> > >> >> "set parea 0.7 7 0.75
10"<br>>> > >> >> Else<br>>> > >>
>> "set parea 0.7 6.5 0.5 8"<br>>> > >> >>
Endif<br>>> > >> >><br>>> > >> >>
"set axlim 0 100"<br>>> > >> >> "set lon 0
100"<br>>> > >> >> "set grid on 1 1"<br>>> >
>> >><br>>> > >> >> "set z "_zmin " "
_zmax<br>>> > >> >> "set lon 0 100"<br>>> >
>> >> "set clevs -900"<br>>> > >> >> "set
gxout contour"<br>>> > >> >><br>>> > >>
>> *-------------------------------------<br>>> > >>
>> * Draw pressure lines<br>>> > >> >>
*-------------------------------------<br>>> > >>
>><br>>> > >> >> If (DrawPLev = 0)<br>>>
> >> >> "set ylab off"<br>>> > >>
>> Else<br>>> > >> >> "set ylab
on"<br>>> > >> >> "set ylopts 1 5
0.18"<br>>> > >> >> "set xlopts 1 3
0.18"<br>>> > >> >> Endif<br>>> > >>
>><br>>> > >> >> "d lon"<br>>> >
>> >><br>>> > >> >>
*--------------------------------------<br>>> > >> >> *
Determine corners of skewt/logp frame<br>>> > >> >>
*--------------------------------------<br>>> > >>
>><br>>> > >> >> "q w2xy 100 "_pmin<br>>>
> >> >> rxloc=subwrd(result,3)<br>>> > >>
>> tyloc=subwrd(result,6)<br>>> > >> >> "q w2xy 0
"_pmax<br>>> > >> >> lxloc=subwrd(result,3)<br>>>
> >> >> byloc=subwrd(result,6)<br>>> > >>
>><br>>> > >> >> If (DrawPLev = 1 & LblAxes =
1)<br>>> > >> >> "set strsiz 0.15"<br>>>
> >> >> "set string 1 c 3 0"<br>>> > >>
>> If (PageType = "Portrait")<br>>> > >>
>> * "draw string 0.5 10.5 hPa."<br>>> > >> >>
Else<br>>> > >> >> "draw
string 0.5 8.35 hPa."<br>>> > >> >>
Endif<br>>> > >> >> Endif<br>>> > >>
>><br>>> > >> >>
*---------------------------------------------------<br>>> >
>> >> * Calculate & draw actual height lines using temp
data<br>>> > >> >>
*---------------------------------------------------<br>>> >
>> >><br>>> > >> >> If (DrawZLev >
0)<br>>> > >> >> say "Calculating observed
height levels from temp/pressure<br>>> data."<br>>> >
>> >> zz=1<br>>> > >> >> "set
gxout stat"<br>>> > >> >> "set x
"_xval<br>>> > >> >> "set y "_yval<br>>>
> >> >> "set t "_tval<br>>> > >>
>> count=0<br>>> > >> >> while (zz
< _zmax)<br>>> > >> >> "set z
"zz<br>>> > >> >>
pp.zz=subwrd(result,4)<br>>> > >> >>
lpp.zz=log(pp.zz)<br>>> > >> >>
"d "sndtemp<br>>> > >> >>
rec=sublin(result,8)<br>>> > >> >>
tt=subwrd(rec,4)<br>>> > >> >>
if (tt > -900)<br>>> > >> >>
tk=tt+273.15<br>>> > >> >>
count=count+1<br>>> > >> >>
zzm=zz-1<br>>> > >> >>
If (count = 1)<br>>> > >> >>
If (DrawZLev = 2)<br>>> >
>> >>
htlb="ASL"<br>>> > >> >>
height.zz=SfcElev<br>>> > >> >>
Else<br>>> > >>
>>
htlb="AGL"<br>>> > >> >>
height.zz=0<br>>> > >> >>
Endif<br>>> > >> >>
sfcz=height.zz<br>>> >
>> >> Else<br>>> > >>
>><br>>> > >> >><br>>> >
DZ=29.2857*(lpp.zzm-lpp.zz)*(lpp.zz*tk+lpp.zzm*tkold)/(lpp.zz+lpp.zzm)<br>>>
> >> >>
height.zz=height.zzm+DZ<br>>> > >> >>
highz=height.zz<br>>> > >>
>> Endif<br>>> > >>
>> else<br>>> > >> >>
height.zz = -9999<br>>> > >> >>
endif<br>>> > >> >>
tkold=tk<br>>> > >> >>
zz=zz+1<br>>> > >> >> endwhile<br>>>
> >> >><br>>> > >> >>
maxht=int(highz/1000)<br>>> > >> >> if
(int(sfcz/1000) = sfcz/1000)<br>>> > >> >>
minht=int(sfcz/1000)<br>>> > >> >>
else<br>>> > >> >>
minht=1+int(sfcz/1000)<br>>> > >> >>
endif<br>>> > >> >><br>>> > >> >>
ht=minht<br>>> > >> >> "set line 1 3
1"<br>>> > >> >> "set strsiz 0.10"<br>>>
> >> >> "set string 1 l 3 0"<br>>> > >>
>> while (ht <= maxht)<br>>> > >> >>
zz=1<br>>> > >> >>
while (height.zz/1000 <= ht)<br>>> > >> >>
zz=zz+1<br>>> > >> >>
endwhile<br>>> > >> >>
zzm=zz-1<br>>> > >> >>
PBelow=pp.zzm<br>>> > >> >>
PAbove=pp.zz<br>>> > >> >>
HBelow=height.zzm<br>>> > >> >>
HAbove=height.zz<br>>> > >> >>
DZ=HAbove-HBelow<br>>> > >> >>
DP=PAbove-PBelow<br>>> > >> >>
Del=ht*1000-HBelow<br>>> > >> >>
Est=PBelow+Del*DP/DZ<br>>> > >> >>
If (Est >= _pmin & Est <= _pmax)<br>>> > >>
>> "q w2xy 1 " Est<br>>> >
>> >>
yloc=subwrd(result,6)<br>>> > >> >>
"draw line " lxloc " " yloc " " rxloc " "
yloc<br>>> > >> >>
"draw string 0.22 "yloc-0.05" "ht<br>>> > >> >>
Endif<br>>> > >> >>
ht=ht+1<br>>> > >> >>
endwhile<br>>> > >> >> "set strsiz
0.10"<br>>> > >> >> "set string 1"<br>>>
> >> >> If (LblAxes = 1)<br>>> > >>
>> If (PageType = "Portrait")<br>>> >
>> >> "draw string 0.25 10.85
km"<br>>> > >> >> "draw
string 0.25 10.75 "htlb<br>>> > >> >>
"draw string 0.25 10.65 OBS"<br>>> > >>
>> Else<br>>> > >> >>
"draw string 0.25 8.35 km"<br>>> > >>
>> "draw string 0.25 8.25
"htlb<br>>> > >> >> "draw
string 0.25 8.15 OBS"<br>>> > >> >>
Endif<br>>> > >> >> Endif<br>>> >
>> >> Endif<br>>> > >> >><br>>> >
>> >><br>>> > >> >>
*---------------------------------------------------<br>>> >
>> >> * Draw height levels (height above MSL using Std
Atm)<br>>> > >> >>
*---------------------------------------------------<br>>> >
>> >><br>>> > >> >> If (DrawZSTD =
1)<br>>> > >> >> "set strsiz 0.10"<br>>>
> >> >>
minht=30.735*(1-pow(_pmax/1013.26,0.287))<br>>> > >> >>
minht=int(minht+0.5)<br>>> > >> >>
maxht=30.735*(1-pow(_pmin/1013.26,0.287))<br>>> > >> >>
maxht=int(maxht)<br>>> > >> >> "set gxout
stat"<br>>> > >> >> zcount=minht<br>>>
> >> >> while (zcount <= maxht)<br>>> >
>> >>
plev=1013.26*pow((1-zcount/30.735),3.4843)<br>>> > >>
>> "q w2xy 0 "plev<br>>> > >>
>> yloc=subwrd(result,6)<br>>> > >>
>> "draw string 0 "yloc-0.05" "zcount<br>>>
> >> >> zcount=zcount+1<br>>> >
>> >> endwhile<br>>> > >> >>
"set strsiz 0.10"<br>>> > >> >> If (LblAxes =
1)<br>>> > >> >> If (PageType =
"Portrait")<br>>> > >> >>
"draw string 0 10.85 km"<br>>> > >> >>
"draw string 0 10.75 ASL"<br>>> > >> >>
"draw string 0 10.65 STD"<br>>> >
>> >> Else<br>>> > >>
>> "draw string 0 8.35 km"<br>>>
> >> >> "draw string 0 8.25
ASL"<br>>> > >> >> "draw
string 0 8.15 STD"<br>>> > >> >>
Endif<br>>> > >> >> Endif<br>>> >
>> >> Endif<br>>> > >> >><br>>> >
>> >><br>>> > >> >>
*-----------------------<br>>> > >> >> * Plot
temperature lines<br>>> > >> >>
*-----------------------<br>>> > >> >><br>>> >
>> >> If (DrawTemp = 1)<br>>> > >> >>
"set strsiz 0.1"<br>>> > >> >> "set z
"_zmin " " _zmax<br>>> > >> >> "set line
"TempCol " " TempLine " "TempThk<br>>> > >> >>
"set string 1 c 3 0"<br>>> > >> >> "set gxout
stat"<br>>> > >> >>
maxtline=GetTemp(100,_pmax)<br>>> > >> >>
mintline=GetTemp(0,_pmin)<br>>> > >> >><br>>>
> >> >>
maxtline=tempint*int(maxtline/tempint)<br>>> > >> >>
mintline=tempint*int(mintline/tempint)<br>>> > >>
>><br>>> > >> >>
tloop=mintline<br>>> > >> >> While (tloop <=
maxtline)<br>>> > >> >>
Botxtemp=GetXLoc(tloop,_pmax)<br>>> > >> >>
"q w2xy "Botxtemp " " _pmax<br>>> > >> >>
Botxloc=subwrd(result,3)<br>>> > >>
>> Botyloc=byloc<br>>> > >>
>> Topxtemp=GetXLoc(tloop,_pmin)<br>>> >
>> >> "q w2xy "Topxtemp " "
_pmin<br>>> > >> >>
Topxloc=subwrd(result,3)<br>>> > >> >>
Topyloc=tyloc<br>>> > >> >>
If (Botxtemp <= 100 | Topxtemp <= 100)<br>>> > >>
>> If (Topxtemp >
100)<br>>> > >> >>
Slope=(Topyloc-Botyloc)/(Topxtemp-Botxtemp)<br>>> > >>
>>
b=Topyloc-Slope*Topxtemp<br>>> > >> >>
Topyloc=Slope*100+b<br>>> > >>
>> Topxloc=rxloc<br>>>
> >> >> Endif<br>>>
> >> >> If (Botxtemp <
0)<br>>> > >> >>
Slope=(Topyloc-Botyloc)/(Topxtemp-Botxtemp)<br>>> > >>
>>
b=Botyloc-Slope*Botxtemp<br>>> > >> >>
Botyloc=b<br>>> > >> >>
Botxloc=lxloc<br>>> >
>> >> Else<br>>> >
>> >> "set strsiz
0.15"<br>>> > >> >>
"draw string " Botxloc-0.05 " " Botyloc-0.15 " " tloop<br>>>
> >> >> Endif<br>>>
> >> >> "draw line "Botxloc
" " Botyloc " " Topxloc " " Topyloc<br>>> > >> >>
Endif<br>>> > >> >>
tloop=tloop+tempint<br>>> > >> >>
EndWhile<br>>> > >> >> If (LblAxes =
1)<br>>> > >> >> "set strsiz
0.15"<br>>> > >> >> "set string 1
c"<br>>> > >> >> If (PageType =
"Portrait")<br>>> > >> >>
"draw string 4.0 0.35 Temperatura (`3.`0C)"<br>>> > >>
>> "draw string 7.7 0.35 Viento m/s
"<br>>> > >> >> Else<br>>>
> >> >> "draw string 3.5 0.15
Temperatura (`3.`0C)"<br>>> > >> >>
Endif<br>>> > >> >> Endif<br>>> >
>> >> Endif<br>>> > >> >><br>>> >
>> >><br>>> > >> >>
*------------------<br>>> > >> >> * Plot dry
adiabats<br>>> > >> >> *------------------<br>>>
> >> >><br>>> > >> >> If (DrawThet =
1)<br>>> > >> >>
temp=GetTemp(100,_pmin)<br>>> > >> >>
maxtheta=GetThet2(temp,-100,_pmin)<br>>> > >> >>
maxtheta=thetaint*int(maxtheta/thetaint)<br>>> > >> >>
temp=GetTemp(0,_pmax)<br>>> > >> >>
mintheta=GetThet2(temp,-100,_pmax)<br>>> > >> >>
mintheta=thetaint*int(mintheta/thetaint)<br>>> > >>
>><br>>> > >> >> "set lon 0
100"<br>>> > >> >> "set y 1"<br>>> >
>> >> "set z 1"<br>>> > >> >>
tloop=mintheta<br>>> > >> >> "set line "ThetCol"
"ThetLine " "ThetThk<br>>> > >> >> While (tloop
<= maxtheta)<br>>> > >> >>
PTemp=LiftDry(tloop,1000,_pmin,1,_pmin,_pmax)<br>>> > >>
>> tloop=tloop+thetaint<br>>> > >>
>> Endwhile<br>>> > >> >>
Endif<br>>> > >> >><br>>> > >> >>
*------------------------<br>>> > >> >> * Plot mixing
ratio lines<br>>> > >> >>
*------------------------<br>>> > >> >><br>>>
> >> >> If (DrawMix = 1)<br>>> > >> >>
If (MixStop < _pmin)<br>>> > >> >>
MixStop = _pmin<br>>> > >> >>
Endif<br>>> > >> >> "set string 1 l"<br>>>
> >> >> "set z "_zmin " " _zmax<br>>> >
>> >> "set cint
1"<br>></div></div></blockquote></div><br>_______________________________________________<br>gradsusr
mailing list<br><a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div><pre>***************************************************
The information contained in this e-mail message
is intended only for the use of the recipient(s)
named above and may contain information that is
privileged, confidential, and/or proprietary.
If you are not the intended recipient, you may not
review, copy or distribute this message. If you have
received this communication in error, please notify
the sender immediately by e-mail, and delete the original message.
***************************************************
</pre></div>
<br>_______________________________________________<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>
<br></blockquote></div><br>
<pre>***************************************************
The information contained in this e-mail message
is intended only for the use of the recipient(s)
named above and may contain information that is
privileged, confidential, and/or proprietary.
If you are not the intended recipient, you may not
review, copy or distribute this message. If you have
received this communication in error, please notify
the sender immediately by e-mail, and delete the original message.
***************************************************