[netcdfgroup] Writing files > 4GB with 16-bit variales

• To: netcdfgroup@xxxxxxxxxxxxxxxx
• Subject: [netcdfgroup] Writing files > 4GB with 16-bit variales
• From: Natarajan CS <csnataraj@xxxxxxxxx>
• Date: Wed, 4 Nov 2009 10:08:50 -0600

Hello all,
     I am trying to write a large file by combining multiple files. When I
try to write data (type short) greater than 4GB to a file, I get the
following error  *"ncx.c:1810: ncx_put_size_t: Assertion `*ulp <=
4294967295U' failed."
*I was wondering if this is an issue with my code or if it is a known and
resolved issue or even if there is a workaround? Appreciate any help!

Cheers.
Code :
...........
short data[1490][1636][1720];
int NX=1720, NY=1636, NZ=1490;
for (k=1; k<=2; k++) {
                sprintf(name,"test%d.nc",k);
                nc_open(name, NC_NOWRITE, &ncid);

                if (nc_inq_ndims(ncid,&ndims) != NC_NOERR)
                        printf("Another error");
                size_t dim_siz[ndims];
                printf("The file has %d dimensions\n",ndims);
                for (dim =0; dim < ndims; dim++)  {
                        nc_inq_dim(ncid, dim, dim_name, &dim_siz[dim]);
                        printf("Name of dimension %d is %s\n",dim,
dim_name);
                        printf("Size of dimension %d is %d\n",dim,
dim_siz[dim]);
                }
                int varid;
                if (nc_inq_varid(ncid, "segmented", &varid) != NC_NOERR)
                        printf("Error 1");
                if (nc_get_var_short(ncid, varid, &data[745*(k-1)) !=
NC_NOERR)
                        printf("Error 2");
                nc_close(ncid);
        }
        printf("Done reading data\n");
        if(retval=nc_create(FILE_NAME, NC_CLOBBER | NC_64BIT_OFFSET,
&ncid2))
                ERR(retval);
        if ((retval = nc_def_dim(ncid, "x", NX, &x_dimid)))
                ERR(retval);
        if ((retval = nc_def_dim(ncid, "y", NY, &y_dimid)))
                ERR(retval);
        if ((retval = nc_def_dim(ncid, "z", NZ, &z_dimid)))
                ERR(retval);


        dimids[0] = x_dimid;
        dimids[1] = y_dimid;
        dimids[2] = z_dimid;

        if ((retval = nc_def_var(ncid, "segmented", NC_SHORT, 3, dimids,
&varid)))
                ERR(retval);

        if ((retval = nc_enddef(ncid)))
                ERR(retval);

        if ((retval = nc_put_var_short(ncid, varid, &data[0][0][0])))
                ERR(retval);

        if ((retval = nc_close(ncid)))
                ERR(retval);
   printf("*** SUCCESS writing example file simple_xy.nc!\n");


   return 0;
}

• Follow-Ups:
  • Re: [netcdfgroup] Writing files > 4GB with 16-bit variales
    • From: Mark Rivers