[gradsusr] looping program

[Muhammad Yunus Ahmad Mazuki]

My advice is that you use different naming in the variable R1 and R2 for
this. I would suggest R1-R10 for 10 days, where
'R1=ave(R,time=13Z10APR2004,time=22Z10APR2004)'
...
'R10=ave(R,time=13Z19APR2004,time=22Z19APR2004)'
You just change the day day part in the date.

Average R1-R10 (eg RAve=(R1+...R10)/10
Change your R2 into something else to avoid confliction, perhaps
'Rsomething=log10 (RAve)'

This however is only short term workaround. If you want to automate this
process, a complicated script might be needed..

Muhammad Yunus

On Fri, Dec 14, 2012 at 2:20 PM, wendi harjupa <wendiharjupa at gmail.com>wrote:

> Dear Jeff
>
> Thank you very much for your solution,
>
> I need to make average for 10 days, how to make loop for the next day
> (same time/13-22)??
>
> Best regards,
> On Fri, Dec 14, 2012 at 1:06 PM, Jeff Duda <jeffduda319 at gmail.com> wrote:
>
>> Change the line
>> 'R1=ave(R, t=1, t=3600)'
>>
>> to
>>
>> 'R1=ave(R,time=13Z10APR2004,time=22Z10APR2004)'
>>
>> Jeff Duda
>>
>>  On Thu, Dec 13, 2012 at 8:42 PM, wendi harjupa <wendiharjupa at gmail.com>wrote:
>>
>>>  Dear Sir/Madam
>>>
>>> Please I need your help.
>>> I have one set data (10 days data)
>>>
>>> How to make looping to make reflectivity average for only 10 hours on
>>> every day (13hrs-22hrs only)??
>>>
>>> Below is my ctl and gs script
>>>
>>> try.ctl
>>>
>>> DSET
>>> /misc/home/wendi/data1/radar/kototabang/XDR/data/%y4%m2%d2/%d2%h2.dat
>>> TITLE  RADAR
>>> OPTIONS  TEMPLATE LITTLE_ENDIAN
>>> UNDEF  -1000
>>> XDEF    321 LINEAR    99.688120996    0.004492369
>>> YDEF    321 LINEAR    -1.083708891    0.004522556
>>> ZDEF    39  LINEAR   1.0 0.5
>>> TDEF    3600 LINEAR 00:00Z10apr2004  4mn
>>> VARS    2
>>> ref 39  99  Reflectivity(dBZ)
>>> vel 39  99  Doppler Velocity(m/s)
>>> ENDVARS
>>>
>>>
>>> script.gs
>>>
>>> 'reinit'
>>> 'open try.ctl'
>>> 'set grads off'
>>> 'set mproj latlon'
>>> 'set mpdset hires'
>>> 'set gxout shaded'
>>> 'set lon 100  100.82'
>>> 'set lat -0.7  0'
>>> 'set z 3600'
>>> 'set parea 1 7.5 2.8 8.3'
>>> 'set clevs 0 5 10 15 20 25 30 35 40 45 50 55 60'
>>> 'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'
>>> 'set t 1 3600'
>>> 'R=pow(10, ref/10)'
>>>
>>> 'R1=ave(R, t=1, t=3600)'
>>> 'set t 1'
>>> 'R2=log10 (R1)'
>>> 'd 10*R2'
>>> 'run /misc/home/wendi/grads-2.0.1/cbarn 0.7 0.7 7.4 5.5'
>>> 'draw string  7.2 8.1 [dBz]'
>>> 'draw xlab longitude'
>>> 'draw ylab latitude'
>>> 'set mproj latlon'
>>> 'set mpdset hires'
>>> 'set grid off'
>>>
>>> 'close 1'
>>>
>>> 'printim try.gif white'
>>> 'disable print'
>>>
>>> 'close 1'
>>>
>>> Thank you very much for your help
>>>
>>> Best regards,
>>> wendi
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> --
>>> -------------
>>> Wendi Harjupa. ST
>>> Shimane University Graduate School of Engineering Sciences
>>> Department of Electronic Control Systems
>>> Remote Sensing Laboratory
>>> s119447
>>> Cellphone : +81-080-4268-6676
>>> email : wendi at rslab.riko.shimane-u.ac.jp
>>>
>>>
>>> _______________________________________________
>>> gradsusr mailing list
>>> gradsusr at gradsusr.org
>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>
>>>
>>
>>
>> --
>> Jeff Duda
>> Graduate research assistant
>> University of Oklahoma School of Meteorology
>> Center for Analysis and Prediction of Storms
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
>
> --
> -------------
> Wendi Harjupa. ST
> Shimane University Graduate School of Engineering Sciences
> Department of Electronic Control Systems
> Remote Sensing Laboratory
> s119447
> Cellphone : +81-080-4268-6676
> email : wendi at rslab.riko.shimane-u.ac.jp
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20121214/932ee42e/attachment-0003.html