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