[gradsusr] RMSE for accumulated Rainfall

[Lyndon Mark Olaguera]

Hello fellow grad user,

i want to find the rmse for the accumulated rainfall,I pasted my sample
script below but it is not working..grads is giving an error.

Any help will be appreciated.



'open obsnudge/observation_nudging.ctl'
'open model/fnl_nonudging.ctl'
Time_val=2
Time_prev=1
'set t 'Time_val
'define rn=(rainc+rainnc)-(rainc(t='Time_prev')-rainnc(t='Time_prev'))'
'd rn'
'define dif=rn.1-rn.2'
 'define sq=dif*dif'
 'define am=amean(sq,lon=115,lon=135,lat=4,lat=25)'
 'define rmse=sqrt(am)'
 'd rmse'
 a=subwrd(result,4)
 output=a
 res=write(rmsewind_v.csv,output,append)
 'undefine dif'
 'undefine sq'
 'undefine am'

Time_val=Time_val+1
Time_prev=Time_prev+1
endwhile
quit



Lyndon
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130515/9289cfdf/attachment-0003.html