<div dir="ltr"><div>the 'close' command does not operate as it does in Matplotlib or MATLAB in closing figures. In Grads, 'close' closes the data file. Since your open statement is outside the loop, that probably explains why you are only getting one map. The output from the script should be cluing you into that fact by saying something like "CAN'T DISPLAY - NO FILES OPEN"<br><br></div>Jeff Duda<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Nov 7, 2015 at 12:15 PM, Ramchandra Karki <span dir="ltr"><<a href="mailto:rammetro@hotmail.com" target="_blank">rammetro@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><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: <a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>> Subject: gradsusr Digest, Vol 69, Issue 13<br>> To: <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>> Date: Sat, 7 Nov 2015 12:00:01 -0500<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. 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 <<a href="mailto:boddapati.anil2@gmail.com" target="_blank">boddapati.anil2@gmail.com</a>><br>> Subject: Re: [gradsusr] netcdf file opening problem<br>> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>> Message-ID:<br>> <CAJdYfFHSGiGSFYLmZtAERp++hM4gTgy0r1d5rLqha2C=_qL=<a href="mailto:pw@mail.gmail.com" target="_blank">pw@mail.gmail.com</a>><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 <<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>> 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 <a href="http://www.iges.org/grads/gadoc/descriptorfile.html#VARS" target="_blank">http://www.iges.org/grads/gadoc/descriptorfile.html#VARS</a> for complete<br>> > details.<br>> ><br>> > Jeff Duda<br>> ><br>> > On Fri, Nov 6, 2015 at 5:29 AM, Boddapati Anil <<a href="mailto:boddapati.anil2@gmail.com" target="_blank">boddapati.anil2@gmail.com</a>><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 ^<a href="http://MOD08_D3.A2011003.051.nc" target="_blank">MOD08_D3.A2011003.051.nc</a><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>> >> <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>> > --<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>> > <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>> -- <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: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20151107/67269e98/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20151107/67269e98/attachment-0001.html</a> <br>> <br>> ------------------------------<br>> <br>> Message: 2<br>> Date: Sat, 7 Nov 2015 14:49:20 +0330<br>> From: leila amin<br>> i <<a href="mailto:labi475@gmail.com" target="_blank">labi475@gmail.com</a>><br>> Subject: Re: [gradsusr] maskout problems with positive and negative<br>> numbers.<br>> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>> Message-ID:<br>> <CAChC=<a href="mailto:xRM57ij5bvio1xHjUeP-aaBe%2BTXCoRXEjnVEW6-XJBmtA@mail.gmail.com" target="_blank">xRM57ij5bvio1xHjUeP-aaBe+TXCoRXEjnVEW6-XJBmtA@mail.gmail.com</a>><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 <<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>> 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] <<a href="mailto:j.c.mikovitz@nasa.gov" target="_blank">j.c.mikovitz@nasa.gov</a>> 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 <<a href="mailto:chris.gilroy@gmail.com" target="_blank">chris.gilroy@gmail.com</a>><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>> >> > <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>> >> 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>> > --<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 <<a href="mailto:chris.gilroy@gmail.com" target="_blank">chris.gilroy@gmail.com</a>><br>> Subject: [gradsusr] g2ctl/gribmap and 6hr/12hr GFS<br>> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>> Message-ID:<br>> <<a href="mailto:CAFCiPemZ48O9Y5zqGNB5y8GRrpuRCnLyftECffCQT9BHAWdVfQ@mail.gmail.com" target="_blank">CAFCiPemZ48O9Y5zqGNB5y8GRrpuRCnLyftECffCQT9BHAWdVfQ@mail.gmail.com</a>><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: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20151107/61568f88/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20151107/61568f88/attachment-0001.html</a> <br>> <br>> ------------------------------<br>> <br>> _______________________________________________<br>> gradsusr mailing list<br>> gradsusr@gra<br>> <a href="http://dsusr.org" target="_blank">dsusr.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 69, Issue 13<br>> ****************************************<br>> <br></div></div> </div></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" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature">Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br></div>
</div>