<div dir="ltr"><div><div><div><div>Hello fellow grad user,<br></div> <br></div>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.<br><br></div>

<div>Any help will be appreciated.<br><br><br><br>&#39;open obsnudge/observation_nudging.ctl&#39;<br>&#39;open model/fnl_nonudging.ctl&#39;<br>Time_val=2<br>Time_prev=1<br>&#39;set t &#39;Time_val<br>&#39;define rn=(rainc+rainnc)-(rainc(t=&#39;Time_prev&#39;)-rainnc(t=&#39;Time_prev&#39;))&#39;<br>

&#39;d rn&#39;<br>&#39;define dif=rn.1-rn.2&#39;<br> &#39;define sq=dif*dif&#39;<br> &#39;define am=amean(sq,lon=115,lon=135,lat=4,lat=25)&#39;<br> &#39;define rmse=sqrt(am)&#39;<br> &#39;d rmse&#39;<br> a=subwrd(result,4)<br>

 output=a<br> res=write(rmsewind_v.csv,output,append)<br> &#39;undefine dif&#39;<br> &#39;undefine sq&#39;<br> &#39;undefine am&#39;<br><br>Time_val=Time_val+1<br>Time_prev=Time_prev+1<br>endwhile<br>quit<br><br></div><br>

<br></div>Lyndon<br></div>