int main(void) { float data[X_LEN * Y_LEN * Z_LEN]; int i, ncmode, ncid, dimids[NDIMS], var; size_t start[NDIMS] = {0, 0, 0}; size_t count[NDIMS] = {X_LEN, Y_LEN, Z_LEN}; ptrdiff_t stride[NDIMS] = {1, 1, 1}; /* Initialize data. */ for (i = 0; i < (X_LEN * Y_LEN * Z_LEN); i++) data[i] = i; printf("*** Testing netcdf-4 writes with compressed data...\n"); ncmode = NC_CLOBBER|NC_NETCDF4; if (nc_create(FILE_NAME, NC_NETCDF4, &ncid)) ERR_RET; if (nc_def_dim(ncid, "time", X_LEN, &dimids[0])) ERR_RET; if (nc_def_dim(ncid, "lat", Y_LEN, &dimids[1])) ERR_RET; if (nc_def_dim(ncid, "lon", Z_LEN, &dimids[2])) ERR_RET; if (nc_def_var(ncid, "test", NC_FLOAT, NDIMS, dimids, &var)) ERR_RET; if (nc_def_var_deflate(ncid, var, 1, 1, 2)) ERR_RET; if (nc_enddef(ncid)) ERR_RET; if (nc_put_vars_float(ncid, var, start, count, stride, data)) ERR_RET; if (nc_close(ncid)) ERR_RET; printf("*** Tests successful!\n"); return 0; }
int dump_file3(const float *data, int docompression, int usedefdim) { int ncmode, ncid, dimids[NDIMS], var; size_t start[NDIMS] = {0, 0, 0}, count[NDIMS] = {X_LEN, Y_LEN, Z_LEN}; ptrdiff_t stride[NDIMS] = {1, 1, 1}; ncmode = NC_CLOBBER|NC_NETCDF4; if (nc_create(FILE_NAME, ncmode, &ncid)) ERR_RET; if (nc_def_dim(ncid, "time", X_LEN, &dimids[0])) ERR_RET; if (nc_def_dim(ncid, "lat", Y_LEN, &dimids[1])) ERR_RET; if (nc_def_dim(ncid, "lon", Z_LEN, &dimids[2])) ERR_RET; if (nc_def_var(ncid, "test", NC_FLOAT, NDIMS, dimids, &var)) ERR_RET; if (docompression) if (nc_def_var_deflate(ncid, var, 1, 1, 1)) ERR_RET; if (nc_enddef(ncid)) ERR_RET; if (nc_put_vars_float(ncid, var, start, count, stride, data)) ERR_RET; if (nc_close(ncid)) ERR_RET; return 0; }
int ex_put_partial_one_attr( int exoid, ex_entity_type obj_type, ex_entity_id obj_id, int64_t start_num, int64_t num_ent, int attrib_index, const void *attrib ) { int status; int attrid, obj_id_ndx, temp; size_t num_entries_this_obj, num_attr; size_t start[2], count[2]; ptrdiff_t stride[2]; char errmsg[MAX_ERR_LENGTH]; const char* dnumobjent; const char* dnumobjatt; const char* vattrbname; exerrval = 0; /* clear error code */ /* Determine index of obj_id in id array */ if (obj_type != EX_NODAL) { obj_id_ndx = ex_id_lkup(exoid,obj_type,obj_id); if (exerrval != 0) { if (exerrval == EX_NULLENTITY) { sprintf(errmsg, "Warning: no attributes allowed for NULL %s %"PRId64" in file id %d", ex_name_of_object(obj_type),obj_id,exoid); ex_err("ex_put_partial_one_attr",errmsg,EX_MSG); return (EX_WARN); /* no attributes for this element block */ } else { sprintf(errmsg, "Error: no %s id %"PRId64" in id array in file id %d", ex_name_of_object(obj_type), obj_id, exoid); ex_err("ex_put_partial_one_attr",errmsg,exerrval); return (EX_FATAL); } } } switch (obj_type) { case EX_SIDE_SET: dnumobjent = DIM_NUM_SIDE_SS(obj_id_ndx); dnumobjatt = DIM_NUM_ATT_IN_SS(obj_id_ndx); vattrbname = VAR_SSATTRIB(obj_id_ndx); break; case EX_NODE_SET: dnumobjent = DIM_NUM_NOD_NS(obj_id_ndx); dnumobjatt = DIM_NUM_ATT_IN_NS(obj_id_ndx); vattrbname = VAR_NSATTRIB(obj_id_ndx); break; case EX_EDGE_SET: dnumobjent = DIM_NUM_EDGE_ES(obj_id_ndx); dnumobjatt = DIM_NUM_ATT_IN_ES(obj_id_ndx); vattrbname = VAR_ESATTRIB(obj_id_ndx); break; case EX_FACE_SET: dnumobjent = DIM_NUM_FACE_FS(obj_id_ndx); dnumobjatt = DIM_NUM_ATT_IN_FS(obj_id_ndx); vattrbname = VAR_FSATTRIB(obj_id_ndx); break; case EX_ELEM_SET: dnumobjent = DIM_NUM_ELE_ELS(obj_id_ndx); dnumobjatt = DIM_NUM_ATT_IN_ELS(obj_id_ndx); vattrbname = VAR_ELSATTRIB(obj_id_ndx); break; case EX_NODAL: dnumobjent = DIM_NUM_NODES; dnumobjatt = DIM_NUM_ATT_IN_NBLK; vattrbname = VAR_NATTRIB; break; case EX_EDGE_BLOCK: dnumobjent = DIM_NUM_ED_IN_EBLK(obj_id_ndx); dnumobjatt = DIM_NUM_ATT_IN_EBLK(obj_id_ndx); vattrbname = VAR_EATTRIB(obj_id_ndx); break; case EX_FACE_BLOCK: dnumobjent = DIM_NUM_FA_IN_FBLK(obj_id_ndx); dnumobjatt = DIM_NUM_ATT_IN_FBLK(obj_id_ndx); vattrbname = VAR_FATTRIB(obj_id_ndx); break; case EX_ELEM_BLOCK: dnumobjent = DIM_NUM_EL_IN_BLK(obj_id_ndx); dnumobjatt = DIM_NUM_ATT_IN_BLK(obj_id_ndx); vattrbname = VAR_ATTRIB(obj_id_ndx); break; default: exerrval = 1005; sprintf(errmsg, "Internal Error: unrecognized object type in switch: %d in file id %d", obj_type,exoid); ex_err("ex_putt_partial_one_attr",errmsg,EX_MSG); return (EX_FATAL); } /* inquire id's of previously defined dimensions */ if (ex_get_dimension(exoid, dnumobjent,"entries", &num_entries_this_obj, &temp, "ex_put_partial_one_attr") != NC_NOERR) return EX_FATAL; if (start_num + num_ent -1 > num_entries_this_obj) { exerrval = EX_BADPARAM; sprintf(errmsg, "Error: start index (%"PRId64") + count (%"PRId64") is larger than total number of entities (%"ST_ZU") in file id %d", start_num, num_ent, num_entries_this_obj, exoid); ex_err("ex_put_partial_one_attr",errmsg,exerrval); return (EX_FATAL); } if (ex_get_dimension(exoid, dnumobjatt,"attributes", &num_attr, &temp, "ex_put_partial_one_attr") != NC_NOERR) return EX_FATAL; if (attrib_index < 1 || attrib_index > (int)num_attr) { exerrval = EX_FATAL; sprintf(errmsg, "Error: Invalid attribute index specified: %d. Valid range is 1 to %"ST_ZU" for %s %"PRId64" in file id %d", attrib_index, num_attr, ex_name_of_object(obj_type), obj_id, exoid); ex_err("ex_put_partial_one_attr",errmsg,exerrval); return (EX_FATAL); } if ((status = nc_inq_varid(exoid, vattrbname, &attrid)) != NC_NOERR) { exerrval = status; sprintf(errmsg, "Error: failed to locate attribute variable for %s %"PRId64" in file id %d", ex_name_of_object(obj_type),obj_id,exoid); ex_err("ex_put_partial_one_attr",errmsg,exerrval); return (EX_FATAL); } /* write out the attributes */ start[0] = start_num-1; start[1] = attrib_index-1; count[0] = num_ent; count[1] = 1; stride[0] = 1; stride[1] = num_attr; if (count[0] == 0) start[0] = 0; if (ex_comp_ws(exoid) == 4) { status = nc_put_vars_float(exoid, attrid, start, count, stride, attrib); } else { status = nc_put_vars_double(exoid, attrid, start, count, stride, attrib); } if (status != NC_NOERR) { exerrval = status; sprintf(errmsg, "Error: failed to put attribute %d for %s %"PRId64" in file id %d", attrib_index, ex_name_of_object(obj_type), obj_id, exoid); ex_err("ex_put_partial_one_attr",errmsg,exerrval); return (EX_FATAL); } return(EX_NOERR); }
int ex_put_one_attr( int exoid, int obj_type, int obj_id, int attrib_index, const void *attrib ) { int numobjentdim, numattrdim, attrid, obj_id_ndx; long num_entries_this_obj, num_attr; size_t start[2], count[2]; ptrdiff_t stride[2]; int error; char errmsg[MAX_ERR_LENGTH]; const char* tname; const char* vobjids; const char* dnumobjent = 0; const char* dnumobjatt = 0; const char* vattrbname = 0; switch (obj_type) { case EX_EDGE_BLOCK: tname = "edge block"; vobjids = VAR_ID_ED_BLK; break; case EX_FACE_BLOCK: tname = "face block"; vobjids = VAR_ID_FA_BLK; break; case EX_ELEM_BLOCK: tname = "element block"; vobjids = VAR_ID_EL_BLK; break; default: exerrval = EX_BADPARAM; sprintf(errmsg, "Error: Bad block type (%d) specified for file id %d", obj_type, exoid ); ex_err("ex_put_attr",errmsg,exerrval); return (EX_FATAL); } exerrval = 0; /* clear error code */ /* Determine index of obj_id in vobjids array */ obj_id_ndx = ex_id_lkup(exoid,vobjids,obj_id); if (exerrval != 0) { if (exerrval == EX_NULLENTITY) { sprintf(errmsg, "Warning: no attributes allowed for NULL %s %d in file id %d", tname,obj_id,exoid); ex_err("ex_put_one_attr",errmsg,EX_MSG); return (EX_WARN); /* no attributes for this element block */ } else { sprintf(errmsg, "Error: no %s id %d in %s array in file id %d", tname, obj_id, vobjids, exoid); ex_err("ex_put_one_attr",errmsg,exerrval); return (EX_FATAL); } } switch (obj_type) { case EX_EDGE_BLOCK: dnumobjent = DIM_NUM_ED_IN_EBLK(obj_id_ndx); dnumobjatt = DIM_NUM_ATT_IN_EBLK(obj_id_ndx); vattrbname = VAR_EATTRIB(obj_id_ndx); break; case EX_FACE_BLOCK: dnumobjent = DIM_NUM_FA_IN_FBLK(obj_id_ndx); dnumobjatt = DIM_NUM_ATT_IN_FBLK(obj_id_ndx); vattrbname = VAR_FATTRIB(obj_id_ndx); break; case EX_ELEM_BLOCK: dnumobjent = DIM_NUM_EL_IN_BLK(obj_id_ndx); dnumobjatt = DIM_NUM_ATT_IN_BLK(obj_id_ndx); vattrbname = VAR_ATTRIB(obj_id_ndx); break; } /* inquire id's of previously defined dimensions */ if ((numobjentdim = ncdimid (exoid, dnumobjent)) == -1) { if (ncerr == NC_EBADDIM) { exerrval = ncerr; sprintf(errmsg, "Error: no %s with id %d in file id %d", tname, obj_id, exoid); ex_err("ex_put_one_attr",errmsg,exerrval); return (EX_FATAL); } else { exerrval = ncerr; sprintf(errmsg, "Error: failed to locate number of entries for %s %d in file id %d", tname, obj_id, exoid); ex_err("ex_put_one_attr",errmsg,exerrval); return (EX_FATAL); } } if (ncdiminq (exoid, numobjentdim, (char *) 0, &num_entries_this_obj) == -1) { exerrval = ncerr; sprintf(errmsg, "Error: failed to get number of entries for %s %d in file id %d", tname,obj_id,exoid); ex_err("ex_put_one_attr",errmsg,exerrval); return (EX_FATAL); } if ((numattrdim = ncdimid(exoid, dnumobjatt)) == -1) { exerrval = ncerr; sprintf(errmsg, "Error: number of attributes not defined for %s %d in file id %d", tname,obj_id,exoid); ex_err("ex_put_one_attr",errmsg,EX_MSG); return (EX_FATAL); /* number of attributes not defined */ } if (ncdiminq (exoid, numattrdim, (char *) 0, &num_attr) == -1) { exerrval = ncerr; sprintf(errmsg, "Error: failed to get number of attributes for block %d in file id %d", obj_id,exoid); ex_err("ex_put_one_attr",errmsg,exerrval); return (EX_FATAL); } if (attrib_index < 1 || attrib_index > num_attr) { exerrval = EX_FATAL; sprintf(errmsg, "Error: Invalid attribute index specified: %d. Valid range is 1 to %ld for %s %d in file id %d", attrib_index, num_attr, tname, obj_id, exoid); ex_err("ex_put_one_attr",errmsg,exerrval); return (EX_FATAL); } if ((attrid = ncvarid (exoid, vattrbname)) == -1) { exerrval = ncerr; sprintf(errmsg, "Error: failed to locate attribute variable for %s %d in file id %d", tname,obj_id,exoid); ex_err("ex_put_one_attr",errmsg,exerrval); return (EX_FATAL); } /* write out the attributes */ start[0] = 0; start[1] = attrib_index-1; count[0] = num_entries_this_obj; count[1] = 1; stride[0] = 1; stride[1] = num_attr; if (nc_flt_code(exoid) == NC_FLOAT) { error = nc_put_vars_float(exoid, attrid, start, count, stride, ex_conv_array(exoid,WRITE_CONVERT,attrib, (int)num_attr*num_entries_this_obj)); } else { error = nc_put_vars_double(exoid, attrid, start, count, stride, ex_conv_array(exoid,WRITE_CONVERT,attrib, (int)num_attr*num_entries_this_obj)); } if (error == -1) { exerrval = ncerr; sprintf(errmsg, "Error: failed to put attribute %d for %s %d in file id %d", attrib_index, tname, obj_id, exoid); ex_err("ex_put_one_attr",errmsg,exerrval); return (EX_FATAL); } return(EX_NOERR); }