Hello Tahereh<div><br></div><div>First of all I made a slight mistake in my last mail. For RMSE calculations the term &quot;A&quot; in the formula should be not the parameter you want to calculate but the error of the parameter. That is A should be error value. Something like model minus observation.</div>
<div><br></div><div>To answer your question, as I am not familiar with you data. I would like you to download outgoing longwave radiation (OLR) data from this website. <a href="http://www.esrl.noaa.gov/psd/data/gridded/data.interp_OLR.html">http://www.esrl.noaa.gov/psd/data/gridded/data.interp_OLR.html</a>. Download the <a href="http://olr.daily.mean.nc">olr.daily.mean.nc</a>. Now write the following lines in the terminal where you have downloaded and stored your OLR data.</div>
<div><br></div><div>sdfopen <a href="http://olr.day.mean.nc">olr.day.mean.nc</a></div><div>set t 1 4</div><div>set lon 85 100</div><div>set lat 0 10</div><div>define RMSE=sqrt(ave((olr*olr),t=1,t=4))</div><div>d RMSE</div>
<div><br></div><div>This should work. I have tried with the same dataset and the output has come. If this is successful, try the same with you dataset and see what the result is.</div><div>Now the mistake in this method is the term &quot;olr&quot; is not error. RMSE formula is sqrt(ave(error*error),t=1,t=4)). So usually error must be added, not the parameter itself. Hence, you write another command before defining RMSE, where you define error by denoting &quot;define error=(forecast-observation)&quot;. You can read this short webpage to understand RMSE. <a href="http://www.weatheroffice.gc.ca/verification/scores/rmse_e.html">http://www.weatheroffice.gc.ca/verification/scores/rmse_e.html</a>.</div>
<div><br></div><div>With regards</div><div>Amulya Chevuturi<br><br><div class="gmail_quote">On Tue, Dec 25, 2012 at 12:42 PM, T. haghroosta <span dir="ltr">&lt;<a href="mailto:haghroosta@hotmail.com" target="_blank">haghroosta@hotmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


<div><div dir="ltr">
Dear Amulya,<br> <br>I wrote these commands as you said in a .gs file (<a href="http://test.gs" target="_blank">test.gs</a>) and after that I tried to run the <a href="http://test.gs" target="_blank">test.gs</a>, I got nothing. I mean when I used the command &quot;run <a href="http://test.gs" target="_blank">test.gs</a>&quot;, grads did nothin, niether error nor output!<br>
 <br> <br>&#39;open tem.ctl&#39;<br>&#39;set t 1 4&#39;<br>&#39;set lon 85 100&#39;<br>&#39;set lat 0 10&#39;<br>&#39;define RMSE=sqrt(ave(T*T,t=1,t=4))&#39;<br>&#39;d RMSE&#39;<br><br> and also in another .gs file :<br> <br>
 <br><div><div>&#39;open tem.ctl&#39;<br>&#39;set t 1 4&#39;<br>&#39;set lon 85 100&#39;<br>&#39;set lat 0 10&#39;<br>&#39;define RMSE=sqrt(aave(lh*lh,lon=85,lon=100,lat=0,lat=10))&#39;<br>&#39;d RMSE&#39;</div><div> </div>
<div>I also couldn&#39;t work with the sample scripts.</div><div> </div><div>Could you please lead me how to do? </div><div> </div><div>Regards,</div><div> </div><div>Tahereh</div><br><br>&gt; ------------------------------<br>
&gt; <br>&gt; Message: 3<br>&gt; Date: Tue, 25 Dec 2012 04:07:21 +0000<br>&gt; From: &quot;T.!
  haghroosta&quot; &lt;<a href="mailto:haghroosta@hotmail.com" target="_blank">haghroosta@hotmail.com</a>&gt;<br>&gt; Subject: [gradsusr] Root Mean Square Error<br>&gt; To: &lt;<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>&gt;<br>
&gt; Message-ID: &lt;DUB121-W43F87358FE9C823BDD5BA2DD3A0@phx.gbl&gt;<br>&gt; Content-Type: text/plain; charset=&quot;windows-1256&quot;<div class="im"><br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; Dear All, Could you please tell me how I can calculate RMSE (root mean square error) in Grads and draw relevant graphs?  Regards, Tahereh                                               <br>
</div>&gt; -------------- next part --------------<br>&gt; An HTML attachment was scrubbed...<br>&gt; URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20121225/e2532241/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20121225/e2532241/attachment-0001.html</a> <br>
&gt; <br>&gt; ------------------------------<br>&gt; <br>&gt; Message: 4<br>&gt; Date: Tue, 25 Dec 2012 09:55:21 +0530<br>&gt; From: Amulya Chevuturi &lt;<a href="mailto:amulya.chevuturi@gmail.com" target="_blank">amulya.chevuturi@gmail.com</a>&gt;<br>
&gt; Subject: Re: [gradsusr] Root Mean Square Error<br>&gt; To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>&gt;<br>&gt; Message-ID:<u></u>&gt;         &lt;CAOV67LkiGd7ayjps3XTMW-3FWq5V=Ao4R4-oz7BUowZN0gkhqA@mail.g

<a href="http://mail.com" target="_blank">mail.com</a>&gt;<br>&gt; Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<div class="im"><br>&gt; <br>&gt; Hello Tahereh,<br>&gt; <br>&gt; This is a website containing quite many scripts for various analysis. There<br>
&gt; is a RMSE calculation script also called <a href="http://ltrend.gs" target="_blank">ltrend.gs</a>. Please look it and see<br>&gt; if it is useful for your work.<br>&gt; <a href="http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html" target="_blank">http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html</a><br>
&gt; <br>&gt; Usually the RMSE calculation is done through the commands:<br>&gt; <br>&gt; &quot;define RMSE=sqrt(ave((A*A),t=t1,t=tn)&quot;<br>&gt; Where,<br>&gt; RMSE=Root Mean Square Error<br>&gt; A=Is the parameter you want to calculate<br>
&gt; t1=Time step 1 for you series<br>&gt; tn=Last time step for your series<br>&gt; This will give you the plot of a time period.<br>&gt; <br>&gt; OR<br>&gt; <br>&gt; &quot;define RMSE=sqrt(aave((A*A),lon=lon1,lon=lonn,lat=lat1,lat=latn)&quot;<br>
&gt; Where<br>&gt; RMSE=Root Mean Square Error<br>&gt; A=Is the parameter you want to calculate<br></div>&gt; Lat1, Latn, Lon1 and Lonn = Tel!
 l the extent of your study<div class="im"><br>&gt; This will give you a time series of RMSE.<br>&gt; <br>&gt; Hopefully this clears your doubt.<br>&gt; <br>&gt; With regards,<br>&gt; Amulya Chevuturi<br>&gt; School of Environmental Sciences<br>
&gt; Jawaharlal Nehru University<br>&gt; New Delhi<br>&gt; <br>&gt; On Tue, Dec 25, 2012 at 9:37 AM, T. haghroosta &lt;<a href="mailto:haghroosta@hotmail.com" target="_blank">haghroosta@hotmail.com</a>&gt;wrote:<br>&gt; <br>
&gt; &gt;  Dear All,<br>&gt; &gt;<br>&gt; &gt; Could you please tell me how I can calculate RMSE (root mean square error)<br>&gt; &gt; in Grads and draw relevant graphs?<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt; Regards,<br>
&gt; &gt;<br>&gt; &gt; Tahereh<br>&gt; &gt;<br>&gt; &gt; _______________________________________________<br>&gt; &gt; gradsusr mailing list<br>&gt; &gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
&gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt; &gt;<br>&gt; &gt;<br></div>&gt; -------------- next part --------------<br>&gt; An HTML attachment was scrubbed...<br>
&gt; URL: <a href="http://grad" target="_blank">http://grad</a>!
 <a href="http://susr.org/pipermail/gradsusr/attachments/20121225/c861c614/attachment-0" target="_blank">susr.org/pipermail/gradsusr/attachments/20121225/c861c614/attachment-0</a>

001.html <br>&gt; <br>&gt; ------------------------------<br>&gt; <br>&gt; Message: 5<br>&gt; Date: Tue, 25 Dec 2012 10:08:26 +0530<br>&gt; From: madhu latha &lt;<a href="mailto:madhulatha11@gmail.com" target="_blank">madhulatha11@gmail.com</a>&gt;<br>
&gt; Subject: Re: [gradsusr] problem for prebufr data<br>&gt; To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>&gt;<br>&gt; Message-ID:<br>&gt;         &lt;<a href="mailto:CAAhFJCSWRsxZMmPfNuZosfDgbzfJETwccC%2BsJxoSgMroxmWmTw@mail.gmail.com" target="_blank">CAAhFJCSWRsxZMmPfNuZosfDgbzfJETwccC+sJxoSgMroxmWmTw@mail.gmail.com</a>&gt;<br>
&gt; Content-Type: text/plain; charset=&quot;windows-1252&quot;<br>&gt; <br>&gt; hi<br>&gt; <br>&gt;   Thank you<br>&gt; First itself i have converted prebufr file to litle endian prebufr file<br>&gt; using ssrc.exe<br>&gt; then i have assimilated data using 3dvar.<br>
&gt; Then mentioned problem came<br>&gt; <br>&gt; thank you.<br>&gt; <br>&gt; On Mon, Dec 24, 2012 at 11:56 PM, Rabah Hachelaf &lt;<a href="mailto:r.hachelaf@gmail.com" target="_blank">r.hachelaf@gmail.com</a>&gt;wrote:<br>
&gt; <br>&gt; &gt; Hello,<br>&gt; &gt;<br>&gt; &gt; try to convert your prepbufr file from big endian to littele endian buy<br>&gt; &gt; the ssrc.c  too!
 l.<br>&gt; &gt;<br>&gt; &gt; ssrc.exe &lt;Big Endian prepbufr file&gt; Little Endian prepbufr file<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt; 2012/12/24 madhu latha &lt;<a href="mailto:madhulatha11@gmail.com" target="_blank">madhulatha11@gmail.com</a>&gt;<br>
&gt; &gt;<br>&gt; &gt;&gt; hi<br>&gt; &gt;&gt; earlier i have tried the suggestions in the link by setting the<br>&gt; &gt;&gt; environment<br>&gt; &gt;&gt; setenv F_UFMTENDIAN &quot;little:94-99&quot;<br>&gt; &gt;&gt; but no use<br>
&gt; &gt;&gt;<br>&gt; &gt;&gt; then itried to do reblock.c mentioned in the same link.<br>&gt; &gt;&gt; but i did not understand how to use the utility reblock.c<br>&gt; &gt;&gt;<br>&gt; &gt;&gt; kindly suggest.<br>&gt; &gt;&gt;<br>
&gt; &gt;&gt;<br>&gt; &gt;&gt; On Mon, Dec 24, 2012 at 10:15 PM, Rabah Hachelaf &lt;<a href="mailto:r.hachelaf@gmail.com" target="_blank">r.hachelaf@gmail.com</a>&gt;wrote:<br>&gt; &gt;&gt;<br>&gt; &gt;&gt;&gt; Hello,<br>
&gt; &gt;&gt;&gt; see this link , it may be help you.<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt; <a href="http://www.mmm.ucar.edu/wrf/users/wrfda/Docs/readBufr.htm" target="_blank">http://www.mmm.ucar.edu/wrf/users/wrfda/Docs/readBufr.htm</a><br>
&gt; &gt;&gt;&gt;<u></u>&gt; &gt;&gt;&gt; --<br>&gt; &gt;&gt;&gt; Cordialement,<br>&gt; &amp;gt

;&gt;&gt; Best regards,<br>&gt; &gt;&gt;&gt; Rabah Hachelaf<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;      ____<br>&gt; &gt;&gt;&gt;     (        )<br>&gt; &gt;&gt;&gt;    (        )<br>&gt; &gt;&gt;&gt;   (___ __)<br>&gt; &gt;&gt;&gt;      /////////<br>
&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt; 2012/12/24 madhu latha &lt;<a href="mailto:madhulatha11@gmail.com" target="_blank">madhulatha11@gmail.com</a>&gt;<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;&gt; Dear all<br>
&gt; &gt;&gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;&gt; Iam trying to run 3dvar on intel compiler<br>&gt; &gt;&gt;&gt;&gt; i have successfully installed wrf3dvar<br>&gt; &gt;&gt;&gt;&gt; but when iam trying to assimilate prepbufr data it is showing the<br>
&gt; &gt;&gt;&gt;&gt; following error<br>&gt; &gt;&gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;&gt; Error message<br>&gt; &gt;&gt;&gt;&gt; **************BUFR ARCHIVE LIBRARY ABORT*****************<br>&gt; &gt;&gt;&gt;&gt; BUFRLIB: OPENBF - ERROR READING INPUT FILE CONNECTED TO !
 UNIT 96 WHEN<br>&gt; &gt;&gt;&gt;&gt; CHECKING<br>&gt; &gt;&gt;&gt;&gt; FOR ?BUFR? IN FIRST 4 BYTES OF RECORD<br>&gt; &gt;&gt;&gt;&gt; **************BUFR ARCHIVE LIBRARY ABORT*****************<br>&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;&gt; can anyone know how to fix the problem<br>&gt; &gt;&gt;&gt;&gt; Any help will be appreciated.<br>&gt; &gt;&gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;&gt; Thank you.<br>&gt; &gt;&gt;&gt;&gt;<br>
&gt; &gt;&gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;&gt; --<br>&gt; &gt;&gt;&gt;&gt; *With Regards<br>&gt; &gt;&gt;&gt;&gt; Madhu Latha Akkisetti<br>&gt; &gt;&gt;&gt;&gt; SRF<br>&gt; &gt;&gt;&gt;&gt; NARL &lt;<a href="http://www.narl.gov.in" target="_blank">http://www.narl.gov.in</a>&gt;<br>
&gt; &gt;&gt;&gt;&gt; Department of space<br>&gt; &gt;&gt;&gt;&gt; Gadanki*<div class="im"><br>&gt; &gt;&gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;&gt; _______________________________________________<br>&gt; &gt;&gt;&gt;&gt; gradsusr mailing list<br>
&gt; &gt;&gt;&gt;&gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br></div>&gt; &gt;&gt;&gt;&gt; <a href="http://gradsusr" target="_blank">http://gradsusr</a>!
 .org/mailman/listinfo/gradsusr<br>&gt; &gt;&gt;&gt;&gt;<br>&gt; &gt;&amp;g

t;&gt;&gt;<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt; --<br>&gt; &gt;&gt;&gt; Cordialement,<br>&gt; &gt;&gt;&gt; Best regards,<br>&gt; &gt;&gt;&gt; Rabah Hachelaf<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;      ____<br>
&gt; &gt;&gt;&gt;     (          )<br>&gt; &gt;&gt;&gt;    (            )<br>&gt; &gt;&gt;&gt;   (___ __)<br>&gt; &gt;&gt;&gt;      /////////<div class="im"><br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt; _______________________________________________<br>
&gt; &gt;&gt;&gt; gradsusr mailing list<br>&gt; &gt;&gt;&gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>&gt; &gt;&gt;&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;&gt;<br>&gt; &gt;&gt;<br>&gt; &gt;&gt;<br></div>&gt; &gt;&gt; --<br>&gt; &gt;&gt; *With Regards<br>&gt; &gt;&gt; Madhu Latha Akkisetti<br>&gt; &gt;&gt; SRF<br>&gt; &gt;&gt; NARL &lt;<a href="http://www.narl.gov.in" target="_blank">http://www.narl.gov.in</a>&gt;<br>
&gt; &gt;&gt; Department of space<br>&gt; &gt;&gt; Gadanki*<br>&gt; &gt;&gt;<br>&gt; &gt;&gt; ________!
 _______________________________________<div class="im"><br>&gt; &gt;&gt; gradsusr mailing list<br>&gt; &gt;&gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>&gt; &gt;&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; &gt;&gt;<br>&gt; &gt;&gt;<br>&gt; &gt;<br>&gt; &gt;<br></div>&gt; &gt; --<br>&gt; &gt; Cordialement,<br>&gt; &gt; Best regards,<br>&gt; &gt; Rabah Hachelaf<br>&gt; &gt;<br>&gt; &gt;      ____<br>&gt; &gt;     (          )<br>
&gt; &gt;    (            )<br>&gt; &gt;   (___ __)<br>&gt; &gt;      /////////<div class="im"><br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt; _______________________________________________<br>&gt; &gt; gradsusr mailing list<br>
&gt; &gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>&gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt; &gt;<br>&gt; &gt;<br>&gt; <br>&gt; <br></div>&gt; -- <br>&gt; *With Regards<br>&gt; Madhu Latha Akkisetti<br>&gt; SRF<br>&gt; NARL &lt;<a href="http://www.narl.gov.in" target="_blank">http://www.narl.gov.in</a>&gt;<br>
&gt; Department of space<br>&gt; Gadanki*<br>&gt; -------------- next part --------------<br>&gt; An HTML attachment !
 was scrubbed...<br>&gt; URL: <a href="http://gradsusr.org/pipermail/gradsusr/at" target="_blank">http://gradsusr.org/pipermail/gradsusr/at</a>

tachments/20121225/0f3c0bbd/attachment.html <br>&gt; <br>&gt; ------------------------------<div class="im"><br>&gt; <br>&gt; _______________________________________________<br>&gt; gradsusr mailing list<br>&gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt; <br>&gt; <br></div>&gt; End of gradsusr Digest, Vol 34, Issue 74<br>&gt; ****************************************<br>
</div>                                               </div></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">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>