[gradsusr] Error in using skew.gs

Wed Dec 19 05:47:39 EST 2012

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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20121219/94f79652/attachment-0003.html 