<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Dear Amulya,<BR> <BR>Thank you for your compelet explanations. Because I don't have predicted data and observated data in one file, I cannot do this right now. I will do it in future. I think this way will work. <BR> <BR>Regards,<BR> <BR>Tahereh<br> <BR><div><div id="SkyDrivePlaceholder"></div>> From: gradsusr-request@gradsusr.org<br>> Subject: gradsusr Digest, Vol 34, Issue 76<br>> To: gradsusr@gradsusr.org<br>> Date: Tue, 25 Dec 2012 23:52:21 -0500<br>> <br>> Send gradsusr mailing list submissions to<br>> gradsusr@gradsusr.org<br>> <br>> To subscribe or unsubscribe via the World Wide Web, visit<br>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> or, via email, send a message with subject or body 'help' to<br>> gradsusr-request@gradsusr.org<br>> <br>> You can reach the person managing the list at<br>> gradsusr-owner@gradsusr.org<br>> <br>> When replying, please edit your Subject line so it is more specific<br>> than "Re: Contents of gradsusr digest..."<br>> <br>> <br>> Today's Topics:<br>> <br>> 1. Re: Root Mean Square Error (Amulya Chevuturi)<br>> 2. Re: Error in using skew.gs (Amulya Chevuturi)<br>> <br>> <br>> ----------------------------------------------------------------------<br>> <br>> Message: 1<br>> Date: Wed, 26 Dec 2012 09:04:21 +0530<br>> From: Amulya Chevuturi <amulya.chevuturi@gmail.com><br>> Subject: Re: [gradsusr] Root Mean Square Error<br>> To: GrADS Users Forum <gradsusr@gradsusr.org><br>> Message-ID:<br>> <CAOV67LmHzuRqbatcdRobwo3ej_H_tH_qgNnQ0o3db8Me3rR7rg@mail.gmail.com><br>> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> Hello Tahereh<br>> <br>> First of all I made a slight mistake in my last mail. For RMSE calculations<br>> the term "A" in the formula should be not the parameter you want to<br>> calculate but the error of the parameter. That is A should be error value.<br>> Something like
model minus observation.<br>> <br>> To answer your question, as I am not familiar with you data. I would like<br>> you to download outgoing longwave radiation (OLR) data from this website.<br>> http://www.esrl.noaa.gov/psd/data/gridded/data.interp_OLR.html. Download<br>> the olr.daily.mean.nc. Now write the following lines in the terminal where<br>> you have downloaded and stored your OLR data.<br>> <br>> sdfopen olr.day.mean.nc<br>> set t 1 4<br>> set lon 85 100<br>> set lat 0 10<br>> define RMSE=sqrt(ave((olr*olr),t=1,t=4))<br>> d RMSE<br>> <br>> This should work. I have tried with the same dataset and the output has<br>> come. If this is successful, try the same with you dataset and see what the<br>> result is.<br>> Now the mistake in this method is the term "olr" is not error. RMSE formula<br>> is sqrt(ave(error*error),t=1,t=4)). So usually error must be added, not the<br>> parameter itself. Hence, you write another command before defining RMSE,<br>> where you define error by denoting "define error=(forecast-observation)".<br>> You can read this short webpage to understand RMSE.<br>> http://www.weatheroffice.gc.ca/verification/scores/rmse_e.html.<br>> <br>> With regards<br>> Amulya Chevuturi<br>> <br>> On Tue, Dec 25, 2012 at 12:42 PM, T. haghroosta <haghroosta@hotmail.com>wrote:<br>> <br>> > Dear Amulya,<br>> ><br>> > I wrote these commands as you said in a .gs file (test.gs) and after that<br>> > I tried to run the test.gs, I got nothing. I mean when I used the command<br>> > "run test.gs", 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>> > 'open tem.ctl'<br>> > 'set t 1 4'<br>> &g
t; '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'<br>> ><br>> > I also couldn't work with the sample scripts.<br>> ><br>> > Could you please lead me how to do?<br>> ><br>> > Regards,<br>> ><br>> > Tahereh<br>> ><br>> ><br>> > > ------------------------------<br>> > ><br>> > > Message: 3<br>> > > Date: Tue, 25 Dec 2012 04:07:21 +0000<br>> > > From: "T.! haghroosta" <haghroosta@hotmail.com><br>> > > Subject: [gradsusr] Root Mean Square Error<br>> > > To: <gradsusr@gradsusr.org><br>> > > Message-ID: <DUB121-W43F87358FE9C823BDD5BA2DD3A0@phx.gbl><br>> > > Content-Type: text/plain; charset="windows-1256"<br>> ><br>> > ><br>> > ><br>> > ><br>> > ><br>> > ><br>> > > Dear All, Could you please tell me how I can calculate RMSE (root mean<br>> > square error) in Grads and draw relevant graphs? Regards, Tahereh<br>> > > -------------- next part --------------<br>> > > An HTML attachment was scrubbed...<br>> > > URL:<br>> > http://gradsusr.org/pipermail/gradsusr/attachments/20121225/e2532241/attachment-0001.html<br>> > ><br>> > > ------------------------------<br>> > ><br>> > > Message: 4<br>> > > Date: Tue, 25 Dec 2012 09:55:21 +0530<br>> > > From: Amulya Chevuturi <amulya.chevuturi@gmail.com><br>> > > Subject: Re: [gradsusr] Root Mean Square Error<br>> > > To: GrADS Users Forum <gradsusr@gradsusr.org><br>> > > Message-ID:**><br>> > <CAOV67LkiGd7ayjps3XTMW-3FWq5V=Ao4R4-oz7BUowZN0gkhqA@mail.g mail.com><br>> > > Content-Type: text/plain; charset="iso-8859-1"<br>> ><br>> > ><br>> > > Hello Tahereh,<br>> > ><br>> > > This is a
website containing quite many scripts for various analysis.<br>> > There<br>> > > is a RMSE calculation script also called ltrend.gs. Please look it and<br>> > see<br>> > > if it is useful for your work.<br>> > > http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html<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>> > > Lat1, Latn, Lon1 and Lonn = Tel! l the extent of your study<br>> ><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 <haghroosta@hotmail.com<br>> > >wrote:<br>> > ><br>> > > > Dear All,<br>> > > ><br>> > > > Could you please tell me how I can calculate RMSE (root mean square<br>> > error)<br>> > > > in Grads and draw relevant graphs?<br>> > > ><br>> > > ><br>> > > > Regards,<br>> > > ><br>> > > > Tahereh<br>> > > ><br>> > > > ________________
_______________________________<br>> > > > gradsusr mailing list<br>> > > > gradsusr@gradsusr.org<br>> > > > http://gradsusr.org/mailman/listinfo/gradsusr<br>> > > ><br>> > > ><br>> > > -------------- next part --------------<br>> > > An HTML attachment was scrubbed...<br>> > > URL: http://grad!<br>> > susr.org/pipermail/gradsusr/attachments/20121225/c861c614/attachment-0001.html<br>> > ><br>> > > ------------------------------<br>> > ><br>> > > Message: 5<br>> > > Date: Tue, 25 Dec 2012 10:08:26 +0530<br>> > > From: madhu latha <madhulatha11@gmail.com><br>> > > Subject: Re: [gradsusr] problem for prebufr data<br>> > > To: GrADS Users Forum <gradsusr@gradsusr.org><br>> > > Message-ID:<br>> > > <CAAhFJCSWRsxZMmPfNuZosfDgbzfJETwccC+sJxoSgMroxmWmTw@mail.gmail.com><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 <r.hachelaf@gmail.com<br>> > >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 <madhulatha11@gmail.com><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 <<br>> > r.hachelaf@gmail.com>wrote:<br>> > > >><br>> > > >>> Hello,<br>> > > >>> see this link , it may be help you.<br>> > > >>><br>> > > >>> http://www.mmm.ucar.edu/wrf/users/wrfda/Docs/readBufr.htm<br>> > > >>>**> >>> --<br>> > > >>> Cordialement,<br>> > > &gt ;>> Best regards,<br>> > > >>> Rabah Hachelaf<br>> > > >>><br>> > > >>> ____<br>> > > >>> ( )<br>> > > >>> ( )<br>> > > >>> (___ __)<br>> > > >>> /////////<br>> > > >>><br>> > > >>><br>> > > >>> 2012/12/24 madhu latha <madhulatha11@gmail.com><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<br>> > 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 <http://www.narl.gov.in><br>> > > >>>> Department of space<br>> > > >>>> Gadanki*<br>> ><br>> > > >>>><br>> > > >>>> _______________________________________________<br>> > > >>>> gradsusr mailing list<br>> > > >>>> gradsusr@gradsusr.org<br>> > > >>>> http://gradsusr! .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>> > > >>> /////////<br>> ><br>> > > >>><br>> > > >
>><br>> > > >>> _______________________________________________<br>> > > >>> gradsusr mailing list<br>> > > >>> gradsusr@gradsusr.org<br>> > > >>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> > > >>><br>> > > >>><br>> > > >><br>> > > >><br>> > > >> --<br>> > > >> *With Regards<br>> > > >> Madhu Latha Akkisetti<br>> > > >> SRF<br>> > > >> NARL <http://www.narl.gov.in><br>> > > >> Department of space<br>> > > >> Gadanki*<br>> > > >><br>> > > >> ________! _______________________________________<br>> ><br>> > > >> gradsusr mailing list<br>> > > >> gradsusr@gradsusr.org<br>> > > >> http://gradsusr.org/mailman/listinfo/gradsusr<br>> > > >><br>> > > >><br>> > > ><br>> > > ><br>> > > > --<br>> > > > Cordialement,<br>> > > > Best regards,<br>> > > > Rabah Hachelaf<br>> > > ><br>> > > > ____<br>> > > > ( )<br>> > > > ( )<br>> > > > (___ __)<br>> > > > /////////<br>> ><br>> > > ><br>> > > ><br>> > > > _______________________________________________<br>> > > > gradsusr mailing list<br>> > > > gradsusr@gradsusr.org<br>> > > > http://gradsusr.org/mailman/listinfo/gradsusr<br>> > > ><br>> > > ><br>> > ><br>> > ><br>> > > --<br>> > > *With Regards<br>> > > Madhu Latha Akkisetti<br>> > > SRF<br>> > > NARL <http://www.narl.gov.in><br>> > > Department of space<br>> > > Gadanki*<br>> > > -------------- next part --------------<br>&
gt; > > An HTML attachment ! was scrubbed...<br>> > > URL: http://gradsusr.org/pipermail/gradsusr/attachments/20121225/0f3c0bbd/attachment.html<br>> > ><br>> > > ------------------------------<br>> ><br>> > ><br>> > > _______________________________________________<br>> > > gradsusr mailing list<br>> > > gradsusr@gradsusr.org<br>> > > http://gradsusr.org/mailman/listinfo/gradsusr<br>> > ><br>> > ><br>> > > End of gradsusr Digest, Vol 34, Issue 74<br>> > > ****************************************<br>> ><br>> > _______________________________________________<br>> > gradsusr mailing list<br>> > gradsusr@gradsusr.org<br>> > http://gradsusr.org/mailman/listinfo/gradsusr<br>> ><br>> ><br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20121226/4061926d/attachment-0001.html <br>> <br>> ------------------------------<br>> <br>> Message: 2<br>> Date: Wed, 26 Dec 2012 10:22:13 +0530<br>> From: Amulya Chevuturi <amulya.chevuturi@gmail.com><br>> Subject: Re: [gradsusr] Error in using skew.gs<br>> To: GrADS Users Forum <gradsusr@gradsusr.org><br>> Message-ID:<br>> <CAOV67LnHv1FsWgwvRaCAageYZx1nikFDb+FotYoXpE2+sM1LqA@mail.gmail.com><br>> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> Hello Mr. Jeff<br>> <br>> I have mailed you the data (.dat) via a website that helps send larger<br>> files. The file sent you was a similar but different .dat file as my<br>> original was too large to be sent. Thus, I have mailed the data to your<br>> email id directly rather than the gradsuser mailing list. I hope this was<br>> not a problem. Please let me know if you have received the mail with the<br>> data or do I need to send it again. It might be sent to the spam folder by<br>> mistake.<
br>> <br>> In this mail I am attaching my skew.gs file, the .ctl file to go with the<br>> .dat file sent to you and error text that is showing in the terminal when I<br>> run skew.gs file with this output file.<br>> <br>> Thank you for all your help.<br>> <br>> With regards,<br>> Amulya Chevuturi<br>> <br>> <br>> <br>> <br>> <br>> On Mon, Dec 24, 2012 at 11:21 PM, Jeff Duda <jeffduda319@gmail.com> wrote:<br>> <br>> > Amulya,<br>> > Can you send me your data (or put it on a server somewhere for me to<br>> > obtain)? I would like to try using yours and my version of the script<br>> > using your data. If I find something that works, I will let you know.<br>> ><br>> > Jeff<br>> ><br>> ><br>> > On Mon, Dec 24, 2012 at 1:50 AM, Amulya Chevuturi <<br>> > amulya.chevuturi@gmail.com> wrote:<br>> ><br>> >> Dear Grads users<br>> >><br>> >> You suggestions are quite accurate. The td values undefined for Lev 1000.<br>> >> This is because of the high terrain area in the region. So model is not<br>> >> simulating this region. The values are undefined even if in the ARWpost I<br>> >> add "extrapolate" as "true". Do I need to add another term "interpolate"? I<br>> >> have not come across this term in namelist. ARWpost.<br>> >> Also to answer Mr. Mahakur question, the options byteswapped is required<br>> >> for the ctl file to run. I have made the temperature profiles and am<br>> >> attaching herewith.<br>> >><br>> >> I have added the Lines Mr. Jeff suggested and this is the result I have<br>> >> obtained.<br>> >> ga-> skew.gs<br>> >> No hardcopy metafile open<br>> >> All files closed; all defined objects released;<br>> >> All GrADS attributes have been reinitialized<br>> >> Initial Time is 00Z08JAN2007<br>> >> Create gif images as w
ell (1=yes ; 0=no)<br>> >> 0<br>> >> Please enter coordinated for sounding<br>> >><br>> >> Latitude values between: 4.96654 and 50.0419<br>> >> Enter Latitude<br>> >> 33<br>> >><br>> >> Longitude values between: 48.3967 and 105.603<br>> >> Enter Longitude<br>> >> 78<br>> >><br>> >> Overlay two time periods 1=yes ; 0=no<br>> >> 0<br>> >><br>> >> Please enter time to plot<br>> >> Available times: 1 to 17<br>> >> 5<br>> >> Drawing temperature sounding.<br>> >> Drawing dewpoint sounding.<br>> >> Printing Grid -- 1 Values -- Undef = -9.99e+08<br>> >> -9.99e+08<br>> >><br>> >> 1000<br>> >> getrh<br>> >> 13.743<br>> >><br>> >> 1000<br>> >> 0.00992721580569856<br>> >> Non-numeric args to numeric operation<br>> >> Error occurred on line 2176<br>> >> In file skew.gs<br>> >> Error occurred on line 2432<br>> >> In file skew.gs<br>> >> Error occurred on line 1109<br>> >> In file skew.gs<br>> >> Error occurred on line 69<br>> >> In file skew.gs<br>> >><br>> >> Thank you so much for all your help. Please let me know how to proceed.<br>> >> As you told Mr. Jeff even if there is no interpolation, Does this problem<br>> >> still occur?<br>> >><br>> >><br>> >><br>> >> On Sat, Dec 22, 2012 at 12:06 AM, Jeff Duda <jeffduda319@gmail.com>wrote:<br>> >><br>> >>> Okay, modify that section of code, which was<br>> >>><br>> >>> Sfctemp=subwrd(rec,4)<br>> >>> "d "snddewp<br>> >>> say result<br>> >>><br>> >>> to the following:<br>> >>><br>> >>> Sfctemp=subwrd(rec,4)<br>> >>> 'set
gxout print'<br>> >>> "d "snddewp<br>> >>> say result<br>> >>><br>> >>> as I'm not sure what that output means either. My suspicion is that the<br>> >>> dewpoint value it's getting is undefined for some reason. Are you trying<br>> >>> to plot a sounding in an area with elevated terrain? When you<br>> >>> postprocessed with ARWpost, did you set the interpolate option to .true.?<br>> >>> If you try to display the td field with lev = 1000 and some areas are<br>> >>> undefined, then you didn't set ARWpost to interpolate, but that shouldn't<br>> >>> be the problem. It sounds like you're trying to run this the exact same<br>> >>> way I did, as I've also run this code with ARWpost processed WRF output.<br>> >>> If we can't find a solution, I can send you a copy of my version of the<br>> >>> script for you to try, although I've already looked and it's not<br>> >>> significantly different from what you have.<br>> >>><br>> >>> Jeff<br>> >>><br>> >>> On Fri, Dec 21, 2012 at 12:13 AM, Amulya Chevuturi <<br>> >>> amulya.chevuturi@gmail.com> wrote:<br>> >>><br>> >>>> Hello Jeff,<br>> >>>><br>> >>>> Thank you so much for the help.<br>> >>>><br>> >>>> I added your suggestions to skew.gs. The SfcPlev is showing 1000,<br>> >>>> which I dont think is wrong. But the result of Sfcdewp is something I am<br>> >>>> not understanding. This is the following result I am getting on running<br>> >>>> skew.gs after additions made of "say result" and "say SfcPlev".<br>> >>>><br>> >>>> ga-> skew.gs<br>> >>>> No hardcopy metafile open<br>> >>>> All files closed; all defined objects released;<br>> >>>
> All GrADS attributes have been reinitialized<br>> >>>> Initial Time is 00Z08JAN2007<br>> >>>> Create gif images as well (1=yes ; 0=no)<br>> >>>> 0<br>> >>>> Please enter coordinated for sounding<br>> >>>><br>> >>>> Latitude values between: 4.96654 and 50.0419<br>> >>>> Enter Latitude<br>> >>>> 33<br>> >>>><br>> >>>> Longitude values between: 48.3967 and 105.603<br>> >>>> Enter Longitude<br>> >>>> 78<br>> >>>><br>> >>>> Overlay two time periods 1=yes ; 0=no<br>> >>>> 0<br>> >>>><br>> >>>> Please enter time to plot<br>> >>>> Available times: 1 to 17<br>> >>>> 5<br>> >>>> Drawing temperature sounding.<br>> >>>> Drawing dewpoint sounding.<br>> >>>> Data Type = grid<br>> >>>> Dimensions = -1 -1<br>> >>>> I Dimension = -999 to -999<br>> >>>> J Dimension = -999 to -999<br>> >>>> Sizes = 1 1 1<br>> >>>> Undef value = -9.99e+08<br>> >>>> Undef count = 1 Valid count = 0<br>> >>>> Min, Max = 1e+30 1e+30<br>> >>>><br>> >>>> 1000<br>> >>>> getrh<br>> >>>> 13.743<br>> >>>> 1e+30<br>> >>>> 1000<br>> >>>> 0.00992721580569856<br>> >>>> 286.893<br>> >>>> 43.385364499893<br>> >>>> 1573660000<br>> >>>> Non-numeric args to numeric operation<br>> >>>> Error occurred on line 2437<br>> >>>> In file skew.gs<br>> >>>> Error occurred on line 1108<br>> >>>> In file skew.gs<br>> >>>> Error occurred on line 69<br>> >>>> In file skew.gs<br>> &
gt;>>><br>> >>>> I have also attached the new skew.gs file along with my ctl file. I<br>> >>>> rechecked my ctl file, the file does have variable td. How am I to proceed<br>> >>>> further.<br>> >>>><br>> >>>> Thank you so much for taking the time to helping me out.<br>> >>>><br>> >>>> With regards,<br>> >>>> Amulya Chevuturi<br>> >>>><br>> >>>><br>> >>>><br>> >>>> On Wed, Dec 19, 2012 at 10:40 PM, Jeff Duda <jeffduda319@gmail.com>wrote:<br>> >>>><br>> >>>>> As you can see from your output, the script has the value 1e+30 stored<br>> >>>>> for dewpoint. That is clearly not correct. Therefore, there is a problem<br>> >>>>> with the variable dewp getting passed into this function, which is named<br>> >>>>> Thetae. The other error output tells you where the call to that function<br>> >>>>> came from (line 1106). That line is<br>> >>>>><br>> >>>>> SfcThee=Thetae(Sfctemp,Sfcdewp,SfcPlev)<br>> >>>>><br>> >>>>> So there must be a problem with the variable Sfcdewp outside of the<br>> >>>>> Thetae function. The variable Sfcdewp is computed in the line above 1106.<br>> >>>>> Clearly something is not right there. The 10 or so lines above this are<br>> >>>>> where the determination is made of how to obtain the surface dewpoint. See<br>> >>>>> the line "d "snddewp? I would insert a 'say result' command directly after<br>> >>>>> that command and see what the output shows. Also insert the command say<br>> >>>>> SfcPlev in that same area. The error could be from either source.<br>> >>>>><br>> >>>>> Lastly, make
sure that you have the variable td in your control file.<br>> >>>>> I'm pretty sure I saw it there when I looked at it a few days ago, but I no<br>> >>>>> longer have that file, so I can't check. If you don't have that, that is<br>> >>>>> the source of the problem.<br>> >>>>><br>> >>>>> Jeff<br>> >>>>><br>> >>>>><br>> >>>>> On Wed, Dec 19, 2012 at 4:47 AM, Amulya Chevuturi <<br>> >>>>> amulya.chevuturi@gmail.com> wrote:<br>> >>>>><br>> >>>>>> Thank you so much for the help that you are providing me MR. Jeff.<br>> >>>>>><br>> >>>>>> So as you suggested, I added the say commands for the line 2424. This<br>> >>>>>> is for the script variable mix. What is the term "getrh" in that? Is it a<br>> >>>>>> function or a variable. As when I input say getrh the output is getrh. So I<br>> >>>>>> added these lines in my gs file before the script for mix variable:<br>> >>>>>><br>> >>>>>> say getrh<br>> >>>>>> say temp<br>> >>>>>> say dewp<br>> >>>>>> say pres<br>> >>>>>> say ws<br>> >>>>>><br>> >>>>>><br>> >>>>>> This is the output I am getting now.<br>> >>>>>><br>> >>>>>> ga-> skew.gs<br>> >>>>>> No hardcopy metafile open<br>> >>>>>> All files closed; all defined objects released;<br>> >>>>>> All GrADS attributes have been reinitialized<br>> >>>>>> Initial Time is 00Z08JAN2007<br>> >>>>>> Create gif images as well (1=yes ; 0=no)<br>> >>>>>> 0<br>> >>>&g
t;>> Please enter coordinated for sounding<br>> >>>>>><br>> >>>>>> Latitude values between: 4.96654 and 50.0419<br>> >>>>>> Enter Latitude<br>> >>>>>> 33<br>> >>>>>><br>> >>>>>> Longitude values between: 48.3967 and 105.603<br>> >>>>>> Enter Longitude<br>> >>>>>> 78<br>> >>>>>><br>> >>>>>> Overlay two time periods 1=yes ; 0=no<br>> >>>>>> 0<br>> >>>>>><br>> >>>>>> Please enter time to plot<br>> >>>>>> Available times: 1 to 17<br>> >>>>>> 5<br>> >>>>>> Drawing temperature sounding.<br>> >>>>>> Drawing dewpoint sounding.<br>> >>>>>> getrh<br>> >>>>>> 13.743<br>> >>>>>> 1e+30<br>> >>>>>> 1000<br>> >>>>>> 0.00992721580569856<br>> >>>>>> 286.893<br>> >>>>>> 43.385364499893<br>> >>>>>> 1573660000<br>> >>>>>> Non-numeric args to numeric operation<br>> >>>>>> Error occurred on line 2435<br>> >>>>>> In file skew.gs<br>> >>>>>> Error occurred on line 1106<br>> >>>>>> In file skew.gs<br>> >>>>>> Error occurred on line 69<br>> >>>>>> In file skew.gs<br>> >>>>>><br>> >>>>>> Any help would be most appreciated.<br>> >>>>>><br>> >>>>>> With regards,<br>> >>>>>> Amulya Chevuturi<br>> >>>>>> Junior Research Fellow (Ph.D. Student)<br>> >>>>>> Jawaharlal Nehru University<br>> >>>>>> New Del
hi, India<br>> >>>>>><br>> >>>>>><br>> >>>>>> On Tue, Dec 18, 2012 at 11:45 PM, Jeff Duda <jeffduda319@gmail.com>wrote:<br>> >>>>>><br>> >>>>>>> Amulya,<br>> >>>>>>> Thank you for adding the new information. I can tell that the value<br>> >>>>>>> of the script variable mix, which is given as 1573660000, cannot be right.<br>> >>>>>>> Check line 2424 of your script code. That is where the variable mix is<br>> >>>>>>> computed. Insert say statements for all the component variables there.<br>> >>>>>>> With this kind of investigation you will eventually be able to discover the<br>> >>>>>>> source of the error. My gut tells me one of the input variables you<br>> >>>>>>> supplied to the script in the main function call plotskew() is not what the<br>> >>>>>>> script is expecting it to be. I'm guessing it is the dewpoint variable.<br>> >>>>>>> But I saw your control file and it appears there should not be a problem<br>> >>>>>>> with that. Nonetheless, insert say statements to find out what the values<br>> >>>>>>> are, then go from there. If you need help knowing what to do next or where<br>> >>>>>>> to look, please ask.<br>> >>>>>>><br>> >>>>>>> Jeff<br>> >>>>>>><br>> >>>>>>><br>> >>>>>>> On Tue, Dec 18, 2012 at 1:21 AM, Amulya Chevuturi <<br>> >>>>>>> amulya.chevuturi@gmail.com> wrote:<br>> >>>>>>><br>> >>>>>>>> Thank you very much for replying Mr. Jeff.<br>> >>>>>>>><br>> >>>>>
;>>> I have added the lines as stated. In the output all of them are<br>> >>>>>>>> showing numeric values only. So I cannot make out what the problem is.<br>> >>>>>>>> Could you or anyone in the grads user list help me out further.<br>> >>>>>>>><br>> >>>>>>>> This the output I am getting at the terminal when I am running the<br>> >>>>>>>> skew.gs with the changes you suggested.<br>> >>>>>>>><br>> >>>>>>>> ga-> skew.gs<br>> >>>>>>>> No hardcopy metafile open<br>> >>>>>>>> All files closed; all defined objects released;<br>> >>>>>>>> All GrADS attributes have been reinitialized<br>> >>>>>>>> Initial Time is 00Z08JAN2007<br>> >>>>>>>> Create gif images as well (1=yes ; 0=no)<br>> >>>>>>>> 0<br>> >>>>>>>> Please enter coordinated for sounding<br>> >>>>>>>><br>> >>>>>>>> Latitude values between: 4.96654 and 50.0419<br>> >>>>>>>> Enter Latitude<br>> >>>>>>>> 33<br>> >>>>>>>><br>> >>>>>>>> Longitude values between: 48.3967 and 105.603<br>> >>>>>>>> Enter Longitude<br>> >>>>>>>> 78<br>> >>>>>>>><br>> >>>>>>>> Overlay two time periods 1=yes ; 0=no<br>> >>>>>>>> 0<br>> >>>>>>>><br>> >>>>>>>> Please enter time to plot<br>> >>>>>>>> Available times: 1 to 17<br>> >>>>>>>> 5<br>> >>>>>>>> Drawing temperature sounding.<b
r>> >>>>>>>> Drawing dewpoint sounding.<br>> >>>>>>>> 286.893<br>> >>>>>>>> 43.385364499893<br>> >>>>>>>> 1573660000<br>> >>>>>>>> Non-numeric args to numeric operation<br>> >>>>>>>> Error occurred on line 2430<br>> >>>>>>>> In file skew.gs<br>> >>>>>>>> Error occurred on line 1106<br>> >>>>>>>> In file skew.gs<br>> >>>>>>>> Error occurred on line 69<br>> >>>>>>>> In file skew.gs<br>> >>>>>>>><br>> >>>>>>>> I am attaching the skew.gs file with the changes added. Please let<br>> >>>>>>>> me know how I can further correct the error. Thank you for the help in<br>> >>>>>>>> advance.<br>> >>>>>>>><br>> >>>>>>>> With regards,<br>> >>>>>>>> Amulya Chevuturi<br>> >>>>>>>> Junior Research Fellow (Ph.D. Student)<br>> >>>>>>>> Jawaharlal Nehru University<br>> >>>>>>>> New Delhi, India<br>> >>>>>>>><br>> >>>>>>>> On Tue, Dec 18, 2012 at 1:06 AM, Jeff Duda <jeffduda319@gmail.com>wrote:<br>> >>>>>>>><br>> >>>>>>>>> Amulya,<br>> >>>>>>>>> I looked through your script and I can't find the cause of the<br>> >>>>>>>>> error. I've used this script on WRF-ARW output myself before (also<br>> >>>>>>>>> postprocessed using ARWpost). I remember having some issues getting the<br>> >>>>>>>>> script to run properly the first time, but I
got around that somehow. The<br>> >>>>>>>>> errors indicate that in this line<br>> >>>>>>>>><br>> >>>>>>>>> thetae=theta*exp((3.376/TLcl-0.00254)*mix*1.0+0.00081*mix)<br>> >>>>>>>>><br>> >>>>>>>>> one of theta, TLcl, or mix is not being computed correctly. I'm<br>> >>>>>>>>> not sure which of those it is, although you could put SAY statements just<br>> >>>>>>>>> above that line for each of the variables to see which it is. What I mean<br>> >>>>>>>>> is, insert these lines<br>> >>>>>>>>> say theta<br>> >>>>>>>>> say TLcl<br>> >>>>>>>>> say mix<br>> >>>>>>>>> Above the line I copied above. Then see what values get spit<br>> >>>>>>>>> out. One of those is either a NaN of some sort or has characters where<br>> >>>>>>>>> there should be numbers.<br>> >>>>>>>>><br>> >>>>>>>>><br>> >>>>>>>>> Jeff Duda<br>> >>>>>>>>><br>> >>>>>>>>> On Mon, Dec 17, 2012 at 12:56 AM, Amulya Chevuturi <<br>> >>>>>>>>> amulya.chevuturi@gmail.com> wrote:<br>> >>>>>>>>><br>> >>>>>>>>>> Dear GrADS users,<br>> >>>>>>>>>><br>> >>>>>>>>>> I have recently tried to use the skew.gs script to plot<br>> >>>>>>>>>> Skew-T/Log-P graphs for my WRF (ARW) output. I use ARWPost as my post<br>> >>>>>>>>>> processing software. This is the problem that is cropping up.<
br>> >>>>>>>>>><br>> >>>>>>>>>> ga-> skew.gs<br>> >>>>>>>>>> No hardcopy metafile open<br>> >>>>>>>>>> All files closed; all defined objects released;<br>> >>>>>>>>>> All GrADS attributes have been reinitialized<br>> >>>>>>>>>> Initial Time is 00Z08JAN2007<br>> >>>>>>>>>> Create gif images as well (1=yes ; 0=no)<br>> >>>>>>>>>> 0<br>> >>>>>>>>>> Please enter coordinated for sounding<br>> >>>>>>>>>><br>> >>>>>>>>>> Latitude values between: 4.96654 and 50.0419<br>> >>>>>>>>>> Enter Latitude<br>> >>>>>>>>>> 33<br>> >>>>>>>>>><br>> >>>>>>>>>> Longitude values between: 48.3967 and 105.603<br>> >>>>>>>>>> Enter Longitude<br>> >>>>>>>>>> 78<br>> >>>>>>>>>><br>> >>>>>>>>>> Overlay two time periods 1=yes ; 0=no<br>> >>>>>>>>>> 0<br>> >>>>>>>>>><br>> >>>>>>>>>> Please enter time to plot<br>> >>>>>>>>>> Available times: 1 to 17<br>> >>>>>>>>>> 5<br>> >>>>>>>>>> Drawing temperature sounding.<br>> >>>>>>>>>> Drawing dewpoint sounding.<br>> >>>>>>>>>> Non-numeric args to numeric operation<br>> >>>>>>>>>> Error occurred on line 2427<br>> >>>>>>>>>> In file skew.gs<br>> >>&
gt;>>>>>>> Error occurred on line 1106<br>> >>>>>>>>>> In file skew.gs<br>> >>>>>>>>>> Error occurred on line 69<br>> >>>>>>>>>> In file skew.gs<br>> >>>>>>>>>><br>> >>>>>>>>>> Though I think it is plotting the temperature and dewpoint. The<br>> >>>>>>>>>> rest of the parameters in the line 69 is not being plotted due to some<br>> >>>>>>>>>> errors in parameters. I am attaching my ctl file and skew.gs.<br>> >>>>>>>>>> Please let me know how to correct the problem.<br>> >>>>>>>>>><br>> >>>>>>>>>> Thanking you in advance.<br>> >>>>>>>>>><br>> >>>>>>>>>> With regards,<br>> >>>>>>>>>> Amulya Chevuturi<br>> >>>>>>>>>> Junior Research Fellow (Ph.D. Student)<br>> >>>>>>>>>> Jawaharlal Nehru University<br>> >>>>>>>>>> New Delhi, India<br>> >>>>>>>>>><br>> >>>>>>>>>> _______________________________________________<br>> >>>>>>>>>> gradsusr mailing list<br>> >>>>>>>>>> gradsusr@gradsusr.org<br>> >>>>>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> >>>>>>>>>><br>> >>>>>>>>>><br>> >>>>>>>>><br>> >>>>>>>>><br>> >>>>>>>>> --<br>> >>>>>>>>> Jeff Duda<br>> >>>>>>>>> Graduate research assistant<br>>
; >>>>>>>>> University of Oklahoma School of Meteorology<br>> >>>>>>>>> Center for Analysis and Prediction of Storms<br>> >>>>>>>>><br>> >>>>>>>>><br>> >>>>>>>>> _______________________________________________<br>> >>>>>>>>> gradsusr mailing list<br>> >>>>>>>>> gradsusr@gradsusr.org<br>> >>>>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> >>>>>>>>><br>> >>>>>>>>><br>> >>>>>>>><br>> >>>>>>>> _______________________________________________<br>> >>>>>>>> gradsusr mailing list<br>> >>>>>>>> gradsusr@gradsusr.org<br>> >>>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> >>>>>>>><br>> >>>>>>>><br>> >>>>>>><br>> >>>>>>><br>> >>>>>>> --<br>> >>>>>>> Jeff Duda<br>> >>>>>>> Graduate research assistant<br>> >>>>>>> University of Oklahoma School of Meteorology<br>> >>>>>>> Center for Analysis and Prediction of Storms<br>> >>>>>>><br>> >>>>>>><br>> >>>>>>> _______________________________________________<br>> >>>>>>> gradsusr mailing list<br>> >>>>>>> gradsusr@gradsusr.org<br>> >>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> >>>>>>><br>> >>>>>>><br>> >>>>>><br>> >>>>>> _______________________________________________<br>> >
>>>>> gradsusr mailing list<br>> >>>>>> gradsusr@gradsusr.org<br>> >>>>>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> >>>>>><br>> >>>>>><br>> >>>>><br>> >>>>><br>> >>>>> --<br>> >>>>> Jeff Duda<br>> >>>>> Graduate research assistant<br>> >>>>> University of Oklahoma School of Meteorology<br>> >>>>> Center for Analysis and Prediction of Storms<br>> >>>>><br>> >>>>><br>> >>>>> _______________________________________________<br>> >>>>> gradsusr mailing list<br>> >>>>> gradsusr@gradsusr.org<br>> >>>>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> >>>>><br>> >>>>><br>> >>>><br>> >>>> _______________________________________________<br>> >>>> gradsusr mailing list<br>> >>>> gradsusr@gradsusr.org<br>> >>>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> >>>><br>> >>>><br>> >>><br>> >>><br>> >>> --<br>> >>> Jeff Duda<br>> >>> Graduate research assistant<br>> >>> University of Oklahoma School of Meteorology<br>> >>> Center for Analysis and Prediction of Storms<br>> >>><br>> >>> _______________________________________________<br>> >>> gradsusr mailing list<br>> >>> gradsusr@gradsusr.org<br>> >>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> >>><br>> >>><br>> >><br>> >> _______________________________________________<br>> >> gradsusr mailing list<br>> >> gradsusr@gradsusr.org<br>> >> http://gradsusr.org/mailman/listinfo/gradsusr<b
r>> >><br>> >><br>> ><br>> ><br>> > --<br>> > Jeff Duda<br>> > Graduate research assistant<br>> > University of Oklahoma School of Meteorology<br>> > Center for Analysis and Prediction of Storms<br>> ><br>> > _______________________________________________<br>> > gradsusr mailing list<br>> > gradsusr@gradsusr.org<br>> > http://gradsusr.org/mailman/listinfo/gradsusr<br>> ><br>> ><br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20121226/bb6baa45/attachment.html <br>> -------------- next part --------------<br>> A non-text attachment was scrubbed...<br>> Name: 20120911_d03.ctl<br>> Type: application/octet-stream<br>> Size: 8196 bytes<br>> Desc: not available<br>> Url : http://gradsusr.org/pipermail/gradsusr/attachments/20121226/bb6baa45/attachment.obj <br>> -------------- next part --------------<br>> A non-text attachment was scrubbed...<br>> Name: Error.rtf<br>> Type: application/rtf<br>> Size: 1417 bytes<br>> Desc: not available<br>> Url : http://gradsusr.org/pipermail/gradsusr/attachments/20121226/bb6baa45/attachment.rtf <br>> -------------- next part --------------<br>> A non-text attachment was scrubbed...<br>> Name: skew.gs<br>> Type: application/octet-stream<br>> Size: 67805 bytes<br>> Desc: not available<br>> Url : http://gradsusr.org/pipermail/gradsusr/attachments/20121226/bb6baa45/attachment-0001.obj <br>> <br>> ------------------------------<br>> <br>> _______________________________________________<br>> gradsusr mailing list<br>> gradsusr@gradsusr.org<br>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> <br>> <br>> End of gradsusr Digest, Vol 34, Issue 76<br>> ****************************************<br></div> </div></body>
</html>