[gradsusr] Error in using skew.gs
Jeff Duda
jeffduda319 at gmail.com
Mon Dec 24 12:51:03 EST 2012
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20121224/ddd5c89f/attachment-0003.html
More information about the gradsusr
mailing list