Hi Dennis:
Thanks for help! It works.
Two years ago, the ncl version I used didn't have skewt_func.ncl in the lib. I got skewt_func_1.ncl offline. I have been using Skewt that way since then. Thanks for the correction.
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>
> Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Date: Tuesday, February 28, 2012, 11:41 AM
> 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:
> [ncl-talk]
> >> 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 Tue Feb 28 13:40:53 2012
This archive was generated by hypermail 2.1.8 : Mon Mar 05 2012 - 14:12:07 MST