[gradsusr] gradsusr Digest, Vol 2, Issue 40

Geeta Geeta geeta124 at hotmail.com
Tue Apr 27 12:09:16 EDT 2010 

Dear Stephen, 
Thank you very much for all the help. 
It is plotting and also printing the TMPprs. 
Can I write this output to a file?
I tried with fwrite but this does not seem to work. 
Kindly suggest. 

geeta



> From: gradsusr-request at gradsusr.org
> Subject: gradsusr Digest, Vol 2, Issue 40
> To: gradsusr at gradsusr.org
> Date: Tue, 27 Apr 2010 11:01:22 -0400
> 
> 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: Gribmap (James Burn)
>    2. Re: error with gr2stn function (Stephen McMillan)
> 
> 
> ----------------------------------------------------------------------
> 
> Message: 1
> Date: Tue, 27 Apr 2010 09:38:50 +0100
> From: James Burn <james at weatherbadger.com>
> Subject: Re: [gradsusr] Gribmap
> To: gradsusr at gradsusr.org
> Message-ID: <201004270938.51071.james at weatherbadger.com>
> Content-Type: text/plain; charset="iso-8859-15"
> 
> Andrew
> In my cron-run shell scripts I change directory to the location of the data 
> first before running gridmap.
> 
> "cd $HOME/grads_etc/data " 
> 
> I've had it like this for some time, and can't remember why I did it like 
> this, but perhaps it was to get around the problem you've got?
> 
> Good luck
> James
> 
> On Tuesday 27 April 2010 06:18:58 Andrew Revering wrote:
> > I'm new to Linux and I'm pulling my hair out trying to figure out why
> > gribmap will not run when executed as a cron job.
> >
> >
> >
> > In crontab I execute a .sh script that has something that looks like:
> >
> >
> >
> > #!/bin/bash
> >
> >
> >
> > /wrf/wrfems/util/grads/bin/grib2ctl.pl -verf
> > /wrf/wrfems/runs/float/emsprd/grib/arw_0000.grib >
> > /wrf/wrfems/runs/float/emsprd/grib/arw_0000.ctl
> >
> >
> >
> > /wrf/wrfems/util/grads/bin/gribmap -I
> > /wrf/wrfems/runs/float/emsprd/grib/arw_0000.ctl
> >
> >
> >
> >
> >
> >
> >
> > The last line is the problem child. Note that the script is running, and
> > the first line in the script (converting GRIB to a CTL file) works, but
> > nothing happens when trying to execute the gribmap executable.
> >
> >
> >
> > Also I should note that if I type that same line in the terminal it DOES
> > work, just not in crontab.
> >
> >
> >
> > Any ideas?
> 
> -- 
> 
> Reply to : james at weatherbadger.com
> Visit http://www.weatherbadger.com for some weather forecasts.
> 
> 
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100427/892ffb92/attachment-0001.html 
> 
> ------------------------------
> 
> Message: 2
> Date: Tue, 27 Apr 2010 11:01:13 -0400
> From: Stephen McMillan <smcmillan at planalytics.com>
> Subject: Re: [gradsusr] error with gr2stn function
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> 	<s2v1f726e2f1004270801u5011e3cbu4b58b50581a43120 at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
> 
> Geeta,
> Try setting or re-setting your prnopts, for example:
> 'set gxout print'
> 'set prnopts %0.2f'
> 'd gr2stn(TMPprs,77.5,12.5)'
> 
> If you still don't get a result...
> If you display a map with x and y varying, t = 1, and lev = 1000, do you get
> a sensible image of the variable tmpprs?  Do you also get a reasonable value
> if you don't interpolate?, for example:
> 'd TMPprs(lon=77.5,lat=12.5)'  ...should give you the same value as 'd
> gr2stn(TMPprs,77.5,12.5)' if both are at the same data gridpoint.
> 
> Good luck!
> Stephen McMillan
> 
> On Tue, Apr 27, 2010 at 1:36 AM, Geeta Geeta <geeta124 at hotmail.com> wrote:
> 
> >  Dear All,
> > I have tried to fix the x and y as suggested. But I have not got the result
> > . I am posting the message given by grads.
> >
> > ga-> reinit
> > No hardcopy metafile open
> > All files closed; all defined objects released
> > All GrADS attributes have been reinitialized
> > ga-> open 12mar2008
> > Scanning description file:  12mar2008
> > Data file gdas2008031200f00 is open as file 1
> > LON set to 0 360
> > LAT set to -90 90
> > LEV set to 1000 1000
> > Time values set: 2008:3:12:0 2008:3:12:0
> > ga-> set t 1
> > Time values set: 2008:3:12:0 2008:3:12:0
> > ga-> set x 1
> > LON set to 0 0
> > ga-> set y 1
> > LAT set to -90 -90
> > ga-> set lev 1000
> > LEV set to 1000 1000
> > ga-> set gxout print
> > ga-> d gr2stn(TMPprs,77.5,12.5)
> > Printing Stations -- 1 Reports -- Undef = 9.999e+20
> > xxxxxxxx 77.5     12.5     1000
> >
> > ga->
> > ga-> d gr2stn(TMPprs,0.0,-90.0)
> > Printing Stations -- 1 Reports -- Undef = 9.999e+20
> > xxxxxxxx 0        -90      1000
> >
> > ga-> d gr2stn(TMPprs,50.0,-45.0)
> > Printing Stations -- 1 Reports -- Undef = 9.999e+20
> > xxxxxxxx 50       -45      1000
> >
> > ga->
> >
> >
> > It is printing the lat, lon and the level only.
> > Kindly suggest.
> > geeta
> >
> >
> >
> > > From: gradsusr-request at gradsusr.org
> > > Subject: gradsusr Digest, Vol 2, Issue 34
> > > To: gradsusr at gradsusr.org
> > > Date: Mon, 26 Apr 2010 02:46:40 -0400
> > >
> > > 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. gr2stn function (Geeta Geeta)
> > > 2. scorr function (Carlos Batista)
> > > 3. Re: gr2stn function (Stephen McMillan)
> > > 4. Re: scorr function (Stephen McMillan)
> > > 5. Unable to read Monthly Northern Hemisphere Sea-Level Pressure
> > > (DSS SLP Grid-binary) (zowiti at ncst.go.ke)
> > > 6. error in gr2stn script (Geeta Geeta)
> > >
> > >
> > > ----------------------------------------------------------------------
> > >
> > > Message: 1
> > > Date: Sat, 24 Apr 2010 07:45:42 +0000
> > > From: Geeta Geeta <geeta124 at hotmail.com>
> > > Subject: [gradsusr] gr2stn function
> > > To: <gradsusr at gradsusr.org>
> > > Message-ID: <SNT117-W12B787EAD3ABD70AA641528A060 at phx.gbl>
> > > Content-Type: text/plain; charset="iso-8859-1"
> > >
> > >
> > > Dear all,
> > > Kindly tell me how to use the "gr2stn " function of grads.
> > > How the station data file will look like?
> > > Any suggestions
> > > geeta
> > >
> > >
> > > _________________________________________________________________
> > > The amazing world in sharp snaps
> > > http://news.in.msn.com/gallery/archive.aspx
> > > -------------- next part --------------
> > > An HTML attachment was scrubbed...
> > > URL:
> > http://gradsusr.org/pipermail/gradsusr/attachments/20100424/61b8f017/attachment-0001.html
> > >
> > > ------------------------------
> > >
> > > Message: 2
> > > Date: Sat, 24 Apr 2010 11:02:31 -0300
> > > From: Carlos Batista <krlosbatist at gmail.com>
> > > Subject: [gradsusr] scorr function
> > > To: gradsusr at gradsusr.org
> > > Message-ID:
> > > <y2gdd4d468f1004240702zaf186fd7ufa21b4b504df48a1 at mail.gmail.com>
> > > Content-Type: text/plain; charset="iso-8859-1"
> > >
> > > Dear all,
> > > Kindly tell me how to use the "scorr " function of grads.
> > > I don't how use its this function for to do correlation between
> > > precipitation and flow rivers. Someone know some gs for to do
> > correlations
> > > between this data, please, useding the scorr function and other functions
> > ?
> > > Any suggestions
> > > Krlos
> > > e-mail; krlosbatist at gmail.com
> > > -------------- next part --------------
> > > An HTML attachment was scrubbed...
> > > URL:
> > http://gradsusr.org/pipermail/gradsusr/attachments/20100424/c4d8f6e8/attachment-0001.html
> > >
> > > ------------------------------
> > >
> > > Message: 3
> > > Date: Sat, 24 Apr 2010 10:42:04 -0400
> > > From: Stephen McMillan <smcmillan at planalytics.com>
> > > Subject: Re: [gradsusr] gr2stn function
> > > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > > Message-ID:
> > > <s2u1f726e2f1004240742j1329b506kcf3e48202ec2924c at mail.gmail.com>
> > > Content-Type: text/plain; charset="iso-8859-1"
> > >
> > > Dear geeta,
> > > Please see this post on the the gradsusr archives:
> > > http://gradsusr.org/pipermail/gradsusr/2010-March/000077.html
> > > You can search for additional posts at
> > > http://gradsusr.org/pipermail/gradsusr/
> > > GrADS documentation at:
> > http://sidata.org/grads/gadoc/gradfuncgr2stn.html
> > >
> > > Good luck--
> > > Stephen McMillan
> > >
> > > On Sat, Apr 24, 2010 at 3:45 AM, Geeta Geeta <geeta124 at hotmail.com>
> > wrote:
> > >
> > > > Dear all,
> > > > Kindly tell me how to use the "gr2stn " function of grads.
> > > > How the station data file will look like?
> > > > Any suggestions
> > > > geeta
> > > >
> > > >
> > > > ------------------------------
> > > > Catch the changing security environment Get it now.<
> > http://news.in.msn.com/internalsecurity/>
> > > >
> > > > _______________________________________________
> > > > 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/20100424/591f216e/attachment-0001.html
> > >
> > > ------------------------------
> > >
> > > Message: 4
> > > Date: Sat, 24 Apr 2010 12:03:04 -0400
> > > From: Stephen McMillan <smcmillan at planalytics.com>
> > > Subject: Re: [gradsusr] scorr function
> > > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > > Message-ID:
> > > <j2m1f726e2f1004240903q63515e10pcba71b2095589713 at mail.gmail.com>
> > > Content-Type: text/plain; charset="iso-8859-1"
> > >
> > > Dear Krlos,
> > > If you haven't already, first review the correct syntax and example in
> > the
> > > GrADS Documentation Index:
> > http://sidata.org/grads/gadoc/gradfuncscorr.html
> > > Other functions can be found in the same Documentation Index.
> > >
> > > The basic idea is that you have two previously defined variables, whose
> > > horizontal dimensions both include the geographical area you wish to
> > > spatially correlate. For example, let's say you have two variables var1
> > and
> > > var2 and you wanted to correlate them over a 10 x 10 degree domain over
> > > somewhere at a fixed time:
> > > 'd scorr(var1,var2,lon=-90,lon=-80,lat=30,lat=40)'
> > > should give you something like this:
> > > result value = -0.2356
> > > You could capture above value in script by including this after display
> > > command:
> > > corr=subwrd(result,4)
> > > You can also create a time loop of correlations, 50 timesteps for
> > example:
> > > 'set t 1 50'
> > > 'd tloop(scorr(var1,var2,....)'
> > > Create a correlation time series of above loop:
> > > 'define corrts=tloop(scorr(var1,var2,...)'
> > > Display the time series:
> > > 'set x 1';'set y 1'
> > > 'd corrts'
> > >
> > > The tcorr function (see http://sidata.org/grads/gadoc/gradfunctcorr.html
> > )
> > > may be more appropriate for correlating your data, especially if they are
> > > time series at a fixed point.
> > >
> > > A search of the GrADS user archive (
> > http://gradsusr.org/pipermail/gradsusr/)
> > > may yield more appropriate answers.
> > >
> > > Best wishes--
> > > Stephen McMillan
> > >
> > > On Sat, Apr 24, 2010 at 10:02 AM, Carlos Batista <krlosbatist at gmail.com
> > >wrote:
> > >
> > > > Dear all,
> > > > Kindly tell me how to use the "scorr " function of grads.
> > > > I don't how use its this function for to do correlation between
> > > > precipitation and flow rivers. Someone know some gs for to do
> > correlations
> > > > between this data, please, useding the scorr function and other
> > functions ?
> > > > Any suggestions
> > > > Krlos
> > > > e-mail; krlosbatist at gmail.com
> > > >
> > > >
> > > > _______________________________________________
> > > > 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/20100424/cadeffd4/attachment-0001.html
> > >
> > > ------------------------------
> > >
> > > Message: 5
> > > Date: Sun, 25 Apr 2010 06:37:16 -0700
> > > From: "zowiti at ncst.go.ke" <zowiti at ncst.go.ke>
> > > Subject: [gradsusr] Unable to read Monthly Northern Hemisphere
> > > Sea-Level Pressure (DSS SLP Grid-binary)
> > > To: <gradsusr at gradsusr.org>
> > > Message-ID: <4363b2a4$4cb7762b$7c5d40f6$@com>
> > > Content-Type: text/plain; charset="iso-8859-1"
> > >
> > > Dear Users,
> > > I am reading the DSS SLP [5? x 5?, 0E to 355E and 15N to 90N (72 x 16
> > Longitude/Latitude)] data in Grads via the ctl below. When I display the
> > slp, I get a distorted plot. I'm unable to figure out the problem in the
> > ctl.
> > >
> > > dset d:\jf\molydata.bin
> > > options sequential
> > > undef -9999
> > > title dss_slp
> > > xdef 72 linear 0 5.0
> > > ydef 15 linear 15 5.0
> > > zdef 1 levels 1000
> > > TDEF 1335 LINEAR jan1899 1mo
> > > vars 1
> > > slp 0 99 pressure
> > > endvars
> > > Additional data information:
> > > 1080 binary integer data values in 15-bit bytes. Data values are to
> > nearest tenth millibar and are the rightmost 14 bits of each byte. The grid
> > of Y(I,J) is dimensioned Y(72,15). As I goes from 1 to 72, longitude goes
> > from 00E-355E; as J goes from 1 to 15, latitude goes from 15N-85N.
> > >
> > > ZABLONE OWITI
> > > GRADUATE STUDENT
> > > Nanjing University of Information, Science and Technology
> > > College of International Education
> > >
> > > ====================================================
> > >
> > > -------------- next part --------------
> > > An HTML attachment was scrubbed...
> > > URL:
> > http://gradsusr.org/pipermail/gradsusr/attachments/20100425/ccb08070/attachment-0001.html
> > >
> > > ------------------------------
> > >
> > > Message: 6
> > > Date: Mon, 26 Apr 2010 06:46:39 +0000
> > > From: Geeta Geeta <geeta124 at hotmail.com>
> > > Subject: [gradsusr] error in gr2stn script
> > > To: grads grads <gradsusr at gradsusr.org>
> > > Message-ID: <SNT117-W28F33B5C383D2CB799A3548A040 at phx.gbl>
> > > Content-Type: text/plain; charset="iso-8859-1"
> > >
> > >
> > > Dear All,
> > > I am trying to use the "gr2stn " function of grads but failed. I am
> > posting the following message from grads. My gridded data file is the gfs
> > grib files.
> > >
> > > ga-> reinit
> > > No hardcopy metafile open
> > > All files closed; all defined objects released
> > > All GrADS attributes have been reinitialized
> > > ga-> open 12mar2008
> > > Scanning description file: 12mar2008
> > > Data file gdas2008031200f00 is open as file 1
> > > LON set to 0 360
> > > LAT set to -90 90
> > > LEV set to 1000 1000
> > > Time values set: 2008:3:12:0 2008:3:12:0
> > > ga-> set t 1
> > > Time values set: 2008:3:12:0 2008:3:12:0
> > > ga-> set lev 1000
> > > LEV set to 1000 1000
> > > ga-> d gr2stn(TMP2M,77.5,12.5)
> > > Error from GR2STN: 2nd argument must be stns
> > > Operation Error: Error from gr2stn function
> > > Error ocurred at column 1
> > > DISPLAY error: Invalid expression
> > > Expression = gr2stn(tmp2m,77.5,12.5)
> > >
> > >
> > >
> > > If I use this function with the WRF output, I get the following message.
> > > ga-> reinit
> > > No hardcopy metafile open
> > > All files closed; all defined objects released
> > > All GrADS attributes have been reinitialized
> > > ga-> open test
> > > Scanning description file: test
> > > Data file test.dat is open as file 1
> > > LON set to 60.2713 89.7287
> > > LAT set to 0.57003 23.9053
> > > LEV set to 1000 1000
> > > Time values set: 2010:4:24:0 2010:4:24:0
> > > ga-> set t 1
> > > Time values set: 2010:4:24:0 2010:4:24:0
> > > ga-> set lev 1000
> > > LEV set to 1000 1000
> > > ga-> d gr2stn(T2, 77.5, 12.5)
> > > Error from GR2STN: Grid must have linear scaling
> > > Operation Error: Error from gr2stn function
> > > Error ocurred at column 1
> > > DISPLAY error: Invalid expression
> > > Expression = gr2stn(t2,77.5,12.5)
> > > ga->
> > >
> > > I Think that since WRF output does not have fixed dx and dy so the error.
> >
> > >
> > > Kindly suggest.
> > > Geeta
> > >
> > >
> > > _________________________________________________________________
> > > Catch the latest in the world of fashion
> > > http://lifestyle.in.msn.com/
> > > -------------- next part --------------
> > > An HTML attachment was scrubbed...
> > > URL:
> > http://gradsusr.org/pipermail/gradsusr/attachments/20100426/690c0685/attachment.html
> > >
> > > ------------------------------
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> > > End of gradsusr Digest, Vol 2, Issue 34
> > > ***************************************
> >
> > ------------------------------
> > Invest your money wisely post Budget Sign up now.<http://news.in.msn.com/moneyspecial/>
> >
> > _______________________________________________
> > 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/20100427/e3848da1/attachment.html 
> 
> ------------------------------
> 
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
> 
> 
> End of gradsusr Digest, Vol 2, Issue 40
> ***************************************
 		 	   		  
_________________________________________________________________
Bollywood This Decade
http://entertainment.in.msn.com/bollywoodthisdecade/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100427/4149d168/attachment-0003.html

More information about the gradsusr mailing list