[gradsusr] make time series data in one plot
wendi harjupa
wendiharjupa at gmail.com
Sun Dec 8 23:22:18 EST 2013
Dear Sir/Madam
I need help in my program...
I have one set data of radar. I made ctl file as follow :
DSET /media/HD-PVRU2/%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 15 LINEAR 13:00Z10apr2004 4mn
VARS 2
ref 39 99 Reflectivity(dBZ)
vel 39 99 Doppler Velocity(m/s)
ENDVARS
As you see, one folder data contain 24 files of data (24 hours)
I made already script program to plot one hour data as follow :
### Reflectivity overlayed topography ###
## Open and Initiating data ##
'reinit'
'open karajo.ctl'
'set grads off'
'set mproj latlon'
'set mpdset hires'
'set gxout shaded'
'set lon 99.65 101.2'
'set lat -1.075 0.35'
'set z 5'
'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 0 9 14 4 11 5 13 3 10 7 12 8 2 6'
## average 60 minute ##
'set t 1 15'
'R=pow(10, ref/10)'
'R1=ave(R, t=1, t=15)'
'set t 1'
'R2=log10 (R1)'
## Plotting ##
'd 10*R2'
'run /usr/local/lib/grads/cbarn 0.7 0.7 7.4 5.5'
'draw string 7.2 8.1 [dBz]'
'draw title CAPPI of Reflectivity at 3 km'
'draw xlab longitude'
'draw ylab latitude'
'set mproj latlon'
'set mpdset hires'
'set grid off'
'close 1'
'printim 1013.jpg white'
'disable print'
'close 1'
*My question is, how to plot 24 hours data directly??*
*I mean. when I execute program, I can get 24 figures in one time...*
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/20131209/c708cf53/attachment.html
More information about the gradsusr
mailing list