[gradsusr] error with gr2stn function

Geeta Geeta geeta124 at hotmail.com
Tue Apr 27 01:36:03 EDT 2010


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
> ***************************************
 		 	   		  
_________________________________________________________________
Climate, controversies and the changing signatures of nature
http://green.in.msn.com/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100427/883cd119/attachment-0003.html 


More information about the gradsusr mailing list