[gradsusr] GrADS help, plotting shaded contour precip/ptype products over a basemap/background
Stephen McMillan
smcmillan at planalytics.com
Tue Nov 9 14:25:12 EST 2010
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.
***************************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101109/ca536ce7/attachment-0003.html
More information about the gradsusr
mailing list