mirror of
https://github.com/Unidata/netcdf-c.git
synced 2024-11-27 07:30:33 +08:00
3db4f013bf
Specific changes:
1. Add dap4 code: libdap4 and dap4_test.
Note that until the d4ts server problem is solved, dap4 is turned off.
2. Modify various files to support dap4 flags:
configure.ac, Makefile.am, CMakeLists.txt, etc.
3. Add nc_test/test_common.sh. This centralizes
the handling of the locations of various
things in the build tree: e.g. where is
ncgen.exe located. See nc_test/test_common.sh
for details.
4. Modify .sh files to use test_common.sh
5. Obsolete separate oc2 by moving it to be part of
netcdf-c. This means replacing code with netcdf-c
equivalents.
5. Add --with-testserver to configure.ac to allow
override of the servers to be used for --enable-dap-remote-tests.
6. There were multiple versions of nctypealignment code. Try to
centralize in libdispatch/doffset.c and include/ncoffsets.h
7. Add a unit test for the ncuri code because of its complexity.
8. Move the findserver code out of libdispatch and into
a separate, self contained program in ncdap_test and dap4_test.
9. Move the dispatch header files (nc{3,4}dispatch.h) to
.../include because they are now shared by modules.
10. Revamp the handling of TOPSRCDIR and TOPBUILDDIR for shell scripts.
11. Make use of MREMAP if available
12. Misc. minor changes e.g.
- #include <config.h> -> #include "config.h"
- Add some no-install headers to /include
- extern -> EXTERNL and vice versa as needed
- misc header cleanup
- clean up checking for misc. unix vs microsoft functions
13. Change copyright decls in some files to point to LICENSE file.
14. Add notes to RELEASENOTES.md
1221 lines
22 KiB
C
1221 lines
22 KiB
C
/** \file
|
|
The V2 API Functions.
|
|
|
|
Copyright 1996, University Corporation for Atmospheric Research
|
|
See \ref copyright file for copying and redistribution conditions.
|
|
*/
|
|
|
|
#ifndef NO_NETCDF_2
|
|
|
|
#include "config.h"
|
|
#include <stdlib.h>
|
|
#include <stdio.h>
|
|
#include <stdarg.h>
|
|
#include "netcdf.h"
|
|
#include "math.h"
|
|
/* The subroutines in error.c emit no messages unless NC_VERBOSE bit
|
|
* is on. They call exit() when NC_FATAL bit is on. */
|
|
int ncopts = (NC_FATAL | NC_VERBOSE) ;
|
|
int ncerr = NC_NOERR ;
|
|
|
|
#if SIZEOF_LONG == SIZEOF_SIZE_T
|
|
/*
|
|
* We don't have to copy the arguments to switch from 'long'
|
|
* to 'size_t' or 'ptrdiff_t'. Use dummy macros.
|
|
*/
|
|
|
|
# define NDIMS_DECL
|
|
# define A_DECL(name, type, ndims, rhs) \
|
|
const type *const name = ((const type *)(rhs))
|
|
|
|
# define A_FREE(name)
|
|
|
|
# define A_INIT(lhs, type, ndims, rhs)
|
|
|
|
#else
|
|
/*
|
|
* We do have to copy the arguments to switch from 'long'
|
|
* to 'size_t' or 'ptrdiff_t'. In my tests on an SGI,
|
|
* any additional cost was lost in measurement variation.
|
|
*
|
|
* This stanza is true on Windows with MinGW-64
|
|
*/
|
|
|
|
# include "onstack.h"
|
|
|
|
static int
|
|
nvdims(int ncid, int varid)
|
|
{
|
|
int ndims=-1, status;
|
|
|
|
if ((status = nc_inq_varndims(ncid, varid, &ndims)))
|
|
{
|
|
nc_advise("ncvdims", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return ndims;
|
|
}
|
|
|
|
/* Used to avoid errors on 64-bit windows related to
|
|
c89 macros and flow control/conditionals. */
|
|
static void* nvmalloc(off_t size) {
|
|
if(size < 0)
|
|
return NULL;
|
|
|
|
return malloc(size);
|
|
|
|
}
|
|
|
|
#define NDIMS_DECL const int ndims = nvdims(ncid, varid); \
|
|
|
|
|
|
# define A_DECL(name, type, ndims, rhs) \
|
|
type *const name = (type*) nvmalloc((ndims) * sizeof(type))
|
|
|
|
|
|
#if 0
|
|
ALLOC_ONSTACK(name, type, ndims)
|
|
#endif
|
|
|
|
# define A_FREE(name) \
|
|
FREE_ONSTACK(name)
|
|
|
|
# define A_INIT(lhs, type, ndims, rhs) \
|
|
{ \
|
|
if((off_t)ndims >= 0) { \
|
|
const long *lp = rhs; \
|
|
type *tp = lhs; \
|
|
type *const end = lhs + ndims; \
|
|
while(tp < end) \
|
|
{ \
|
|
*tp++ = (type) *lp++; \
|
|
} \
|
|
} \
|
|
} \
|
|
\
|
|
if ((off_t)ndims < 0) {nc_advise("nvdims",NC_EMAXDIMS,"ndims %d",ndims); return -1;}
|
|
|
|
|
|
#endif
|
|
|
|
typedef signed char schar;
|
|
|
|
/*
|
|
* Computes number of record variables in an open netCDF file, and an array of
|
|
* the record variable ids, if the array parameter is non-null.
|
|
*/
|
|
static int
|
|
numrecvars(int ncid, int* nrecvarsp, int *recvarids)
|
|
{
|
|
int status = NC_NOERR;
|
|
int nvars = 0;
|
|
int ndims = 0;
|
|
int nrecvars = 0;
|
|
int varid;
|
|
int recdimid;
|
|
int dimids[MAX_NC_DIMS];
|
|
|
|
status = nc_inq_nvars(ncid, &nvars);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
|
|
status = nc_inq_unlimdim(ncid, &recdimid);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
|
|
if (recdimid == -1) {
|
|
*nrecvarsp = 0;
|
|
return NC_NOERR;
|
|
}
|
|
nrecvars = 0;
|
|
for (varid = 0; varid < nvars; varid++) {
|
|
status = nc_inq_varndims(ncid, varid, &ndims);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
status = nc_inq_vardimid(ncid, varid, dimids);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
if (ndims > 0 && dimids[0] == recdimid) {
|
|
if (recvarids != NULL)
|
|
recvarids[nrecvars] = varid;
|
|
nrecvars++;
|
|
}
|
|
}
|
|
*nrecvarsp = nrecvars;
|
|
return NC_NOERR;
|
|
}
|
|
|
|
|
|
/*
|
|
* Computes record size (in bytes) of the record variable with a specified
|
|
* variable id. Returns size as 0 if not a record variable.
|
|
*/
|
|
static int
|
|
ncrecsize(int ncid, int varid, size_t *recsizep)
|
|
{
|
|
int status = NC_NOERR;
|
|
int recdimid;
|
|
nc_type type;
|
|
int ndims;
|
|
int dimids[MAX_NC_DIMS];
|
|
int id;
|
|
int size;
|
|
|
|
*recsizep = 0;
|
|
status = nc_inq_unlimdim(ncid, &recdimid);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
status = nc_inq_vartype(ncid, varid, &type);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
status = nc_inq_varndims(ncid, varid, &ndims);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
status = nc_inq_vardimid(ncid, varid, dimids);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
if (ndims == 0 || dimids[0] != recdimid) {
|
|
return NC_NOERR;
|
|
}
|
|
size = nctypelen(type);
|
|
for (id = 1; id < ndims; id++) {
|
|
size_t len;
|
|
status = nc_inq_dimlen(ncid, dimids[id], &len);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
size *= (int)len;
|
|
}
|
|
*recsizep = (size_t)size;
|
|
return NC_NOERR;
|
|
}
|
|
|
|
|
|
/*
|
|
* Retrieves the dimension sizes of a variable with a specified variable id in
|
|
* an open netCDF file. Returns -1 on error.
|
|
*/
|
|
static int
|
|
dimsizes(int ncid, int varid, size_t *sizes)
|
|
{
|
|
int status = NC_NOERR;
|
|
int ndims;
|
|
int id;
|
|
int dimids[MAX_NC_DIMS];
|
|
|
|
status = nc_inq_varndims(ncid, varid, &ndims);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
status = nc_inq_vardimid(ncid, varid, dimids);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
if (ndims == 0 || sizes == NULL)
|
|
return NC_NOERR;
|
|
for (id = 0; id < ndims; id++) {
|
|
size_t len;
|
|
status = nc_inq_dimlen(ncid, dimids[id], &len);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
sizes[id] = len;
|
|
}
|
|
return NC_NOERR;
|
|
}
|
|
|
|
|
|
/*
|
|
* Retrieves the number of record variables, the record variable ids, and the
|
|
* record size of each record variable. If any pointer to info to be returned
|
|
* is null, the associated information is not returned. Returns -1 on error.
|
|
*/
|
|
int
|
|
nc_inq_rec(
|
|
int ncid,
|
|
size_t *nrecvarsp,
|
|
int *recvarids,
|
|
size_t *recsizes)
|
|
{
|
|
int status = NC_NOERR;
|
|
int nvars = 0;
|
|
int recdimid;
|
|
int varid;
|
|
int rvarids[MAX_NC_VARS];
|
|
int nrvars = 0;
|
|
|
|
status = nc_inq_nvars(ncid, &nvars);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
|
|
status = nc_inq_unlimdim(ncid, &recdimid);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
|
|
if (recdimid == -1)
|
|
return NC_NOERR;
|
|
|
|
status = numrecvars(ncid, &nrvars, rvarids);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
|
|
if (nrecvarsp != NULL)
|
|
*nrecvarsp = (size_t)nrvars;
|
|
|
|
if (recvarids != NULL)
|
|
for (varid = 0; varid < nrvars; varid++)
|
|
recvarids[varid] = rvarids[varid];
|
|
|
|
if (recsizes != NULL)
|
|
for (varid = 0; varid < nrvars; varid++) {
|
|
size_t rsize;
|
|
status = ncrecsize(ncid, rvarids[varid], &rsize);
|
|
if (status != NC_NOERR)
|
|
return status;
|
|
recsizes[varid] = rsize;
|
|
}
|
|
return NC_NOERR;
|
|
}
|
|
|
|
|
|
/*
|
|
* Write one record's worth of data, except don't write to variables for which
|
|
* the address of the data to be written is NULL. Return -1 on error. This is
|
|
* the same as the ncrecput() in the library, except that can handle errors
|
|
* better.
|
|
*/
|
|
int
|
|
nc_put_rec(
|
|
int ncid,
|
|
size_t recnum,
|
|
void* const* datap)
|
|
{
|
|
int status = NC_NOERR;
|
|
int varid;
|
|
int rvarids[MAX_NC_VARS];
|
|
int nrvars;
|
|
size_t start[MAX_NC_DIMS];
|
|
size_t edges[MAX_NC_DIMS];
|
|
|
|
status = numrecvars(ncid, &nrvars, rvarids);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
|
|
if (nrvars == 0)
|
|
return NC_NOERR;
|
|
|
|
start[0] = recnum;
|
|
for (varid = 1; varid < nrvars; varid++)
|
|
start[varid] = 0;
|
|
|
|
for (varid = 0; varid < nrvars; varid++) {
|
|
if (datap[varid] != NULL) {
|
|
status = dimsizes(ncid, rvarids[varid], edges);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
|
|
edges[0] = 1; /* only 1 record's worth */
|
|
status = nc_put_vara(ncid, rvarids[varid], start, edges, datap[varid]);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
}
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
/*
|
|
* Read one record's worth of data, except don't read from variables for which
|
|
* the address of the data to be read is null. Return -1 on error. This is
|
|
* the same as the ncrecget() in the library, except that can handle errors
|
|
* better.
|
|
*/
|
|
int
|
|
nc_get_rec(
|
|
int ncid,
|
|
size_t recnum,
|
|
void **datap)
|
|
{
|
|
int status = NC_NOERR;
|
|
int varid;
|
|
int rvarids[MAX_NC_VARS];
|
|
int nrvars;
|
|
size_t start[MAX_NC_DIMS];
|
|
size_t edges[MAX_NC_DIMS];
|
|
|
|
status = numrecvars(ncid, &nrvars, rvarids);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
|
|
if (nrvars == 0)
|
|
return NC_NOERR;
|
|
|
|
start[0] = recnum;
|
|
for (varid = 1; varid < nrvars; varid++)
|
|
start[varid] = 0;
|
|
|
|
for (varid = 0; varid < nrvars; varid++) {
|
|
if (datap[varid] != NULL) {
|
|
status = dimsizes(ncid, rvarids[varid], edges);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
edges[0] = 1; /* only 1 record's worth */
|
|
status = nc_get_vara(ncid, rvarids[varid], start, edges, datap[varid]);
|
|
if(status != NC_NOERR)
|
|
return status;
|
|
}
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
/*
|
|
*/
|
|
void
|
|
nc_advise(const char *routine_name, int err, const char *fmt,...)
|
|
{
|
|
va_list args;
|
|
|
|
if(NC_ISSYSERR(err))
|
|
ncerr = NC_SYSERR;
|
|
else
|
|
ncerr = err;
|
|
|
|
if( ncopts & NC_VERBOSE )
|
|
{
|
|
(void) fprintf(stderr,"%s: ", routine_name);
|
|
va_start(args ,fmt);
|
|
(void) vfprintf(stderr,fmt,args);
|
|
va_end(args);
|
|
if(err != NC_NOERR)
|
|
{
|
|
(void) fprintf(stderr,": %s",
|
|
nc_strerror(err));
|
|
}
|
|
(void) fputc('\n',stderr);
|
|
(void) fflush(stderr); /* to ensure log files are current */
|
|
}
|
|
|
|
if( (ncopts & NC_FATAL) && err != NC_NOERR )
|
|
{
|
|
exit(ncopts);
|
|
}
|
|
}
|
|
|
|
/* End error handling */
|
|
|
|
int
|
|
nccreate(const char* path, int cmode)
|
|
{
|
|
int ncid;
|
|
const int status = nc_create(path, cmode, &ncid);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("nccreate", status, "filename \"%s\"", path);
|
|
return -1;
|
|
}
|
|
return ncid;
|
|
}
|
|
|
|
|
|
int
|
|
ncopen(const char *path, int mode)
|
|
{
|
|
int ncid;
|
|
const int status = nc_open(path, mode, &ncid);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncopen", status, "filename \"%s\"", path);
|
|
return -1;
|
|
}
|
|
return ncid;
|
|
}
|
|
|
|
|
|
int
|
|
ncredef(int ncid)
|
|
{
|
|
const int status = nc_redef(ncid);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncredef", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
int
|
|
ncendef(int ncid)
|
|
{
|
|
const int status = nc_enddef(ncid);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncendef", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
int
|
|
ncclose(int ncid)
|
|
{
|
|
const int status = nc_close(ncid);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncclose", status, "ncid %d", ncid);
|
|
return -1;
|
|
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
int
|
|
ncinquire(
|
|
int ncid,
|
|
int* ndims,
|
|
int* nvars,
|
|
int* natts,
|
|
int* recdim
|
|
)
|
|
{
|
|
int nd, nv, na;
|
|
const int status = nc_inq(ncid, &nd, &nv, &na, recdim);
|
|
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncinquire", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
/* else */
|
|
|
|
if(ndims != NULL)
|
|
*ndims = (int) nd;
|
|
|
|
if(nvars != NULL)
|
|
*nvars = (int) nv;
|
|
|
|
if(natts != NULL)
|
|
*natts = (int) na;
|
|
|
|
return ncid;
|
|
}
|
|
|
|
|
|
int
|
|
ncsync(int ncid)
|
|
{
|
|
const int status = nc_sync(ncid);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncsync", status, "ncid %d", ncid);
|
|
return -1;
|
|
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
int
|
|
ncabort(int ncid)
|
|
{
|
|
const int status = nc_abort(ncid);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncabort", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
int
|
|
ncdimdef(
|
|
int ncid,
|
|
const char* name,
|
|
long length
|
|
)
|
|
{
|
|
int dimid;
|
|
int status = NC_NOERR;
|
|
if(length < 0) {
|
|
status = NC_EDIMSIZE;
|
|
nc_advise("ncdimdef", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
status = nc_def_dim(ncid, name, (size_t)length, &dimid);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncdimdef", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return dimid;
|
|
}
|
|
|
|
|
|
int
|
|
ncdimid(int ncid, const char* name)
|
|
{
|
|
int dimid;
|
|
const int status = nc_inq_dimid(ncid, name, &dimid);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncdimid", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return dimid;
|
|
}
|
|
|
|
|
|
int
|
|
ncdiminq(
|
|
int ncid,
|
|
int dimid,
|
|
char* name,
|
|
long* length
|
|
)
|
|
{
|
|
size_t ll;
|
|
const int status = nc_inq_dim(ncid, dimid, name, &ll);
|
|
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncdiminq", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
/* else */
|
|
|
|
if(length != NULL)
|
|
*length = (int) ll;
|
|
|
|
return dimid;
|
|
}
|
|
|
|
|
|
int
|
|
ncdimrename(
|
|
int ncid,
|
|
int dimid,
|
|
const char* name
|
|
)
|
|
{
|
|
const int status = nc_rename_dim(ncid, dimid, name);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncdimrename", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return dimid;
|
|
}
|
|
|
|
|
|
int
|
|
ncvardef(
|
|
int ncid,
|
|
const char* name,
|
|
nc_type datatype,
|
|
int ndims,
|
|
const int* dim
|
|
)
|
|
{
|
|
int varid = -1;
|
|
const int status = nc_def_var(ncid, name, datatype, ndims, dim, &varid);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncvardef", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return varid;
|
|
}
|
|
|
|
|
|
int
|
|
ncvarid(
|
|
int ncid,
|
|
const char* name
|
|
)
|
|
{
|
|
int varid = -1;
|
|
const int status = nc_inq_varid(ncid, name, &varid);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncvarid", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return varid;
|
|
}
|
|
|
|
|
|
int
|
|
ncvarinq(
|
|
int ncid,
|
|
int varid,
|
|
char* name,
|
|
nc_type* datatype,
|
|
int* ndims,
|
|
int* dim,
|
|
int* natts
|
|
)
|
|
{
|
|
int nd, na;
|
|
const int status = nc_inq_var(ncid, varid, name, datatype,
|
|
&nd, dim, &na);
|
|
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncvarinq", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
/* else */
|
|
|
|
if(ndims != NULL)
|
|
*ndims = (int) nd;
|
|
|
|
if(natts != NULL)
|
|
*natts = (int) na;
|
|
|
|
return varid;
|
|
}
|
|
|
|
|
|
int
|
|
ncvarput1(
|
|
int ncid,
|
|
int varid,
|
|
const long* index,
|
|
const void* value
|
|
)
|
|
{
|
|
NDIMS_DECL
|
|
A_DECL(coordp, size_t, (size_t)ndims, index);
|
|
A_INIT(coordp, size_t, (size_t)ndims, index);
|
|
{
|
|
const int status = nc_put_var1(ncid, varid, coordp, value);
|
|
A_FREE(coordp);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncvarput1", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
int
|
|
ncvarget1(
|
|
int ncid,
|
|
int varid,
|
|
const long* index,
|
|
void* value
|
|
)
|
|
{
|
|
NDIMS_DECL
|
|
A_DECL(coordp, size_t, ndims, index);
|
|
A_INIT(coordp, size_t, ndims, index);
|
|
{
|
|
const int status = nc_get_var1(ncid, varid, coordp, value);
|
|
A_FREE(coordp);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncdimid", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
int
|
|
ncvarput(
|
|
int ncid,
|
|
int varid,
|
|
const long* start,
|
|
const long* count,
|
|
const void* value
|
|
)
|
|
{
|
|
NDIMS_DECL
|
|
A_DECL(stp, size_t, ndims, start);
|
|
A_DECL(cntp, size_t, ndims, count);
|
|
A_INIT(stp, size_t, ndims, start);
|
|
A_INIT(cntp, size_t, ndims, count);
|
|
{
|
|
const int status = nc_put_vara(ncid, varid, stp, cntp, value);
|
|
A_FREE(cntp);
|
|
A_FREE(stp);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncvarput", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
int
|
|
ncvarget(
|
|
int ncid,
|
|
int varid,
|
|
const long* start,
|
|
const long* count,
|
|
void* value
|
|
)
|
|
{
|
|
NDIMS_DECL
|
|
A_DECL(stp, size_t, ndims, start);
|
|
A_DECL(cntp, size_t, ndims, count);
|
|
A_INIT(stp, size_t, ndims, start);
|
|
A_INIT(cntp, size_t, ndims, count);
|
|
{
|
|
const int status = nc_get_vara(ncid, varid, stp, cntp, value);
|
|
A_FREE(cntp);
|
|
A_FREE(stp);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncvarget", status, "ncid %d; varid %d", ncid, varid);
|
|
return -1;
|
|
}
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
int
|
|
ncvarputs(
|
|
int ncid,
|
|
int varid,
|
|
const long* start,
|
|
const long* count,
|
|
const long* stride,
|
|
const void* value
|
|
)
|
|
{
|
|
if(stride == NULL)
|
|
return ncvarput(ncid, varid, start, count, value);
|
|
/* else */
|
|
{
|
|
|
|
NDIMS_DECL
|
|
A_DECL(stp, size_t, ndims, start);
|
|
A_DECL(cntp, size_t, ndims, count);
|
|
A_DECL(strdp, ptrdiff_t, ndims, stride);
|
|
A_INIT(stp, size_t, ndims, start);
|
|
A_INIT(cntp, size_t, ndims, count);
|
|
A_INIT(strdp, ptrdiff_t, ndims, stride);
|
|
{
|
|
const int status = nc_put_vars(ncid, varid, stp, cntp, strdp, value);
|
|
A_FREE(strdp);
|
|
A_FREE(cntp);
|
|
A_FREE(stp);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncvarputs", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
}
|
|
return 0;
|
|
}
|
|
}
|
|
|
|
|
|
int
|
|
ncvargets(
|
|
int ncid,
|
|
int varid,
|
|
const long* start,
|
|
const long* count,
|
|
const long* stride,
|
|
void* value
|
|
)
|
|
{
|
|
if(stride == NULL)
|
|
return ncvarget(ncid, varid, start, count, value);
|
|
/* else */
|
|
{
|
|
NDIMS_DECL
|
|
A_DECL(stp, size_t, ndims, start);
|
|
A_DECL(cntp, size_t, ndims, count);
|
|
A_DECL(strdp, ptrdiff_t, ndims, stride);
|
|
A_INIT(stp, size_t, ndims, start);
|
|
A_INIT(cntp, size_t, ndims, count);
|
|
A_INIT(strdp, ptrdiff_t, ndims, stride);
|
|
{
|
|
const int status = nc_get_vars(ncid, varid, stp, cntp, strdp, value);
|
|
A_FREE(strdp);
|
|
A_FREE(cntp);
|
|
A_FREE(stp);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncvargets", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
}
|
|
return 0;
|
|
}
|
|
}
|
|
|
|
|
|
int
|
|
ncvarputg(
|
|
int ncid,
|
|
int varid,
|
|
const long* start,
|
|
const long* count,
|
|
const long* stride,
|
|
const long* map,
|
|
const void* value
|
|
)
|
|
{
|
|
int ndims = 0;
|
|
if(map == NULL)
|
|
return ncvarputs(ncid, varid, start, count, stride, value);
|
|
/* else */
|
|
{
|
|
ptrdiff_t *imp=NULL;
|
|
if (map != NULL) {
|
|
int ret = NC_NOERR;
|
|
/* make map[ndims-1] number of elements instead of bytes */
|
|
int i, el_size;
|
|
nc_type type;
|
|
ret = nc_inq_varndims(ncid, varid, &ndims);
|
|
if(ret) return ret;
|
|
ret = nc_inq_vartype(ncid, varid, &type);
|
|
if(ret) return ret;
|
|
el_size = nctypelen(type);
|
|
imp = (ptrdiff_t*) malloc(ndims * sizeof(ptrdiff_t));
|
|
for (i=0; i<ndims; i++) imp[i] = map[i] / el_size;
|
|
}
|
|
|
|
{
|
|
A_DECL(stp, size_t, ndims, start);
|
|
A_DECL(cntp, size_t, ndims, count);
|
|
A_DECL(strdp, ptrdiff_t, ndims, stride);
|
|
A_INIT(stp, size_t, ndims, start);
|
|
A_INIT(cntp, size_t, ndims, count);
|
|
A_INIT(strdp, ptrdiff_t, ndims, stride);
|
|
{
|
|
const int status = nc_put_varm(ncid, varid,
|
|
stp, cntp, strdp, imp, value);
|
|
if (imp!=NULL) free(imp);
|
|
A_FREE(strdp);
|
|
A_FREE(cntp);
|
|
A_FREE(stp);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncvarputg", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
}
|
|
return 0;
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
int
|
|
ncvargetg(
|
|
int ncid,
|
|
int varid,
|
|
const long* start,
|
|
const long* count,
|
|
const long* stride,
|
|
const long* map,
|
|
void* value
|
|
)
|
|
{
|
|
int ndims = 0;
|
|
if(map == NULL)
|
|
return ncvargets(ncid, varid, start, count, stride, value);
|
|
/* else */
|
|
{
|
|
ptrdiff_t *imp=NULL;
|
|
if (map != NULL) {
|
|
int ret = NC_NOERR;
|
|
/* make map[ndims-1] number of elements instead of bytes */
|
|
int i, el_size;
|
|
nc_type type;
|
|
ret = nc_inq_varndims(ncid, varid, &ndims);
|
|
if(ret) return ret;
|
|
ret = nc_inq_vartype(ncid, varid, &type);
|
|
if(ret) return ret;
|
|
el_size = nctypelen(type);
|
|
imp = (ptrdiff_t*) malloc(ndims * sizeof(ptrdiff_t));
|
|
for (i=0; i<ndims; i++) imp[i] = map[i] / el_size;
|
|
}
|
|
|
|
{
|
|
A_DECL(stp, size_t, ndims, start);
|
|
A_DECL(cntp, size_t, ndims, count);
|
|
A_DECL(strdp, ptrdiff_t, ndims, stride);
|
|
A_INIT(stp, size_t, ndims, start);
|
|
A_INIT(cntp, size_t, ndims, count);
|
|
A_INIT(strdp, ptrdiff_t, ndims, stride);
|
|
{
|
|
const int status = nc_get_varm(ncid, varid,
|
|
stp, cntp, strdp, imp, value);
|
|
if (imp!=NULL) free(imp);
|
|
A_FREE(strdp);
|
|
A_FREE(cntp);
|
|
A_FREE(stp);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncvargetg", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
}
|
|
return 0;
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
int
|
|
ncvarrename(
|
|
int ncid,
|
|
int varid,
|
|
const char* name
|
|
)
|
|
{
|
|
const int status = nc_rename_var(ncid, varid, name);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncvarrename", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return varid;
|
|
}
|
|
|
|
|
|
int
|
|
ncattput(
|
|
int ncid,
|
|
int varid,
|
|
const char* name,
|
|
nc_type datatype,
|
|
int len,
|
|
const void* value
|
|
)
|
|
{
|
|
const int status = nc_put_att(ncid, varid, name, datatype, len, value);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncattput", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
int
|
|
ncattinq(
|
|
int ncid,
|
|
int varid,
|
|
const char* name,
|
|
nc_type* datatype,
|
|
int* len
|
|
)
|
|
{
|
|
size_t ll;
|
|
const int status = nc_inq_att(ncid, varid, name, datatype, &ll);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncattinq", status,
|
|
"ncid %d; varid %d; attname \"%s\"",
|
|
ncid, varid, name);
|
|
return -1;
|
|
}
|
|
|
|
if(len != NULL)
|
|
*len = (int) ll;
|
|
|
|
return 1;
|
|
|
|
}
|
|
|
|
|
|
int
|
|
ncattget(
|
|
int ncid,
|
|
int varid,
|
|
const char* name,
|
|
void* value
|
|
)
|
|
{
|
|
const int status = nc_get_att(ncid, varid, name, value);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncattget", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return 1;
|
|
}
|
|
|
|
|
|
int
|
|
ncattcopy(
|
|
int ncid_in,
|
|
int varid_in,
|
|
const char* name,
|
|
int ncid_out,
|
|
int varid_out
|
|
)
|
|
{
|
|
const int status = nc_copy_att(ncid_in, varid_in, name, ncid_out, varid_out);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncattcopy", status, "%s", name);
|
|
return -1;
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
int
|
|
ncattname(
|
|
int ncid,
|
|
int varid,
|
|
int attnum,
|
|
char* name
|
|
)
|
|
{
|
|
const int status = nc_inq_attname(ncid, varid, attnum, name);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncattname", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return attnum;
|
|
}
|
|
|
|
|
|
int
|
|
ncattrename(
|
|
int ncid,
|
|
int varid,
|
|
const char* name,
|
|
const char* newname
|
|
)
|
|
{
|
|
const int status = nc_rename_att(ncid, varid, name, newname);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncattrename", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return 1;
|
|
}
|
|
|
|
|
|
int
|
|
ncattdel(
|
|
int ncid,
|
|
int varid,
|
|
const char* name
|
|
)
|
|
{
|
|
const int status = nc_del_att(ncid, varid, name);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncattdel", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return 1;
|
|
}
|
|
|
|
#endif /* NO_NETCDF_2 */
|
|
|
|
#ifndef NO_NETCDF_2
|
|
|
|
int
|
|
ncsetfill(
|
|
int ncid,
|
|
int fillmode
|
|
)
|
|
{
|
|
int oldmode = -1;
|
|
const int status = nc_set_fill(ncid, fillmode, &oldmode);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncsetfill", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return oldmode;
|
|
}
|
|
|
|
|
|
int
|
|
ncrecinq(
|
|
int ncid,
|
|
int* nrecvars,
|
|
int* recvarids,
|
|
long* recsizes
|
|
)
|
|
{
|
|
size_t nrv = 0;
|
|
size_t *rs = NULL;
|
|
int status = NC_NOERR;
|
|
|
|
rs = (size_t*)malloc(sizeof(size_t)*NC_MAX_VARS);
|
|
if(rs == NULL)
|
|
return NC_ENOMEM;
|
|
|
|
status = nc_inq_rec(ncid, &nrv, recvarids, rs);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncrecinq", status, "ncid %d", ncid);
|
|
if(rs != NULL) free(rs);
|
|
return -1;
|
|
}
|
|
|
|
if(nrecvars != NULL)
|
|
*nrecvars = (int) nrv;
|
|
|
|
if(recsizes != NULL)
|
|
{
|
|
size_t ii;
|
|
for(ii = 0; ii < nrv; ii++)
|
|
{
|
|
recsizes[ii] = (long) rs[ii];
|
|
}
|
|
}
|
|
|
|
if(rs != NULL) free(rs);
|
|
|
|
return (int) nrv;
|
|
}
|
|
|
|
|
|
int
|
|
ncrecget(
|
|
int ncid,
|
|
long recnum,
|
|
void** datap
|
|
)
|
|
{
|
|
const int status = nc_get_rec(ncid, (size_t)recnum, datap);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncrecget", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
|
|
int
|
|
ncrecput(
|
|
int ncid,
|
|
long recnum,
|
|
void* const* datap
|
|
)
|
|
{
|
|
const int status = nc_put_rec(ncid, (size_t)recnum, datap);
|
|
if(status != NC_NOERR)
|
|
{
|
|
nc_advise("ncrecput", status, "ncid %d", ncid);
|
|
return -1;
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
#endif /* NO_NETCDF_2 */
|