[gradsusr] error with gr2stn function

Tue Apr 27 11:01:13 EDT 2010

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-0003.html 