<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'><div>Hi,</div>I have a monthly data set in one binary file and I want to get annual map but only one map is plotted with this script. can any body help me to correct this script. I want to read those data in ascii too.<div><br></div><div><br></div><div><br><div>indata = "TMAX20KM.CTL"</div><div>datname = "MAXT"</div><div>'open 'indata''</div><div>'set lat 26 31'</div><div>'set lon 80 89'</div><div>'set lev 1'</div><div>'set gxout shaded'</div><div>y= 1</div><div>while (y <=29)</div><div> yr = 70+y</div><div> b=y*12</div><div> a=12*y-11</div><div> 'define yavg=ave(tt, t='a', t='b', 1)'</div><div> 'set mpdset hires'</div><div> 'set display color white'</div><div> 'clear'</div><div> 'd yavg'</div><div> 'draw title Annual 'datname': 19'yr</div><div> 'printim 'yr'.png'</div><div> 'close 1'</div><div>y=y+1</div><div>endwhile</div><div>'quit'</div><i><div><i><br></i></div><div><i><br></i></div>Regards<br>Ramchandra Karki<br></i><br><br><div>> From: gradsusr-request@gradsusr.org<br>> Subject: gradsusr Digest, Vol 69, Issue 13<br>> To: gradsusr@gradsusr.org<br>> Date: Sat, 7 Nov 2015 12:00:01 -0500<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: netcdf file opening problem (Boddapati Anil)<br>> 2. Re: maskout problems with positive and negative numbers.<br>> (leila amini)<br>> 3. g2ctl/gribmap and 6hr/12hr GFS (Christopher Gilroy)<br>> <br>> <br>> ----------------------------------------------------------------------<br>> <br>> Message: 1<br>> Date: Sat, 7 Nov 2015 16:45:54 +0530<br>> From: Boddapati Anil <boddapati.anil2@gmail.com><br>> Subject: Re: [gradsusr] netcdf file opening problem<br>> To: GrADS Users Forum <gradsusr@gradsusr.org><br>> Message-ID:<br>> <CAJdYfFHSGiGSFYLmZtAERp++hM4gTgy0r1d5rLqha2C=_qL=pw@mail.gmail.com><br>> Content-Type: text/plain; charset="utf-8"<br>> <br>> Thank you very much for your support in my job<br>> <br>> On Fri, Nov 6, 2015 at 10:12 PM, Jeff Duda <jeffduda319@gmail.com> wrote:<br>> <br>> > Grads variable names have requirements. I believe those are a limit of 16<br>> > characters, and there can be only alphanumeric characters in the name. In a<br>> > control file, you can rename a netcdf array to meet the grads variable name<br>> > requirements using the => syntax. For example:<br>> ><br>> > Atmospheric_Water_Vapor_Mean=>meanwvap<br>> ><br>> > See http://www.iges.org/grads/gadoc/descriptorfile.html#VARS for complete<br>> > details.<br>> ><br>> > Jeff Duda<br>> ><br>> > On Fri, Nov 6, 2015 at 5:29 AM, Boddapati Anil <boddapati.anil2@gmail.com><br>> > wrote:<br>> ><br>> >> Dear Users,<br>> >><br>> >> I'm using one MODIS file it was converted from hdf to netcdf<br>> >><br>> >> for that i wrote one discription file ,using that it was opening but<br>> >> cann't display the data<br>> >><br>> >><br>> >><br>> >> *Discription File:*dset ^MOD08_D3.A2011003.051.nc<br>> >> dtype netcdf<br>> >> undef -999<br>> >> TITLE MODIS/Terra Aerosol Cloud Water Vapor Ozone Daily L3 Global 1Deg CMG<br>> >> xdef 360 linear -179.5 1<br>> >> ydef 180 linear -89.5 1<br>> >> tdef 1 linear 03Jan2011 1dy<br>> >> zdef 1 linear 1 1<br>> >> vars 1<br>> >> Atmospheric_Water_Vapor_Mean 0 99999 y,x Precipitable Water Vapor (IR<br>> >> Retrieval) Total Column: Mean<br>> >> endvars<br>> >><br>> >> *ERROR*<br>> >> ga-> d atmospheric_wat<br>> >> Error: Variable atmospheric_wat not in netcdf file<br>> >> Data Request Error: Error for variable 'atmospheric_wat'<br>> >> Error ocurred at column 1<br>> >> DISPLAY error: Invalid expression<br>> >> Expression = atmospheric_wat<br>> >><br>> >> --<br>> >><br>> >><br>> >> *Boddapati Anil(PF-822)CSIR Fourth Paradigm Institute*<br>> >><br>> >> *Wind Tunnel Road,Belur campous*<br>> >><br>> >> *Bengalore-37*<br>> >><br>> >><br>> >> *contact:+91-8867766688 +91-7022665333*<br>> >><br>> >><br>> >> _______________________________________________<br>> >> gradsusr mailing list<br>> >> gradsusr@gradsusr.org<br>> >> http://gradsusr.org/mailman/listinfo/gradsusr<br>> >><br>> >><br>> ><br>> ><br>> > --<br>> > Jeff Duda<br>> > Graduate research assistant<br>> > University of Oklahoma School of Meteorology<br>> > Center for Analysis and Prediction of Storms<br>> ><br>> > _______________________________________________<br>> > gradsusr mailing list<br>> > gradsusr@gradsusr.org<br>> > http://gradsusr.org/mailman/listinfo/gradsusr<br>> ><br>> ><br>> <br>> <br>> -- <br>> <br>> <br>> *Boddapati Anil(PF-822)CSIR Fourth Paradigm Institute*<br>> <br>> *Wind Tunnel Road,Belur campous*<br>> <br>> *Bengalore-37*<br>> <br>> <br>> *contact:+91-8867766688 +91-7022665333*<br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20151107/67269e98/attachment-0001.html <br>> <br>> ------------------------------<br>> <br>> Message: 2<br>> Date: Sat, 7 Nov 2015 14:49:20 +0330<br>> From: leila amin<br>> i <labi475@gmail.com><br>> Subject: Re: [gradsusr] maskout problems with positive and negative<br>> numbers.<br>> To: GrADS Users Forum <gradsusr@gradsusr.org><br>> Message-ID:<br>> <CAChC=xRM57ij5bvio1xHjUeP-aaBe+TXCoRXEjnVEW6-XJBmtA@mail.gmail.com><br>> Content-Type: text/plain; charset=UTF-8<br>> <br>> hi<br>> please help me ,how can i draw sum of rain in 24 hours on map<br>> what is my script?<br>> please send for me command<br>> thanks<br>> leila<br>> <br>> On 10/29/15, Jeff Duda <jeffduda319@gmail.com> wrote:<br>> > It would probably help if you left the skip function for the very end. I<br>> > don't see a need for using the skip command to merely define a new<br>> > variable. skip() is a display type command.<br>> ><br>> > define negto10 = const(maskout(zcomp,zcomp-1),10,-u)<br>> ><br>> > define mor4140 = const(maskout(negto10,41-negto10),40,-u)<br>> ><br>> > d skip(mor4140,8,8)<br>> ><br>> > Jeff Duda<br>> ><br>> > On Thu, Oct 29, 2015 at 6:43 AM, Mikovitz, Colleen (LARC-E302)[SCIENCE<br>> > SYSTEMS AND APPLICATIONS, INC] <j.c.mikovitz@nasa.gov> wrote:<br>> ><br>> >> You can make a negative version of your variable and do a similar<br>> >> process. With a field I have on, limiting values between 200 and 300:<br>> >><br>> >> 'negto200=const(maskout(toaup,toaup-200),200,-u)'<br>> >> 'revnegto200=-negto200'<br>> >> 'keep=revnegto200+300'<br>> >> 'backtogood=-1*(maskout(keep,keep)-300)'<br>> >> 'final=const(back<br>> togood,300,-u)'<br>> >> 'd final?<br>> >><br>> >><br>> >> On Oct 28, 2015, at 10:27 PM, Christopher Gilroy <chris.gilroy@gmail.com><br>> >> wrote:<br>> >><br>> >> > Alright, so I'm unsure if maskout is really the option I need for both<br>> >> of these, but is there anyway to take any value less than 1 (even<br>> >> negatives) and replace them with 10? I seem to be able to do that with:<br>> >> ><br>> >> > 'define negto10= skip(const(maskout(zcomp,zcomp-0), 10, -u), 8, 8)'<br>> >> ><br>> >> > That gives me all positive numbers, great so far. Now, I'd like to also<br>> >> change the results of negto10 and turn all numbers larger than 41 to 40.<br>> >> I<br>> >> can't seem to figure that part out. I've tried all sorts of variants of:<br>> >> ><br>> >> > 'define mor4140 = skip(const(maskout(negto10,negto10+41), 40, -u), 8,<br>> >> > 8)'<br>> >> ><br>> >> > all with no luck. Any thoughts?<br>> >> > _______________________________________________<br>> >> > gradsusr mailing list<br>> >> > gradsusr@gradsusr.org<br>> >> > http://gradsusr.org/mailman/listinfo/gradsusr<br>> >><br>> >><br>> >> _______________________________________________<br>> >> gradsusr mailing list<br>> >> gradsusr@gradsusr.org<br>> >> http://gradsusr.org/mailman/listinfo/gradsusr<br>> >><br>> ><br>> ><br>> ><br>> > --<br>> > Jeff Duda<br>> > Graduate research assistant<br>> > University of Oklahoma School of Meteorology<br>> > Center for Analysis and Prediction of Storms<br>> ><br>> <br>> <br>> <br>> ------------------------------<br>> <br>> Message: 3<br>> Date: Sat, 7 Nov 2015 09:52:47 -0500<br>> From: Christ<br>> opher Gilroy <chris.gilroy@gmail.com><br>> Subject: [gradsusr] g2ctl/gribmap and 6hr/12hr GFS<br>> To: GrADS Users Forum <gradsusr@gradsusr.org><br>> Message-ID:<br>> <CAFCiPemZ48O9Y5zqGNB5y8GRrpuRCnLyftECffCQT9BHAWdVfQ@mail.gmail.com><br>> Content-Type: text/plain; charset="utf-8"<br>> <br>> I've always been curious for those who run their plots using scripting<br>> approaches and download in 3 or 6-hour increments, how do you handle past<br>> hour 240?<br>> <br>> For example, downloading 6hr files, starting at hr 0, t41=240 and then<br>> t43=252 and I'm curious if there's a way to make t42=252, t43=264, etc. I'm<br>> just looking for a way to make every t = a working timestep.<br>> <br>> So as an example, if you made a looping grads script for all 384 hours,<br>> using 6h files the only way I knew of to "get around" the 2 t-step<br>> increments after 240 was doing something like:<br>> <br>> i=1<br>> while (i <= 65)<br>> <br>> if (i <= 41)<br>> trecord=i<br>> endif<br>> if (i=42)<br>> trecord = 43<br>> endif<br>> if (i=43)<br>> trecord = 45<br>> endif<br>> <br>> 'd tmp2m'<br>> <br>> The only reason I cared to do that even was because otherwise every other<br>> image past hr240 would be blank since there is no t=42, t=44, etc.<br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20151107/61568f88/attachment-0001.html <br>> <br>> ------------------------------<br>> <br>> _______________________________________________<br>> gradsusr mailing list<br>> gradsusr@gra<br>> dsusr.org<br>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> <br>> <br>> End of gradsusr Digest, Vol 69, Issue 13<br>> ****************************************<br>> <br></div></div> </div></body>
</html>