geeta,<div><br></div><div>Yes, you can write this output (result of &#39;d gr2stn....&#39; using gxout print) to a file, but you did not say what type of output.  A simple way to output to a text file, for example:</div><div>
<div>outfile=&#39;c:/temp/tmpprs.dat&#39;</div><div>&#39;set t 1&#39;;&#39;set x 1&#39;;&#39;set y 1&#39;</div><div>&#39;set gxout print&#39;</div><div>&#39;set prnopts %0.2f&#39;</div><div>ts=1</div><div>while(ts&lt;=10)</div>
<div>  &#39;d gr2stn(TMPprs(t=&#39;ts&#39;),77.5,12.5)&#39;</div><div>  tmp=sublin(result,3)</div><div>  rc=write(outfile,ts&#39;,&#39;tmp)</div><div>  ts=ts+1</div><div>endwhile</div><div>rc=close(outfile)</div><div><br>
</div><div>Above should produce a file (.../tmpprs.dat, in this case) with an output similar to this, where first value is the timestep, and second value is the interpolated value for that timestep:</div><div><div><div>1,308.80</div>
<div>2,311.90</div><div>3,311.30</div><div>4,307.60</div><div>5,303.90</div><div>6,301.30</div><div>7,300.50</div><div>8,302.90</div><div>9,308.20</div><div>10,311.50</div><div><br></div></div><div><br></div></div></div><div>
<div>You may also find this OpenGrADS Cookbook reference useful, depending on the type of output you want:</div><div><a href="http://cookbooks.opengrads.org/index.php?title=Recipe-002:_Saving_GrADS_variable_data_to_a_text_file">http://cookbooks.opengrads.org/index.php?title=Recipe-002:_Saving_GrADS_variable_data_to_a_text_file</a></div>
<div><br></div><div>Using a Google or other web search, you should be able to come across other useful help by specifying something like &quot;ascii output site:<a href="http://gradsusr.org">gradsusr.org</a>&quot; to search the User archives.</div>
<div><br></div><div>Stephen McMillan<br><br><div class="gmail_quote">On Tue, Apr 27, 2010 at 12:09 PM, Geeta Geeta <span dir="ltr">&lt;<a href="mailto:geeta124@hotmail.com">geeta124@hotmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">




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