[gradsusr] Root Mean Square Error

Amulya Chevuturi amulya.chevuturi at gmail.com
Wed Dec 26 06:45:23 EST 2012 

Hello Tahereh

It is good that you understood. Even if you dont have predicted data and
observed data in one file, this would work by the following method.

open forecast.ctl
set t 1 4
set lon 85 100
set lat 0 10
define model=Model parameter
close 1
open observation.ctl
set t 1 4
set lon 85 100
set lat 0 10
define observation=Observation parameter
define erro=model-observation
define RMSE=sqrt(ave((error*error),t=1,t=4))
d RMSE

I have just sent you a tentative script. Change it according to your needs
and see if it works. Best of luck in your endeavor.

With regards,
Amulya


On Wed, Dec 26, 2012 at 5:02 PM, T. haghroosta <haghroosta at hotmail.com>wrote:

>  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 an! other 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 :
> > >
> &gt! ; >
>
> > > 'open tem.ctl'
> > > 'set t 1 4'
> > &g t; '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"
> > >
> > > >
> > &! gt; > 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>
> > > &gt! ; 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
> > > > >>>>
> > > > &gt! ;>>>
>
> > > > >>>> 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
> > > > >>>&! gt; 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
> > > > G! adanki*
> > > > -------------- next part --------------
> & gt; > > 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/grads!
> usr/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 s! end it again. It might be sent to the spam
> folder by
> > mistake.< br>>
> > 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
> > >>
> > &g! t;> 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
> > >&gt! ; Initial Time is 00Z08JAN2007
> > >> Create gif images as w ell (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
> > >> E! rror 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:
> > >>>
> &g! t; >>> 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, alt! hough 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
> > >>>&gt! ; 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.
> > >>>> D! ata 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 oc! curred on line 69
> > >>>> In file skew.gs
> > & gt;>>>
> > >>>> 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
> > >>>&g! t;> 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
> > >>>&gt! ;>> 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
> > >>>&g t;>> 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
> > >>>&g! t;>> 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)
> > >>>>>&gt! ; Jawaharlal Nehru University
> > >>>>>> New Del hi, 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! .
> > >>>>>>>>
> > >>>>&gt ;>>> 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 im! ages 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
> > >>&gt! ;>>>>>> 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
> > >>>>>>>>&! gt;> processing software. This is the problem that is
> cropping up.< br>> >>>>>>>>>>
> > >>>>>>>>>> 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
> > >&g! t;>>>>>>>> In file skew.gs
> > >>& gt;>>>>>>> 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 r! egards,
> > >>>>>>>>>> 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
> > &! gt;>>>>>>>> Graduate research assistant
> &gt ; >>>>>>>>> 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
> > >>>>>>>
> > >>>>>>>
> > >>>>>>
> > >>>&gt! ;>> _______________________________________________
>
> > > >>>>> 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
> > >>>>>
> > >>>>>
> > >>>>
> > >>&! gt;> _______________________________________________
>
> > >>>> 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 gradsus! r.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
> > ****************************************
> *
> **
>
> _______________________________________________
> 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/6a1f349c/attachment-0003.html

More information about the gradsusr mailing list