<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Dear Jeff,<BR> <BR>I tried different commeds in C:\grads, when I saw ga-> I wrote these commends: <BR><br>1) open E:\WRF\first\test.ctl<BR>2) open <E>:\<WRF><first>\<test.ctl><BR>3) open <E>:\<WRF\first>\<test.ctl><BR><br>and some others, but unfortunately no one worked.<BR><br>Could you please write me the commend completely ?<BR> <BR>Regards,<BR> <BR>Tahereh<BR> <BR> <BR><br><br> <BR><div>> From: gradsusr-request@gradsusr.org<br>> Subject: gradsusr Digest, Vol 34, Issue 55<br>> To: gradsusr@gradsusr.org<br>> Date: Tue, 18 Dec 2012 13:16:45 -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: Error in using skew.gs (Jeff Duda)<br>> 2. Re: How to open a .ctl file in a different path (Jeff Duda)<br>> <br>> <br>> ----------------------------------------------------------------------<br>> <br>> Message: 1<br>> Date: Tue, 18 Dec 2012 12:15:58 -0600<br>> From: Jeff Duda <jeffduda319@gmail.com><br>> Subject: Re: [gradsusr] Error in using skew.gs<br>> To: GrADS Users Forum <gradsusr@gradsusr.org><br>> Message-ID:<br>> <CAAig09B36=1fktE9zZ=hTwFP-EojmUBM=Vug6k9nVmBhbL-afw@mail.gmail.com><br>> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> Amulya,<br>> Thank you for adding the new information. I can tell that the value of the<br>> script variable mix, which is given as 1573660000
, cannot be right. Check<br>> line 2424 of your script code. That is where the variable mix is<br>> computed. Insert say statements for all the component variables there.<br>> With this kind of investigation you will eventually be able to discover the<br>> source of the error. My gut tells me one of the input variables you<br>> supplied to the script in the main function call plotskew() is not what the<br>> script is expecting it to be. I'm guessing it is the dewpoint variable.<br>> But I saw your control file and it appears there should not be a problem<br>> with that. Nonetheless, insert say statements to find out what the values<br>> are, then go from there. If you need help knowing what to do next or where<br>> to look, please ask.<br>> <br>> Jeff<br>> <br>> On Tue, Dec 18, 2012 at 1:21 AM, Amulya Chevuturi <<br>> amulya.chevuturi@gmail.com> wrote:<br>> <br>> > Thank you very much for replying Mr. Jeff.<br>> ><br>> > I have added the lines as stated. In the output all of them are showing<br>> > numeric values only. So I cannot make out what the problem is. Could you or<br>> > anyone in the grads user list help me out further.<br>> ><br>> > This the output I am getting at the terminal when I am running the skew.gswith the changes you suggested.<br>> ><br>> > ga-> skew.gs<br>> > No hardcopy metafile open<br>> > All files closed; all defined objects released;<br>> > All GrADS attributes have been reinitialized<br>> > Initial Time is 00Z08JAN2007<br>> > Create gif images as well (1=yes ; 0=no)<br>> > 0<br>> > Please enter coordinated for sounding<br>> ><br>> > Latitude values between: 4.96654 and 50.0419<br>> > Enter Latitude<br>> > 33<br>> ><br>> > Longitude values between: 48.3967 and 105.603<br>> > Enter Longitude<br>> > 78<br>> ><br>> > Overlay two time periods 1=yes ; 0=no<br>> > 0<br>
> ><br>> > Please enter time to plot<br>> > Available times: 1 to 17<br>> > 5<br>> > Drawing temperature sounding.<br>> > Drawing dewpoint sounding.<br>> > 286.893<br>> > 43.385364499893<br>> > 1573660000<br>> > Non-numeric args to numeric operation<br>> > Error occurred on line 2430<br>> > In file skew.gs<br>> > Error occurred on line 1106<br>> > In file skew.gs<br>> > Error occurred on line 69<br>> > In file skew.gs<br>> ><br>> > I am attaching the skew.gs file with the changes added. Please let me<br>> > know how I can further correct the error. Thank you for the help in advance.<br>> ><br>> > With regards,<br>> > Amulya Chevuturi<br>> > Junior Research Fellow (Ph.D. Student)<br>> > Jawaharlal Nehru University<br>> > New Delhi, India<br>> ><br>> > On Tue, Dec 18, 2012 at 1:06 AM, Jeff Duda <jeffduda319@gmail.com> wrote:<br>> ><br>> >> Amulya,<br>> >> I looked through your script and I can't find the cause of the error.<br>> >> I've used this script on WRF-ARW output myself before (also postprocessed<br>> >> using ARWpost). I remember having some issues getting the script to run<br>> >> properly the first time, but I got around that somehow. The errors<br>> >> indicate that in this line<br>> >><br>> >> thetae=theta*exp((3.376/TLcl-0.00254)*mix*1.0+0.00081*mix)<br>> >><br>> >> one of theta, TLcl, or mix is not being computed correctly. I'm not sure<br>> >> which of those it is, although you could put SAY statements just above that<br>> >> line for each of the variables to see which it is. What I mean is, insert<br>> >> these lines<br>> >> say theta<br>> >> say TLcl<br>> >> say mix<br>> >> Above the line I copied above. Then see what values get spit out. One<br>> >>
; of those is either a NaN of some sort or has characters where there should<br>> >> be numbers.<br>> >><br>> >><br>> >> Jeff Duda<br>> >><br>> >> On Mon, Dec 17, 2012 at 12:56 AM, Amulya Chevuturi <<br>> >> amulya.chevuturi@gmail.com> wrote:<br>> >><br>> >>> Dear GrADS users,<br>> >>><br>> >>> I have recently tried to use the skew.gs script to plot Skew-T/Log-P<br>> >>> graphs for my WRF (ARW) output. I use ARWPost as my post processing<br>> >>> software. This is the problem that is cropping up.<br>> >>><br>> >>> ga-> skew.gs<br>> >>> No hardcopy metafile open<br>> >>> All files closed; all defined objects released;<br>> >>> All GrADS attributes have been reinitialized<br>> >>> Initial Time is 00Z08JAN2007<br>> >>> Create gif images as well (1=yes ; 0=no)<br>> >>> 0<br>> >>> Please enter coordinated for sounding<br>> >>><br>> >>> Latitude values between: 4.96654 and 50.0419<br>> >>> Enter Latitude<br>> >>> 33<br>> >>><br>> >>> Longitude values between: 48.3967 and 105.603<br>> >>> Enter Longitude<br>> >>> 78<br>> >>><br>> >>> Overlay two time periods 1=yes ; 0=no<br>> >>> 0<br>> >>><br>> >>> Please enter time to plot<br>> >>> Available times: 1 to 17<br>> >>> 5<br>> >>> Drawing temperature sounding.<br>> >>> Drawing dewpoint sounding.<br>> >>> Non-numeric args to numeric operation<br>> >>> Error occurred on line 2427<br>> >>> In file skew.gs<br>> >>> Error occurred on line 1106<br>> >>> In file skew.gs<br>> >>> Error occurred on line 69<br>> >>> In file skew.gs<br>> >
;>><br>> >>> Though I think it is plotting the temperature and dewpoint. The rest of<br>> >>> the parameters in the line 69 is not being plotted due to some errors in<br>> >>> parameters. I am attaching my ctl file and skew.gs. Please let me know<br>> >>> how to correct the problem.<br>> >>><br>> >>> Thanking you in advance.<br>> >>><br>> >>> With regards,<br>> >>> Amulya Chevuturi<br>> >>> Junior Research Fellow (Ph.D. Student)<br>> >>> Jawaharlal Nehru University<br>> >>> New Delhi, India<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>> >> gradsusr mailing list<br>> >> gradsusr@gradsusr.org<br>> >> http://gradsusr.org/mailman/listinfo/gradsusr<br>> >><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>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20121218/b950e8f2/
attachment-0001.html <br>> <br>> ------------------------------<br>> <br>> Message: 2<br>> Date: Tue, 18 Dec 2012 12:16:44 -0600<br>> From: Jeff Duda <jeffduda319@gmail.com><br>> Subject: Re: [gradsusr] How to open a .ctl file in a different path<br>> To: GrADS Users Forum <gradsusr@gradsusr.org><br>> Message-ID:<br>> <CAAig09AW042aDjXSThtWzaV-84SAjcjMrUUzPHUN3X+VYgh=_w@mail.gmail.com><br>> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> What is the exact command you entered, and from what directory?<br>> <br>> Jeff<br>> <br>> On Tue, Dec 18, 2012 at 1:52 AM, T. haghroosta <haghroosta@hotmail.com>wrote:<br>> <br>> > I couldn't open it again. I tries but It send me an open error : can't<br>> > open description file<br>> ><br>> > my files path is :<br>> ><br>> > E:\WRF\first\test.ctl<br>> ><br>> > the test.dat file is also in that path. My grads can read data from own<br>> > path (C:\grads\data\). How can I write the path of my file in open<br>> > command?<br>> ><br>> > Tahereh<br>> ><br>> > > From: gradsusr-request@gradsusr.org<br>> > > Subject: gradsusr Digest, Vol 34, Issue 40<br>> > > To: gradsusr@gradsusr.org<br>> > > Date: Mon, 17 Dec 2012 13:13:24 -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>> > > W! hen 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: How to open a .ctl file in a different path (Jeff Duda)<br>> > > 2. Re: looping program (Jeff Duda)<br>> > ><br>> > ><br>> > > ----------------------------------------------------------------------<br>> > ><br>> > > Message: 1<br>> > > Date: Mon, 17 Dec 2012 12:09:27 -0600<br>> > > From: Jeff Duda <jeffduda319@gmail.com><br>> > > Subject: Re: [gradsusr] How to open a .ctl file in a different path<br>> > > To: GrADS Users Forum <gradsusr@gradsusr.org><br>> > > Message-ID:<br>> > > <CAAig09Bbh4ir8-Jz7W5iDcYxjgEGpwmiab-nsnkVY_tfm6DWSw@mail.gmail.com><br>> > > Content-Type: text/plain; charset="iso-8859-1"<br>> > ><br>> > > The full path name of your file can be expressed as<br>> > ><br>> > > <Drive name>:\<directories containing file>/<file name><br>> > ><br>> > > So ! in your case, you'd want to reference the control file as<br>> > > E:\w hatever\whatever\whatever\control_file.ctl.<br>> > ><br>> > > Make sure the dset entry in your control file is also changed<br>> > accordingly.<br>> > > If the incorrect pathname is specified there, you won't be able to open<br>> > the<br>> > > data.<br>> > ><br>> > > Jeff Duda<br>> > ><br>> > > On Mon, Dec 17, 2012 at 6:30 AM, T. haghroosta <haghroosta@hotmail.com<br>> > >wrote:<br>> > ><br>> > > > Dear All,<br>> > > ><br>> > > > I have a .ctl file on my drive E and my grads software is installed on<br>> > my<br>> > > > drive c (windows 7). Unfortunately I don't know how can I write "open"<br>> > > > command to address that file in drive E to open it! Could
you please<br>> > learn<br>> > > > me?<br>> > > ><br>> > > > Regards,<br>> > > ><br>> > > > Tahereh<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>> > > -------------- next part --------------<br>> > > An HTML attachment was scrubbed...<br>> > > URL:<br>> > http://gradsusr.org/pipermail/gradsusr/attachments/20121217/8aa433c3/attachment-0001.html<br>> > ><br>> > > ------------------------------<br>> > ><br>> > > Message: 2<br>> > > Date: Mon, 17 Dec 2012 12:13:22 -0600<br>> > > From: Jeff Duda <jeffduda319@gmail.com><br>> > > Subject: Re: [gradsusr] looping program<br>> > > To: GrADS Users Forum <gradsusr@gradsusr.org><br>> > > Message-ID:<br>> > > <CAAig09BS5=RxxN+830O6F=uo3-FpS7OgwV5=E=Jr5K2gqNFRwA@mail.gmail.com><br>> > > Content-Type: text/plain; charset="iso-8859-1"<br>> > ><br>> > > Wendi,<br>> > > Remove LITTLE_ENDIAN from your options line. There's no such option.<br>> > ><br>> > > Try displaying ref or vel upon first opening the! file. If your file<br>> > names<br>> > > don't match the template strings yo u've set in the DSET line of your<br>> > > control file, then you will get the "Cannot contour grid" error. Other<br>> > > than that, your script looks okay. Just make sure you can display all the<br>> >
> base data first before trying to manipulate it.<br>> > ><br>> > > Jeff<br>> > ><br>> > > On Mon, Dec 17, 2012 at 5:13 AM, wendi harjupa <wendiharjupa@gmail.com<br>> > >wrote:<br>> > ><br>> > > > Dear Mr. Jeff<br>> > > ><br>> > > > I have tried the script that you sent but, still I could not do this<br>> > > > would you please help me more.<br>> > > ><br>> > > > Here some question<br>> > > > 1. Is my control file ok?<br>> > > > I want to use single control file<br>> > > > 2. base on my gs and control file below, I got error message<br>> > > > 'Cannot contour grid - all undefined values'<br>> > > > 'Cannot plot color bar: No shading information'<br>> > > > 'No hardcopy metafile open'<br>> > > > what the meaning of error message??<br>> > > ! ><br>> > > > below is control file (*try.ctl)* file<br>> > > ><br>> > > > DSET /misc/home/wendi/data1/radar/kototabang/XDR/data/*<br>> > > > %y4%m2%d2/%d2%h2.dat*<br>> > > > TITLE RADAR<br>> > > > OPTIONS TEMPLATE LITTLE_ENDIAN<br>> > > > UNDEF -1000<br>> > > > XDEF 321 LINEAR 99.688120996 0.004492369<br>> > > > YDEF 321 LINEAR -1.083708891 0.004522556<br>> > > > ZDEF 39 LINEAR 1.0 0.5<br>> > > > TDEF 300 LINEAR 13:00Z10apr2004 4mn<br>> > > > VARS 2<br>> > > > ref 39 99 Reflectivity(dBZ)<br>> > > > vel 39 99 Doppler Velocity(m/s)<br>> > > > ENDVARS<br>> > > ><br>> > > > this is my gs file (averaging for 2 days only)<br>> > > ><br>> > > > 'reinit'<br>> > > > 'open *try.ctl*'<br>> > > > 'set grads off'<br>> > > > 'set mproj latlon'<br>> > > > 'set mpdset hires'<br>> > > >
'set gxout shaded'<br>> > > > 'set lon 100 100.82'<br>> > > > 'set lat -0.7 0'<br>> > > > 'set z 5'<br>> > > &g! t; 'set parea 1 7.5 2.8 8.3'<br>> > > > 'set clevs 0 5 10 15 20 25 3 0 35 40 45 50 55 60'<br>> > > > 'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'<br>> > > > 'set t 1 300'<br>> > > > 'R=pow(10, ref/10)'<br>> > > ><br>> > > ><br>> > 'R1=ave(R,time=13Z10apr2004,time=22Z10apr2004)+ave(R,time=13Z11apr2004,time=22Z11apr2004)'<br>> > > > 'R2=R1/2'<br>> > > > 'R3=log10 (R2)'<br>> > > > 'd 10*R3'<br>> > > > 'run /misc/home/wendi/grads-2.0.1/cbarn 0.7 0.7 7.4 5.5'<br>> > > > 'draw string 7.2 8.1 [dBz]'<br>> > > > 'draw xlab longitude'<br>> > > > 'draw ylab latitude'<br>> > > > 'set mproj latlon'<br>> > > > 'set mpdset hires'<br>> > > > 'set grid off'<br>> > > ><br>> > > > 'close 1'<br>> > > ><br>> > > > 'printim baru.gif white'<br>> > > > 'disable print'<br>> > > ><br>> > > > 'close 1'<br>> > > ><br>> > > > Thank you very much for your help<br>> > > > best regards,<br>> > > > wendi<br>> > > ><br>> > > ><br>> > > ><br>> > > ><br>> > > ><br>> > > > On Sun, Dec 16, 2012 at 7:03 PM, wendi harjupa <!<br>> > wendiharjupa@gmail.com>wrote:<br>> > > ><br>> > > >> Dear Mr. Jeff<br>> > > >><br>> > > >> Thank you very much for your detail explanation<br>> > > >> I will try it.<br>> > > >><br>> > > >> Best regards,<br>> > > >> wendi<br>> > > >><br>> > > >><br>> > > >> On Sat, Dec 15, 2012 at 2:24 AM, Jeff Duda &
lt;jeffduda319@gmail.com><br>> > wrote:<br>> > > >><br>> > > >>> Wendi,<br>> > > >>> Assuming your radar data are available at the same time every day,<br>> > > >>> create an overall average and add daily averages to it. Then take the<br>> > > >>> average of the averages at the end. If you have differing numbers of<br>> > radar<br>> > > >>> valid times each day, you should take a weighted average based on the<br>> > > >>> number of data slices available each day.<br>> > > >>><br>> > > >>> Do you plan to access each day's radar data using a separate control<br>> > > ! >>> file or using a single control file? This distinction ma kes a<br>> > difference<br>> > > >>> in how best to code your script. You could either do it one day at a<br>> > time,<br>> > > >>> or sum all days in one line. That is, you could do it like this:<br>> > > >>><br>> > > >>> 'reinit'<br>> > > >>> 'sumaverage = 0'<br>> > > >>> dd = 1<br>> > > >>> while (dd <= 10)<br>> > > >>> 'open [DAY dd CONTROL FILE]'<br>> > > >>> 'set grads off'<br>> > > >>> 'set mproj latlon'<br>> > > >>> 'set mpdset hires'<br>> > > >>> 'set gxout shaded'<br>> > > >>> 'set lon 100 100.82'<br>> > > >>> 'set lat -0.7 0'<br>> > > >>> 'set z 3600'<br>> > > >>> 'set parea 1 7.5 2.8 8.3'<br>> > > >>> 'set clevs 0 5 10 15 20 25 30 35 40 45 50 55 60'<br>> > > >>> 'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'<br>> > > >>> 'set t 1 3600'<br>> > > >>> 'R=pow(10, ref/10)'<br>> > > >>><br>> > > >>&
gt; day = d! d + 9<br>> > > >>> 'R1=ave(R,time=13Z'dd'APR2004,time=22Z'dd'APR2004)'<br>> > > >>> 'sumaverage = sumaverage + R1'<br>> > > >>> dd = dd + 1<br>> > > >>> endwhile<br>> > > >>> 'average = sumaverage / 10'<br>> > > >>> ...color table options...<br>> > > >>> ...map options...<br>> > > >>> 'd average'<br>> > > >>> 'printim try.gif white'<br>> > > >>><br>> > > >>> OR like this<br>> > > >>><br>> > > >>> 'reinit'<br>> > > >>> 'open [ALL INCLUSIVE CONTROL FILE]'<br>> > > >>> 'set grads off'<br>> > > >>> 'set mproj latlon'<br>> > > >>> 'set mpdset hires'<br>> > > >>> 'set gxout shaded'<br>> > > >>> 'set lon 100 100.82'<br>> > > >>> 'set lat -0.7 0'<br>> > > >>> 'set z 3600'<br>> > > >>> 'set parea 1 7.5 2.8 8.3'<br>> > > >>> 'set clevs 0 5 10 15 20 25 30 35 40 45 50 55 60'<br>> > > >>> 'set c! cols 9 14 4 11 5 13 3 10 7 12 8 2 6'<br>> > > >>> 'set t 1 360 0'<br>> > > >>> 'R=pow(10, ref/10)'<br>> > > >>><br>> > > >>><br>> > 'R1=ave(R,time=13Z10APR2004,time=22Z10APR2004)+ave(R,time=13Z11APR2004,time=22Z11APR2004)+ave(R,time=13Z12APR2004,time=22Z12APR2004)+ave(R,time=13Z13APR2004,time=22Z13APR2004)+ave(R,time=13Z14APR2004,time=22Z14APR2004)+ave(R,time=13Z15APR2004,time=22Z15APR2004)+ave(R,time=13Z16APR2004,time=22Z16APR2004)+ave(R,time=13Z17APR2004,time=22Z17APR2004)+ave(R,time=13Z18APR2004,time=22Z18APR2004)+ave(R,time=13Z19APR2004,time=22Z19APR2004)'<br>> > > >>> ...color table options...<br>> > > >>> ...map options...<br>> > > >>> 'd R1/10'<br>> > &g
t; >>> 'printim try.gif white'<br>> > > >>><br>> > > >>> Jeff<br>> > > >>><br>> > > >>><br>> > > >>> On Fri, Dec 14, 2012 at 12:20 AM, wendi harjupa <<br>> > wendiharjupa@gmail.com>wrote:<br>> > > >>><br>> > > >>>> Dear Jeff<br>> > > >>>><br>> > > >>>> Thank! you very much for your solution,<br>> > > >>>><br>> > > >>>> I need to make average for 10 days, how to make loop for the next<br>> > day<br>> > > >>>> (same time/13-22)??<br>> > > >>>><br>> > > >>>> Best regards,<br>> > > >>>> On Fri, Dec 14, 2012 at 1:06 PM, Jeff Duda <jeffduda319@gmail.com<br>> > >wrote:<br>> > > >>>><br>> > > >>>>> Change the line<br>> > > >>>>> 'R1=ave(R, t=1, t=3600)'<br>> > > >>>>><br>> > > >>>>> to<br>> > > >>>>><br>> > > >>>>> 'R1=ave(R,time=13Z10APR2004,time=22Z10APR2004)'<br>> > > >>>>><br>> > > >>>>> Jeff Duda<br>> > > >>>>><br>> > > >>>>> On Thu, Dec 13, 2012 at 8:42 PM, wendi harjupa <<br>> > > >>>>> wendiharjupa@gmail.com> wrote:<br>> > > >>>>><br>> > > >>>>>> Dear Sir/Madam<br>> > > &g! t;>>>>><br>> > > >>>>>> Please I need your help.<br>> > > >>>>>> I have one set data (10 days data)<br>> > > >>>>>><br>> > > >>>>>> How to make looping to make reflectivity average for only 10<br>> > hours on<br>> > > >>>
;>>> every day (13hrs-22hrs only)??<br>> > > >>>>>><br>> > > >>>>>> Below is my ctl and gs script<br>> > > >>>>>><br>> > > >>>>>> try.ctl<br>> > > >>>>>><br>> > > >>>>>> DSET<br>> > > >>>>>><br>> > /misc/home/wendi/data1/radar/kototabang/XDR/data/%y4%m2%d2/%d2%h2.dat<br>> > > >>>>>> TITLE RADAR<br>> > > >>>>>> OPTIONS TEMPLATE LITTLE_ENDIAN<br>> > > >>>>>> UNDEF -1000<br>> > > >>>>>> XDEF 321 LINEAR 99.688120996 0.004492369<br>> > > >>>>>> YDEF 321 LINEAR -1.083708891 0.004522556<br>> > > >>>>>> ZDEF 39! LINEAR 1.0 0.5<br>> > > >>>>>> TDEF 3600 LINEAR 00:00Z10apr2004 4mn<br>> > > >>>>>> VARS 2<br>> > > >>>>>> ref 39 99 Reflectivity(dBZ)<br>> > > >>>>>> vel 39 99 Doppler Velocity(m/s)<br>> > > >>>>>> ENDVARS<br>> > > >>>>>><br>> > > >>>>>><br>> > > >>>>>> script.gs<br>> > > >>>>>><br>> > > >>>>>> 'reinit'<br>> > > >>>>>> 'open try.ctl'<br>> > > >>>>>> 'set grads off'<br>> > > >>>>>> 'set mproj latlon'<br>> > > >>>>>> 'set mpdset hires'<br>> > > >>>>>> 'set gxout shaded'<br>> > > >>>>>> 'set lon 100 100.82'<br>> > > >>>>>> 'set lat -0.7 0'<br>> > > >>>>>> 'set z 3600'<br>> > > >>>>>> 'set parea 1 7.5 2.8 8.3'<br>> > > >>>>>>
! 'set clevs 0 5 10 15 20 25 30 35 40 45 50 55 60'<br>> > > >>>& gt;>> 'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'<br>> > > >>>>>> 'set t 1 3600'<br>> > > >>>>>> 'R=pow(10, ref/10)'<br>> > > >>>>>><br>> > > >>>>>> 'R1=ave(R, t=1, t=3600)'<br>> > > >>>>>> 'set t 1'<br>> > > >>>>>> 'R2=log10 (R1)'<br>> > > >>>>>> 'd 10*R2'<br>> > > >>>>>> 'run /misc/home/wendi/grads-2.0.1/cbarn 0.7 0.7 7.4 5.5'<br>> > > >>>>>> 'draw string 7.2 8.1 [dBz]'<br>> > > >>>>>> 'draw xlab longitude'<br>> > > >>>>>> 'draw ylab latitude'<br>> > > >>>>>> 'set mproj latlon'<br>> > > >>>>>> 'set mpdset hires'<br>> > > >>>>>> 'set grid off'<br>> > > >>>>>><br>> > > >>>>>> 'close 1'<br>> > > >>>>>><br>> > > >>>>>> 'printim try.gif white'<br>> > > >>>>>! > 'disable print'<br>> > > >>>>>><br>> > > >>>>>> 'close 1'<br>> > > >>>>>><br>> > > >>>>>> Thank you very much for your help<br>> > > >>>>>><br>> > > >>>>>> Best regards,<br>> > > >>>>>> wendi<br>> > > >>>>>><br>> > > >>>>>><br>> > > >>>>>><br>> > > >>>>>><br>> > > >>>>>><br>> > > >>>>>><br>> > > >>>>>><br>> > > >>>>>><br>> > > >>>>>> --<br>&
gt; > > >>>>>> -------------<br>> > > >>>>>> Wendi Harjupa. ST<br>> > > >>>>>> Shimane University Graduate School of Engineering Sciences<br>> > > >>>>>> Department of Electronic Control Systems<br>> > > >>>>>> Remote Sensing Laboratory<br>> > > >>>>>> s119447<br>> > > ! >>>>>> Cellphone : +81-080-4268-6676<br>> > > >&gt ;>>>> email : wendi@rslab.riko.shimane-u.ac.jp<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>> > > >>>>> gradsusr mailing list<br>> > > ! >>>>> gradsusr@gradsusr.org<br>> > > >>>>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> > > >>>>><br>> > > >>>>><br>> > > >>>><br>> > > >>>><br>> > > >>>> --<br>> > >
>>>> -------------<br>> > > >>>> Wendi Harjupa. ST<br>> > > >>>> Shimane University Graduate School of Engineering Sciences<br>> > > >>>> Department of Electronic Control Systems<br>> > > >>>> Remote Sensing Laboratory<br>> > > >>>> s119447<br>> > > >>>> Cellphone : +81-080-4268-6676<br>> > > >>>> email : wendi@rslab.riko.shimane-u.ac.jp<br>> > > >>>><br>> > > >>>><br>> > > >>>> _______________________________________________<br>> > > >>>> gradsusr mailing list<br>> > > >>>> gradsusr@gradsusr.org<br>> > > >>>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> > > &! gt;>>><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>> > > >>> gradsusr mailing list<br>> > > >>> gradsusr@gradsusr.org<br>> > > >>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> > > >>><br>> > > >>><br>> > > >><br>> > > >><br>> > > >> --<br>> > > >> -------------<br>> > > >> Wendi Harjupa. ST<br>> > > >> Shimane University Graduate School of Engineering Sciences<br>> > > >> Department of Electronic Control Systems<br>> > > >> Remote Sensing Laboratory<br>> >
> >> s119447<br>> > > >> Cellphone : +81-080-4268-6676<br>> > > >> email : wendi@rslab.riko.shimane-u.ac.jp<br>> > >! >><br>> > > >><br>> > > ><br>> > > ><br>> > > > --<br>> > > > -------------<br>> > > > Wendi Harjupa. ST<br>> > > > Shimane University Graduate School of Engineering Sciences<br>> > > > Department of Electronic Control Systems<br>> > > > Remote Sensing Laboratory<br>> > > > s119447<br>> > > > Cellphone : +81-080-4268-6676<br>> > > > email : wendi@rslab.riko.shimane-u.ac.jp<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>> > > -------------- next part --------------<br>> > > An HTML attachment was scrubbed...<br>> > > URL: http://gradsusr.org/pipermail/gradsusr/att!<br>> > achments/20121217/19251fc8/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 34, Issue 40<br>> > > ****************************************<br>> ><br>> > _______________________________________________<br>> > g
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/20121218/ddfd8694/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 34, Issue 55<br>> ****************************************<br></div> </div></body>
</html>