[gradsusr] Accumulated Rain (6 hourly data)

[Ricardo Hallak]

Hi,
just try:
'd rainc+rainnc-(rainc(t='Time_val'-1)-rainnc(t='Time_val'-1))'
once Time_val is a variable, and the number 1 is not. 
So, Time_val have to be put between quotes, and the number 1, not.
Hope this helps,
Ricardo

On Tue, 7 May 2013 17:20:42 +0800, Lyndon Mark Olaguera wrote
> Hello fellow GRAD user,
> 
> I have a problem in plotting the accumulated rain in grads.,
> my data comes in 6 hourly...
> my sample script is shown below.
> the dates are:
> 08/06/2012/00
> 08/06/2012/06
> 
> 08/06/2012/12
> 08/06/2012/08
> 08/07/2012/00
> 08/07/2012/06
> 08/07/2012/12
> 08/07/2012/18
> 08/08/2012/00
> 08/08/2012/06
> 08/08/2012/12
> 08/08/2012/18
> 08/09/2012/00
> 08/09/2012/06
> 08/09/2012/12
> 08/09/2012/18
> 08/10/2012/00
> 
> I want to plot the accumulated rain for each hour but my script does not ouput it correctly.Could anyone help me correct my script?
> 
> i'll appreciate any help.
> 
> thanks
> 
> 'open surfnudge/surf.ctl'
> 'set display color white'
> 'clear'
> 'set gxout shaded'
> 'set mpdset hires'
> 'set lat  4   24'
> 'set lon  115 134'  
> 
> prompt 'Max Time: '
> pull Max_time
> prompt 'Title Header: '
> pull Title_header
> 
> Time_val=2
> *Time_val=8
> while(Time_val<=Max_time)
> 'set t 'Time_val
> DateStr = subwrd(result,4)
> 'c'
> *Every 3 hours for model output
> 'd rainc+rainnc-(rainc(t='Time_val-1')-rainnc(t='Time_val-1'))'
> 'draw title 'Title_header'\'DateStr
> 'cbarn'
> 'printim 'Title_header'_'DateStr'.png'
> 'c'
> Time_val=Time_val+1
> endwhile
> 'quit'

Ricardo Hallak

Departamento de Ciências Atmosféricas (ACA) 
Instituto de Astronomia, Geofísica e Ciências Atmosféricas (IAG) 
Universidade de São Paulo (USP) 
Rua do Matão, 1226 - Cidade Universitária - São Paulo SP 05508-900 
Phone.: +55 (11) 3091-2850; 3091-4703 
Fax: +55 (11) 3091-4714

 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130507/b465c0f3/attachment-0003.html