Re: skewt

Why do you have

;**************************************************
load "$NCARG_ROOT/share/ncarg/nclscripts/csm/gsn_code.ncl"
load "./skewt_func_1.ncl" ; <???????????
load "$NCARG_ROOT/share/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/share/ncarg/nclscripts/niwa/mapareas.ncl"
load "$NCARG_ROOT/share/ncarg/nclscripts/niwa/maplabels.ncl"
load "$NCARG_ROOT/share/ncarg/nclscripts/niwa/EcoConnectUtils.ncl"
;**************************************************

I think you are executing 'skewt_func_1.ncl' and **not** your script!

Just loading

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.ncl"

yields

Variable: skewt_bkgd
Type: obj
Total Size: 4 bytes
             1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
Number Of Attributes: 7
   Tlcl : 9.313728
   Plcl : 919.7243
   Shox : 1.125401
   Pwat : 2.222433
   Cape : 204.3857
   dataspec : skewt_test_new_data xyDataSpecClass 31
   data : skewt_test_new_data coordArraysClass 9

Variable: skewt_data
Type: obj
Total Size: 4 bytes
             1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
Number Of Attributes: 7
   data : skewt_test_new_data coordArraysClass 9
   dataspec : skewt_test_new_data xyDataSpecClass 31
   Cape : 204.3857
   Pwat : 2.222433
   Shox : 1.125401
   Plcl : 919.7243
   Tlcl : 9.313728
(0) 204.3857
(0) 204.3857

On 2/27/12 5:08 PM, Yang Yang wrote:
> Please run test.ncl with all the data files at the same directory.
>
> Thanks for help!
>
> Ed
> --- On Tue, 2/28/12, Dennis Shea<shea@ucar.edu> wrote:
>
>> From: Dennis Shea<shea@ucar.edu>
>> Subject: Re: skewt
>> To: "Yang Yang"<yang123yang@yahoo.com>
>> Date: Tuesday, February 28, 2012, 7:56 AM
>> yes
>>
>> On 2/27/12 4:54 PM, Yang Yang wrote:
>>> My data file is very small in size (only a few Ks). Can
>> I just attach them on a email to send them to you? Thanks!
>>>
>>> Ed
>>>
>>> --- On Tue, 2/28/12, Dennis Shea<shea@ucar.edu>
>> wrote:
>>>
>>>> From: Dennis Shea<shea@ucar.edu>
>>>> Subject: Re: skewt
>>>> To: "Yang Yang"<yang123yang@yahoo.com>
>>>> Date: Tuesday, February 28, 2012, 7:35 AM
>>>> http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP
>>>>
>>>> On 2/27/12 4:31 PM, Yang Yang wrote:
>>>>> Can you tell me how to upload the data and
>> files to
>>>> you? Thanks!
>>>>> Ed
>>>>>
>>>>> --- On Tue, 2/28/12, Dennis Shea<shea@ucar.edu>
>>>> wrote:
>>>>>
>>>>>> From: Dennis Shea<shea@ucar.edu>
>>>>>> Subject: Re: skewt
>>>>>> To: "Yang Yang"<yang123yang@yahoo.com>
>>>>>> Date: Tuesday, February 28, 2012, 7:19 AM
>>>>>> I wrote thye skewT script and donated
>>>>>> it to NCL.
>>>>>> No one else is going to answer.
>>>>>> D
>>>>>>
>>>>>> On 02/27/2012 04:17 PM, Yang Yang wrote:
>>>>>>> Thanks Dennis. I think that it is
>> better for me
>>>> to
>>>>>> repost my problem on ncl-talk to see if
>> anyone is
>>>> able to
>>>>>> help me.
>>>>>>>
>>>>>>> Ed
>>>>>>>
>>>>>>> --- On Tue, 2/28/12, Dennis Shea<shea@ucar.edu>
>>>>>> wrote:
>>>>>>>
>>>>>>>> From: Dennis Shea<shea@ucar.edu>
>>>>>>>> Subject: Re: skewt
>>>>>>>> To: "Yang Yang"<yang123yang@yahoo.com>
>>>>>>>> Date: Tuesday, February 28, 2012,
>> 7:09 AM
>>>>>>>> This is offline.
>>>>>>>>
>>>>>>>> You can send me your data and a
>> clean
>>>> scrip.
>>>>>>>> I will run here. We have 6.0.0
>> installed.
>>>>>>>>
>>>>>>>> I am very busy with my own work so
>> I'm not
>>>> sure
>>>>>>>> when I will be able to look
>>>>>>>>
>>>>>>>>
>>>>>>>> D
>>>>>>>>
>>>>>>>> On 2/27/12 4:02 PM, Yang Yang
>> wrote:
>>>>>>>>> Hi Dennis:
>>>>>>>>>
>>>>>>>>> Thank you very much for your
>> help!
>>>>>>>>>
>>>>>>>>> Sorry for my last email with a
>> wrong
>>>> subject. I
>>>>>> clicked
>>>>>>>> the wrong one to reply.
>>>>>>>>>
>>>>>>>>> I tried your suggestion. It
>> still does
>>>> not
>>>>>> work. Please
>>>>>>>> see my script below. My ncl version
>> is
>>>> V.5.2. The
>>>>>>>> information after the run is:
>>>>>>>>>
>>>>>>>>> warning:Attempt to reference
>> attribute
>>>> (Cape)
>>>>>> which is
>>>>>>>> undefined
>>>>>>>>> (0)
>> -999
>>>>>>>>> warning:Attempt to reference
>> attribute
>>>> (Cape)
>>>>>> which is
>>>>>>>> undefined
>>>>>>>>> (0)
>> -999
>>>>>>>>>
>>>>>>>>> Regards
>>>>>>>>> Ed
>>>>>>>>>
>>>>>>>>> ;*************************
>>>>>>>>> ;ncl script
>>>>>>>>> ;**************************
>>>>>>>>>
>> wks =
>>>> gsn_open_wks
>>>>>> ("x11",
>>>>>>>> "skewt_test_new")
>>>>>>>>>
>>>>>>>>> ; --- Create background skew-T;
>> plot
>>>> sounding
>>>>>>>> ---------------------
>>>>>>>>>
>>>>>>>>>
>> skewtOpts
>>>>>>>>
>>>> =
>>>>>> True
>>>>>>>>>
>>>> skewtOpts@DrawColAreaFill =
>>>>>> True
>>>>>>>> ;
>> default is False
>>>>>>>>>
>>>>>> skewtOpts@tiMainString =
>>>>>>>> "Raob Data; No Winds"
>>>>>>>>>
>>>>>>>>>
>>>> dataOpts =
>>>>>> True
>>>>>>>>
>> ; no
>>>> options
>>>>>>>>>
>>>>>>>>>
>>>> dataOpts@PrintZ = True
>>>>>>>>>
>> skewt_bkgd
>>>> =
>>>>>> skewT_BackGround (wks,
>>>>>>>> skewtOpts)
>>>>>>>>>
>> skewt_data
>>>> =
>>>>>>>>
>> skewT_PlotData (wks,
>>>> skewt_bkgd,
>>>>>> pp,tt,tc,h
>>>>>>>> \
>>>>>>>>>
>>>>>>>>
>>>>>>>>
>> ,
>>>> wsp,dir, dataOpts)
>>>>>>>>>
>> draw
>>>> (skewt_bkgd)
>>>>>>>>>
>> draw
>>>> (skewt_data)
>>>>>>>>>
>> frame(wks)
>>>>>>>>>
>>>>>>>>> print(skewt_bkgd@Cape)
>>>>>>>>> print(skewt_data@Cape)
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> --- On Mon, 2/27/12, Dennis
>> Shea<shea@ucar.edu>
>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>>> From: Dennis Shea<shea@ucar.edu>
>>>>>>>>>> Subject: Re:
>> skewt
>>>>>>>>>> To: "Yang Yang"<yang123yang@yahoo.com>
>>>>>>>>>> Date: Monday, February 27,
>> 2012,
>>>> 1:24 PM
>>>>>>>>>> Use
>>>>>>>>>>
>>>> skewt_bkgd@Cape
>>>>>>>>>>
>>>>>>>>>> On 2/26/12 6:15 PM, Yang
>> Yang
>>>> wrote:
>>>>>>>>>>> Hello:
>>>>>>>>>>>
>>>>>>>>>>> In the description of
>> skewt, it
>>>> shows:
>>>>>>>>>>> Example: Let skewT be
>> the
>>>> returned
>>>>>> object. The
>>>>>>>> data mat
>>>>>>>>>> be retrieved via:
>>>>>>>>>>>
>>>>>> cape =
>>>>>>>> skewT@Cape
>>>>>>>>>>>
>>>>>> pwat =
>>>>>>>> skewT@Pwat
>>>>>>>>>>>
>>>>>>>>>>> I wrote the following
>> in my
>>>> script as:
>>>>>>>>>>> ;...............
>>>>>>>>>>>
>>>> skewt_bkgd
>>>>>> =
>>>>>>>> skewT_BackGround (wks,
>>>>>>>>>> skewtOpts)
>>>>>>>>>>>
>>>> skewt_data
>>>>>> =
>>>>>>>>>>
>>>> skewT_PlotData (wks,
>>>>>> skewt_bkgd,
>>>>>>>> pp,tt,tc,h
>>>>>>>>>> \
>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>> ,
>>>>>> wsp,dir, dataOpts)
>>>>>>>>>>>
>>>> draw
>>>>>> (skewt_bkgd)
>>>>>>>>>>>
>>>> draw
>>>>>> (skewt_data)
>>>>>>>>>>>
>>>> frame(wks)
>>>>>>>>>>> print(skewt_data@Cape)
>>>>>>>>>>>
>> ;.........................
>>>>>>>>>>>
>>>>>>>>>>> However, when I run it,
>> I got
>>>> the
>>>>>> following
>>>>>>>>>> information:
>>>>>>>>>>> warning:Attempt to
>> reference
>>>> attribute
>>>>>> (Cape)
>>>>>>>> which is
>>>>>>>>>> undefined
>>>>>>>>>>> (0)
>>>> -999
>>>>>>>>>>>
>>>>>>>>>>> But the CAPE value was
>> shown in
>>>> the
>>>>>> diagram
>>>>>>>> generated
>>>>>>>>>> by my script.
>>>>>>>>>>>
>>>>>>>>>>> Can you tell me what is
>> wrong?
>>>> I use
>>>>>> NCL
>>>>>>>>>>>
>>>>>> Copyright
>> (C) 1995-2010
>>>>>>>> - All Rights
>>>>>>>>>> Reserved
>>>>>>>>>>>
>>>>>> University
>> Corporation
>>>>>>>> for Atmospheric
>>>>>>>>>> Research
>>>>>>>>>>>
>>>> NCAR
>>>>>> Command Language
>>>>>>>> Version 5.2.0
>>>>>>>>>>>
>>>> The use
>>>>>> of this
>>>>>>>> software is governed
>>>>>>>>>> by a License Agreement.
>>>>>>>>>>>
>>>> See http://www.ncl.ucar.edu/ for more
>>>>>>>>>> details.
>>>>>>>>>>>
>>>>>>>>>>> Thank you very much!
>>>>>>>>>>>
>>>>>>>>>>> Ed
>>>>>>>>>>>
>>>>>>>>
>>>> _______________________________________________
>>>>>>>>>>> ncl-talk mailing list
>>>>>>>>>>> List instructions,
>> subscriber
>>>> options,
>>>>>>>> unsubscribe:
>>>>>>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>>>>>
>>>>>>>>
>>>>>>
>>>>>>
>>>>
>>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Feb 27 20:41:38 2012