<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
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>> From: gradsusr-request@gradsusr.org<br>> Subject: gradsusr Digest, Vol 2, Issue 40<br>> To: gradsusr@gradsusr.org<br>> Date: Tue, 27 Apr 2010 11:01:22 -0400<br>> <br>> Send gradsusr mailing list submissions to<br>> gradsusr@gradsusr.org<br>> <br>> To subscribe or unsubscribe via the World Wide Web, visit<br>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> or, via email, send a message with subject or body 'help' to<br>> gradsusr-request@gradsusr.org<br>> <br>> You can reach the person managing the list at<br>> gradsusr-owner@gradsusr.org<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. Re: Gribmap (James Burn)<br>> 2. Re: error with gr2stn function (Stephen McMillan)<br>> <br>> <br>> ----------------------------------------------------------------------<br>> <br>> Message: 1<br>> Date: Tue, 27 Apr 2010 09:38:50 +0100<br>> From: James Burn <james@weatherbadger.com><br>> Subject: Re: [gradsusr] Gribmap<br>> To: gradsusr@gradsusr.org<br>> Message-ID: <201004270938.51071.james@weatherbadger.com><br>> Content-Type: text/plain; charset="iso-8859-15"<br>> <br>> Andrew<br>> In my cron-run shell scripts I change directory to the location of the data <br>> first before running gridmap.<br>> <br>> "cd $HOME/grads_etc/data " <br>> <br>> I've had it like this for some time, and can't remember why I did it like <br>> this, but perhaps it was to get around the problem you've got?<br>> <br>> Good luck<br>> James<br>> <br>> On Tuesday 27 April 2010 06:18:58 Andrew Revering wrote:<br>> > I'm new to Linux and I'm pulling my hair out trying to figure out why<br>> > gribmap will not run when executed as a cron job.<br>> ><br>> ><br>> ><br>> > In crontab I execute a .sh script that has something that looks like:<br>> ><br>> ><br>> ><br>> > #!/bin/bash<br>> ><br>> ><br>> ><br>> > /wrf/wrfems/util/grads/bin/grib2ctl.pl -verf<br>> > /wrf/wrfems/runs/float/emsprd/grib/arw_0000.grib ><br>> > /wrf/wrfems/runs/float/emsprd/grib/arw_0000.ctl<br>> ><br>> ><br>> ><br>> > /wrf/wrfems/util/grads/bin/gribmap -I<br>> > /wrf/wrfems/runs/float/emsprd/grib/arw_0000.ctl<br>> ><br>> ><br>> ><br>> ><br>> ><br>> ><br>> ><br>> > The last line is the problem child. Note that the script is running, and<br>> > the first line in the script (converting GRIB to a CTL file) works, but<br>> > nothing happens when trying to execute the gribmap executable.<br>> ><br>> ><br>> ><br>> > Also I should note that if I type that same line in the terminal it DOES<br>> > work, just not in crontab.<br>> ><br>> ><br>> ><br>> > Any ideas?<br>> <br>> -- <br>> <br>> Reply to : james@weatherbadger.com<br>> Visit http://www.weatherbadger.com for some weather forecasts.<br>> <br>> <br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100427/892ffb92/attachment-0001.html <br>> <br>> ------------------------------<br>> <br>> Message: 2<br>> Date: Tue, 27 Apr 2010 11:01:13 -0400<br>> From: Stephen McMillan <smcmillan@planalytics.com><br>> Subject: Re: [gradsusr] error with gr2stn function<br>> To: GrADS Users Forum <gradsusr@gradsusr.org><br>> Message-ID:<br>> <s2v1f726e2f1004270801u5011e3cbu4b58b50581a43120@mail.gmail.com><br>> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> Geeta,<br>> Try setting or re-setting your prnopts, for example:<br>> 'set gxout print'<br>> 'set prnopts %0.2f'<br>> 'd gr2stn(TMPprs,77.5,12.5)'<br>> <br>> If you still don't get a result...<br>> If you display a map with x and y varying, t = 1, and lev = 1000, do you get<br>> a sensible image of the variable tmpprs? Do you also get a reasonable value<br>> if you don't interpolate?, for example:<br>> 'd TMPprs(lon=77.5,lat=12.5)' ...should give you the same value as 'd<br>> gr2stn(TMPprs,77.5,12.5)' if both are at the same data gridpoint.<br>> <br>> Good luck!<br>> Stephen McMillan<br>> <br>> On Tue, Apr 27, 2010 at 1:36 AM, Geeta Geeta <geeta124@hotmail.com> wrote:<br>> <br>> > Dear All,<br>> > I have tried to fix the x and y as suggested. But I have not got the result<br>> > . 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: gradsusr-request@gradsusr.org<br>> > > Subject: gradsusr Digest, Vol 2, Issue 34<br>> > > To: gradsusr@gradsusr.org<br>> > > Date: Mon, 26 Apr 2010 02:46:40 -0400<br>> > ><br>> > > Send gradsusr mailing list submissions to<br>> > > gradsusr@gradsusr.org<br>> > ><br>> > > To subscribe or unsubscribe via the World Wide Web, visit<br>> > > http://gradsusr.org/mailman/listinfo/gradsusr<br>> > > or, via email, send a message with subject or body 'help' to<br>> > > gradsusr-request@gradsusr.org<br>> > ><br>> > > You can reach the person managing the list at<br>> > > gradsusr-owner@gradsusr.org<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) (zowiti@ncst.go.ke)<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 <geeta124@hotmail.com><br>> > > Subject: [gradsusr] gr2stn function<br>> > > To: <gradsusr@gradsusr.org><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>> > > http://news.in.msn.com/gallery/archive.aspx<br>> > > -------------- next part --------------<br>> > > An HTML attachment was scrubbed...<br>> > > URL:<br>> > http://gradsusr.org/pipermail/gradsusr/attachments/20100424/61b8f017/attachment-0001.html<br>> > ><br>> > > ------------------------------<br>> > ><br>> > > Message: 2<br>> > > Date: Sat, 24 Apr 2010 11:02:31 -0300<br>> > > From: Carlos Batista <krlosbatist@gmail.com><br>> > > Subject: [gradsusr] scorr function<br>> > > To: gradsusr@gradsusr.org<br>> > > Message-ID:<br>> > > <y2gdd4d468f1004240702zaf186fd7ufa21b4b504df48a1@mail.gmail.com><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<br>> > correlations<br>> > > between this data, please, useding the scorr function and other functions<br>> > ?<br>> > > Any suggestions<br>> > > Krlos<br>> > > e-mail; krlosbatist@gmail.com<br>> > > -------------- next part --------------<br>> > > An HTML attachment was scrubbed...<br>> > > URL:<br>> > http://gradsusr.org/pipermail/gradsusr/attachments/20100424/c4d8f6e8/attachment-0001.html<br>> > ><br>> > > ------------------------------<br>> > ><br>> > > Message: 3<br>> > > Date: Sat, 24 Apr 2010 10:42:04 -0400<br>> > > From: Stephen McMillan <smcmillan@planalytics.com><br>> > > Subject: Re: [gradsusr] gr2stn function<br>> > > To: GrADS Users Forum <gradsusr@gradsusr.org><br>> > > Message-ID:<br>> > > <s2u1f726e2f1004240742j1329b506kcf3e48202ec2924c@mail.gmail.com><br>> > > Content-Type: text/plain; charset="iso-8859-1"<br>> > ><br>> > > Dear geeta,<br>> > > Please see this post on the the gradsusr archives:<br>> > > http://gradsusr.org/pipermail/gradsusr/2010-March/000077.html<br>> > > You can search for additional posts at<br>> > > http://gradsusr.org/pipermail/gradsusr/<br>> > > GrADS documentation at:<br>> > http://sidata.org/grads/gadoc/gradfuncgr2stn.html<br>> > ><br>> > > Good luck--<br>> > > Stephen McMillan<br>> > ><br>> > > On Sat, Apr 24, 2010 at 3:45 AM, Geeta Geeta <geeta124@hotmail.com><br>> > 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.<<br>> > http://news.in.msn.com/internalsecurity/><br>> > > ><br>> > > > _______________________________________________<br>> > > > gradsusr mailing list<br>> > > > gradsusr@gradsusr.org<br>> > > > http://gradsusr.org/mailman/listinfo/gradsusr<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:<br>> > http://gradsusr.org/pipermail/gradsusr/attachments/20100424/591f216e/attachment-0001.html<br>> > ><br>> > > ------------------------------<br>> > ><br>> > > Message: 4<br>> > > Date: Sat, 24 Apr 2010 12:03:04 -0400<br>> > > From: Stephen McMillan <smcmillan@planalytics.com><br>> > > Subject: Re: [gradsusr] scorr function<br>> > > To: GrADS Users Forum <gradsusr@gradsusr.org><br>> > > Message-ID:<br>> > > <j2m1f726e2f1004240903q63515e10pcba71b2095589713@mail.gmail.com><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<br>> > the<br>> > > GrADS Documentation Index:<br>> > http://sidata.org/grads/gadoc/gradfuncscorr.html<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<br>> > 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<br>> > 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 http://sidata.org/grads/gadoc/gradfunctcorr.html<br>> > )<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 (<br>> > http://gradsusr.org/pipermail/gradsusr/)<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 <krlosbatist@gmail.com<br>> > >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<br>> > correlations<br>> > > > between this data, please, useding the scorr function and other<br>> > functions ?<br>> > > > Any suggestions<br>> > > > Krlos<br>> > > > e-mail; krlosbatist@gmail.com<br>> > > ><br>> > > ><br>> > > > _______________________________________________<br>> > > > gradsusr mailing list<br>> > > > gradsusr@gradsusr.org<br>> > > > http://gradsusr.org/mailman/listinfo/gradsusr<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:<br>> > http://gradsusr.org/pipermail/gradsusr/attachments/20100424/cadeffd4/attachment-0001.html<br>> > ><br>> > > ------------------------------<br>> > ><br>> > > Message: 5<br>> > > Date: Sun, 25 Apr 2010 06:37:16 -0700<br>> > > From: "zowiti@ncst.go.ke" <zowiti@ncst.go.ke><br>> > > Subject: [gradsusr] Unable to read Monthly Northern Hemisphere<br>> > > Sea-Level Pressure (DSS SLP Grid-binary)<br>> > > To: <gradsusr@gradsusr.org><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<br>> > Longitude/Latitude)] data in Grads via the ctl below. When I display the<br>> > slp, I get a distorted plot. I'm unable to figure out the problem in the<br>> > 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<br>> > nearest tenth millibar and are the rightmost 14 bits of each byte. The grid<br>> > of Y(I,J) is dimensioned Y(72,15). As I goes from 1 to 72, longitude goes<br>> > 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:<br>> > http://gradsusr.org/pipermail/gradsusr/attachments/20100425/ccb08070/attachment-0001.html<br>> > ><br>> > > ------------------------------<br>> > ><br>> > > Message: 6<br>> > > Date: Mon, 26 Apr 2010 06:46:39 +0000<br>> > > From: Geeta Geeta <geeta124@hotmail.com><br>> > > Subject: [gradsusr] error in gr2stn script<br>> > > To: grads grads <gradsusr@gradsusr.org><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<br>> > posting the following message from grads. My gridded data file is the gfs<br>> > 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>> > ><br>> > > Kindly suggest.<br>> > > Geeta<br>> > ><br>> > ><br>> > > _________________________________________________________________<br>> > > Catch the latest in the world of fashion<br>> > > http://lifestyle.in.msn.com/<br>> > > -------------- next part --------------<br>> > > An HTML attachment was scrubbed...<br>> > > URL:<br>> > http://gradsusr.org/pipermail/gradsusr/attachments/20100426/690c0685/attachment.html<br>> > ><br>> > > ------------------------------<br>> > ><br>> > > _______________________________________________<br>> > > gradsusr mailing list<br>> > > gradsusr@gradsusr.org<br>> > > http://gradsusr.org/mailman/listinfo/gradsusr<br>> > ><br>> > ><br>> > > End of gradsusr Digest, Vol 2, Issue 34<br>> > > ***************************************<br>> ><br>> > ------------------------------<br>> > Invest your money wisely post Budget Sign up now.<http://news.in.msn.com/moneyspecial/><br>> ><br>> > _______________________________________________<br>> > gradsusr mailing list<br>> > gradsusr@gradsusr.org<br>> > http://gradsusr.org/mailman/listinfo/gradsusr<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: http://gradsusr.org/pipermail/gradsusr/attachments/20100427/e3848da1/attachment.html <br>> <br>> ------------------------------<br>> <br>> _______________________________________________<br>> gradsusr mailing list<br>> gradsusr@gradsusr.org<br>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> <br>> <br>> End of gradsusr Digest, Vol 2, Issue 40<br>> ***************************************<br> <br /><hr />The battle for the FIH Hockey World Cup <a href='http://specials.msn.co.in/sp10/hockey/index.aspx' target='_new'>Drag n' drop</a></body>
</html>