<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'>Hi Jeff<div>Thank you. It worked out. Now the map is different. I think its correct too but need to verify further and check.<br><br><i>Regards<br>Ramchandra Karki<br></i><br><div>> From: gradsusr-request@gradsusr.org<br>> Subject: gradsusr Digest, Vol 69, Issue 14<br>> To: gradsusr@gradsusr.org<br>> Date: Sat, 7 Nov 2015 15:13:00 -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. monthly to yearly script (Ramchandra Karki)<br>> 2. Re: monthly to yearly script (Jeff Duda)<br>> <br>> <br>> ----------------------------------------------------------------------<br>> <br>> Message: 1<br>> Date: Sat, 7 Nov 2015 23:45:18 +0530<br>> From: Ramchandra Karki <rammetro@hotmail.com><br>> Subject: [gradsusr] monthly to yearly script<br>> To: "gradsusr@gradsusr.org" <gradsusr@gradsusr.org><br>> Message-ID: <SNT148-W790CC6E477256724A9E228D0170@phx.gbl><br>> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> Hi,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.<br>> <br>> <br>> indata = "TMAX20KM.CTL"datname = "MAXT"'open 'indata'''set lat 26 31''set lon 80 89''set lev 1''set gxout shaded'y= 1while (y <=29) yr = 70+y b=y*12 a=12*y-11 'define yavg=ave(tt, t='a', t='b', 1)' 'set mpdset hires' 'set display color white' 'clear' 'd yavg' 'draw title Annual 'datname': 19'yr 'printim 'yr'.png' 'close 1'y=y+1endwhile'quit'<br>> <br>> Regards<br>> Ramchandra Karki<br>> <br>> <br>> > 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>> <br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20151107/75ea8337/attachment-0001.html <br>> <br>> ------------------------------<br>> <br>> Message: 2<br>> Date: Sat, 7 Nov 2015 14:13:36 -0600<br>> From: Jeff Duda <jeffduda319@gmail.com><br>> Subject: Re: [gradsusr] monthly to yearly script<br>> To: GrADS Users Forum <gradsusr@gradsusr.org><br>> Message-ID:<br>> <CAAig09CTaDMM08inRi=invPwoike-_xpaTPswPzO-QYn50sm3g@mail.gmail.com><br>> Content-Type: text/plain; charset="utf-8"<br>> <br>> the 'close' command does not operate as it does in Matplotlib or MATLAB in<br>> closing figures. In Grads, 'close' closes the data file. Since your open<br>> statement is outside the loop, that probably explains why you are only<br>> getting one map. The output from the script should be cluing you into that<br>> fact by saying something like "CAN'T DISPLAY - NO FILES OPEN"<br>> <br>> Jeff Duda<br>> <br>> On Sat, Nov 7, 2015 at 12:15 PM, Ramchandra Karki <rammetro@hotmail.com><br>> wrote:<br>> <br>> > Hi,<br>> > I have a monthly data set in one binary file and I want to get annual map<br>> > but only one map is plotted with this script. can any body help me to<br>> > correct this script. I want to read those data in ascii too.<br>> ><br>> ><br>> ><br>> > indata = "TMAX20KM.CTL"<br>> > datname = "MAXT"<br>> > 'open 'indata''<br>> > 'set lat 26 31'<br>> > 'set lon 80 89'<br>> > 'set lev 1'<br>> > 'set gxout shaded'<br>> > y= 1<br>> > while (y <=29)<br>> > yr = 70+y<br>> > b=y*12<br>> > a=12*y-11<br>> > 'define yavg=ave(tt, t='a', t='b', 1)'<br>> > 'set mpdset hires'<br>> > 'set display color white'<br>> > 'clear'<br>> > 'd yavg'<br>> > 'draw title Annual 'datname': 19'yr<br>> > 'printim 'yr'.png'<br>> > 'close 1'<br>> > y=y+1<br>> > endwhile<br>> > 'quit'<br>> ><br>> ><br>> ><br>> ><br>> > *RegardsRamchandra Karki*<br>> ><br>> > > 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><br>> > wrote:<br>> > ><br>> > > > Grads variable names have requirements. I believe those are a limit of<br>> > 16<br>> > > > characters, and there can be only alphanumeric characters in the name.<br>> > In a<br>> > > > control file, you can rename a netcdf array to meet the grads variable<br>> > 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<br>> > complete<br>> > > > details.<br>> > > ><br>> > > > Jeff Duda<br>> > > ><br>> > > > On Fri, Nov 6, 2015 at 5:29 AM, Boddapati Anil <<br>> > 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<br>> > 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:<br>> > 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.<br>> > 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,z<br>> comp-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 <<br>> > chris.gilroy@gmail.com><br>> > > >> wrote:<br>> > > >><br>> > > >> > Alright, so I'm unsure if maskout is really the option I need for<br>> > 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<br>> > 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<br>> > also<br>> > > >> change the results of negto10 and turn all numbers larger than 41 to<br>> > 40.<br>> > > >> I<br>> > > >> can't seem to figure that part out. I've tried all <br>> sorts of variants<br>> > of:<br>> > > >> ><br>> > > >> > 'define mor4140 = skip(const(maskout(negto10,negto10+41), 40, -u),<br>> > 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 <br>> 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.<br>> > 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:<br>> > 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>> ><br>> > _______________________________________________<br>> > g<br>> radsusr 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>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20151107/b032dc28/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 69, Issue 14<br>> ****************************************<br>> <br></div></div> </div></body>
</html>