<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: arial,helvetica,sans-serif; font-size: 12pt; color: #000000'><div style="color: rgb(0, 0, 0); ">There is still a sign error in your calculation of rn1 and rn2. It should be:</div><div style="color: rgb(0, 0, 0); "><br></div><div>'define rn1=(rainc+rainnc)-(rainc(t='Time_prev')<font class="Apple-style-span" color="#ff0000">+</font>rainnc(t='Time_prev'))'<br>'define rn2=(rainc.2+rainnc.2)-(rainc.2(t='Time_prev')<font class="Apple-style-span" color="#ff0000">+</font>rainnc.2(t='Time_prev'))'<br><br></div><div><hr id="zwchr" style="color: rgb(0, 0, 0); "><b style="color: rgb(0, 0, 0); ">From: </b>"Lyndon Mark Olaguera" <olagueralyndonmark429@gmail.com><br><b style="color: rgb(0, 0, 0); ">To: </b>"GrADS Users Forum" <gradsusr@gradsusr.org><br><b style="color: rgb(0, 0, 0); ">Sent: </b>Wednesday, May 15, 2013 9:47:48 PM<br><b style="color: rgb(0, 0, 0); ">Subject: </b>Re: [gradsusr] RMSE for accumulated Rainfall<br><br><div dir="ltr" style="color: rgb(0, 0, 0); "><div><div>thanks for the kindest reply..This is a big help..<br></div>Here's my modified script:<br><br>'open 1D_nudging/surfnudge/surface_grid_nudging.ctl'<br>'open 1D_nudging/model/fnl_nonudging.ctl'<br>
<br>prompt 'Max Time: '<br>pull Max_time<br>prompt 'Title Header: '<br>pull Title_header<br>Time_val=2<br>Time_prev=1<br>while(Time_val<=Max_time)<br><br>'set t 'Time_val<br>'define rn1=(rainc+rainnc)-(rainc(t='Time_prev')-rainnc(t='Time_prev'))'<br>
'define rn2=(rainc.2+rainnc.2)-(rainc.2(t='Time_prev')-rainnc.2(t='Time_prev'))'<br>'define dif=rn1-rn2'<br>'define sq=dif*dif'<br>'define am=amean(sq,lon=115,lon=135,lat=4,lat=25)'<br>
'define rmse=sqrt(am)'<br>'d rmse'<br><br>rmse=subwrd(result,4)<br>output=rmse<br>res=write(rmse_rain.csv,output,append)<br><br>Time_val=Time_val+1<br>Time_prev=Time_prev+1<br>endwhile<br>quit<br><br><br>
</div>
This is working now..Thanks for your suggestions.<br><div><div><br><br></div></div></div><div class="gmail_extra" style="color: rgb(0, 0, 0); "><br><br><div class="gmail_quote">On Thu, May 16, 2013 at 12:20 AM, Jeff Duda <span dir="ltr"><<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">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.<span class="HOEnZb"><font color="#888888"><div>
<br></div><div>Jeff</div></font></span></div><div class="gmail_extra">
<br><br><div class="gmail_quote"><div><div class="h5">On Wed, May 15, 2013 at 2:53 AM, Lyndon Mark Olaguera <span dir="ltr"><<a href="mailto:olagueralyndonmark429@gmail.com" target="_blank">olagueralyndonmark429@gmail.com</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, May 15, 2013 at 1:06 PM, Jeff Duda <span dir="ltr"><<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:;border-left-style:solid;padding-left:1ex"><div>'define fcstrain = rainc.1 + rainnc.1 - rainc.1(t-1) - rainnc.1(t-1)'</div>
<div>'define obsrain = rainc.2 + rainnc.2 - rainc.2(t-1) - rainnc.2(t-1)'</div>
<div>'define rainerror = fcstrain - obsrain'</div></blockquote></div><br>Hello,</div><div class="gmail_extra">Thanks for the reply..</div><div class="gmail_extra">I have another question on this line..</div>
<div class="gmail_extra"><br></div><div class="gmail_extra"><div style="font-family:arial,sans-serif;font-size:13px">'define fcstrain = rainc.1 + rainnc.1 - rainc.1(t-1) - rainnc.1(t-1)'</div><div style="font-family:arial,sans-serif;font-size:13px">
'define obsrain = rainc.2 + rainnc.2 - rainc.2(t-1) - rainnc.2(t-1)'</div><div style="font-family:arial,sans-serif;font-size:13px">'define rainerror = fcstrain - obsrain'</div><div style="font-family:arial,sans-serif;font-size:13px">
<br></div><div style="font-family:arial,sans-serif;font-size:13px">should this be </div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px"><div>'define fcstrain = (rainc.1 + rainnc.1) - (rainc.1(t-1) + rainnc.1(t-1))'</div>
<div>'define obsrain =( rainc.2 + rainnc.2) - (rainc.2(t-1) + rainnc.2(t-1))'</div><div><br></div><div>plus not minus?</div><div><br></div><div><br></div><div><br></div><div>I will try this later</div><span><font color="#888888">
<div>Lyndon</div></font></span></div><div style="font-family:arial,sans-serif;font-size:13px"><br></div></div></div>
<br></div></div><div class="im">_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></div></blockquote></div><div class="im"><br><br clear="all"><div><br></div>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br>
</div></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div>
<br>_______________________________________________<br>gradsusr mailing list<br>gradsusr@gradsusr.org<br>http://gradsusr.org/mailman/listinfo/gradsusr<br></div></div></body></html>