<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'>
I couldn't open it again. I tries but It send me an open error : can't 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 path (C:\grads\data\). How can I write the path of my file in open command? <BR> <BR>Tahereh<br><BR><div>> 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>> 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: 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 accordingly.<br>> If the incorrect pathname is specified there, you won't be able to open the<br>> data.<br>> <br>> Jeff Duda<br>> <br>> On Mon, Dec 17, 2012 at 6:30 AM, T. haghroosta <haghroosta@hotmail.com>wrote:<br>> <br>> > Dear All,<br>> ><br>> > I have a .ctl file on my drive E and my grads software is installed on 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 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: 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 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>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>> > '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>> > '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 <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 <jeffduda319@gmail.com> 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 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 difference<br>> >>> in how best to code your script. You could either do it one day at a 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>> >>> day = dd + 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 ccols 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>> >>> '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>> >>> 'printim try.gif white'<br>> >>><br>> >>> Jeff<br>> >>><br>> >>><br>> >>> On Fri, Dec 14, 2012 at 12:20 AM, wendi harjupa <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 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>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>> >>>>>><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 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>> >>>>>> /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>> >>>>>> -------------<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>> >>>>><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>> >>>><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/attachments/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></div> </div></body>
</html>