[gradsusr] Error in using skew.gs

Amulya Chevuturi amulya.chevuturi at gmail.com
Tue Dec 25 23:52:13 EST 2012 

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-0003.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-0006.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-0003.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-0007.obj

More information about the gradsusr mailing list