[gradsusr] Error in using skew.gs
Amulya Chevuturi
amulya.chevuturi at gmail.com
Fri Dec 21 01:13:39 EST 2012
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.gsafter 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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20121221/1a28246b/attachment-0003.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 20070108_d01_a1.ctl
Type: application/octet-stream
Size: 10153 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20121221/1a28246b/attachment-0006.obj
-------------- next part --------------
A non-text attachment was scrubbed...
Name: skew.gs
Type: application/octet-stream
Size: 67787 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20121221/1a28246b/attachment-0007.obj
More information about the gradsusr
mailing list