Geeta,<div>Try setting or re-setting your prnopts, for example:</div><div>'set gxout print'</div><div>'set prnopts %0.2f'</div><div>'d gr2stn(TMPprs,77.5,12.5)'</div><div><br></div><div>If you still don't get a result...</div>
<div>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:</div><div>'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.</div>
<div><br></div><div>Good luck!</div><div>Stephen McMillan<br><br><div class="gmail_quote">On Tue, Apr 27, 2010 at 1:36 AM, Geeta Geeta <span dir="ltr"><<a href="mailto:geeta124@hotmail.com">geeta124@hotmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div>
Dear All, <br>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. <br><br>ga-> reinit<br>No hardcopy metafile open<br>All files closed; all defined objects released<br>
All GrADS attributes have been reinitialized<br>ga-> open 12mar2008<br>Scanning description file: 12mar2008<br>Data file gdas2008031200f00 is open as file 1<br>LON set to 0 360 <br>LAT set to -90 90 <br>LEV set to 1000 1000 <br>
Time values set: 2008:3:12:0 2008:3:12:0 <br>ga-> set t 1<br>Time values set: 2008:3:12:0 2008:3:12:0 <br>ga-> set x 1<br>LON set to 0 0 <br>ga-> set y 1<br>LAT set to -90 -90 <br>ga-> set lev 1000<br>LEV set to 1000 1000 <br>
ga-> set gxout print<br>ga-> d gr2stn(TMPprs,77.5,12.5)<br>Printing Stations -- 1 Reports -- Undef = 9.999e+20<br>xxxxxxxx 77.5 12.5 1000 <br><br>ga-> <br>ga-> d gr2stn(TMPprs,0.0,-90.0)<br>Printing Stations -- 1 Reports -- Undef = 9.999e+20<br>
xxxxxxxx 0 -90 1000 <br><br>ga-> d gr2stn(TMPprs,50.0,-45.0)<br>Printing Stations -- 1 Reports -- Undef = 9.999e+20<br>xxxxxxxx 50 -45 1000 <br><br>ga-> <br><br><br>It is printing the lat, lon and the level only. <br>
Kindly suggest. <br>geeta<br><br><br><br>> From: <a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>> Subject: gradsusr Digest, Vol 2, Issue 34<br>> To: <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
> Date: Mon, 26 Apr 2010 02:46:40 -0400<br>> <br>> Send gradsusr mailing list submissions to<br>> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>> <br>> To subscribe or unsubscribe via the World Wide Web, visit<br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>> or, via email, send a message with subject or body 'help' to<br>> <a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>
> <br>> You can reach the person managing the list at<br>> <a href="mailto:gradsusr-owner@gradsusr.org" target="_blank">gradsusr-owner@gradsusr.org</a><br>> <br>> When replying, please edit your Subject line so it is more specific<br>
> than "Re: Contents of gradsusr digest..."<br>> <br>> <br>> Today's Topics:<br>> <br>> 1. gr2stn function (Geeta Geeta)<br>> 2. scorr function (Carlos Batista)<br>> 3. Re: gr2stn function (Stephen McMillan)<br>
> 4. Re: scorr function (Stephen McMillan)<br>> 5. Unable to read Monthly Northern Hemisphere Sea-Level Pressure<br>> (DSS SLP Grid-binary) (<a href="mailto:zowiti@ncst.go.ke" target="_blank">zowiti@ncst.go.ke</a>)<br>
> 6. error in gr2stn script (Geeta Geeta)<br>> <br>> <br>> ----------------------------------------------------------------------<br>> <br>> Message: 1<br>> Date: Sat, 24 Apr 2010 07:45:42 +0000<br>
> From: Geeta Geeta <<a href="mailto:geeta124@hotmail.com" target="_blank">geeta124@hotmail.com</a>><br>> Subject: [gradsusr] gr2stn function<br>> To: <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
> Message-ID: <SNT117-W12B787EAD3ABD70AA641528A060@phx.gbl><br>> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> <br>> Dear all, <br>> Kindly tell me how to use the "gr2stn " function of grads. <br>
> How the station data file will look like?<br>> Any suggestions<br>> geeta<br>> <br>> <br>> _________________________________________________________________<br>> The amazing world in sharp snaps<br>
> <a href="http://news.in.msn.com/gallery/archive.aspx" target="_blank">http://news.in.msn.com/gallery/archive.aspx</a><br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20100424/61b8f017/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20100424/61b8f017/attachment-0001.html</a> <br>
> <br>> ------------------------------<br>> <br>> Message: 2<br>> Date: Sat, 24 Apr 2010 11:02:31 -0300<br>> From: Carlos Batista <<a href="mailto:krlosbatist@gmail.com" target="_blank">krlosbatist@gmail.com</a>><br>
> Subject: [gradsusr] scorr function<br>> To: <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>> Message-ID:<br>> <<a href="mailto:y2gdd4d468f1004240702zaf186fd7ufa21b4b504df48a1@mail.gmail.com" target="_blank">y2gdd4d468f1004240702zaf186fd7ufa21b4b504df48a1@mail.gmail.com</a>><br>
> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> Dear all,<br>> Kindly tell me how to use the "scorr " function of grads.<br>> I don't how use its this function for to do correlation between<br>
> precipitation and flow rivers. Someone know some gs for to do correlations<br>> between this data, please, useding the scorr function and other functions ?<br>> Any suggestions<br>> Krlos<br>> e-mail; <a href="mailto:krlosbatist@gmail.com" target="_blank">krlosbatist@gmail.com</a><br>
> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20100424/c4d8f6e8/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20100424/c4d8f6e8/attachment-0001.html</a> <br>
> <br>> ------------------------------<br>> <br>> Message: 3<br>> Date: Sat, 24 Apr 2010 10:42:04 -0400<br>> From: Stephen McMillan <<a href="mailto:smcmillan@planalytics.com" target="_blank">smcmillan@planalytics.com</a>><br>
> Subject: Re: [gradsusr] gr2stn function<br>> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>> Message-ID:<br>> <<a href="mailto:s2u1f726e2f1004240742j1329b506kcf3e48202ec2924c@mail.gmail.com" target="_blank">s2u1f726e2f1004240742j1329b506kcf3e48202ec2924c@mail.gmail.com</a>><br>
> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> Dear geeta,<br>> Please see this post on the the gradsusr archives:<br>> <a href="http://gradsusr.org/pipermail/gradsusr/2010-March/000077.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/2010-March/000077.html</a><br>
> You can search for additional posts at<br>> <a href="http://gradsusr.org/pipermail/gradsusr/" target="_blank">http://gradsusr.org/pipermail/gradsusr/</a><br>> GrADS documentation at: <a href="http://sidata.org/grads/gadoc/gradfuncgr2stn.html" target="_blank">http://sidata.org/grads/gadoc/gradfuncgr2stn.html</a><br>
> <br>> Good luck--<br>> Stephen McMillan<br>> <br>> On Sat, Apr 24, 2010 at 3:45 AM, Geeta Geeta <<a href="mailto:geeta124@hotmail.com" target="_blank">geeta124@hotmail.com</a>> wrote:<br>> <br>> > Dear all,<br>
> > Kindly tell me how to use the "gr2stn " function of grads.<br>> > How the station data file will look like?<br>> > Any suggestions<br>> > geeta<br>> ><br>> ><br>> > ------------------------------<br>
> > Catch the changing security environment Get it now.<<a href="http://news.in.msn.com/internalsecurity/" target="_blank">http://news.in.msn.com/internalsecurity/</a>><br>> ><br>> > _______________________________________________<br>
> > gradsusr mailing list<br>> > <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> ><br>> ><br>> <br>> ***************************************************<br>> The information contained in this e-mail message <br>> is intended only for the use of the recipient(s) <br>> named above and may contain information that is <br>
> privileged, confidential, and/or proprietary. <br>> If you are not the intended recipient, you may not<br>> review, copy or distribute this message. If you have<br>> received this communication in error, please notify <br>
> the sender immediately by e-mail, and delete the original message.<br>> ***************************************************<br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>
> URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20100424/591f216e/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20100424/591f216e/attachment-0001.html</a> <br>
> <br>> ------------------------------<br>> <br>> Message: 4<br>> Date: Sat, 24 Apr 2010 12:03:04 -0400<br>> From: Stephen McMillan <<a href="mailto:smcmillan@planalytics.com" target="_blank">smcmillan@planalytics.com</a>><br>
> Subject: Re: [gradsusr] scorr function<br>> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>> Message-ID:<br>> <<a href="mailto:j2m1f726e2f1004240903q63515e10pcba71b2095589713@mail.gmail.com" target="_blank">j2m1f726e2f1004240903q63515e10pcba71b2095589713@mail.gmail.com</a>><br>
> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> Dear Krlos,<br>> If you haven't already, first review the correct syntax and example in the<br>> GrADS Documentation Index: <a href="http://sidata.org/grads/gadoc/gradfuncscorr.html" target="_blank">http://sidata.org/grads/gadoc/gradfuncscorr.html</a><br>
> Other functions can be found in the same Documentation Index.<br>> <br>> The basic idea is that you have two previously defined variables, whose<br>> horizontal dimensions both include the geographical area you wish to<br>
> spatially correlate. For example, let's say you have two variables var1 and<br>> var2 and you wanted to correlate them over a 10 x 10 degree domain over<br>> somewhere at a fixed time:<br>> 'd scorr(var1,var2,lon=-90,lon=-80,lat=30,lat=40)'<br>
> should give you something like this:<br>> result value = -0.2356<br>> You could capture above value in script by including this after display<br>> command:<br>> corr=subwrd(result,4)<br>> You can also create a time loop of correlations, 50 timesteps for example:<br>
> 'set t 1 50'<br>> 'd tloop(scorr(var1,var2,....)'<br>> Create a correlation time series of above loop:<br>> 'define corrts=tloop(scorr(var1,var2,...)'<br>> Display the time series:<br>
> 'set x 1';'set y 1'<br>> 'd corrts'<br>> <br>> The tcorr function (see <a href="http://sidata.org/grads/gadoc/gradfunctcorr.html" target="_blank">http://sidata.org/grads/gadoc/gradfunctcorr.html</a>)<br>
> may be more appropriate for correlating your data, especially if they are<br>> time series at a fixed point.<br>> <br>> A search of the GrADS user archive (<a href="http://gradsusr.org/pipermail/gradsusr/" target="_blank">http://gradsusr.org/pipermail/gradsusr/</a>)<br>
> may yield more appropriate answers.<br>> <br>> Best wishes--<br>> Stephen McMillan<br>> <br>> On Sat, Apr 24, 2010 at 10:02 AM, Carlos Batista <<a href="mailto:krlosbatist@gmail.com" target="_blank">krlosbatist@gmail.com</a>>wrote:<br>
> <br>> > Dear all,<br>> > Kindly tell me how to use the "scorr " function of grads.<br>> > I don't how use its this function for to do correlation between<br>> > precipitation and flow rivers. Someone know some gs for to do correlations<br>
> > between this data, please, useding the scorr function and other functions ?<br>> > Any suggestions<br>> > Krlos<br>> > e-mail; <a href="mailto:krlosbatist@gmail.com" target="_blank">krlosbatist@gmail.com</a><br>
> ><br>> ><br>> > _______________________________________________<br>> > gradsusr mailing list<br>> > <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> ><br>> ><br>> <br>> ***************************************************<br>> The information contained in this e-mail message <br>> is intended only for the use of the recipient(s) <br>> named above and may contain information that is <br>
> privileged, confidential, and/or proprietary. <br>> If you are not the intended recipient, you may not<br>> review, copy or distribute this message. If you have<br>> received this communication in error, please notify <br>
> the sender immediately by e-mail, and delete the original message.<br>> ***************************************************<br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>
> URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20100424/cadeffd4/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20100424/cadeffd4/attachment-0001.html</a> <br>
> <br>> ------------------------------<br>> <br>> Message: 5<br>> Date: Sun, 25 Apr 2010 06:37:16 -0700<br>> From: "<a href="mailto:zowiti@ncst.go.ke" target="_blank">zowiti@ncst.go.ke</a>" <<a href="mailto:zowiti@ncst.go.ke" target="_blank">zowiti@ncst.go.ke</a>><br>
> Subject: [gradsusr] Unable to read Monthly Northern Hemisphere<br>> Sea-Level Pressure (DSS SLP Grid-binary)<br>> To: <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
> Message-ID: <4363b2a4$4cb7762b$7c5d40f6$@com><br>> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> Dear Users,<br>> 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. <br>
> <br>> dset d:\jf\molydata.bin <br>> options sequential <br>> undef -9999 <br>> title dss_slp <br>> xdef 72 linear 0 5.0 <br>> ydef 15 linear 15 5.0 <br>> zdef 1 levels 1000 <br>
> TDEF 1335 LINEAR jan1899 1mo <br>> vars 1 <br>> slp 0 99 pressure <br>> endvars <br>> Additional data information:<br>> 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.<br>
> <br>> ZABLONE OWITI<br>> GRADUATE STUDENT<br>> Nanjing University of Information, Science and Technology<br>> College of International Education<br>> <br>> ====================================================<br>
> <br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20100425/ccb08070/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20100425/ccb08070/attachment-0001.html</a> <br>
> <br>> ------------------------------<br>> <br>> Message: 6<br>> Date: Mon, 26 Apr 2010 06:46:39 +0000<br>> From: Geeta Geeta <<a href="mailto:geeta124@hotmail.com" target="_blank">geeta124@hotmail.com</a>><br>
> Subject: [gradsusr] error in gr2stn script<br>> To: grads grads <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>> Message-ID: <SNT117-W28F33B5C383D2CB799A3548A040@phx.gbl><br>
> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> <br>> Dear All, <br>> 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. <br>
> <br>> ga-> reinit<br>> No hardcopy metafile open<br>> All files closed; all defined objects released<br>> All GrADS attributes have been reinitialized<br>> ga-> open 12mar2008<br>> Scanning description file: 12mar2008<br>
> Data file gdas2008031200f00 is open as file 1<br>> LON set to 0 360 <br>> LAT set to -90 90 <br>> LEV set to 1000 1000 <br>> Time values set: 2008:3:12:0 2008:3:12:0 <br>> ga-> set t 1<br>> Time values set: 2008:3:12:0 2008:3:12:0 <br>
> ga-> set lev 1000<br>> LEV set to 1000 1000 <br>> ga-> d gr2stn(TMP2M,77.5,12.5)<br>> Error from GR2STN: 2nd argument must be stns<br>> Operation Error: Error from gr2stn function<br>> Error ocurred at column 1<br>
> DISPLAY error: Invalid expression <br>> Expression = gr2stn(tmp2m,77.5,12.5)<br>> <br>> <br>> <br>> If I use this function with the WRF output, I get the following message. <br>> ga-> reinit<br>
> No hardcopy metafile open<br>> All files closed; all defined objects released<br>> All GrADS attributes have been reinitialized<br>> ga-> open test<br>> Scanning description file: test<br>> Data file test.dat is open as file 1<br>
> LON set to 60.2713 89.7287 <br>> LAT set to 0.57003 23.9053 <br>> LEV set to 1000 1000 <br>> Time values set: 2010:4:24:0 2010:4:24:0 <br>> ga-> set t 1<br>> Time values set: 2010:4:24:0 2010:4:24:0 <br>
> ga-> set lev 1000<br>> LEV set to 1000 1000 <br>> ga-> d gr2stn(T2, 77.5, 12.5)<br>> Error from GR2STN: Grid must have linear scaling<br>> Operation Error: Error from gr2stn function<br>> Error ocurred at column 1<br>
> DISPLAY error: Invalid expression <br>> Expression = gr2stn(t2,77.5,12.5)<br>> ga-> <br>> <br>> I Think that since WRF output does not have fixed dx and dy so the error. <br>> <br>> Kindly suggest. <br>
> Geeta<br>> <br>> <br>> _________________________________________________________________<br>> Catch the latest in the world of fashion<br>> <a href="http://lifestyle.in.msn.com/" target="_blank">http://lifestyle.in.msn.com/</a><br>
> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20100426/690c0685/attachment.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20100426/690c0685/attachment.html</a> <br>
> <br>> ------------------------------<br>> <br>> _______________________________________________<br>> gradsusr mailing list<br>> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>> <br>> <br>> End of gradsusr Digest, Vol 2, Issue 34<br>> ***************************************<br>
<br><hr>Invest your money wisely post Budget <a href="http://news.in.msn.com/moneyspecial/" target="_blank">Sign up now.</a></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div>
<pre>***************************************************
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.
***************************************************