[gradsusr] looping program

[wendi harjupa]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20121214/c3f26290/attachment-0003.html