NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.

To learn about what's going on, see About the Archive Site.

Re: beta1 bug

NOTE: The netcdf-hdf mailing list is no longer active. The list archives are made available for historical reasons.

Ed Hartnett wrote:

Jeff Whitaker <jswhit@xxxxxxxxxxx> writes:


Ed:  The attached C program creates a netcdf4_classic file that is
unreadable (ncdump'ing it gives an 'HDF Error').

-Jeff

--
Jeffrey S. Whitaker         Phone  : (303)497-6313
Meteorologist               FAX    : (303)497-6449
NOAA/OAR/PSD  R/PSD1        Email  : Jeffrey.S.Whitaker@xxxxxxxx
325 Broadway                Office : Skaggs Research Cntr 1D-124
Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg


#include <stdio.h>
#include <string.h>
#include <netcdf.h>
#define FILENAME "testfile.nc"
#define ATT1 "foo"
#define ATT2 "bar"
#define LONG_NAME "attribute"
#define ATT3 "long_name"
#define ERR(e) {printf("Error: %s\n", nc_strerror(e)); return 2;}
int
main()
{
       int ncid, retval, att;

       if ((retval = nc_set_default_format(NC_FORMAT_NETCDF4_CLASSIC, NULL))) 
ERR(retval);

       if ((retval = nc_create(FILENAME, NC_CLOBBER, &ncid))) ERR(retval);
       if ((retval = nc_enddef(ncid))) ERR(retval);

       att = 1;

       if ((retval = nc_redef(ncid))) ERR(retval);
       if ((retval = nc_put_att(ncid, NC_GLOBAL, ATT1, NC_INT, 1, &att))) 
ERR(retval);
       if ((retval = nc_enddef(ncid))) ERR(retval);

       att = 2;

       if ((retval = nc_redef(ncid))) ERR(retval);
       if ((retval = nc_put_att(ncid, NC_GLOBAL, ATT2, NC_INT, 1, &att))) 
ERR(retval);
       if ((retval = nc_enddef(ncid))) ERR(retval);

       if ((retval = nc_redef(ncid))) ERR(retval);
       if ((retval = nc_put_att_text(ncid, NC_GLOBAL, ATT3, strlen(LONG_NAME), 
LONG_NAME))) ERR(retval);
       if ((retval = nc_enddef(ncid))) ERR(retval);

       if ((retval = nc_close(ncid))) ERR(retval);
}



Thanks Jeff!

This is not an ncdump bug; nc_open fails on this file!

I have added it to the netCDF-4 tests and will let you know when I get
it solved...

Thanks,

Ed
OK Ed - thanks for looking into it. I'm running into another potential bug, probably in the HDF5 lib. I get 

H5Oalloc.c:1350: failed assertion `oh->chunk[curr_msg->chunkno].gap == 0'
Abort
when trying to rename a variable in my python script. Haven't been able to trigger the error in a equivalent C program yet, but if I do I'll forward it. 

-Jeff

--
Jeffrey S. Whitaker         Phone  : (303)497-6313
Meteorologist               FAX    : (303)497-6449
NOAA/OAR/PSD  R/PSD1        Email  : Jeffrey.S.Whitaker@xxxxxxxx
325 Broadway                Office : Skaggs Research Cntr 1D-124
Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg

==============================================================================
To unsubscribe netcdf-hdf, visit:
http://www.unidata.ucar.edu/mailing-list-delete-form.html
==============================================================================
  • 2007 messages navigation, sorted by:
    1. Thread
    2. Subject
    3. Author
    4. Date
    5. ↑ Table Of Contents
  • Search the netcdf-hdf archives: