[gradsusr] Time series of hourly accumulated rainfall plotting problem

[Sahidul]

Dear Help,
I would like plot time series of WRF hourly rainfall plot for 48hour by
finxing lat & lone. Iam able plot time series of
Accumulated rainfall But I would like to plot hour accumulated
rainfall,Please help me to solve the error
******************************
**********************************************
'open file.ctl'
'set lat 13.1'
'set lon 80.00'
'set t 1 49'
'set grads off'
'set grid off'
'set tlsupp year'
'set ccolor 2'
'set cthick 7'
'd rainc+rainnc'
'draw title Preciption in mm/hour'
'draw ylab Accumilated Rainfall(mm)'
'printim CHENNAI.gif white'
'quit'
*****************************************************************************
*****But I would like to plot hourly accumulated rainfall (each hour)
*******
*****************************************************************************
Using the following script I am trying to plot hourly rainfall time series
'open file.ctl'
tt=1
'set lat 13.1'
'set lon 80.00'
while ( tt <= 49)
'set gxout line'
'set t ' tt
'set grads off'
'set grid off'
'set ccolor 2'
'set cthick 7'
'define
rn=(rainc(t='tt+1')+rainnc(t='tt+1'))-(rainc(t='tt')+rainnc(t='tt'))'
'd rn'
fmt='%6.2f'
tmpval=subwrd(result,4)
drain=math_format(fmt,tmpval)
say drain
tt=tt+1
endwhile
'set gxout line'
'd drain'
'draw ylab Hourly Rainfall(mm)'
'printim CHENNAI.gif white'
'quit'
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110830/5b7441f1/attachment-0003.html