[ncl-talk] ncl-talk Digest, Vol 176, Issue 22
Tanya Peevey
trpeevey at gmail.com
Thu Jul 19 11:11:16 MDT 2018
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180719/5dbd481c/attachment.html>
More information about the ncl-talk
mailing list