[ncl-talk] ncl-talk Digest, Vol 176, Issue 22

David Adams dave.k.adams at gmail.com
Thu Jul 19 11:28:29 MDT 2018


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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180719/a1ad4f15/attachment-0001.html>


More information about the ncl-talk mailing list