Hello Tahereh<div><br></div><div>First of all I made a slight mistake in my last mail. For RMSE calculations the term "A" 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 "olr" 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 "define error=(forecast-observation)". 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"><<a href="mailto:haghroosta@hotmail.com" target="_blank">haghroosta@hotmail.com</a>></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 "run <a href="http://test.gs" target="_blank">test.gs</a>", grads did nothin, niether error nor output!<br>
<br> <br>'open tem.ctl'<br>'set t 1 4'<br>'set lon 85 100'<br>'set lat 0 10'<br>'define RMSE=sqrt(ave(T*T,t=1,t=4))'<br>'d RMSE'<br><br> and also in another .gs file :<br> <br>
<br><div><div>'open tem.ctl'<br>'set t 1 4'<br>'set lon 85 100'<br>'set lat 0 10'<br>'define RMSE=sqrt(aave(lh*lh,lon=85,lon=100,lat=0,lat=10))'<br>'d RMSE'</div><div> </div>
<div>I also couldn'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>> ------------------------------<br>
> <br>> Message: 3<br>> Date: Tue, 25 Dec 2012 04:07:21 +0000<br>> From: "T.!
haghroosta" <<a href="mailto:haghroosta@hotmail.com" target="_blank">haghroosta@hotmail.com</a>><br>> Subject: [gradsusr] Root Mean Square Error<br>> To: <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
> Message-ID: <DUB121-W43F87358FE9C823BDD5BA2DD3A0@phx.gbl><br>> Content-Type: text/plain; charset="windows-1256"<div class="im"><br>> <br>> <br>> <br>> <br>> <br>> 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>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> 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>
> <br>> ------------------------------<br>> <br>> Message: 4<br>> Date: Tue, 25 Dec 2012 09:55:21 +0530<br>> From: Amulya Chevuturi <<a href="mailto:amulya.chevuturi@gmail.com" target="_blank">amulya.chevuturi@gmail.com</a>><br>
> Subject: Re: [gradsusr] Root Mean Square Error<br>> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>> Message-ID:<u></u>> <CAOV67LkiGd7ayjps3XTMW-3FWq5V=Ao4R4-oz7BUowZN0gkhqA@mail.g
<a href="http://mail.com" target="_blank">mail.com</a>><br>> Content-Type: text/plain; charset="iso-8859-1"<div class="im"><br>> <br>> Hello Tahereh,<br>> <br>> This is a website containing quite many scripts for various analysis. There<br>
> is a RMSE calculation script also called <a href="http://ltrend.gs" target="_blank">ltrend.gs</a>. Please look it and see<br>> if it is useful for your work.<br>> <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>
> <br>> Usually the RMSE calculation is done through the commands:<br>> <br>> "define RMSE=sqrt(ave((A*A),t=t1,t=tn)"<br>> Where,<br>> RMSE=Root Mean Square Error<br>> A=Is the parameter you want to calculate<br>
> t1=Time step 1 for you series<br>> tn=Last time step for your series<br>> This will give you the plot of a time period.<br>> <br>> OR<br>> <br>> "define RMSE=sqrt(aave((A*A),lon=lon1,lon=lonn,lat=lat1,lat=latn)"<br>
> Where<br>> RMSE=Root Mean Square Error<br>> A=Is the parameter you want to calculate<br></div>> Lat1, Latn, Lon1 and Lonn = Tel!
l the extent of your study<div class="im"><br>> This will give you a time series of RMSE.<br>> <br>> Hopefully this clears your doubt.<br>> <br>> With regards,<br>> Amulya Chevuturi<br>> School of Environmental Sciences<br>
> Jawaharlal Nehru University<br>> New Delhi<br>> <br>> On Tue, Dec 25, 2012 at 9:37 AM, T. haghroosta <<a href="mailto:haghroosta@hotmail.com" target="_blank">haghroosta@hotmail.com</a>>wrote:<br>> <br>
> > Dear All,<br>> ><br>> > Could you please tell me how I can calculate RMSE (root mean square error)<br>> > in Grads and draw relevant graphs?<br>> ><br>> ><br>> > Regards,<br>
> ><br>> > Tahereh<br>> ><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>> ><br></div>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>
> 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>> <br>> ------------------------------<br>> <br>> Message: 5<br>> Date: Tue, 25 Dec 2012 10:08:26 +0530<br>> From: madhu latha <<a href="mailto:madhulatha11@gmail.com" target="_blank">madhulatha11@gmail.com</a>><br>
> Subject: Re: [gradsusr] problem for prebufr data<br>> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>> Message-ID:<br>> <<a href="mailto:CAAhFJCSWRsxZMmPfNuZosfDgbzfJETwccC%2BsJxoSgMroxmWmTw@mail.gmail.com" target="_blank">CAAhFJCSWRsxZMmPfNuZosfDgbzfJETwccC+sJxoSgMroxmWmTw@mail.gmail.com</a>><br>
> Content-Type: text/plain; charset="windows-1252"<br>> <br>> hi<br>> <br>> Thank you<br>> First itself i have converted prebufr file to litle endian prebufr file<br>> using ssrc.exe<br>> then i have assimilated data using 3dvar.<br>
> Then mentioned problem came<br>> <br>> thank you.<br>> <br>> On Mon, Dec 24, 2012 at 11:56 PM, Rabah Hachelaf <<a href="mailto:r.hachelaf@gmail.com" target="_blank">r.hachelaf@gmail.com</a>>wrote:<br>
> <br>> > Hello,<br>> ><br>> > try to convert your prepbufr file from big endian to littele endian buy<br>> > the ssrc.c too!
l.<br>> ><br>> > ssrc.exe <Big Endian prepbufr file> Little Endian prepbufr file<br>> ><br>> ><br>> > 2012/12/24 madhu latha <<a href="mailto:madhulatha11@gmail.com" target="_blank">madhulatha11@gmail.com</a>><br>
> ><br>> >> hi<br>> >> earlier i have tried the suggestions in the link by setting the<br>> >> environment<br>> >> setenv F_UFMTENDIAN "little:94-99"<br>> >> but no use<br>
> >><br>> >> then itried to do reblock.c mentioned in the same link.<br>> >> but i did not understand how to use the utility reblock.c<br>> >><br>> >> kindly suggest.<br>> >><br>
> >><br>> >> On Mon, Dec 24, 2012 at 10:15 PM, Rabah Hachelaf <<a href="mailto:r.hachelaf@gmail.com" target="_blank">r.hachelaf@gmail.com</a>>wrote:<br>> >><br>> >>> Hello,<br>
> >>> see this link , it may be help you.<br>> >>><br>> >>> <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>
> >>><u></u>> >>> --<br>> >>> Cordialement,<br>> &gt
;>> Best regards,<br>> >>> Rabah Hachelaf<br>> >>><br>> >>> ____<br>> >>> ( )<br>> >>> ( )<br>> >>> (___ __)<br>> >>> /////////<br>
> >>><br>> >>><br>> >>> 2012/12/24 madhu latha <<a href="mailto:madhulatha11@gmail.com" target="_blank">madhulatha11@gmail.com</a>><br>> >>><br>> >>>> Dear all<br>
> >>>><br>> >>>><br>> >>>> Iam trying to run 3dvar on intel compiler<br>> >>>> i have successfully installed wrf3dvar<br>> >>>> but when iam trying to assimilate prepbufr data it is showing the<br>
> >>>> following error<br>> >>>><br>> >>>><br>> >>>> Error message<br>> >>>> **************BUFR ARCHIVE LIBRARY ABORT*****************<br>> >>>> BUFRLIB: OPENBF - ERROR READING INPUT FILE CONNECTED TO !
UNIT 96 WHEN<br>> >>>> CHECKING<br>> >>>> FOR ?BUFR? IN FIRST 4 BYTES OF RECORD<br>> >>>> **************BUFR ARCHIVE LIBRARY ABORT*****************<br>> >>>><br>
> >>>><br>> >>>> can anyone know how to fix the problem<br>> >>>> Any help will be appreciated.<br>> >>>><br>> >>>> Thank you.<br>> >>>><br>
> >>>><br>> >>>> --<br>> >>>> *With Regards<br>> >>>> Madhu Latha Akkisetti<br>> >>>> SRF<br>> >>>> NARL <<a href="http://www.narl.gov.in" target="_blank">http://www.narl.gov.in</a>><br>
> >>>> Department of space<br>> >>>> Gadanki*<div class="im"><br>> >>>><br>> >>>> _______________________________________________<br>> >>>> gradsusr mailing list<br>
> >>>> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br></div>> >>>> <a href="http://gradsusr" target="_blank">http://gradsusr</a>!
.org/mailman/listinfo/gradsusr<br>> >>>><br>> >&g
t;>><br>> >>><br>> >>><br>> >>> --<br>> >>> Cordialement,<br>> >>> Best regards,<br>> >>> Rabah Hachelaf<br>> >>><br>> >>> ____<br>
> >>> ( )<br>> >>> ( )<br>> >>> (___ __)<br>> >>> /////////<div class="im"><br>> >>><br>> >>><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>> >>><br>> >><br>> >><br></div>> >> --<br>> >> *With Regards<br>> >> Madhu Latha Akkisetti<br>> >> SRF<br>> >> NARL <<a href="http://www.narl.gov.in" target="_blank">http://www.narl.gov.in</a>><br>
> >> Department of space<br>> >> Gadanki*<br>> >><br>> >> ________!
_______________________________________<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>> >><br>> ><br>> ><br></div>> > --<br>> > Cordialement,<br>> > Best regards,<br>> > Rabah Hachelaf<br>> ><br>> > ____<br>> > ( )<br>
> > ( )<br>> > (___ __)<br>> > /////////<div class="im"><br>> ><br>> ><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>> ><br>> <br>> <br></div>> -- <br>> *With Regards<br>> Madhu Latha Akkisetti<br>> SRF<br>> NARL <<a href="http://www.narl.gov.in" target="_blank">http://www.narl.gov.in</a>><br>
> Department of space<br>> Gadanki*<br>> -------------- next part --------------<br>> An HTML attachment !
was scrubbed...<br>> URL: <a href="http://gradsusr.org/pipermail/gradsusr/at" target="_blank">http://gradsusr.org/pipermail/gradsusr/at</a>
tachments/20121225/0f3c0bbd/attachment.html <br>> <br>> ------------------------------<div class="im"><br>> <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>> <br></div>> End of gradsusr Digest, Vol 34, Issue 74<br>> ****************************************<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>