[ncl-talk] ncl-talk Digest, Vol 176, Issue 22
Adam Phillips
asphilli at ucar.edu
Thu Jul 19 11:47:45 MDT 2018
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/1e1011a1/attachment.html>
More information about the ncl-talk
mailing list