[gradsusr] Root Mean Square Error

T. haghroosta haghroosta at hotmail.com
Wed Dec 26 06:32:40 EST 2012 

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

More information about the gradsusr mailing list