I don't see a call to end definition mode - nc_enddef - needs to be
done after defining dimensions and attributes and before writing data.
-Roy
On May 25, 2009, at 8:52 AM, Marston Johnston wrote:
Hi,
I'm trying to write a series of arrays to a single netcdf file that
will be about 1/2 Gb when done. The datasets are 1D, 2D, and 3D
arrays of packed data.
I'm having trouble writing the data correctly using the examples
given in the manual and the netcdf example directory.
My problem is that when I check the data before writing, everything
looks fine but when I check the data written to netcdf file the
arrays are all screwed up with the max and mins for each level being
incorrect and some slices of the 3D array being garbled.
any help would be greatly appreciated!
/Marston
Here's how I setup the file for writing:
// create netCDF dataset with perfomance options: enter define mode
status = nc_create(outfile,NC_64BIT_OFFSET,&ncid);
if(status != NC_NOERR) handle_error(status);
// Define a title for the file
status = nc_put_att_text(ncid,NC_GLOBAL,"Title",strlen(title),title);
// define dimensions: from name and length
status = nc_def_dim(ncid,"Lat",latlen,&lat_dimid);
if (status != NC_NOERR) handle_error(status);
status = nc_def_dim(ncid,"Lon",lonlen,&lon_dimid);
if (status != NC_NOERR) handle_error(status);
status = nc_def_dim(ncid,"Lev",levlen,&lev_dimid);
if (status != NC_NOERR) handle_error(status);
dim[0] = lev_dimid;
dim2d[0] = lat_dimid;
dim2d[1] = lon_dimid;
dim3d[0] = lev_dimid;
dim3d[1] = lat_dimid;
dim3d[2] = lon_dimid;
Print("Define variables: from name, type");
status = nc_def_var(ncid,"Latitude",NC_FLOAT,2,dim2d,&laid);
if(status != NC_NOERR) handle_error(status);
status = nc_def_var(ncid,"Longitude",NC_FLOAT,2,dim2d,&loid);
if(status != NC_NOERR) handle_error(status);
status = nc_def_var(ncid,"Levels",NC_SHORT,1,dim,&leid);
if(status != NC_NOERR) handle_error(status);
// Setting up the temperature variable
status = nc_def_var(ncid,"Temperature",NC_SHORT,3,dim3d,&tid);
if(status != NC_NOERR) handle_error(status);
status =
nc_put_att_text(ncid,tid,"Units",strlen(TEMP_UNITS),TEMP_UNITS);
if (status != NC_NOERR) handle_error(status);
pos = RetPos(sptr,130);
status = nc_put_att_double(ncid,tid,"scale",NC_DOUBLE,1,&sptr[pos]-
>scale);
if (status != NC_NOERR) handle_error(status);
status = nc_put_att_double(ncid,tid,"offset",NC_DOUBLE,1,&sptr[pos]-
>offset);
if (status != NC_NOERR) handle_error(status);
// Fetching the data, packing it, checking it before writing to file
pos = RetPos(sptr,130);
swd->t = GetNlevsP(f,wk->array_size,sptr[pos]->nlev,sptr[pos]-
>pnum,sptr[pos]->scale,sptr[pos]->offset);
// Checking the Maxs and mins for each level
int lev;
for(lev = 0; lev < wk->nlev; lev++) {
printf("Lev: %d\t",lev);
MMSI(swd->t[lev],wk->nlat,wk->nlon,sptr[pos]->scale,sptr[pos]-
>offset);
}
status = nc_put_var_short(ncid,tid,&swd->t[0][0][0]);
if(status != NC_NOERR) handle_error(status);
Marston Johnston, PhD Student
Radio and Space Science
Chalmers University of Technology
Gothenburg Sweden
marston at chalmers.se
cell: +46 736 637 909
Skype ID: marston.johnston
_______________________________________________
netcdfgroup mailing list
netcdfgroup@xxxxxxxxxxxxxxxx
For list information or to unsubscribe, visit:
http://www.unidata.ucar.edu/mailing_lists/
**********************
"The contents of this message do not reflect any position of the U.S.
Government or NOAA."
**********************
Roy Mendelssohn
Supervisory Operations Research Analyst
NOAA/NMFS
Environmental Research Division
Southwest Fisheries Science Center
1352 Lighthouse Avenue
Pacific Grove, CA 93950-2097
e-mail: Roy.Mendelssohn@xxxxxxxx (Note new e-mail address)
voice: (831)-648-9029
fax: (831)-648-8440
www: http://www.pfeg.noaa.gov/
"Old age and treachery will overcome youth and skill."
"From those who have been given much, much will be expected"
