[ncl-talk] ncl-talk Digest, Vol 176, Issue 22
David Adams
dave.k.adams at gmail.com
Thu Jul 19 13:14:49 MDT 2018
Thanks Adam,
I have always used Kerry Emanuel´s code for my CAPE work, so I will check
this out.
saludos,
Dave
On Thu, Jul 19, 2018 at 12:48 PM Adam Phillips <asphilli at ucar.edu> wrote:
> Hi Dave,
> See the attachment containing capethermo at the bottom of a previous
> ncl-talk thread here:
> http://www.ncl.ucar.edu/Support/talk_archives/2007/0535.html
> If you have further questions please post them to the ncl-talk email list.
> Adam
>
> On Thu, Jul 19, 2018 at 11:29 AM David Adams <dave.k.adams at gmail.com>
> wrote:
>
>> Hi Tanya et al.
>> related to Soma´s question. Where can I find the actual code that ncl
>> uses to calculate CAPE, LCL, etc?
>> I use the ncl skewT plot all of the time, but Cape depends critically on
>> which parcel is lifted and whether it´s a reversible or pseudoadiabatic
>> process etc.... So it would be good to know how it´s being calcuated.
>>
>> thanks,
>> Dave
>>
>> On Thu, Jul 19, 2018 at 12:11 PM Tanya Peevey <trpeevey at gmail.com> wrote:
>>
>>> Soma,
>>>
>>> There is a function in NCL to plot a skew-T diagram and it is
>>> skewT_PlotData.
>>>
>>> Tanya
>>>
>>>
>>> On Wed, Jul 18, 2018 at 12:00 PM <ncl-talk-request at ucar.edu> wrote:
>>>
>>>> Send ncl-talk mailing list submissions to
>>>> ncl-talk at ucar.edu
>>>>
>>>> To subscribe or unsubscribe via the World Wide Web, visit
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>> or, via email, send a message with subject or body 'help' to
>>>> ncl-talk-request at ucar.edu
>>>>
>>>> You can reach the person managing the list at
>>>> ncl-talk-owner at ucar.edu
>>>>
>>>> When replying, please edit your Subject line so it is more specific
>>>> than "Re: Contents of ncl-talk digest..."
>>>>
>>>>
>>>> Today's Topics:
>>>>
>>>> 1. Re: Save in netcdf (Dennis Shea)
>>>> 2. Re: Save in netcdf (Soma Roy)
>>>> 3. Re: Save in netcdf (Be?ta Szab?-Tak?cs)
>>>> 4. Re: Save in netcdf (Soma Roy)
>>>> 5. Vertical wind shear calculation within 400 km radius of TC
>>>> (Md. Jalal Uddin)
>>>> 6. Re: Vertical wind shear calculation within 400 km radius of
>>>> TC (Rick Brownrigg)
>>>> 7. Skew-T diagram (Soma Roy)
>>>>
>>>>
>>>> ----------------------------------------------------------------------
>>>>
>>>> Message: 1
>>>> Date: Wed, 18 Jul 2018 00:44:24 -0600
>>>> From: Dennis Shea <shea at ucar.edu>
>>>> To: Soma Roy <somaroy892 at gmail.com>
>>>> Cc: "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
>>>> Subject: Re: [ncl-talk] Save in netcdf
>>>> Message-ID:
>>>> <CAOF1d_4DBES=
>>>> Ddz_VgcO+9yd3ykWYzUJ6MftxpKuyeZ4eQTr1Q at mail.gmail.com>
>>>> Content-Type: text/plain; charset="utf-8"
>>>>
>>>> ;---Read needed variables from file
>>>>
>>>> itStrt = 9 ;
>>>> itLast = 12
>>>> lonStrt = 20
>>>> lonLast = 100
>>>>
>>>> u = short2flt(in->u(itStrt:itLast,:,::-1,lonStrt:lonLast))
>>>> SPH = short2flt(in1->q(itStrt:itLast,:,::-1,lonStrt:lonLast))
>>>> printVarSummary(u)
>>>>
>>>> ;===> NO: lat=in->latitude(::-1) Do not explicitly read
>>>> from
>>>> file
>>>> :===> NO: lon=in->longitude(20:100) NCL does the book keeping
>>>> ;===> NO: time=in->time(09:12) as part of data
>>>> structure/object
>>>> input
>>>>
>>>> u&level = u&level*100 ; over write with new values
>>>> u&level at units = "Pa"
>>>> printVarSummary(u) ; LOOK at this
>>>>
>>>> ;--- Zonal moisture transport
>>>> ---------------------------------------------------------------------
>>>> ZMT=u*SPH ; not
>>>> appropriate to use monthly means for
>>>>
>>>> ; for a non-linear quantity
>>>> ZMT at long_name = "zonal moisture transport"
>>>> ZMT at units="(m/s)(kg/kg)"
>>>> copy_VarCoords(u, ZMT)
>>>> printVarSummary(ZMT)
>>>>
>>>> ---
>>>> Use Method 1 to write netCDF ... easier
>>>>
>>>> filo = "example.nc" ; Output file
>>>>
>>>> * system* <
>>>> http://www.ncl.ucar.edu/Document/Functions/Built-in/system.shtml
>>>> >("/bin/rm
>>>> -f example.nc") ; remove any pre-existing file
>>>>
>>>> fout = addfile (filo, "c") ; open output file
>>>>
>>>> ;===================================================================
>>>> ; create global attributes of the file
>>>>
>>>> fAtt = True ; assign file attributes
>>>> fAtt at title = "NCL Simple Approach to netCDF Creation"
>>>> fAtt at source_file = "uvwtz_Dec2016.nc"
>>>> fAtt at Conventions = "None"
>>>> fAtt at creation_date = systemfunc ("date")
>>>> fileattdef( fout, fAtt ) ; copy file attributes
>>>>
>>>> ;===================================================================
>>>> ; make time an UNLIMITED dimension; recommended for most
>>>> applications
>>>> ;===================================================================
>>>> *filedimdef*
>>>> <http://www.ncl.ucar.edu/Document/Functions/Built-in/filedimdef.shtml
>>>> >(ncdf,"time",-1,True)
>>>>
>>>> fout->ZMT = ZMT ; 4D
>>>>
>>>>
>>>> On Tue, Jul 17, 2018 at 11:05 AM, Soma Roy <somaroy892 at gmail.com>
>>>> wrote:
>>>>
>>>> > It again shows the following error message;
>>>> >
>>>> > fatal:Variable (lev) is undefined
>>>> > fatal:["Execute.c":8640]:Execute: Error occurred at or near line 44 in
>>>> > file TEST1.ncl
>>>> >
>>>> > 'latitude', 'longitude', 'level' and 'time' are named as variable's
>>>> > dimensions in input file which is very clear from printVarSummary of
>>>> the
>>>> > data. I modified the code as per your hints. Modified code is attached
>>>> > below. Kindly inform me where I am doing wrong again.
>>>> >
>>>> >
>>>> > Thanks,
>>>> > Soma
>>>> >
>>>> > On Tue, Jul 17, 2018 at 10:09 PM, Karin Meier-Fleischer <
>>>> > meier-fleischer at dkrz.de> wrote:
>>>> >
>>>> >> Hi Soma,
>>>> >>
>>>> >> please, reply always to ncl-talk that everybody is able to follow.
>>>> >>
>>>> >> The problem is, that you changed the name of the level dimension to
>>>> lev.
>>>> >> ZMT has the dimension name level. Change lev to level and it should
>>>> work or
>>>> >>
>>>> >> ZMT!1 = "lev"
>>>> >> ZMT&lev = lev
>>>> >>
>>>> >> -Karin
>>>> >>
>>>> >> Am 17.07.18 um 14:06 schrieb Soma Roy:
>>>> >>
>>>> >> In the directory there is no file named as example.nc.
>>>> >> It is created at the run time.
>>>> >>
>>>> >> Now I am getting the following error massage;
>>>> >>
>>>> >> FileAddVar, in file: NclFile.c, line: 412
>>>> >> fatal:FileAddVar: Dimension (level) is not currently defined, can't
>>>> add
>>>> >> variable
>>>> >> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 77
>>>> in
>>>> >> file TEST1.ncl
>>>> >>
>>>> >> Kindly inform me what I am doing wrong in this.
>>>> >>
>>>> >> Thanking you,
>>>> >> Soma
>>>> >>
>>>> >>
>>>> >> On Tue, Jul 17, 2018, 16:31 Karin Meier-Fleischer <
>>>> >> meier-fleischer at dkrz.de> wrote:
>>>> >>
>>>> >>> Hi Soma,
>>>> >>>
>>>> >>> if the file example.nc already exists NCL can't create a new one
>>>> with
>>>> >>> the same name. You have to remove the file, first.
>>>> >>>
>>>> >>> filo = "example.nc" ; Output file
>>>> >>> * system("rm -rf "+filo)*
>>>> >>> fout = addfile (filo, "c") ; open output file
>>>> >>>
>>>> >>> -Karin
>>>> >>>
>>>> >>> Am 17.07.18 um 12:22 schrieb Soma Roy:
>>>> >>>
>>>> >>> I am trying to save the variable ZMT in NetCDF file following the
>>>> >>> instructions at the below-mentioned link;
>>>> >>>
>>>> >>> http://www.ncl.ucar.edu/Applications/method_2.shtml
>>>> >>>
>>>> >>> Summary of the dataset is as follows;
>>>> >>>
>>>> >>> Variable: lat
>>>> >>> Type: float
>>>> >>> Total Size: 244 bytes
>>>> >>> 61 values
>>>> >>> Number of Dimensions: 1
>>>> >>> Dimensions and sizes: [latitude | 61]
>>>> >>> Coordinates:
>>>> >>> latitude: [ 5..20]
>>>> >>> Number Of Attributes: 2
>>>> >>> units : degrees_north
>>>> >>> long_name : latitude
>>>> >>>
>>>> >>> Variable: lon
>>>> >>> Type: float
>>>> >>> Total Size: 324 bytes
>>>> >>> 81 values
>>>> >>> Number of Dimensions: 1
>>>> >>> Dimensions and sizes: [longitude | 81]
>>>> >>> Coordinates:
>>>> >>> longitude: [70..90]
>>>> >>> Number Of Attributes: 2
>>>> >>> units : degrees_east
>>>> >>> long_name : longitude
>>>> >>>
>>>> >>> Variable: time
>>>> >>> Type: integer
>>>> >>> Total Size: 16 bytes
>>>> >>> 4 values
>>>> >>> Number of Dimensions: 1
>>>> >>> Dimensions and sizes: [time | 4]
>>>> >>> Coordinates:
>>>> >>> time: [1025088..1025160]
>>>> >>> Number Of Attributes: 3
>>>> >>> units : hours since 1900-01-01 00:00:0.0
>>>> >>> long_name : time
>>>> >>> calendar : gregorian
>>>> >>>
>>>> >>>
>>>> >>> Variable: pnew
>>>> >>> Type: integer
>>>> >>> Total Size: 112 bytes
>>>> >>> 28 values
>>>> >>> Number of Dimensions: 1
>>>> >>> Dimensions and sizes: [level | 28]
>>>> >>> Coordinates:
>>>> >>> level: [70..1000]
>>>> >>> Number Of Attributes: 2
>>>> >>> units : Pa
>>>> >>> long_name : pressure_level
>>>> >>> (0) 7000
>>>> >>> (1) 10000
>>>> >>> (2) 12500
>>>> >>> (3) 15000
>>>> >>> ..........
>>>> >>> (26) 97500
>>>> >>> (27) 100000
>>>> >>>
>>>> >>> Variable: ZMT
>>>> >>> Type: float
>>>> >>> Total Size: 2213568 bytes
>>>> >>> 553392 values
>>>> >>> Number of Dimensions: 4
>>>> >>> Dimensions and sizes: [time | 4] x [level | 28] x [latitude | 61]
>>>> x
>>>> >>> [longitude | 81]
>>>> >>> Coordinates:
>>>> >>> time: [1025088..1025160]
>>>> >>> level: [70..1000]
>>>> >>> latitude: [ 5..20]
>>>> >>> longitude: [70..90]
>>>> >>> Number Of Attributes: 3
>>>> >>> units : (m/s)(kg/kg)
>>>> >>> long_name : zonal moisture transport
>>>> >>> _FillValue : -32767
>>>> >>>
>>>> >>> At run-time, I am getting the following error message;
>>>> >>>
>>>> >>> fatal:Could not create (example.nc)
>>>> >>> warning:FileSetFileOption: invalid file or format
>>>> >>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 65
>>>> in
>>>> >>> file TEST1.ncl
>>>> >>>
>>>> >>> My script is attached below.
>>>> >>>
>>>> >>> Please kindly inform me what is wrong with this code and how to fix
>>>> it.
>>>> >>>
>>>> >>> Thanking you,
>>>> >>> Soma
>>>> >>>
>>>> >>>
>>>> >>> _______________________________________________
>>>> >>> ncl-talk mailing listncl-talk at ucar.edu
>>>> >>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>> >>>
>>>> >>>
>>>> >>> _______________________________________________
>>>> >>> ncl-talk mailing list
>>>> >>> ncl-talk at ucar.edu
>>>> >>> List instructions, subscriber options, unsubscribe:
>>>> >>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>> >>>
>>>> >>
>>>> >
>>>> > _______________________________________________
>>>> > ncl-talk mailing list
>>>> > ncl-talk at ucar.edu
>>>> > List instructions, subscriber options, unsubscribe:
>>>> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>> >
>>>> >
>>>> -------------- next part --------------
>>>> An HTML attachment was scrubbed...
>>>> URL: <
>>>> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180718/32cfd5f9/attachment-0001.html
>>>> >
>>>>
>>>> ------------------------------
>>>>
>>>> Message: 2
>>>> Date: Wed, 18 Jul 2018 13:23:56 +0530
>>>> From: Soma Roy <somaroy892 at gmail.com>
>>>> To: Dennis Shea <shea at ucar.edu>
>>>> Cc: ncl-talk <ncl-talk at ucar.edu>
>>>> Subject: Re: [ncl-talk] Save in netcdf
>>>> Message-ID:
>>>> <
>>>> CAAJRSQd4RiU7Ro9m+ar56R7eKcdtKWTVWqRNdp_eJdAcyBt2Yg at mail.gmail.com>
>>>> Content-Type: text/plain; charset="utf-8"
>>>>
>>>> I am getting the following error message after running the attached
>>>> code.
>>>>
>>>> fatal:Variable (ncdf) is undefined
>>>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 50 in
>>>> file
>>>> Test2.ncl
>>>>
>>>> Modified code as per easier method is attached below.
>>>>
>>>> Let me know what is wrong in this code.
>>>>
>>>> Thanks,
>>>> Soma
>>>>
>>>> On Wed, Jul 18, 2018, 12:14 Dennis Shea <shea at ucar.edu> wrote:
>>>>
>>>> > ;---Read needed variables from file
>>>> >
>>>> > itStrt = 9 ;
>>>> > itLast = 12
>>>> > lonStrt = 20
>>>> > lonLast = 100
>>>> >
>>>> > u = short2flt(in->u(itStrt:itLast,:,::-1,lonStrt:lonLast))
>>>> > SPH = short2flt(in1->q(itStrt:itLast,:,::-1,lonStrt:lonLast))
>>>> > printVarSummary(u)
>>>> >
>>>> > ;===> NO: lat=in->latitude(::-1) Do not explicitly read
>>>> from
>>>> > file
>>>> > :===> NO: lon=in->longitude(20:100) NCL does the book keeping
>>>> > ;===> NO: time=in->time(09:12) as part of data
>>>> > structure/object input
>>>> >
>>>> > u&level = u&level*100 ; over write with new values
>>>> > u&level at units = "Pa"
>>>> > printVarSummary(u) ; LOOK at this
>>>> >
>>>> > ;--- Zonal moisture transport
>>>> > ---------------------------------------------------------------------
>>>> > ZMT=u*SPH ; not
>>>> > appropriate to use monthly means for
>>>> >
>>>> > ; for a non-linear quantity
>>>> > ZMT at long_name = "zonal moisture transport"
>>>> > ZMT at units="(m/s)(kg/kg)"
>>>> > copy_VarCoords(u, ZMT)
>>>> > printVarSummary(ZMT)
>>>> >
>>>> > ---
>>>> > Use Method 1 to write netCDF ... easier
>>>> >
>>>> > filo = "example.nc" ; Output file
>>>> >
>>>> > * system* <
>>>> http://www.ncl.ucar.edu/Document/Functions/Built-in/system.shtml>("/bin/rm
>>>> -f example.nc") ; remove any pre-existing file
>>>> >
>>>> > fout = addfile (filo, "c") ; open output file
>>>> >
>>>> > ;===================================================================
>>>> > ; create global attributes of the file
>>>> >
>>>> > fAtt = True ; assign file attributes
>>>> > fAtt at title = "NCL Simple Approach to netCDF Creation"
>>>> > fAtt at source_file = "uvwtz_Dec2016.nc"
>>>> > fAtt at Conventions = "None"
>>>> > fAtt at creation_date = systemfunc ("date")
>>>> > fileattdef( fout, fAtt ) ; copy file attributes
>>>> >
>>>> >
>>>> ;===================================================================
>>>> > ; make time an UNLIMITED dimension; recommended for most
>>>> applications
>>>> >
>>>> ;===================================================================
>>>> > *filedimdef* <
>>>> http://www.ncl.ucar.edu/Document/Functions/Built-in/filedimdef.shtml
>>>> >(ncdf,"time",-1,True)
>>>> >
>>>> > fout->ZMT = ZMT ; 4D
>>>> >
>>>> >
>>>> > On Tue, Jul 17, 2018 at 11:05 AM, Soma Roy <somaroy892 at gmail.com>
>>>> wrote:
>>>> >
>>>> >> It again shows the following error message;
>>>> >>
>>>> >> fatal:Variable (lev) is undefined
>>>> >> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 44
>>>> in
>>>> >> file TEST1.ncl
>>>> >>
>>>> >> 'latitude', 'longitude', 'level' and 'time' are named as variable's
>>>> >> dimensions in input file which is very clear from printVarSummary of
>>>> the
>>>> >> data. I modified the code as per your hints. Modified code is
>>>> attached
>>>> >> below. Kindly inform me where I am doing wrong again.
>>>> >>
>>>> >>
>>>> >> Thanks,
>>>> >> Soma
>>>> >>
>>>> >> On Tue, Jul 17, 2018 at 10:09 PM, Karin Meier-Fleischer <
>>>> >> meier-fleischer at dkrz.de> wrote:
>>>> >>
>>>> >>> Hi Soma,
>>>> >>>
>>>> >>> please, reply always to ncl-talk that everybody is able to follow.
>>>> >>>
>>>> >>> The problem is, that you changed the name of the level dimension to
>>>> lev.
>>>> >>> ZMT has the dimension name level. Change lev to level and it should
>>>> work or
>>>> >>>
>>>> >>> ZMT!1 = "lev"
>>>> >>> ZMT&lev = lev
>>>> >>>
>>>> >>> -Karin
>>>> >>>
>>>> >>> Am 17.07.18 um 14:06 schrieb Soma Roy:
>>>> >>>
>>>> >>> In the directory there is no file named as example.nc.
>>>> >>> It is created at the run time.
>>>> >>>
>>>> >>> Now I am getting the following error massage;
>>>> >>>
>>>> >>> FileAddVar, in file: NclFile.c, line: 412
>>>> >>> fatal:FileAddVar: Dimension (level) is not currently defined, can't
>>>> add
>>>> >>> variable
>>>> >>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 77
>>>> in
>>>> >>> file TEST1.ncl
>>>> >>>
>>>> >>> Kindly inform me what I am doing wrong in this.
>>>> >>>
>>>> >>> Thanking you,
>>>> >>> Soma
>>>> >>>
>>>> >>>
>>>> >>> On Tue, Jul 17, 2018, 16:31 Karin Meier-Fleischer <
>>>> >>> meier-fleischer at dkrz.de> wrote:
>>>> >>>
>>>> >>>> Hi Soma,
>>>> >>>>
>>>> >>>> if the file example.nc already exists NCL can't create a new one
>>>> with
>>>> >>>> the same name. You have to remove the file, first.
>>>> >>>>
>>>> >>>> filo = "example.nc" ; Output file
>>>> >>>> * system("rm -rf "+filo)*
>>>> >>>> fout = addfile (filo, "c") ; open output file
>>>> >>>>
>>>> >>>> -Karin
>>>> >>>>
>>>> >>>> Am 17.07.18 um 12:22 schrieb Soma Roy:
>>>> >>>>
>>>> >>>> I am trying to save the variable ZMT in NetCDF file following the
>>>> >>>> instructions at the below-mentioned link;
>>>> >>>>
>>>> >>>> http://www.ncl.ucar.edu/Applications/method_2.shtml
>>>> >>>>
>>>> >>>> Summary of the dataset is as follows;
>>>> >>>>
>>>> >>>> Variable: lat
>>>> >>>> Type: float
>>>> >>>> Total Size: 244 bytes
>>>> >>>> 61 values
>>>> >>>> Number of Dimensions: 1
>>>> >>>> Dimensions and sizes: [latitude | 61]
>>>> >>>> Coordinates:
>>>> >>>> latitude: [ 5..20]
>>>> >>>> Number Of Attributes: 2
>>>> >>>> units : degrees_north
>>>> >>>> long_name : latitude
>>>> >>>>
>>>> >>>> Variable: lon
>>>> >>>> Type: float
>>>> >>>> Total Size: 324 bytes
>>>> >>>> 81 values
>>>> >>>> Number of Dimensions: 1
>>>> >>>> Dimensions and sizes: [longitude | 81]
>>>> >>>> Coordinates:
>>>> >>>> longitude: [70..90]
>>>> >>>> Number Of Attributes: 2
>>>> >>>> units : degrees_east
>>>> >>>> long_name : longitude
>>>> >>>>
>>>> >>>> Variable: time
>>>> >>>> Type: integer
>>>> >>>> Total Size: 16 bytes
>>>> >>>> 4 values
>>>> >>>> Number of Dimensions: 1
>>>> >>>> Dimensions and sizes: [time | 4]
>>>> >>>> Coordinates:
>>>> >>>> time: [1025088..1025160]
>>>> >>>> Number Of Attributes: 3
>>>> >>>> units : hours since 1900-01-01 00:00:0.0
>>>> >>>> long_name : time
>>>> >>>> calendar : gregorian
>>>> >>>>
>>>> >>>>
>>>> >>>> Variable: pnew
>>>> >>>> Type: integer
>>>> >>>> Total Size: 112 bytes
>>>> >>>> 28 values
>>>> >>>> Number of Dimensions: 1
>>>> >>>> Dimensions and sizes: [level | 28]
>>>> >>>> Coordinates:
>>>> >>>> level: [70..1000]
>>>> >>>> Number Of Attributes: 2
>>>> >>>> units : Pa
>>>> >>>> long_name : pressure_level
>>>> >>>> (0) 7000
>>>> >>>> (1) 10000
>>>> >>>> (2) 12500
>>>> >>>> (3) 15000
>>>> >>>> ..........
>>>> >>>> (26) 97500
>>>> >>>> (27) 100000
>>>> >>>>
>>>> >>>> Variable: ZMT
>>>> >>>> Type: float
>>>> >>>> Total Size: 2213568 bytes
>>>> >>>> 553392 values
>>>> >>>> Number of Dimensions: 4
>>>> >>>> Dimensions and sizes: [time | 4] x [level | 28] x [latitude |
>>>> 61] x
>>>> >>>> [longitude | 81]
>>>> >>>> Coordinates:
>>>> >>>> time: [1025088..1025160]
>>>> >>>> level: [70..1000]
>>>> >>>> latitude: [ 5..20]
>>>> >>>> longitude: [70..90]
>>>> >>>> Number Of Attributes: 3
>>>> >>>> units : (m/s)(kg/kg)
>>>> >>>> long_name : zonal moisture transport
>>>> >>>> _FillValue : -32767
>>>> >>>>
>>>> >>>> At run-time, I am getting the following error message;
>>>> >>>>
>>>> >>>> fatal:Could not create (example.nc)
>>>> >>>> warning:FileSetFileOption: invalid file or format
>>>> >>>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line
>>>> 65 in
>>>> >>>> file TEST1.ncl
>>>> >>>>
>>>> >>>> My script is attached below.
>>>> >>>>
>>>> >>>> Please kindly inform me what is wrong with this code and how to
>>>> fix it.
>>>> >>>>
>>>> >>>> Thanking you,
>>>> >>>> Soma
>>>> >>>>
>>>> >>>>
>>>> >>>> _______________________________________________
>>>> >>>> ncl-talk mailing listncl-talk at ucar.edu
>>>> >>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>> >>>>
>>>> >>>>
>>>> >>>> _______________________________________________
>>>> >>>> ncl-talk mailing list
>>>> >>>> ncl-talk at ucar.edu
>>>> >>>> List instructions, subscriber options, unsubscribe:
>>>> >>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>> >>>>
>>>> >>>
>>>> >>
>>>> >> _______________________________________________
>>>> >> ncl-talk mailing list
>>>> >> ncl-talk at ucar.edu
>>>> >> List instructions, subscriber options, unsubscribe:
>>>> >> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>> >>
>>>> >>
>>>> >
>>>> -------------- next part --------------
>>>> An HTML attachment was scrubbed...
>>>> URL: <
>>>> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180718/4bc152ee/attachment-0001.html
>>>> >
>>>> -------------- next part --------------
>>>> A non-text attachment was scrubbed...
>>>> Name: Test2.ncl
>>>> Type: application/octet-stream
>>>> Size: 2282 bytes
>>>> Desc: not available
>>>> URL: <
>>>> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180718/4bc152ee/attachment-0001.obj
>>>> >
>>>>
>>>> ------------------------------
>>>>
>>>> Message: 3
>>>> Date: Wed, 18 Jul 2018 09:16:14 +0000
>>>> From: Be?ta Szab?-Tak?cs <szabo.b at czechglobe.cz>
>>>> To: Soma Roy <somaroy892 at gmail.com>, Dennis Shea <shea at ucar.edu>
>>>> Cc: ncl-talk <ncl-talk at ucar.edu>
>>>> Subject: Re: [ncl-talk] Save in netcdf
>>>> Message-ID:
>>>> <
>>>> HE1PR08MB2665F224A8E77049AE5A37A3F0530 at HE1PR08MB2665.eurprd08.prod.outlook.com
>>>> >
>>>>
>>>> Content-Type: text/plain; charset="utf-8"
>>>>
>>>> Dear Soma,
>>>> I think the matter is in the line:
>>>>
>>>> filedimdef(ncdf,"time",-1,True)
>>>>
>>>> You should change to
>>>>
>>>> filedimdef(fout,"time",-1,True)
>>>>
>>>> I hope I could help you!
>>>> Best regards,
>>>> Beata
>>>>
>>>> From: ncl-talk [mailto:ncl-talk-bounces at ucar.edu] On Behalf Of Soma Roy
>>>> Sent: Wednesday, July 18, 2018 9:54 AM
>>>> To: Dennis Shea <shea at ucar.edu>
>>>> Cc: ncl-talk <ncl-talk at ucar.edu>
>>>> Subject: Re: [ncl-talk] Save in netcdf
>>>>
>>>> I am getting the following error message after running the attached
>>>> code.
>>>>
>>>> fatal:Variable (ncdf) is undefined
>>>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 50 in
>>>> file Test2.ncl
>>>>
>>>> Modified code as per easier method is attached below.
>>>>
>>>> Let me know what is wrong in this code.
>>>>
>>>> Thanks,
>>>> Soma
>>>>
>>>> On Wed, Jul 18, 2018, 12:14 Dennis Shea <shea at ucar.edu<mailto:
>>>> shea at ucar.edu>> wrote:
>>>> ;---Read needed variables from file
>>>> itStrt = 9 ;
>>>> itLast = 12
>>>> lonStrt = 20
>>>> lonLast = 100
>>>>
>>>> u = short2flt(in->u(itStrt:itLast,:,::-1,lonStrt:lonLast))
>>>> SPH = short2flt(in1->q(itStrt:itLast,:,::-1,lonStrt:lonLast))
>>>> printVarSummary(u)
>>>>
>>>> ;===> NO: lat=in->latitude(::-1) Do not explicitly read
>>>> from file
>>>> :===> NO: lon=in->longitude(20:100) NCL does the book keeping
>>>> ;===> NO: time=in->time(09:12) as part of data
>>>> structure/object input
>>>>
>>>> u&level = u&level*100 ; over write with new values
>>>> u&level at units = "Pa"
>>>> printVarSummary(u) ; LOOK at this
>>>>
>>>> ;--- Zonal moisture transport
>>>> ---------------------------------------------------------------------
>>>> ZMT=u*SPH ; not
>>>> appropriate to use monthly means for
>>>>
>>>> ; for a non-linear quantity
>>>> ZMT at long_name = "zonal moisture transport"
>>>> ZMT at units="(m/s)(kg/kg)"
>>>> copy_VarCoords(u, ZMT)
>>>> printVarSummary(ZMT)
>>>>
>>>> ---
>>>> Use Method 1 to write netCDF ... easier
>>>>
>>>> filo = "example.nc<http://example.nc>" ; Output
>>>> file
>>>>
>>>> system<
>>>> http://www.ncl.ucar.edu/Document/Functions/Built-in/system.shtml>("/bin/rm
>>>> -f example.nc<http://example.nc>") ; remove any pre-existing file
>>>> fout = addfile (filo, "c") ; open output file
>>>>
>>>> ;===================================================================
>>>> ; create global attributes of the file
>>>>
>>>> fAtt = True ; assign file attributes
>>>> fAtt at title = "NCL Simple Approach to netCDF Creation"
>>>> fAtt at source_file = "uvwtz_Dec2016.nc"
>>>> fAtt at Conventions = "None"
>>>> fAtt at creation_date = systemfunc ("date")
>>>> fileattdef( fout, fAtt ) ; copy file attributes
>>>>
>>>> ;===================================================================
>>>>
>>>> ; make time an UNLIMITED dimension; recommended for most
>>>> applications
>>>>
>>>> ;===================================================================
>>>>
>>>> filedimdef<
>>>> http://www.ncl.ucar.edu/Document/Functions/Built-in/filedimdef.shtml
>>>> >(ncdf,"time",-1,True)
>>>>
>>>> fout->ZMT = ZMT ; 4D
>>>>
>>>> On Tue, Jul 17, 2018 at 11:05 AM, Soma Roy <somaroy892 at gmail.com
>>>> <mailto:somaroy892 at gmail.com>> wrote:
>>>> It again shows the following error message;
>>>>
>>>> fatal:Variable (lev) is undefined
>>>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 44 in
>>>> file TEST1.ncl
>>>>
>>>> 'latitude', 'longitude', 'level' and 'time' are named as variable's
>>>> dimensions in input file which is very clear from printVarSummary of the
>>>> data. I modified the code as per your hints. Modified code is attached
>>>> below. Kindly inform me where I am doing wrong again.
>>>>
>>>> Thanks,
>>>> Soma
>>>>
>>>> On Tue, Jul 17, 2018 at 10:09 PM, Karin Meier-Fleischer <
>>>> meier-fleischer at dkrz.de<mailto:meier-fleischer at dkrz.de>> wrote:
>>>> Hi Soma,
>>>>
>>>> please, reply always to ncl-talk that everybody is able to follow.
>>>>
>>>> The problem is, that you changed the name of the level dimension to
>>>> lev. ZMT has the dimension name level. Change lev to level and it should
>>>> work or
>>>>
>>>> ZMT!1 = "lev"
>>>> ZMT&lev = lev
>>>>
>>>> -Karin
>>>> Am 17.07.18 um 14:06 schrieb Soma Roy:
>>>> In the directory there is no file named as example.nc<http://example.nc
>>>> >.
>>>> It is created at the run time.
>>>>
>>>> Now I am getting the following error massage;
>>>>
>>>> FileAddVar, in file: NclFile.c, line: 412
>>>> fatal:FileAddVar: Dimension (level) is not currently defined, can't add
>>>> variable
>>>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 77 in
>>>> file TEST1.ncl
>>>>
>>>> Kindly inform me what I am doing wrong in this.
>>>>
>>>> Thanking you,
>>>> Soma
>>>>
>>>>
>>>> On Tue, Jul 17, 2018, 16:31 Karin Meier-Fleischer <
>>>> meier-fleischer at dkrz.de<mailto:meier-fleischer at dkrz.de>> wrote:
>>>> Hi Soma,
>>>>
>>>> if the file example.nc<http://example.nc> already exists NCL can't
>>>> create a new one with the same name. You have to remove the file, first.
>>>>
>>>> filo = "example.nc<http://example.nc>" ; Output
>>>> file
>>>> system("rm -rf "+filo)
>>>> fout = addfile (filo, "c") ; open output file
>>>>
>>>> -Karin
>>>> Am 17.07.18 um 12:22 schrieb Soma Roy:
>>>> I am trying to save the variable ZMT in NetCDF file following the
>>>> instructions at the below-mentioned link;
>>>>
>>>> http://www.ncl.ucar.edu/Applications/method_2.shtml
>>>>
>>>> Summary of the dataset is as follows;
>>>>
>>>> Variable: lat
>>>> Type: float
>>>> Total Size: 244 bytes
>>>> 61 values
>>>> Number of Dimensions: 1
>>>> Dimensions and sizes: [latitude | 61]
>>>> Coordinates:
>>>> latitude: [ 5..20]
>>>> Number Of Attributes: 2
>>>> units : degrees_north
>>>> long_name : latitude
>>>>
>>>> Variable: lon
>>>> Type: float
>>>> Total Size: 324 bytes
>>>> 81 values
>>>> Number of Dimensions: 1
>>>> Dimensions and sizes: [longitude | 81]
>>>> Coordinates:
>>>> longitude: [70..90]
>>>> Number Of Attributes: 2
>>>> units : degrees_east
>>>> long_name : longitude
>>>>
>>>> Variable: time
>>>> Type: integer
>>>> Total Size: 16 bytes
>>>> 4 values
>>>> Number of Dimensions: 1
>>>> Dimensions and sizes: [time | 4]
>>>> Coordinates:
>>>> time: [1025088..1025160]
>>>> Number Of Attributes: 3
>>>> units : hours since 1900-01-01 00:00:0.0
>>>> long_name : time
>>>> calendar : gregorian
>>>>
>>>>
>>>> Variable: pnew
>>>> Type: integer
>>>> Total Size: 112 bytes
>>>> 28 values
>>>> Number of Dimensions: 1
>>>> Dimensions and sizes: [level | 28]
>>>> Coordinates:
>>>> level: [70..1000]
>>>> Number Of Attributes: 2
>>>> units : Pa
>>>> long_name : pressure_level
>>>> (0) 7000
>>>> (1) 10000
>>>> (2) 12500
>>>> (3) 15000
>>>> ..........
>>>> (26) 97500
>>>> (27) 100000
>>>>
>>>> Variable: ZMT
>>>> Type: float
>>>> Total Size: 2213568 bytes
>>>> 553392 values
>>>> Number of Dimensions: 4
>>>> Dimensions and sizes: [time | 4] x [level | 28] x [latitude | 61] x
>>>> [longitude | 81]
>>>> Coordinates:
>>>> time: [1025088..1025160]
>>>> level: [70..1000]
>>>> latitude: [ 5..20]
>>>> longitude: [70..90]
>>>> Number Of Attributes: 3
>>>> units : (m/s)(kg/kg)
>>>> long_name : zonal moisture transport
>>>> _FillValue : -32767
>>>>
>>>> At run-time, I am getting the following error message;
>>>>
>>>> fatal:Could not create (example.nc<http://example.nc/>)
>>>> warning:FileSetFileOption: invalid file or format
>>>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 65 in
>>>> file TEST1.ncl
>>>>
>>>> My script is attached below.
>>>>
>>>> Please kindly inform me what is wrong with this code and how to fix it.
>>>>
>>>> Thanking you,
>>>> Soma
>>>>
>>>>
>>>> _______________________________________________
>>>>
>>>> ncl-talk mailing list
>>>>
>>>> ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
>>>>
>>>> List instructions, subscriber options, unsubscribe:
>>>>
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>> -------------- next part --------------
>>>> An HTML attachment was scrubbed...
>>>> URL: <
>>>> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180718/93da231e/attachment-0001.html
>>>> >
>>>>
>>>> ------------------------------
>>>>
>>>> Message: 4
>>>> Date: Wed, 18 Jul 2018 15:09:59 +0530
>>>> From: Soma Roy <somaroy892 at gmail.com>
>>>> To: Be?ta Szab?-Tak?cs <szabo.b at czechglobe.cz>
>>>> Cc: ncl-talk <ncl-talk at ucar.edu>
>>>> Subject: Re: [ncl-talk] Save in netcdf
>>>> Message-ID:
>>>> <
>>>> CAAJRSQcMwxKMJxe-jQOXXSuf8D38hSDMOnSefvywc2cEWjHu0Q at mail.gmail.com>
>>>> Content-Type: text/plain; charset="utf-8"
>>>>
>>>> It worked. Thanks a lot Dennis and Beata.
>>>>
>>>> Best,
>>>> Soma
>>>>
>>>> On Wed, Jul 18, 2018, 14:46 Be?ta Szab?-Tak?cs <szabo.b at czechglobe.cz>
>>>> wrote:
>>>>
>>>> > Dear Soma,
>>>> >
>>>> > I think the matter is in the line:
>>>> >
>>>> >
>>>> >
>>>> > filedimdef(ncdf,"time",-1,True)
>>>> >
>>>> >
>>>> >
>>>> > You should change to
>>>> >
>>>> >
>>>> >
>>>> > filedimdef(fout,"time",-1,True)
>>>> >
>>>> >
>>>> >
>>>> > I hope I could help you!
>>>> >
>>>> > Best regards,
>>>> >
>>>> > Beata
>>>> >
>>>> >
>>>> >
>>>> > *From:* ncl-talk [mailto:ncl-talk-bounces at ucar.edu] *On Behalf Of
>>>> *Soma
>>>> > Roy
>>>> > *Sent:* Wednesday, July 18, 2018 9:54 AM
>>>> > *To:* Dennis Shea <shea at ucar.edu>
>>>> > *Cc:* ncl-talk <ncl-talk at ucar.edu>
>>>> > *Subject:* Re: [ncl-talk] Save in netcdf
>>>> >
>>>> >
>>>> >
>>>> > I am getting the following error message after running the attached
>>>> code.
>>>> >
>>>> >
>>>> >
>>>> > fatal:Variable (ncdf) is undefined
>>>> >
>>>> > fatal:["Execute.c":8640]:Execute: Error occurred at or near line 50 in
>>>> > file Test2.ncl
>>>> >
>>>> >
>>>> >
>>>> > Modified code as per easier method is attached below.
>>>> >
>>>> >
>>>> >
>>>> > Let me know what is wrong in this code.
>>>> >
>>>> >
>>>> >
>>>> > Thanks,
>>>> >
>>>> > Soma
>>>> >
>>>> >
>>>> >
>>>> > On Wed, Jul 18, 2018, 12:14 Dennis Shea <shea at ucar.edu> wrote:
>>>> >
>>>> > ;---Read needed variables from file
>>>> >
>>>> > itStrt = 9 ;
>>>> >
>>>> > itLast = 12
>>>> >
>>>> > lonStrt = 20
>>>> >
>>>> > lonLast = 100
>>>> >
>>>> >
>>>> > u = short2flt(in->u(itStrt:itLast,:,::-1,lonStrt:lonLast))
>>>> >
>>>> > SPH = short2flt(in1->q(itStrt:itLast,:,::-1,lonStrt:lonLast))
>>>> > printVarSummary(u)
>>>> >
>>>> > ;===> NO: lat=in->latitude(::-1) Do not explicitly read
>>>> from
>>>> > file
>>>> > :===> NO: lon=in->longitude(20:100) NCL does the book keeping
>>>> > ;===> NO: time=in->time(09:12) as part of data
>>>> > structure/object input
>>>> >
>>>> >
>>>> > u&level = u&level*100 ; over write with new values
>>>> >
>>>> > u&level at units = "Pa"
>>>> >
>>>> > printVarSummary(u) ; LOOK at this
>>>> >
>>>> >
>>>> > ;--- Zonal moisture transport
>>>> > ---------------------------------------------------------------------
>>>> > ZMT=u*SPH ; not
>>>> > appropriate to use monthly means for
>>>> >
>>>> >
>>>> > ; for a non-linear quantity
>>>> >
>>>> > ZMT at long_name = "zonal moisture transport"
>>>> > ZMT at units="(m/s)(kg/kg)"
>>>> > copy_VarCoords(u, ZMT)
>>>> > printVarSummary(ZMT)
>>>> >
>>>> > ---
>>>> >
>>>> > Use Method 1 to write netCDF ... easier
>>>> >
>>>> > filo = "example.nc" ; Output file
>>>> >
>>>> > * system* <
>>>> http://www.ncl.ucar.edu/Document/Functions/Built-in/system.shtml>("/bin/rm
>>>> -f example.nc") ; remove any pre-existing file
>>>> >
>>>> > fout = addfile (filo, "c") ; open output file
>>>> >
>>>> > ;===================================================================
>>>> > ; create global attributes of the file
>>>> >
>>>> > fAtt = True ; assign file attributes
>>>> > fAtt at title = "NCL Simple Approach to netCDF Creation"
>>>> > fAtt at source_file = "uvwtz_Dec2016.nc"
>>>> > fAtt at Conventions = "None"
>>>> > fAtt at creation_date = systemfunc ("date")
>>>> > fileattdef( fout, fAtt ) ; copy file attributes
>>>> >
>>>> >
>>>> ;===================================================================
>>>> >
>>>> > ; make time an UNLIMITED dimension; recommended for most
>>>> applications
>>>> >
>>>> >
>>>> ;===================================================================
>>>> >
>>>> > *filedimdef* <
>>>> http://www.ncl.ucar.edu/Document/Functions/Built-in/filedimdef.shtml
>>>> >(ncdf,"time",-1,True)
>>>> >
>>>> > fout->ZMT = ZMT ; 4D
>>>> >
>>>> >
>>>> >
>>>> > On Tue, Jul 17, 2018 at 11:05 AM, Soma Roy <somaroy892 at gmail.com>
>>>> wrote:
>>>> >
>>>> > It again shows the following error message;
>>>> >
>>>> >
>>>> >
>>>> > fatal:Variable (lev) is undefined
>>>> >
>>>> > fatal:["Execute.c":8640]:Execute: Error occurred at or near line 44 in
>>>> > file TEST1.ncl
>>>> >
>>>> >
>>>> >
>>>> > 'latitude', 'longitude', 'level' and 'time' are named as variable's
>>>> > dimensions in input file which is very clear from printVarSummary of
>>>> the
>>>> > data. I modified the code as per your hints. Modified code is attached
>>>> > below. Kindly inform me where I am doing wrong again.
>>>> >
>>>> >
>>>> >
>>>> > Thanks,
>>>> >
>>>> > Soma
>>>> >
>>>> >
>>>> >
>>>> > On Tue, Jul 17, 2018 at 10:09 PM, Karin Meier-Fleischer <
>>>> > meier-fleischer at dkrz.de> wrote:
>>>> >
>>>> > Hi Soma,
>>>> >
>>>> > please, reply always to ncl-talk that everybody is able to follow.
>>>> >
>>>> > The problem is, that you changed the name of the level dimension to
>>>> lev.
>>>> > ZMT has the dimension name level. Change lev to level and it should
>>>> work or
>>>> >
>>>> > ZMT!1 = "lev"
>>>> > ZMT&lev = lev
>>>> >
>>>> > -Karin
>>>> >
>>>> > Am 17.07.18 um 14:06 schrieb Soma Roy:
>>>> >
>>>> > In the directory there is no file named as example.nc.
>>>> >
>>>> > It is created at the run time.
>>>> >
>>>> >
>>>> >
>>>> > Now I am getting the following error massage;
>>>> >
>>>> >
>>>> >
>>>> > FileAddVar, in file: NclFile.c, line: 412
>>>> >
>>>> > fatal:FileAddVar: Dimension (level) is not currently defined, can't
>>>> add
>>>> > variable
>>>> >
>>>> > fatal:["Execute.c":8640]:Execute: Error occurred at or near line 77 in
>>>> > file TEST1.ncl
>>>> >
>>>> >
>>>> >
>>>> > Kindly inform me what I am doing wrong in this.
>>>> >
>>>> >
>>>> >
>>>> > Thanking you,
>>>> >
>>>> > Soma
>>>> >
>>>> >
>>>> >
>>>> >
>>>> >
>>>> > On Tue, Jul 17, 2018, 16:31 Karin Meier-Fleischer <
>>>> meier-fleischer at dkrz.de>
>>>> > wrote:
>>>> >
>>>> > Hi Soma,
>>>> >
>>>> > if the file example.nc already exists NCL can't create a new one
>>>> with the
>>>> > same name. You have to remove the file, first.
>>>> >
>>>> > filo = "example.nc" ; Output file
>>>> > * system("rm -rf "+filo)*
>>>> > fout = addfile (filo, "c") ; open output file
>>>> >
>>>> > -Karin
>>>> >
>>>> > Am 17.07.18 um 12:22 schrieb Soma Roy:
>>>> >
>>>> > I am trying to save the variable ZMT in NetCDF file following the
>>>> > instructions at the below-mentioned link;
>>>> >
>>>> >
>>>> >
>>>> > http://www.ncl.ucar.edu/Applications/method_2.shtml
>>>> >
>>>> >
>>>> >
>>>> > Summary of the dataset is as follows;
>>>> >
>>>> >
>>>> >
>>>> > Variable: lat
>>>> >
>>>> > Type: float
>>>> >
>>>> > Total Size: 244 bytes
>>>> >
>>>> > 61 values
>>>> >
>>>> > Number of Dimensions: 1
>>>> >
>>>> > Dimensions and sizes: [latitude | 61]
>>>> >
>>>> > Coordinates:
>>>> >
>>>> > latitude: [ 5..20]
>>>> >
>>>> > Number Of Attributes: 2
>>>> >
>>>> > units : degrees_north
>>>> >
>>>> > long_name : latitude
>>>> >
>>>> >
>>>> >
>>>> > Variable: lon
>>>> >
>>>> > Type: float
>>>> >
>>>> > Total Size: 324 bytes
>>>> >
>>>> > 81 values
>>>> >
>>>> > Number of Dimensions: 1
>>>> >
>>>> > Dimensions and sizes: [longitude | 81]
>>>> >
>>>> > Coordinates:
>>>> >
>>>> > longitude: [70..90]
>>>> >
>>>> > Number Of Attributes: 2
>>>> >
>>>> > units : degrees_east
>>>> >
>>>> > long_name : longitude
>>>> >
>>>> >
>>>> >
>>>> > Variable: time
>>>> >
>>>> > Type: integer
>>>> >
>>>> > Total Size: 16 bytes
>>>> >
>>>> > 4 values
>>>> >
>>>> > Number of Dimensions: 1
>>>> >
>>>> > Dimensions and sizes: [time | 4]
>>>> >
>>>> > Coordinates:
>>>> >
>>>> > time: [1025088..1025160]
>>>> >
>>>> > Number Of Attributes: 3
>>>> >
>>>> > units : hours since 1900-01-01 00:00:0.0
>>>> >
>>>> > long_name : time
>>>> >
>>>> > calendar : gregorian
>>>> >
>>>> >
>>>> >
>>>> >
>>>> >
>>>> > Variable: pnew
>>>> >
>>>> > Type: integer
>>>> >
>>>> > Total Size: 112 bytes
>>>> >
>>>> > 28 values
>>>> >
>>>> > Number of Dimensions: 1
>>>> >
>>>> > Dimensions and sizes: [level | 28]
>>>> >
>>>> > Coordinates:
>>>> >
>>>> > level: [70..1000]
>>>> >
>>>> > Number Of Attributes: 2
>>>> >
>>>> > units : Pa
>>>> >
>>>> > long_name : pressure_level
>>>> >
>>>> > (0) 7000
>>>> >
>>>> > (1) 10000
>>>> >
>>>> > (2) 12500
>>>> >
>>>> > (3) 15000
>>>> >
>>>> > ..........
>>>> >
>>>> > (26) 97500
>>>> >
>>>> > (27) 100000
>>>> >
>>>> >
>>>> >
>>>> > Variable: ZMT
>>>> >
>>>> > Type: float
>>>> >
>>>> > Total Size: 2213568 bytes
>>>> >
>>>> > 553392 values
>>>> >
>>>> > Number of Dimensions: 4
>>>> >
>>>> > Dimensions and sizes: [time | 4] x [level | 28] x [latitude | 61] x
>>>> > [longitude | 81]
>>>> >
>>>> > Coordinates:
>>>> >
>>>> > time: [1025088..1025160]
>>>> >
>>>> > level: [70..1000]
>>>> >
>>>> > latitude: [ 5..20]
>>>> >
>>>> > longitude: [70..90]
>>>> >
>>>> > Number Of Attributes: 3
>>>> >
>>>> > units : (m/s)(kg/kg)
>>>> >
>>>> > long_name : zonal moisture transport
>>>> >
>>>> > _FillValue : -32767
>>>> >
>>>> >
>>>> >
>>>> > At run-time, I am getting the following error message;
>>>> >
>>>> >
>>>> >
>>>> > fatal:Could not create (example.nc)
>>>> >
>>>> > warning:FileSetFileOption: invalid file or format
>>>> >
>>>> > fatal:["Execute.c":8640]:Execute: Error occurred at or near line 65 in
>>>> > file TEST1.ncl
>>>> >
>>>> >
>>>> >
>>>> > My script is attached below.
>>>> >
>>>> >
>>>> >
>>>> > Please kindly inform me what is wrong with this code and how to fix
>>>> it.
>>>> >
>>>> >
>>>> >
>>>> > Thanking you,
>>>> >
>>>> > Soma
>>>> >
>>>> >
>>>> >
>>>> > _______________________________________________
>>>> >
>>>> > ncl-talk mailing list
>>>> >
>>>> > ncl-talk at ucar.edu
>>>> >
>>>> > List instructions, subscriber options, unsubscribe:
>>>> >
>>>> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>> >
>>>> >
>>>> >
>>>> > _______________________________________________
>>>> > ncl-talk mailing list
>>>> > ncl-talk at ucar.edu
>>>> > List instructions, subscriber options, unsubscribe:
>>>> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>> >
>>>> >
>>>> >
>>>> >
>>>> > _______________________________________________
>>>> > ncl-talk mailing list
>>>> > ncl-talk at ucar.edu
>>>> > List instructions, subscriber options, unsubscribe:
>>>> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>> >
>>>> >
>>>> >
>>>> >
>>>> -------------- next part --------------
>>>> An HTML attachment was scrubbed...
>>>> URL: <
>>>> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180718/62f5dceb/attachment-0001.html
>>>> >
>>>>
>>>> ------------------------------
>>>>
>>>> Message: 5
>>>> Date: Wed, 18 Jul 2018 07:14:34 +0800
>>>> From: "Md. Jalal Uddin" <dmjalal90 at gmail.com>
>>>> To: Ncl-talk <ncl-talk at ucar.edu>
>>>> Subject: [ncl-talk] Vertical wind shear calculation within 400 km
>>>> radius of TC
>>>> Message-ID:
>>>> <
>>>> CAN5hBbsik39xNZ_vnSAWZPLpVe4g_7suzVnSWHEcRk0Lm56h2w at mail.gmail.com>
>>>> Content-Type: text/plain; charset="utf-8"
>>>>
>>>> Hi all,
>>>>
>>>> I wanted to calculate mean vertical wind shear for each pressure level
>>>> (950-100mb). I used u and v wind component from the NCEP Climate
>>>> Forecast
>>>> System Reanalysis (CFSR) 6-hourly Products.
>>>>
>>>> I am getting the following errors:
>>>> fatal:Number of subscripts do not match number of dimensions of
>>>> variable,(2) Subscripts used, (3) Subscripts expected
>>>> fatal:["Execute.c":8575]:Execute: Error occurred at or near line 41 in
>>>> file
>>>> test_vws_eachlevel.ncl
>>>>
>>>> Number of subscripts do not match number of dimensions of variable,(2)
>>>> Subscripts used, (3) Subscripts expected
>>>> fatal:["Execute.c":8575]:Execute: Error occurred at or near line 47 in
>>>> file
>>>> test_vws_eachlevel.ncl
>>>>
>>>> Please help.
>>>>
>>>> Cheers,
>>>> Jalal
>>>> --
>>>> *Md. Jalal Uddin*
>>>> MSc in Applied Meteorology (English Language)
>>>> Nanjing University of Information, Science and Technology, China
>>>> Jasmine Jiangsu Government Scholar
>>>> Cell: +8613260859092
>>>> Office: WMO RMTC building in middle campus
>>>>
>>>> *AND *
>>>> B.Sc. in Disaster Management (Hons.)
>>>> Patuakhali Science and Technology University, Bangladesh.
>>>> Cell: +8801792052662, +8801838613203
>>>> Email: dmjalal90 at gmail.com
>>>> Web: www.dmjalal90.weebly.com
>>>> Facebook: jalal.hossen.39
>>>> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
>>>> Twitter: dmjalal90
>>>> Skype: dmjalal90
>>>> -------------- next part --------------
>>>> An HTML attachment was scrubbed...
>>>> URL: <
>>>> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180718/5ca1ddec/attachment-0001.html
>>>> >
>>>> -------------- next part --------------
>>>> A non-text attachment was scrubbed...
>>>> Name: test_vws_eachlevel.ncl
>>>> Type: application/octet-stream
>>>> Size: 3004 bytes
>>>> Desc: not available
>>>> URL: <
>>>> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180718/5ca1ddec/attachment-0001.obj
>>>> >
>>>>
>>>> ------------------------------
>>>>
>>>> Message: 6
>>>> Date: Wed, 18 Jul 2018 09:11:01 -0600
>>>> From: Rick Brownrigg <brownrig at ucar.edu>
>>>> To: "Md. Jalal Uddin" <dmjalal90 at gmail.com>
>>>> Cc: Ncl-talk <ncl-talk at ucar.edu>
>>>> Subject: Re: [ncl-talk] Vertical wind shear calculation within 400 km
>>>> radius of TC
>>>> Message-ID:
>>>> <CAGKRhbE0_mW=
>>>> jD4Y4oT5fyKGEfd2uO3B0JZqCJqi6DDE+RUbCA at mail.gmail.com>
>>>> Content-Type: text/plain; charset="utf-8"
>>>>
>>>> Hi,
>>>>
>>>> Unfortunately, the script you sent does not align with the line numbers
>>>> of
>>>> the error messages -- line 47 is a blank line. In any case, the error
>>>> message is pretty self-explanatory -- you are making reference to a 3D
>>>> variable as though it were a 2D variable. You might place a few
>>>> strategic
>>>> "printVarSummary()" statements for the variables involved to make sure
>>>> your
>>>> variables are what you think they are.
>>>>
>>>> Hope that helps...
>>>> Rick
>>>>
>>>> On Tue, Jul 17, 2018 at 5:14 PM, Md. Jalal Uddin <dmjalal90 at gmail.com>
>>>> wrote:
>>>>
>>>> > Hi all,
>>>> >
>>>> > I wanted to calculate mean vertical wind shear for each pressure level
>>>> > (950-100mb). I used u and v wind component from the NCEP Climate
>>>> Forecast
>>>> > System Reanalysis (CFSR) 6-hourly Products.
>>>> >
>>>> > I am getting the following errors:
>>>> > fatal:Number of subscripts do not match number of dimensions of
>>>> > variable,(2) Subscripts used, (3) Subscripts expected
>>>> > fatal:["Execute.c":8575]:Execute: Error occurred at or near line 41 in
>>>> > file test_vws_eachlevel.ncl
>>>> >
>>>> > Number of subscripts do not match number of dimensions of variable,(2)
>>>> > Subscripts used, (3) Subscripts expected
>>>> > fatal:["Execute.c":8575]:Execute: Error occurred at or near line 47 in
>>>> > file test_vws_eachlevel.ncl
>>>> >
>>>> > Please help.
>>>> >
>>>> > Cheers,
>>>> > Jalal
>>>> > --
>>>> > *Md. Jalal Uddin*
>>>> > MSc in Applied Meteorology (English Language)
>>>> > Nanjing University of Information, Science and Technology, China
>>>> > Jasmine Jiangsu Government Scholar
>>>> > Cell: +8613260859092
>>>> > Office: WMO RMTC building in middle campus
>>>> >
>>>> > *AND *
>>>> > B.Sc. in Disaster Management (Hons.)
>>>> > Patuakhali Science and Technology University, Bangladesh.
>>>> > Cell: +8801792052662, +8801838613203
>>>> > Email: dmjalal90 at gmail.com
>>>> > Web: www.dmjalal90.weebly.com
>>>> > Facebook: jalal.hossen.39
>>>> > LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
>>>> > Twitter: dmjalal90
>>>> > Skype: dmjalal90
>>>> >
>>>> > _______________________________________________
>>>> > ncl-talk mailing list
>>>> > ncl-talk at ucar.edu
>>>> > List instructions, subscriber options, unsubscribe:
>>>> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>> >
>>>> >
>>>> -------------- next part --------------
>>>> An HTML attachment was scrubbed...
>>>> URL: <
>>>> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180718/9832e9a7/attachment-0001.html
>>>> >
>>>>
>>>> ------------------------------
>>>>
>>>> Message: 7
>>>> Date: Wed, 18 Jul 2018 20:54:35 +0530
>>>> From: Soma Roy <somaroy892 at gmail.com>
>>>> To: ncl-talk <ncl-talk at ucar.edu>
>>>> Subject: [ncl-talk] Skew-T diagram
>>>> Message-ID:
>>>> <CAAJRSQdZQV2Y5vCOxv4vsZgMaxC=
>>>> 7d_p9Et979qePL_dFB7rww at mail.gmail.com>
>>>> Content-Type: text/plain; charset="utf-8"
>>>>
>>>> Hello,
>>>>
>>>> Is there any function/script available to calculate LCL height and draw
>>>> the
>>>> region of LCL, CAPE ,CINE and LFC on skew-T diagram along with wind
>>>> speed
>>>> and direction.
>>>>
>>>> Kindly inform me.
>>>>
>>>> Thanks,
>>>> Soma
>>>> -------------- next part --------------
>>>> An HTML attachment was scrubbed...
>>>> URL: <
>>>> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180718/3710c8d4/attachment-0001.html
>>>> >
>>>>
>>>> ------------------------------
>>>>
>>>> Subject: Digest Footer
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>> ------------------------------
>>>>
>>>> End of ncl-talk Digest, Vol 176, Issue 22
>>>> *****************************************
>>>>
>>>
>>>
>>> --
>>> Tanya R. Peevey, PhD
>>> Research Scientist II, Global Observing Systems Analysis (GOSA) Group
>>> NOAA ESRL Global Systems Division
>>> 325 Broadway, Boulder, CO 80305
>>> (303) 497-5847
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>
> --
> Adam Phillips
> Associate Scientist, Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/ 303-497-1726
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180719/fd344aff/attachment.html>
More information about the ncl-talk
mailing list