[gradsusr] gradsusr Digest, Vol 59, Issue 57
Geeta Geeta
geetag54 at yahoo.com
Thu Jan 29 01:32:35 EST 2015
Thanks for the help Chuck. It worked. Geeta.
On Thursday, 29 January 2015 10:13 AM, "gradsusr-request at gradsusr.org" <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: Fwrite contd. (Charles Seman - NOAA Federal)
2. Sub: mask file for US land/Ocean (Adv)
3. Great Lakes Wave upgrade (Bill Reilly)
4. Fwrite contd. (Geeta Geeta)
5. Re: gradsusr Digest, Vol 59, Issue 56 (Geeta Geeta)
----------------------------------------------------------------------
Message: 1
Date: Wed, 28 Jan 2015 14:41:46 -0500
From: Charles Seman - NOAA Federal <charles.seman at noaa.gov>
Subject: Re: [gradsusr] Fwrite contd.
To: GrADS Users Forum <gradsusr at gradsusr.org>
Message-ID: <54C93B7A.2090500 at noaa.gov>
Content-Type: text/plain; charset=windows-1252; format=flowed
Geeta,
It may be because the longitude range is (-179.875,179.875); try to use
lon, lat to define obs variables within longitude range (0.,360.).
In your GrADS script, try this:
'open rfaccum.ctl'
'set lon 0.125 359.875'
'set lat -49.875 49.875'
'set z 1' ;* just to be sure
'set t 1' ;* just to be sure
'define rf1=data(t=1)'
...
'printim figure1.gif'
'set gxout fwrite'
'set lon 0.125 359.875' ;* just to be sure
'set lat -49.875 49.875' ;* just to be sure
'set z 1' ;* just to be sure
'set t 1' ;* just to be sure
'd rf'
'disable fwrite'
Then, make a GrADS ctl file Obs22apr2011.ctl for Obs22apr2011.bin which
uses the parameters:
XDEF 1440 LINEAR 0.125 0.25
YDEF 400 LINEAR -49.875 0.25
ZDEF 1 LEVELS 1000
TDEF 1 LINEAR 00z22apr2011 1dy
now for Obs22apr2011.ctl you can use x, y coordinates...
'open Obs22apr2011.ctl'
'set gxout shaded'
'set x 1 1440'
'set y 1 400'
'set z 1'
'set t 1'
'd data'
Hope this helps,
Chuck
On 01/28/2015 02:49 AM, Geeta Geeta wrote:
> Hi Davide and all other users.
> This message is in continuation to my earlier e-mail in the previous
> Volume. I have to sum up all the 8 3 hourly TRMM binary files.
>
> Following was my CTL file.
>
> DSET ^3B42.1104%d2.%h2z.7.precipitation.bin
> options template
> TITLE Precipitation estimates
> UNDEF -9999.9
> XDEF 1440 LINEAR -179.875 0.25
> YDEF 400 LINEAR -49.875 0.25
> ZDEF 1 LEVELS 1000
> TDEF 8 LINEAR 06z22apr2011 3hr
> VARS 1
> data 0 99 3B42 TRMM precipitation estimates [mm/3hr]
> ENDVARS
>
> NOW in MY gs file, I forced grads to write the data.
> 'reinit'
> 'open rfaccum.ctl'
> 'define rf1=data(t=1)'
> 'define rf2=data(t=2)'
> 'define rf3=data(t=3)'
> 'define rf4=data(t=4)'
> 'define rf5=data(t=5)'
> 'define rf6=data(t=6)'
> 'define rf7=data(t=7)'
> 'define rf8=data(t=8)'
> 'define rf=rf1+rf2+rf3+rf4+rf5+rf6+rf7+rf8'
> 'd rf'
> 'printim figure1.gif'
> *_/'set x 1 1440'/_*
> *_/'set y 1 400 '/_*
> *_/'set fwrite Obs22apr2011.bin'/_*
> 'set gxout fwrite'
> 'd rf'
> 'disable fwrite'
> I got same size of O/P file *_/Obs22apr2011.bin. /_*
> *_/
> 2304000 Jan 27 14:17 3B42.110423.00z.7.precipitation.bin
> -rw-r--r-- 1 model staff 2304000 Jan 27 14:18
> 3B42.110423.03z.7.precipitation.bin
> -rw-r--r-- 1 model staff 2304000 Jan 27 14:22
> 3B42_daily.2011.04.22.7.bin
> -rw-r--r-- 1 model staff 2304000 Jan 28 13:12 Obs22apr2011.bin
>
> Now when I PLOT Obs22apr2011.bin file, It is only plotting the eastern
> hemisphere. (Figures attached).
> Following is the CTL file------
> SET Obs22apr2011.bin
> options big_endian
> TITLE Precipitation estimates
> UNDEF -9999.9
> XDEF 1440 LINEAR -179.875 0.25
> YDEF 400 LINEAR -49.875 0.25
> ZDEF 1 LEVELS 1000
> TDEF 1 LINEAR 00z22apr2011 1dy
> VARS 1
> data 0 99 3B42 TRMM precipitation estimates [mm/3hr]
> ENDVARS
> ~
> ~
> ~
> ~
> ~
> ~
> ~
> ~
> ~
> ~
> ~
> ~
> ~
> ~
> ~
> "dailyrf.ctl" 11 lines, 305 characters
>
>
> Kindly suggest.
> /_*
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
--
Please note that Charles.Seman at noaa.gov should be considered my NOAA
email address, not cjs at gfdl.noaa.gov.
********************************************************************
Charles Seman Charles.Seman at noaa.gov
U.S. Department of Commerce / NOAA / OAR
Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547
201 Forrestal Road fax: (609) 987-5063
Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/
********************************************************************
"The contents of this message are mine personally and do not reflect any
official or unofficial position of the United States Federal Government,
the United States Department of Commerce, or NOAA."
------------------------------
Message: 2
Date: Wed, 28 Jan 2015 16:13:13 -0600
From: Adv <advita6 at gmail.com>
Subject: [gradsusr] Sub: mask file for US land/Ocean
To: GrADS Users Forum <gradsusr at gradsusr.org>
Message-ID:
<CAHfbfHtZhn7+kjG9Sh=sp9NBGf90O1-OPtS9zFEe_FEtis1+Ow at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"
Hi,
I would like to do the same calculation shown below, only over US land
region. Could someone help me with mask file ? or else explain me how to
use basemap shape file to mask the ocean. Is that possible to do?
An example of taking an area average of data only over land, given a mask
grid:
ave(aave(maskout(p,mask.3(t=1)),x=1,x=72,y=1,y=46),t=1,t=30)
say result
Thanks,
Adv
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150128/0f53efbd/attachment-0001.html
------------------------------
Message: 3
Date: Thu, 29 Jan 2015 02:06:18 +0100
From: Bill Reilly <bill_reilly at compuserve.com>
Subject: [gradsusr] Great Lakes Wave upgrade
To: gradsusr at gradsusr.org
Message-ID: <54C9878A.8040100 at compuserve.com>
Content-Type: text/plain; charset=windows-1252; format=flowed
Hello all:
I have a question:
I downloaded the new Great Lakes Wave GRIB files today the same way I
always do, and then I used wgrib2 to prepare the data for using in GrADS
(like I've been doing every day for years), which all went fine. But
when trying to open the files in GrADS I get the error "SDF file has no
discernable X coordinate". I've never had this problem with the old GLW
data files, what's different now?... I know the file names and
resolution is different, but what about the data itself?
If anybody can help me I'd really appreciate it.
Regards,
Bill
--
Bill Reilly
(+34) 686-110-748
bill at passageweather.com
www.passageweather.com
------------------------------
Message: 4
Date: Thu, 29 Jan 2015 04:34:17 +0000 (UTC)
From: Geeta Geeta <geetag54 at yahoo.com>
Subject: [gradsusr] Fwrite contd.
To: "gradsusr at gradsusr.org" <gradsusr at gradsusr.org>
Message-ID:
<235282296.1031293.1422506057419.JavaMail.yahoo at mail.yahoo.com>
Content-Type: text/plain; charset="utf-8"
Dear?pyarimohan?
ga-> reinit
No hardcopy metafile openAll files closed; all defined objects releasedAll GrADS attributes have been reinitializedga-> run trmm3hr-accum.gsNo hardcopy metafile openAll files closed; all defined objects releasedAll GrADS attributes have been reinitializedga-> q fwriteFWrite file is closedFWrite file name is: grads.fwriteFWrite byte order is big_endian ; machine byte order is big_endianga->
the CTL file looks like...............SET Obs22apr2011.binoptions big_endianTITLE ?Precipitation estimatesUNDEF ?-9999.9XDEF ?1440 LINEAR -179.875 0.25YDEF ?400 LINEAR ?-49.875 ?0.25ZDEF ? 1 LEVELS 1000TDEF 1 LINEAR ?00z22apr2011 ?1dyVARS 1data ?0 ? 99 ?3B42 TRMM precipitation estimates [mm/3hr]ENDVARS
But still the data is plotted only in the EASTERN hemisphere.?Geeta.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150129/4316f9c4/attachment-0001.html
------------------------------
Message: 5
Date: Thu, 29 Jan 2015 10:04:59 +0530
From: Geeta Geeta <geeta124 at hotmail.com>
Subject: Re: [gradsusr] gradsusr Digest, Vol 59, Issue 56
To: grads grads <gradsusr at gradsusr.org>
Message-ID: <BLU178-W42411E3E565988F9E252A68A300 at phx.gbl>
Content-Type: text/plain; charset="iso-8859-1"
Dear pyarimohan
ga-> reinit
No hardcopy metafile openAll files closed; all defined objects releasedAll GrADS attributes have been reinitializedga-> run trmm3hr-accum.gsNo hardcopy metafile openAll files closed; all defined objects releasedAll GrADS attributes have been reinitializedga-> q fwriteFWrite file is closedFWrite file name is: grads.fwriteFWrite byte order is big_endian ; machine byte order is big_endianga->
the CTL file looks like...............SET Obs22apr2011.binoptions big_endianTITLE Precipitation estimatesUNDEF -9999.9XDEF 1440 LINEAR -179.875 0.25YDEF 400 LINEAR -49.875 0.25ZDEF 1 LEVELS 1000TDEF 1 LINEAR 00z22apr2011 1dyVARS 1data 0 99 3B42 TRMM precipitation estimates [mm/3hr]ENDVARS
But still the data is plotted only in the EASTERN hemisphere. Geeta.
> From: gradsusr-request at gradsusr.org
> Subject: gradsusr Digest, Vol 59, Issue 56
> To: gradsusr at gradsusr.org
> Date: Wed, 28 Jan 2015 12:00:02 -0500
>
> 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: Fwrite contd. (pyarimohan maharana)
> 2. New meteogram scripts (Jennifer Adams)
> 3. GrADS question re: file downloads (webmaster at solargrafix.com)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Wed, 28 Jan 2015 09:32:35 +0100
> From: pyarimohan maharana <chiku777rock at gmail.com>
> Subject: Re: [gradsusr] Fwrite contd.
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <CAGdodSdD4HmdB38Ng96-dzXmcBj8Knac_-5O8-b=e4Km=CDY0A at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Dear Geeta,
>
> use the "q fwrite" command on the GrADS.
>
> This gives you the machine byte order and the undef value.
> Edit the .CTL file and put the numbers and try to plot.
>
> This will work.
> All the best.
>
> Regards
> Pyarimohan
>
>
> On Wed, Jan 28, 2015 at 8:49 AM, Geeta Geeta <geeta124 at hotmail.com> wrote:
>
> > Hi Davide and all other users.
> > This message is in continuation to my earlier e-mail in the previous
> > Volume. I have to sum up all the 8 3 hourly TRMM binary files.
> >
> > Following was my CTL file.
> >
> > DSET ^3B42.1104%d2.%h2z.7.precipitation.bin
> > options template
> > TITLE Precipitation estimates
> > UNDEF -9999.9
> > XDEF 1440 LINEAR -179.875 0.25
> > YDEF 400 LINEAR -49.875 0.25
> > ZDEF 1 LEVELS 1000
> > TDEF 8 LINEAR 06z22apr2011 3hr
> > VARS 1
> > data 0 99 3B42 TRMM precipitation estimates [mm/3hr]
> > ENDVARS
> >
> > NOW in MY gs file, I forced grads to write the data.
> > 'reinit'
> > 'open rfaccum.ctl'
> > 'define rf1=data(t=1)'
> > 'define rf2=data(t=2)'
> > 'define rf3=data(t=3)'
> > 'define rf4=data(t=4)'
> > 'define rf5=data(t=5)'
> > 'define rf6=data(t=6)'
> > 'define rf7=data(t=7)'
> > 'define rf8=data(t=8)'
> > 'define rf=rf1+rf2+rf3+rf4+rf5+rf6+rf7+rf8'
> > 'd rf'
> > 'printim figure1.gif'
> > *'set x 1 1440'*
> > *'set y 1 400 '*
> > *'set fwrite Obs22apr2011.bin'*
> > 'set gxout fwrite'
> > 'd rf'
> > 'disable fwrite'
> > I got same size of O/P file *Obs22apr2011.bin. *
> >
> >
> >
> > *2304000 Jan 27 14:17 3B42.110423.00z.7.precipitation.bin-rw-r--r-- 1
> > model staff 2304000 Jan 27 14:18
> > 3B42.110423.03z.7.precipitation.bin-rw-r--r-- 1 model staff
> > 2304000 Jan 27 14:22 3B42_daily.2011.04.22.7.bin-rw-r--r-- 1 model
> > staff 2304000 Jan 28 13:12 Obs22apr2011.binNow when I PLOT
> > Obs22apr2011.bin file, It is only plotting the eastern hemisphere. (Figures
> > attached). Following is the CTL file------SET Obs22apr2011.binoptions
> > big_endianTITLE Precipitation estimatesUNDEF -9999.9XDEF 1440 LINEAR
> > -179.875 0.25YDEF 400 LINEAR -49.875 0.25ZDEF 1 LEVELS 1000TDEF 1
> > LINEAR 00z22apr2011 1dyVARS 1data 0 99 3B42 TRMM precipitation
> > estimates [mm/3hr]ENDVARS~~~~~~~~~~~~~~~"dailyrf.ctl" 11 lines, 305
> > charactersKindly suggest. *
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150128/b3dd0576/attachment-0001.html
>
> ------------------------------
>
> Message: 2
> Date: Wed, 28 Jan 2015 09:27:55 -0500
> From: Jennifer Adams <jma at cola.iges.org>
> Subject: [gradsusr] New meteogram scripts
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID: <BDB74825-E437-4D2D-A42C-62C515D2D551 at cola.iges.org>
> Content-Type: text/plain; charset="windows-1252"
>
> Dear All,
> I have posted my new scripts to download the required data and draw a 10-day meteogram. Some comments appear at the top of each script with instructions on how to use them. In my testing, it took about 10-12 minutes to download the data to a local file.
>
> ftp://cola.gmu.edu/grads/scripts/meteogram_subset_GDS.gs
> ftp://cola.gmu.edu/grads/scripts/meteogram_GDS.gs
>
> Links in the GrADS script library http://iges.org/grads/gadoc/library.html has also been updated.
> For best results, use GrADS version 2.1.
> ?Jennifer
>
> --
> Jennifer M. Adams
> Center for Ocean-Land-Atmosphere Studies (COLA)
> 111 Research Hall, Mail Stop 2B3
> George Mason University
> 4400 University Drive
> Fairfax, VA 22030
>
>
>
>
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150128/fab17869/attachment-0001.html
>
> ------------------------------
>
> Message: 3
> Date: Wed, 28 Jan 2015 10:38:49 -0500
> From: <webmaster at solargrafix.com>
> Subject: [gradsusr] GrADS question re: file downloads
> To: <gradsusr at gradsusr.org>
> Message-ID: <E98D6BC6ACF04CF3971ACA1F531E5862 at JasonPC>
> Content-Type: text/plain; charset="us-ascii"
>
> Hello
>
>
>
> Is there a way to get GrADS to download the most recent file from the GrADS
> Data Server. I see all of the files on that server are saved in directories
> according to time. I can download files using sdfopen, but I would like to
> automate the process, and have GrADS download the most recent file.
>
>
>
> Thanks,
>
> Jason Warren
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150128/c5ffbde9/attachment-0001.html
>
> ------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> End of gradsusr Digest, Vol 59, Issue 56
> ****************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150129/a73becce/attachment.html
------------------------------
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr
End of gradsusr Digest, Vol 59, Issue 57
****************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150129/2166623c/attachment-0001.html
More information about the gradsusr
mailing list