[gradsusr] GrADS help, plotting shaded contour precip/ptype products over a basemap/background

Stephen McMillan smcmillan at planalytics.com
Wed Nov 10 22:05:50 EST 2010


Meredith,
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.

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.

Best wishes...
Stephen Mc

On Wed, Nov 10, 2010 at 4:26 PM, Meredith Croke <mcroke at airdat.com> wrote:

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

***************************************************
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.
***************************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101110/aaae8f02/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: backmap.gs
Type: application/octet-stream
Size: 1432 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20101110/aaae8f02/attachment-0003.obj 


More information about the gradsusr mailing list