Re: cannot read a netCDF 4 file with group

From: Wei Huang <huangwei_at_nyahnyahspammersnyahnyah>
Date: Fri Oct 05 2012 - 15:35:20 MDT

Wenshan,

We have some issues to read netCDF4 file with group in 6.0.0, and 6.1.0-beta.
We have fixed those issues and will have the 6.1.0 release soon.
If you need the new NCL version sooner, please send us your machine info:
uname -a
and we can try to have a binary for you.

Thanks,

Wei

huangwei@ucar.edu
VETS/CISL
National Center for Atmospheric Research
P.O. Box 3000 (1850 Table Mesa Dr.)
Boulder, CO 80307-3000 USA
(303) 497-8924

On Oct 5, 2012, at 3:23 PM, Wenshan Wang wrote:

> Dear all,
>
> I have a data file in netCDF 4, in_grp.nc (see the attachment). And there are several groups in the file.
>
> However, I cannot read it by addfile or ncl_filedump. I am using NCL 6.0.0 on Ubuntu 12.04 (32bit).
>
> The error message is:
> ncdiminq: ncid 65536: NetCDF: Invalid dimension ID or name
> ncdiminq: ncid 65536: NetCDF: Invalid dimension ID or name
> ncdiminq: ncid 65536: NetCDF: Invalid dimension ID or name
> ncdiminq: ncid 65536: NetCDF: Invalid dimension ID or name
> Segmentation fault (core dumped)
>
> Here is the structure of this file given by ncdump:
> netcdf in_grp {
> dimensions:
> lat = 2 ;
> lev = 3 ;
> lon = 4 ;
> time = UNLIMITED ; // (10 currently)
> variables:
> float lat(lat) ;
> lat:long_name = "Latitude (typically midpoints)" ;
> lat:units = "degrees_north" ;
> float lev(lev) ;
> lev:purpose = "Monotonically increasing coordinate pressure" ;
> lev:units = "hybrid_sigma_pressure" ;
> lev:positive = "down" ;
> lev:A_var = "hyam" ;
> lev:B_var = "hybm" ;
> lev:P0_var = "P0" ;
> lev:PS_var = "PS" ;
> lev:bounds = "ilev" ;
> float lon(lon) ;
> lon:long_name = "Longitude (typically midpoints)" ;
> lon:units = "degrees_east" ;
> double time(time) ;
> float scl ;
>
> // global attributes:
> :Conventions = "CF-1.0" ;
> :history = "History global attribute.\n" ;
> :julian_day = 200000.04 ;
> :RCS_Header = "$Header$" ;
> data:
>
> lat = -90, 90 ;
>
> lev = 100, 500, 1000 ;
>
> lon = 0, 90, 180, 270 ;
>
> time = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 ;
>
> scl = 1 ;
>
> group: level_1_group_1 {
> variables:
> double time(time) ;
> float lat(lat) ;
> float lev(lev) ;
> float lon(lon) ;
> lon:long_name = "level_1_group_1 Longitude (typically midpoints)" ;
> lon:units = "level_1_group_1 degrees_east" ;
> float scl ;
> data:
>
> time = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 ;
>
> lat = -90, 90 ;
>
> lev = 100, 500, 1000 ;
>
> lon = 0, 90, 180, 270 ;
>
> scl = 1.1 ;
>
> group: level_2_group_1 {
> variables:
> float scl ;
> data:
>
> scl = 1.1 ;
> } // group level_2_group_1
> } // group level_1_group_1
>
> group: level_1_group_2 {
> variables:
> double time(time) ;
> float scl ;
> data:
>
> time = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 ;
>
> scl = 1.2 ;
> } // group level_1_group_2
>
> group: level_1_group_3 {
> dimensions:
> rlev = 3 ;
> time2 = UNLIMITED ; // (20 currently)
> variables:
> float rz(rlev) ;
> double time2(time2) ;
> float scl ;
>
> // group attributes:
> :level_1_group_3_global_
> attribute = "level_1_group_3_global_attribute" ;
> data:
>
> rz = 0, 5000, 17000 ;
>
> time2 = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 ;
>
> scl = 1.3 ;
> } // group level_1_group_3
>
> group: level_1_group_4 {
>
> group: level_2_group_4 {
> dimensions:
> rlev = 3 ;
> variables:
> float rz(rlev) ;
> data:
>
> rz = 0, 5000, 17000 ;
>
> group: level_3_group_4 {
> dimensions:
> time2 = UNLIMITED ; // (5 currently)
> time3 = UNLIMITED ; // (6 currently)
> variables:
> double time2(time2) ;
> double time3(time3) ;
> data:
>
> time2 = 1, 2, 3, 4, 5 ;
>
> time3 = 1, 2, 3, 4, 5, 6 ;
> } // group level_3_group_4
> } // group level_2_group_4
> } // group level_1_group_4
>
> group: level_1_group_5 {
>
> group: level_2_group_5 {
>
> group: level_3_group_5 {
> } // group level_3_group_5
> } // group level_2_group_5
> } // group level_1_group_5
> }
>
> Thank you!
>
> Best regards,
> Wenshan
> <in_grp.nc>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Oct 5 15:36:24 2012

This archive was generated by hypermail 2.1.8 : Mon Oct 08 2012 - 15:54:16 MDT