[gradsusr] RMSE for accumulated Rainfall

[Eric Altshuler]

There is still a sign error in your calculation of rn1 and rn2. It should be: 


'define rn1=(rainc+rainnc)-(rainc(t='Time_prev') + rainnc(t='Time_prev'))' 
'define rn2=(rainc.2+rainnc.2)-(rainc.2(t='Time_prev') + rainnc.2(t='Time_prev'))' 


----- Original Message -----
From: "Lyndon Mark Olaguera" <olagueralyndonmark429 at gmail.com> 
To: "GrADS Users Forum" <gradsusr at gradsusr.org> 
Sent: Wednesday, May 15, 2013 9:47:48 PM 
Subject: Re: [gradsusr] RMSE for accumulated Rainfall 




thanks for the kindest reply..This is a big help.. 
Here's my modified script: 

'open 1D_nudging/surfnudge/surface_grid_nudging.ctl' 
'open 1D_nudging/model/fnl_nonudging.ctl' 

prompt 'Max Time: ' 
pull Max_time 
prompt 'Title Header: ' 
pull Title_header 
Time_val=2 
Time_prev=1 
while(Time_val<=Max_time) 

'set t 'Time_val 
'define rn1=(rainc+rainnc)-(rainc(t='Time_prev')-rainnc(t='Time_prev'))' 
'define rn2=(rainc.2+rainnc.2)-(rainc.2(t='Time_prev')-rainnc.2(t='Time_prev'))' 
'define dif=rn1-rn2' 
'define sq=dif*dif' 
'define am=amean(sq,lon=115,lon=135,lat=4,lat=25)' 
'define rmse=sqrt(am)' 
'd rmse' 

rmse=subwrd(result,4) 
output=rmse 
res=write(rmse_rain.csv,output,append) 

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


This is working now..Thanks for your suggestions. 








On Thu, May 16, 2013 at 12:20 AM, Jeff Duda < jeffduda319 at gmail.com > wrote: 



Yes, you can choose to leave the expression in that form. What I have done is simply to distribute the minus sign. The expressions are mathematically identical. 


Jeff 





On Wed, May 15, 2013 at 2:53 AM, Lyndon Mark Olaguera < olagueralyndonmark429 at gmail.com > wrote: 

<blockquote>






On Wed, May 15, 2013 at 1:06 PM, Jeff Duda < jeffduda319 at gmail.com > wrote: 

<blockquote>

'define fcstrain = rainc.1 + rainnc.1 - rainc.1(t-1) - rainnc.1(t-1)' 
'define obsrain = rainc.2 + rainnc.2 - rainc.2(t-1) - rainnc.2(t-1)' 
'define rainerror = fcstrain - obsrain' 


Hello, 
Thanks for the reply.. 
I have another question on this line.. 



'define fcstrain = rainc.1 + rainnc.1 - rainc.1(t-1) - rainnc.1(t-1)' 
'define obsrain = rainc.2 + rainnc.2 - rainc.2(t-1) - rainnc.2(t-1)' 
'define rainerror = fcstrain - obsrain' 


should this be 



'define fcstrain = (rainc.1 + rainnc.1) - (rainc.1(t-1) + rainnc.1(t-1))' 
'define obsrain =( rainc.2 + rainnc.2) - (rainc.2(t-1) + rainnc.2(t-1))' 


plus not minus? 






I will try this later 
Lyndon 



_______________________________________________ 
gradsusr mailing list 
gradsusr at gradsusr.org 
http://gradsusr.org/mailman/listinfo/gradsusr 


</blockquote>





-- 
Jeff Duda 
Graduate research assistant 
University of Oklahoma School of Meteorology 
Center for Analysis and Prediction of Storms 

_______________________________________________ 
gradsusr mailing list 
gradsusr at gradsusr.org 
http://gradsusr.org/mailman/listinfo/gradsusr 


</blockquote>


_______________________________________________ 
gradsusr mailing list 
gradsusr at gradsusr.org 
http://gradsusr.org/mailman/listinfo/gradsusr 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130516/4c183c62/attachment-0003.html