Index: env-vars
===================================================================
--- env-vars	(.../trunk)	(revision 0)
+++ env-vars	(.../branches/fileinf_improvement)	(revision 531)
@@ -0,0 +1,2 @@
+# to incorporate SXP (NC_FILEINF) support:
+setenv CPPFLAGS "-DNDEBUG -DUSE_SXP"
Index: README-fileinf
===================================================================
--- README-fileinf	(.../trunk)	(revision 0)
+++ README-fileinf	(.../branches/fileinf_improvement)	(revision 531)
@@ -0,0 +1,20 @@
+    NetCDF performance metrics
+    --------------------------
+
+    This version of netCDF uses profiling hooks within netCDF calls
+to monitor I/O performance. To activate these hooks, set the environment 
+variable NC_FILEINF to:
+ 
+    BASIC    minimal information, just file details and a summary of
+             performance for "access" (file open, define, close etc),
+             "info" (file/data enquiry), "put" and "get" operations
+
+    YES or NORMAL  the above plus summaries of buffer use and data
+             type conversion performance
+
+    DETAIL   all of the above but with detailed info about _which_
+             "put" and "get" netcdf calls had what performance
+             (grouped by call, that is, nc_put_vara and nc_put_varm
+             are reported separately)
+
+
Index: libsrc/putget.c
===================================================================
--- libsrc/putget.c	(.../trunk)	(revision 531)
+++ libsrc/putget.c	(.../branches/fileinf_improvement)	(revision 531)
@@ -1,11556 +0,0 @@
-/* Do not edit this file. It is produced from the corresponding .m4 source */
-/*
- *	Copyright 1996, University Corporation for Atmospheric Research
- *      See netcdf/COPYRIGHT file for copying and redistribution conditions.
- */
-/* $Id: putget.m4,v 2.60 2006/01/03 20:20:07 russ Exp $ */
-
-#include "nc.h"
-#include <string.h>
-#include <stdlib.h>
-#include <assert.h>
-#include "ncx.h"
-#include "fbits.h"
-#include "onstack.h"
-#ifdef LOCKNUMREC
-#  include <mpp/shmem.h>	/* for SGI/Cray SHMEM routines */
-#  ifdef LN_TEST
-#    include <stdio.h>
-#  endif
-#endif
-
-#undef MIN  /* system may define MIN somewhere and complain */
-#define MIN(mm,nn) (((mm) < (nn)) ? (mm) : (nn))
-
-/* #define ODEBUG 1 */
-
-#if ODEBUG
-#include <stdio.h>
-/*
- * Print the values of an array of size_t
- */
-void
-arrayp(const char *label, size_t count, const size_t *array)
-{
-	(void) fprintf(stderr, "%s", label);
-	(void) fputc('\t',stderr);	
-	for(; count > 0; count--, array++)
-		(void) fprintf(stderr," %lu", (unsigned long)*array);
-	(void) fputc('\n',stderr);	
-}
-#endif /* ODEBUG */
-
-/*
- *  This is how much space is required by the user, as in
- *
- *   vals = malloc(nel * nctypelen(var.type));
- *   ncvarget(cdfid, varid, cor, edg, vals);
- */
-int
-nctypelen(nc_type type) 
-{
-	switch(type){
-	case NC_BYTE :
-	case NC_CHAR :
-		return((int)sizeof(char));
-	case NC_SHORT :
-		return(int)(sizeof(short));
-	case NC_INT :
-		return((int)sizeof(int));
-	case NC_FLOAT :
-		return((int)sizeof(float));
-	case NC_DOUBLE : 
-		return((int)sizeof(double));
-	}
-
-	return -1;
-}
-
-
-/* Begin fill */
-/*
- * This is tunable parameter.
- * It essentially controls the tradeoff between the number of times
- * memcpy() gets called to copy the external data to fill 
- * a large buffer vs the number of times its called to
- * prepare the external data.
- */
-#if	_SX
-/* NEC SX specific optimization */
-#define	NFILL	2048
-#else
-#define	NFILL	16
-#endif
-
-
-
-/*
- * Next 6 type specific functions
- * Fill a some memory with the default special value.
- * Formerly
-NC_arrayfill()
- */
-static int
-NC_fill_schar(
-	void **xpp,
-	size_t nelems)	/* how many */
-{
-	schar fillp[NFILL * sizeof(double)/X_SIZEOF_CHAR];
-
-	assert(nelems <= sizeof(fillp)/sizeof(fillp[0]));
-
-	{
-		schar *vp = fillp;	/* lower bound of area to be filled */
-		const schar *const end = vp + nelems;
-		while(vp < end)
-		{
-			*vp++ = NC_FILL_BYTE;
-		}
-	}
-	return ncx_putn_schar_schar(xpp, nelems, fillp);
-}
-
-static int
-NC_fill_char(
-	void **xpp,
-	size_t nelems)	/* how many */
-{
-	char fillp[NFILL * sizeof(double)/X_SIZEOF_CHAR];
-
-	assert(nelems <= sizeof(fillp)/sizeof(fillp[0]));
-
-	{
-		char *vp = fillp;	/* lower bound of area to be filled */
-		const char *const end = vp + nelems;
-		while(vp < end)
-		{
-			*vp++ = NC_FILL_CHAR;
-		}
-	}
-	return ncx_putn_char_char(xpp, nelems, fillp);
-}
-
-static int
-NC_fill_short(
-	void **xpp,
-	size_t nelems)	/* how many */
-{
-	short fillp[NFILL * sizeof(double)/X_SIZEOF_SHORT];
-
-	assert(nelems <= sizeof(fillp)/sizeof(fillp[0]));
-
-	{
-		short *vp = fillp;	/* lower bound of area to be filled */
-		const short *const end = vp + nelems;
-		while(vp < end)
-		{
-			*vp++ = NC_FILL_SHORT;
-		}
-	}
-	return ncx_putn_short_short(xpp, nelems, fillp);
-}
-
-
-#if (SIZEOF_INT >= X_SIZEOF_INT)
-static int
-NC_fill_int(
-	void **xpp,
-	size_t nelems)	/* how many */
-{
-	int fillp[NFILL * sizeof(double)/X_SIZEOF_INT];
-
-	assert(nelems <= sizeof(fillp)/sizeof(fillp[0]));
-
-	{
-		int *vp = fillp;	/* lower bound of area to be filled */
-		const int *const end = vp + nelems;
-		while(vp < end)
-		{
-			*vp++ = NC_FILL_INT;
-		}
-	}
-	return ncx_putn_int_int(xpp, nelems, fillp);
-}
-
-#elif SIZEOF_LONG == X_SIZEOF_INT
-static int
-NC_fill_int(
-	void **xpp,
-	size_t nelems)	/* how many */
-{
-	long fillp[NFILL * sizeof(double)/X_SIZEOF_INT];
-
-	assert(nelems <= sizeof(fillp)/sizeof(fillp[0]));
-
-	{
-		long *vp = fillp;	/* lower bound of area to be filled */
-		const long *const end = vp + nelems;
-		while(vp < end)
-		{
-			*vp++ = NC_FILL_INT;
-		}
-	}
-	return ncx_putn_int_long(xpp, nelems, fillp);
-}
-
-#else
-#error "NC_fill_int implementation"
-#endif
-
-static int
-NC_fill_float(
-	void **xpp,
-	size_t nelems)	/* how many */
-{
-	float fillp[NFILL * sizeof(double)/X_SIZEOF_FLOAT];
-
-	assert(nelems <= sizeof(fillp)/sizeof(fillp[0]));
-
-	{
-		float *vp = fillp;	/* lower bound of area to be filled */
-		const float *const end = vp + nelems;
-		while(vp < end)
-		{
-			*vp++ = NC_FILL_FLOAT;
-		}
-	}
-	return ncx_putn_float_float(xpp, nelems, fillp);
-}
-
-static int
-NC_fill_double(
-	void **xpp,
-	size_t nelems)	/* how many */
-{
-	double fillp[NFILL * sizeof(double)/X_SIZEOF_DOUBLE];
-
-	assert(nelems <= sizeof(fillp)/sizeof(fillp[0]));
-
-	{
-		double *vp = fillp;	/* lower bound of area to be filled */
-		const double *const end = vp + nelems;
-		while(vp < end)
-		{
-			*vp++ = NC_FILL_DOUBLE;
-		}
-	}
-	return ncx_putn_double_double(xpp, nelems, fillp);
-}
-
-
-
-/*
- * Fill the external space for variable 'varp' values at 'recno'
- * with the appropriate value. If 'varp' is not a record
- * variable, fill the whole thing.
- * Formerly
-xdr_NC_fill()
- */
-int
-fill_NC_var(NC *ncp, const NC_var *varp, size_t recno)
-{
-	char xfillp[NFILL * X_SIZEOF_DOUBLE];
-	const size_t step = varp->xsz;
-	const size_t nelems = sizeof(xfillp)/step;
-	const size_t xsz = varp->xsz * nelems;
-	NC_attr **attrpp = NULL;
-	off_t offset;
-	size_t remaining = varp->len;
-
-	void *xp;
-	int status = NC_NOERR;
-
-	/*
-	 * Set up fill value
-	 */
-	attrpp = NC_findattr(&varp->attrs, _FillValue);
-	if( attrpp != NULL )
-	{
-		/* User defined fill value */
-		if( (*attrpp)->type != varp->type || (*attrpp)->nelems != 1 )
-		{
-			return NC_EBADTYPE;
-		}
-		else
-		{
-			/* Use the user defined value */
-			char *cp = xfillp;
-			const char *const end = &xfillp[sizeof(xfillp)];
-
-			assert(step <= (*attrpp)->xsz);
-
-			for( /*NADA*/; cp < end; cp += step)
-			{
-				(void) memcpy(cp, (*attrpp)->xvalue, step);
-			}
-		}
-	}
-	else
-	{
-		/* use the default */
-		
-		assert(xsz % X_ALIGN == 0);
-		assert(xsz <= sizeof(xfillp));
-	
-		xp = xfillp;
-	
-		switch(varp->type){
-		case NC_BYTE :
-			status = NC_fill_schar(&xp, nelems);
-			break;
-		case NC_CHAR :
-			status = NC_fill_char(&xp, nelems);
-			break;
-		case NC_SHORT :
-			status = NC_fill_short(&xp, nelems);
-			break;
-		case NC_INT :
-			status = NC_fill_int(&xp, nelems);
-			break;
-		case NC_FLOAT :
-			status = NC_fill_float(&xp, nelems);
-			break;
-		case NC_DOUBLE : 
-			status = NC_fill_double(&xp, nelems);
-			break;
-		default :
-			assert("fill_NC_var invalid type" == 0);
-			status = NC_EBADTYPE;
-			break;
-		}
-		if(status != NC_NOERR)
-			return status;
-	
-		assert(xp == xfillp + xsz);
-	}
-
-	/*
-	 * copyout:
-	 * xfillp now contains 'nelems' elements of the fill value
-	 * in external representation.
-	 */
-
-	/*
-	 * Copy it out.
-	 */
-
-	offset = varp->begin;
-	if(IS_RECVAR(varp))
-	{
-		offset += (off_t)ncp->recsize * recno;
-	}
-
-	assert(remaining > 0);
-	for(;;)
-	{
-		const size_t chunksz = MIN(remaining, ncp->chunk);
-		size_t ii;
-		assert(chunksz % X_ALIGN == 0);
-
-		status = ncp->nciop->get(ncp->nciop, offset, chunksz,
-				 RGN_WRITE, &xp);	
-		if(status != NC_NOERR)
-		{
-			return status;
-		}
-
-		/*
-		 * fill the chunksz buffer in units  of xsz
-		 */
-		for(ii = 0; ii < chunksz/xsz; ii++)
-		{
-			(void) memcpy(xp, xfillp, xsz);
-			xp = (char *)xp + xsz;
-		}
-		/*
-		 * Deal with any remainder
-		 */
-		{
-			const size_t rem = chunksz % xsz;
-			if(rem != 0)
-			{
-				(void) memcpy(xp, xfillp, rem);
-				/* xp = (char *)xp + xsz; */
-			}
-
-		}
-
-		status = ncp->nciop->rel(ncp->nciop, offset, RGN_MODIFIED);
-
-		if(status != NC_NOERR)
-		{
-			break;
-		}
-
-		remaining -= chunksz;
-		if(remaining == 0)
-			break;	/* normal loop exit */
-		offset += chunksz;
-
-	}
-
-	return status;
-}
-/* End fill */
-
-
-/*
- * Add a record containing the fill values.
- */
-static int
-NCfillrecord(NC *ncp, const NC_var *const *varpp, size_t recno)
-{
-	size_t ii = 0;
-	for(; ii < ncp->vars.nelems; ii++, varpp++)
-	{
-		if( !IS_RECVAR(*varpp) )
-		{
-			continue;	/* skip non-record variables */
-		}
-		{
-		const int status = fill_NC_var(ncp, *varpp, recno);
-		if(status != NC_NOERR)
-			return status;
-		}
-	}
-	return NC_NOERR;
-}
-
-/*
- * It is advantageous to
- * #define TOUCH_LAST
- * when using memory mapped io.
- */
-#if TOUCH_LAST
-/*
- * Grow the file to a size which can contain recno
- */
-static int
-NCtouchlast(NC *ncp, const NC_var *const *varpp, size_t recno)
-{
-	int status = NC_NOERR;
-	const NC_var *varp = NULL;
-	
-	{
-	size_t ii = 0;
-	for(; ii < ncp->vars.nelems; ii++, varpp++)
-	{
-		if( !IS_RECVAR(*varpp) )
-		{
-			continue;	/* skip non-record variables */
-		}
-		varp = *varpp;
-	}
-	}
-	assert(varp != NULL);
-	assert( IS_RECVAR(varp) );
-	{
-		const off_t offset = varp->begin
-				+ (off_t)(recno-1) * (off_t)ncp->recsize
-				+ (off_t)(varp->len - varp->xsz);
-		void *xp;
-
-
-		status = ncp->nciop->get(ncp->nciop, offset, varp->xsz,
-				 RGN_WRITE, &xp);	
-		if(status != NC_NOERR)
-			return status;
-		(void)memset(xp, 0, varp->xsz);
-		status = ncp->nciop->rel(ncp->nciop, offset, RGN_MODIFIED);
-	}
-	return status;
-}
-#endif /* TOUCH_LAST */
-
-
-/*
- * Ensure that the netcdf file has 'numrecs' records,
- * add records and fill as neccessary.
- */
-static int
-NCvnrecs(NC *ncp, size_t numrecs)
-{
-	int status = NC_NOERR;
-#ifdef LOCKNUMREC
-	ushmem_t myticket = 0, nowserving = 0;
-	ushmem_t numpe = (ushmem_t) _num_pes();
-
-	/* get ticket and wait */
-	myticket = shmem_short_finc((shmem_t *) ncp->lock + LOCKNUMREC_LOCK,
-		ncp->lock[LOCKNUMREC_BASEPE]);
-#ifdef LN_TEST
-		fprintf(stderr,"%d of %d : ticket = %hu\n",
-			_my_pe(), _num_pes(), myticket);
-#endif
-	do {
-		shmem_short_get((shmem_t *) &nowserving,
-			(shmem_t *) ncp->lock + LOCKNUMREC_SERVING, 1,
-			ncp->lock[LOCKNUMREC_BASEPE]);
-#ifdef LN_TEST
-		fprintf(stderr,"%d of %d : serving = %hu\n",
-			_my_pe(), _num_pes(), nowserving);
-#endif
-		/* work-around for non-unique tickets */
-		if (nowserving > myticket && nowserving < myticket + numpe ) {
-			/* get a new ticket ... you've been bypassed */ 
-			/* and handle the unlikely wrap-around effect */
-			myticket = shmem_short_finc(
-				(shmem_t *) ncp->lock + LOCKNUMREC_LOCK,
-				ncp->lock[LOCKNUMREC_BASEPE]);
-#ifdef LN_TEST
-				fprintf(stderr,"%d of %d : new ticket = %hu\n",
-					_my_pe(), _num_pes(), myticket);
-#endif
-		}
-	} while(nowserving != myticket);
-	/* now our turn to check & update value */
-#endif
-
-	if(numrecs > NC_get_numrecs(ncp))
-	{
-
-
-#if TOUCH_LAST
-		status = NCtouchlast(ncp,
-			(const NC_var *const*)ncp->vars.value,
-			numrecs);
-		if(status != NC_NOERR)
-			goto common_return;
-#endif /* TOUCH_LAST */
-
-		set_NC_ndirty(ncp);
-
-		if(!NC_dofill(ncp))
-		{
-			/* Simply set the new numrecs value */
-			NC_set_numrecs(ncp, numrecs);
-		}
-		else
-		{
-			/* Fill each record out to numrecs */
-			size_t cur_nrecs;
-			while((cur_nrecs = NC_get_numrecs(ncp)) < numrecs)
-			{
-				status = NCfillrecord(ncp,
-					(const NC_var *const*)ncp->vars.value,
-					cur_nrecs);
-				if(status != NC_NOERR)
-				{
-					break;
-				}
-				NC_increase_numrecs(ncp, cur_nrecs +1);
-			}
-			if(status != NC_NOERR)
-				goto common_return;
-		}
-
-		if(NC_doNsync(ncp))
-		{
-			status = write_numrecs(ncp);
-		}
-
-	}
-common_return:
-#ifdef LOCKNUMREC
-	/* finished with our lock - increment serving number */
-	(void) shmem_short_finc((shmem_t *) ncp->lock + LOCKNUMREC_SERVING,
-		ncp->lock[LOCKNUMREC_BASEPE]);
-#endif
-	return status;
-}
-
-
-/* 
- * Check whether 'coord' values are valid for the variable.
- */
-static int
-NCcoordck(NC *ncp, const NC_var *varp, const size_t *coord)
-{
-	const size_t *ip;
-	size_t *up;
-
-	if(varp->ndims == 0)
-		return NC_NOERR;	/* 'scalar' variable */
-
-	if(IS_RECVAR(varp))
-	{
-		if(*coord > X_INT_MAX)
-			return NC_EINVALCOORDS; /* sanity check */
-		if(NC_readonly(ncp) && *coord >= NC_get_numrecs(ncp))
-		{
-			if(!NC_doNsync(ncp))
-				return NC_EINVALCOORDS;
-			/* else */
-			{
-				/* Update from disk and check again */
-				const int status = read_numrecs(ncp);
-				if(status != NC_NOERR)
-					return status;
-				if(*coord >= NC_get_numrecs(ncp))
-					return NC_EINVALCOORDS;
-			}
-		}
-		ip = coord + 1;
-		up = varp->shape + 1;
-	}
-	else
-	{
-		ip = coord;
-		up = varp->shape;
-	}
-	
-#ifdef CDEBUG
-fprintf(stderr,"	NCcoordck: coord %ld, count %d, ip %ld\n",
-		coord, varp->ndims, ip );
-#endif /* CDEBUG */
-
-	for(; ip < coord + varp->ndims; ip++, up++)
-	{
-
-#ifdef CDEBUG
-fprintf(stderr,"	NCcoordck: ip %p, *ip %ld, up %p, *up %lu\n",
-			ip, *ip, up, *up );
-#endif /* CDEBUG */
-
-		/* cast needed for braindead systems with signed size_t */
-		if((unsigned long) *ip >= (unsigned long) *up )
-			return NC_EINVALCOORDS;
-	}
-
-	return NC_NOERR;
-}
-
-
-/* 
- * Check whether 'edges' are valid for the variable and 'start'
- */
-/*ARGSUSED*/
-static int
-NCedgeck(const NC *ncp, const NC_var *varp,
-	 const size_t *start, const size_t *edges)
-{
-	const size_t *const end = start + varp->ndims;
-	const size_t *shp = varp->shape;
-
-	if(varp->ndims == 0)
-		return NC_NOERR;	/* 'scalar' variable */
-
-	if(IS_RECVAR(varp))
-	{
-		start++;
-		edges++;
-		shp++;
-	}
-
-	for(; start < end; start++, edges++, shp++)
-	{
-		/* cast needed for braindead systems with signed size_t */
-		if((unsigned long) *edges > *shp ||
-			(unsigned long) *start + (unsigned long) *edges > *shp)
-		{
-			return(NC_EEDGE);
-		}
-	}
-	return NC_NOERR;
-}
-
-
-/* 
- * Translate the (variable, coord) pair into a seek index
- */
-static off_t
-NC_varoffset(const NC *ncp, const NC_var *varp, const size_t *coord)
-{
-	if(varp->ndims == 0) /* 'scalar' variable */
-		return varp->begin;
-
-	if(varp->ndims == 1)
-	{
-		if(IS_RECVAR(varp))
-			return varp->begin +
-				 (off_t)(*coord) * (off_t)ncp->recsize;
-		/* else */
-		return varp->begin + (off_t)(*coord) * (off_t)varp->xsz;
-	}
-	/* else */
-	{
-		off_t lcoord = (off_t)coord[varp->ndims -1];
-
-		size_t *up = varp->dsizes +1;
-		const size_t *ip = coord;
-		const size_t *const end = varp->dsizes + varp->ndims;
-		
-		if(IS_RECVAR(varp))
-			up++, ip++;
-
-		for(; up < end; up++, ip++)
-			lcoord += *up * *ip;
-
-		lcoord *= varp->xsz;
-		
-		if(IS_RECVAR(varp))
-			lcoord += (off_t)(*coord) * ncp->recsize;
-		
-		lcoord += varp->begin;
-		return lcoord;
-	}
-}
-
-
-
-static int
-putNCvx_char_char(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const char *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_char_char(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-
-static int
-putNCvx_schar_schar(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const schar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_schar_schar(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_schar_uchar(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const uchar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_schar_uchar(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_schar_short(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const short *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_schar_short(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_schar_int(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const int *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_schar_int(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_schar_long(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const long *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_schar_long(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_schar_float(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const float *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_schar_float(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_schar_double(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const double *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_schar_double(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-
-static int
-putNCvx_short_schar(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const schar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_short_schar(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_short_uchar(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const uchar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_short_uchar(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_short_short(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const short *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_short_short(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_short_int(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const int *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_short_int(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_short_long(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const long *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_short_long(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_short_float(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const float *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_short_float(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_short_double(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const double *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_short_double(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-
-static int
-putNCvx_int_schar(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const schar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_int_schar(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_int_uchar(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const uchar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_int_uchar(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_int_short(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const short *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_int_short(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_int_int(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const int *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_int_int(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_int_long(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const long *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_int_long(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_int_float(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const float *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_int_float(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_int_double(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const double *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_int_double(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-
-static int
-putNCvx_float_schar(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const schar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_float_schar(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_float_uchar(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const uchar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_float_uchar(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_float_short(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const short *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_float_short(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_float_int(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const int *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_float_int(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_float_long(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const long *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_float_long(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_float_float(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const float *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_float_float(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_float_double(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const double *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_float_double(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-
-static int
-putNCvx_double_schar(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const schar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_double_schar(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_double_uchar(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const uchar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_double_uchar(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_double_short(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const short *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_double_short(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_double_int(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const int *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_double_int(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_double_long(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const long *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_double_long(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_double_float(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const float *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_double_float(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-static int
-putNCvx_double_double(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const double *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nput = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 RGN_WRITE, &xp);	
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_putn_double_double(&xp, nput, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-		{
-			/* not fatal to the loop */
-			status = lstatus;
-		}
-
-		(void) ncp->nciop->rel(ncp->nciop, offset,
-				 RGN_MODIFIED);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nput;
-
-	}
-
-	return status;
-}
-
-
-
-
-static int
-putNCv_text(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const char *value)
-{
-	if(varp->type != NC_CHAR)
-		return NC_ECHAR;
-	return putNCvx_char_char(ncp, varp, start, nelems, value);
-}
-
-static int
-putNCv_schar(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const schar *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return putNCvx_schar_schar(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return putNCvx_short_schar(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return putNCvx_int_schar(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return putNCvx_float_schar(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return putNCvx_double_schar(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-static int
-putNCv_uchar(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const uchar *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return putNCvx_schar_uchar(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return putNCvx_short_uchar(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return putNCvx_int_uchar(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return putNCvx_float_uchar(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return putNCvx_double_uchar(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-static int
-putNCv_short(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const short *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return putNCvx_schar_short(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return putNCvx_short_short(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return putNCvx_int_short(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return putNCvx_float_short(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return putNCvx_double_short(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-static int
-putNCv_int(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const int *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return putNCvx_schar_int(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return putNCvx_short_int(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return putNCvx_int_int(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return putNCvx_float_int(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return putNCvx_double_int(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-static int
-putNCv_long(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const long *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return putNCvx_schar_long(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return putNCvx_short_long(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return putNCvx_int_long(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return putNCvx_float_long(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return putNCvx_double_long(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-static int
-putNCv_float(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const float *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return putNCvx_schar_float(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return putNCvx_short_float(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return putNCvx_int_float(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return putNCvx_float_float(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return putNCvx_double_float(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-static int
-putNCv_double(NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, const double *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return putNCvx_schar_double(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return putNCvx_short_double(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return putNCvx_int_double(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return putNCvx_float_double(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return putNCvx_double_double(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-
-
-
-static int
-getNCvx_char_char(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, char *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_char_char(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-
-static int
-getNCvx_schar_schar(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, schar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_schar_schar(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_schar_uchar(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, uchar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_schar_uchar(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_schar_short(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, short *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_schar_short(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_schar_int(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, int *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_schar_int(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_schar_long(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, long *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_schar_long(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_schar_float(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, float *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_schar_float(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_schar_double(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, double *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_schar_double(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-
-static int
-getNCvx_short_schar(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, schar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_short_schar(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_short_uchar(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, uchar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_short_uchar(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_short_short(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, short *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_short_short(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_short_int(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, int *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_short_int(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_short_long(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, long *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_short_long(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_short_float(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, float *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_short_float(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_short_double(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, double *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_short_double(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-
-static int
-getNCvx_int_schar(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, schar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_int_schar(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_int_uchar(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, uchar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_int_uchar(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_int_short(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, short *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_int_short(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_int_int(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, int *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_int_int(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_int_long(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, long *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_int_long(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_int_float(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, float *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_int_float(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_int_double(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, double *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_int_double(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-
-static int
-getNCvx_float_schar(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, schar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_float_schar(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_float_uchar(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, uchar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_float_uchar(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_float_short(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, short *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_float_short(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_float_int(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, int *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_float_int(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_float_long(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, long *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_float_long(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_float_float(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, float *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_float_float(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_float_double(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, double *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_float_double(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-
-static int
-getNCvx_double_schar(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, schar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_double_schar(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_double_uchar(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, uchar *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_double_uchar(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_double_short(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, short *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_double_short(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_double_int(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, int *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_double_int(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_double_long(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, long *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_double_long(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_double_float(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, float *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_double_float(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-static int
-getNCvx_double_double(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, double *value)
-{
-	off_t offset = NC_varoffset(ncp, varp, start);
-	size_t remaining = varp->xsz * nelems;
-	int status = NC_NOERR;
-	const void *xp;
-
-	if(nelems == 0)
-		return NC_NOERR;
-
-	assert(value != NULL);
-
-	for(;;)
-	{
-		size_t extent = MIN(remaining, ncp->chunk);
-		size_t nget = ncx_howmany(varp->type, extent);
-
-		int lstatus = ncp->nciop->get(ncp->nciop, offset, extent,
-				 0, (void **)&xp);	/* cast away const */
-		if(lstatus != NC_NOERR)
-			return lstatus;
-		
-		lstatus = ncx_getn_double_double(&xp, nget, value);
-		if(lstatus != NC_NOERR && status == NC_NOERR)
-			status = lstatus;
-
-		(void) ncp->nciop->rel(ncp->nciop, offset, 0);	
-
-		remaining -= extent;
-		if(remaining == 0)
-			break; /* normal loop exit */
-		offset += extent;
-		value += nget;
-	}
-
-	return status;
-}
-
-
-
-
-static int
-getNCv_schar(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, schar *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return getNCvx_schar_schar(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return getNCvx_short_schar(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return getNCvx_int_schar(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return getNCvx_float_schar(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return getNCvx_double_schar(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-static int
-getNCv_uchar(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, uchar *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return getNCvx_schar_uchar(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return getNCvx_short_uchar(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return getNCvx_int_uchar(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return getNCvx_float_uchar(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return getNCvx_double_uchar(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-static int
-getNCv_short(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, short *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return getNCvx_schar_short(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return getNCvx_short_short(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return getNCvx_int_short(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return getNCvx_float_short(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return getNCvx_double_short(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-static int
-getNCv_int(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, int *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return getNCvx_schar_int(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return getNCvx_short_int(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return getNCvx_int_int(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return getNCvx_float_int(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return getNCvx_double_int(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-static int
-getNCv_long(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, long *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return getNCvx_schar_long(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return getNCvx_short_long(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return getNCvx_int_long(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return getNCvx_float_long(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return getNCvx_double_long(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-static int
-getNCv_float(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, float *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return getNCvx_schar_float(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return getNCvx_short_float(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return getNCvx_int_float(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return getNCvx_float_float(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return getNCvx_double_float(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-static int
-getNCv_double(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, double *value)
-{
-	switch(varp->type){
-	case NC_CHAR:
-		return NC_ECHAR;
-	case NC_BYTE:
-		return getNCvx_schar_double(ncp, varp, start, nelems,
-			value);
-	case NC_SHORT:
-		return getNCvx_short_double(ncp, varp, start, nelems,
-			value);
-	case NC_INT:
-		return getNCvx_int_double(ncp, varp, start, nelems,
-			value);
-	case NC_FLOAT:
-		return getNCvx_float_double(ncp, varp, start, nelems,
-			value);
-	case NC_DOUBLE: 
-		return getNCvx_double_double(ncp, varp, start, nelems,
-			value);
-	}
-	return NC_EBADTYPE;
-}
-
-
-
-static int
-getNCv_text(const NC *ncp, const NC_var *varp,
-		 const size_t *start, size_t nelems, char *value)
-{
-	if(varp->type != NC_CHAR)
-		return NC_ECHAR;
-	return getNCvx_char_char(ncp, varp, start, nelems, value);
-}
-
-
-/*
- * Copy 'nbytes' contiguous external values
- * from ('inncp', invp', inncoord')
- * to   ('outncp', 'outvp', 'outcoord')
- * 'inncp' shouldn't be the same as 'outncp'.
- * Used only by ncvarcopy()
- */
-static int
-NCxvarcpy(NC *inncp, NC_var *invp, size_t *incoord,
-	NC *outncp, NC_var *outvp, size_t *outcoord, size_t nbytes)
-{
-	int status;
-	off_t inoffset = NC_varoffset(inncp, invp, incoord);
-	off_t outoffset = NC_varoffset(outncp, outvp, outcoord);
-	void *inxp;
-	void *outxp;
-	const size_t chunk = MIN(inncp->chunk, outncp->chunk);
-
-	do {
-		const size_t extent = MIN(nbytes, chunk);
-
-		status = inncp->nciop->get(inncp->nciop, inoffset, extent,
-				 0, &inxp);	
-		if(status != NC_NOERR)
-			return status;
-
-		status = outncp->nciop->get(outncp->nciop, outoffset, extent,
-				 RGN_WRITE, &outxp);	
-		if(status != NC_NOERR)
-		{
-			(void) inncp->nciop->rel(inncp->nciop, inoffset, 0);	
-			break;
-		}
-
-		(void) memcpy(outxp, inxp, extent);
-
-		status = outncp->nciop->rel(outncp->nciop, outoffset,
-			 RGN_MODIFIED);
-		(void) inncp->nciop->rel(inncp->nciop, inoffset, 0);	
-
-		nbytes -= extent;
-		if(nbytes == 0)
-			break; /* normal loop exit */
-		inoffset += extent;
-		outoffset += extent;
-		
-	} while (status == NC_NOERR);
-
-	return status;
-}
-
-
-/*
- *  For ncvar{put,get},
- *  find the largest contiguous block from within 'edges'.
- *  returns the index to the left of this (which may be -1).
- *  Compute the number of contiguous elements and return
- *  that in *iocountp.
- *  The presence of "record" variables makes this routine
- *  overly subtle.
- */
-static int
-NCiocount(const NC *const ncp, const NC_var *const varp,
-	const size_t *const edges,
-	size_t *const iocountp)
-{
-	const size_t *edp0 = edges;
-	const size_t *edp = edges + varp->ndims;
-	const size_t *shp = varp->shape + varp->ndims;
-
-	if(IS_RECVAR(varp))
-	{
-		if(varp->ndims == 1 && ncp->recsize <= varp->len)
-		{
-			/* one dimensional && the only 'record' variable */
-			*iocountp = *edges;
-			return(0);
-		}
-		/* else */
-		edp0++;
-	}
-
-	assert(edges != NULL);
-
-	/* find max contiguous */
-	while(edp > edp0)
-	{
-		shp--; edp--;
-		if(*edp < *shp )
-		{
-			const size_t *zedp = edp;
-			while(zedp >= edp0)
-			{
-				if(*zedp == 0)
-				{
-					*iocountp = 0;
-					goto done;
-				}
-				/* Tip of the hat to segmented architectures */
-				if(zedp == edp0)
-					break;
-				zedp--;
-			}
-			break;
-		}
-		assert(*edp == *shp);
-	}
-
-	/*
-	 * edp, shp reference rightmost index s.t. *(edp +1) == *(shp +1)
-	 *
-	 * Or there is only one dimension.
-	 * If there is only one dimension and it is 'non record' dimension,
-	 * 	edp is &edges[0] and we will return -1.
-	 * If there is only one dimension and and it is a "record dimension",
-	 *	edp is &edges[1] (out of bounds) and we will return 0;
-	 */
-	assert(shp >= varp->shape + varp->ndims -1 
-		|| *(edp +1) == *(shp +1));
-
-	/* now accumulate max count for a single io operation */
-	for(*iocountp = 1, edp0 = edp;
-		 	edp0 < edges + varp->ndims;
-			edp0++)
-	{
-		*iocountp *= *edp0;
-	}
-
-done:
-	return((int)(edp - edges) - 1);
-}
-
-
-/*
- * Set the elements of the array 'upp' to
- * the sum of the corresponding elements of
- * 'stp' and 'edp'. 'end' should be &stp[nelems].
- */
-static void
-set_upper(size_t *upp, /* modified on return */
-	const size_t *stp,
-	const size_t *edp,
-	const size_t *const end)
-{
-	while(upp < end) {
-		*upp++ = *stp++ + *edp++;
-	}
-}
-
-
-/*
- * The infamous and oft-discussed odometer code.
- *
- * 'start[]' is the starting coordinate.
- * 'upper[]' is the upper bound s.t. start[ii] < upper[ii].
- * 'coord[]' is the register, the current coordinate value.
- * For some ii,
- * upp == &upper[ii]
- * cdp == &coord[ii]
- * 
- * Running this routine increments *cdp.
- *
- * If after the increment, *cdp is equal to *upp
- * (and cdp is not the leftmost dimension),
- * *cdp is "zeroed" to the starting value and
- * we need to "carry", eg, increment one place to
- * the left.
- * 
- * TODO: Some architectures hate recursion?
- * 	Reimplement non-recursively.
- */
-static void
-odo1(const size_t *const start, const size_t *const upper,
-	size_t *const coord, /* modified on return */
-	const size_t *upp,
-	size_t *cdp)
-{
-	assert(coord <= cdp && cdp <= coord + NC_MAX_VAR_DIMS);
-	assert(upper <= upp && upp <= upper + NC_MAX_VAR_DIMS);
-	assert(upp - upper == cdp - coord);
-	
-	assert(*cdp <= *upp);
-
-	(*cdp)++;
-	if(cdp != coord && *cdp >= *upp)
-	{
-		*cdp = start[cdp - coord];
-		odo1(start, upper, coord, upp -1, cdp -1);
-	}
-}
-#ifdef _CRAYC
-#pragma _CRI noinline odo1
-#endif
-
-
-
-
-/* Public */
-
-
-int
-nc_put_var1_text(int ncid, int varid, const size_t *coord,
-	const char *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type != NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *coord +1);
-		if(status != NC_NOERR)
-			return status;
-	}
-
-	return putNCv_text(ncp, varp, coord, 1, value);
-}
-
-
-int
-nc_put_var1_uchar(int ncid, int varid, const size_t *coord,
-	const uchar *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *coord +1);
-		if(status != NC_NOERR)
-			return status;
-	}
-
-	return putNCv_uchar(ncp, varp, coord, 1, value);
-}
-
-int
-nc_put_var1_schar(int ncid, int varid, const size_t *coord,
-	const schar *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *coord +1);
-		if(status != NC_NOERR)
-			return status;
-	}
-
-	return putNCv_schar(ncp, varp, coord, 1, value);
-}
-
-int
-nc_put_var1_short(int ncid, int varid, const size_t *coord,
-	const short *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *coord +1);
-		if(status != NC_NOERR)
-			return status;
-	}
-
-	return putNCv_short(ncp, varp, coord, 1, value);
-}
-
-int
-nc_put_var1_int(int ncid, int varid, const size_t *coord,
-	const int *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *coord +1);
-		if(status != NC_NOERR)
-			return status;
-	}
-
-	return putNCv_int(ncp, varp, coord, 1, value);
-}
-
-int
-nc_put_var1_long(int ncid, int varid, const size_t *coord,
-	const long *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *coord +1);
-		if(status != NC_NOERR)
-			return status;
-	}
-
-	return putNCv_long(ncp, varp, coord, 1, value);
-}
-
-int
-nc_put_var1_float(int ncid, int varid, const size_t *coord,
-	const float *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *coord +1);
-		if(status != NC_NOERR)
-			return status;
-	}
-
-	return putNCv_float(ncp, varp, coord, 1, value);
-}
-
-int
-nc_put_var1_double(int ncid, int varid, const size_t *coord,
-	const double *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *coord +1);
-		if(status != NC_NOERR)
-			return status;
-	}
-
-	return putNCv_double(ncp, varp, coord, 1, value);
-}
-
-
-
-int
-nc_get_var1_text(int ncid, int varid, const size_t *coord, char *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type != NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	return getNCv_text(ncp, varp, coord, 1, value);
-}
-
-
-int
-nc_get_var1_uchar(int ncid, int varid, const size_t *coord, uchar *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	return getNCv_uchar(ncp, varp, coord, 1, value);
-}
-
-int
-nc_get_var1_schar(int ncid, int varid, const size_t *coord, schar *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	return getNCv_schar(ncp, varp, coord, 1, value);
-}
-
-int
-nc_get_var1_short(int ncid, int varid, const size_t *coord, short *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	return getNCv_short(ncp, varp, coord, 1, value);
-}
-
-int
-nc_get_var1_int(int ncid, int varid, const size_t *coord, int *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	return getNCv_int(ncp, varp, coord, 1, value);
-}
-
-int
-nc_get_var1_long(int ncid, int varid, const size_t *coord, long *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	return getNCv_long(ncp, varp, coord, 1, value);
-}
-
-int
-nc_get_var1_float(int ncid, int varid, const size_t *coord, float *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	return getNCv_float(ncp, varp, coord, 1, value);
-}
-
-int
-nc_get_var1_double(int ncid, int varid, const size_t *coord, double *value)
-{
-	int status;
-	NC *ncp;
-	const NC_var *varp;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, coord);
-	if(status != NC_NOERR)
-		return status;
-
-	return getNCv_double(ncp, varp, coord, 1, value);
-}
-
-
-
-int
-nc_put_vara_text(int ncid, int varid,
-	 const size_t *start, const size_t *edges, const char *value)
-{
-	int status = NC_NOERR;
-	NC *ncp;
-	const NC_var *varp;
-	int ii;
-	size_t iocount;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type != NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, start);
-	if(status != NC_NOERR)
-		return status;
-	status = NCedgeck(ncp, varp, start, edges);
-	if(status != NC_NOERR)
-		return status;
-
-	if(varp->ndims == 0) /* scalar variable */
-	{
-		return( putNCv_text(ncp, varp, start, 1, value) );
-	}
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *start + *edges);
-		if(status != NC_NOERR)
-			return status;
-
-		if(varp->ndims == 1
-			&& ncp->recsize <= varp->len)
-		{
-			/* one dimensional && the only record variable  */
-			return( putNCv_text(ncp, varp, start, *edges, value) );
-		}
-	}
-
-	/*
-	 * find max contiguous
-	 *   and accumulate max count for a single io operation
-	 */
-	ii = NCiocount(ncp, varp, edges, &iocount);
-
-	if(ii == -1)
-	{
-		return( putNCv_text(ncp, varp, start, iocount, value) );
-	}
-
-	assert(ii >= 0);
-
-
-	{ /* inline */
-	ALLOC_ONSTACK(coord, size_t, varp->ndims);
-	ALLOC_ONSTACK(upper, size_t, varp->ndims);
-	const size_t index = ii;
-
-	/* copy in starting indices */
-	(void) memcpy(coord, start, varp->ndims * sizeof(size_t));
-
-	/* set up in maximum indices */
-	set_upper(upper, start, edges, &upper[varp->ndims]);
-
-	/* ripple counter */
-	while(*coord < *upper)
-	{
-		const int lstatus = putNCv_text(ncp, varp, coord, iocount,
-				 value);
-		if(lstatus != NC_NOERR)
-		{
-			if(lstatus != NC_ERANGE)
-			{
-				status = lstatus;
-				/* fatal for the loop */
-				break;
-			}
-			/* else NC_ERANGE, not fatal for the loop */
-			if(status == NC_NOERR)
-				status = lstatus;
-		}
-		value += iocount;
-		odo1(start, upper, coord, &upper[index], &coord[index]);
-	}
-
-	FREE_ONSTACK(upper);
-	FREE_ONSTACK(coord);
-	} /* end inline */
-
-	return status;
-}
-
-
-int
-nc_put_vara_uchar(int ncid, int varid,
-	 const size_t *start, const size_t *edges, const uchar *value)
-{
-	int status = NC_NOERR;
-	NC *ncp;
-	const NC_var *varp;
-	int ii;
-	size_t iocount;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, start);
-	if(status != NC_NOERR)
-		return status;
-	status = NCedgeck(ncp, varp, start, edges);
-	if(status != NC_NOERR)
-		return status;
-
-	if(varp->ndims == 0) /* scalar variable */
-	{
-		return( putNCv_uchar(ncp, varp, start, 1, value) );
-	}
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *start + *edges);
-		if(status != NC_NOERR)
-			return status;
-
-		if(varp->ndims == 1
-			&& ncp->recsize <= varp->len)
-		{
-			/* one dimensional && the only record variable  */
-			return( putNCv_uchar(ncp, varp, start, *edges, value) );
-		}
-	}
-
-	/*
-	 * find max contiguous
-	 *   and accumulate max count for a single io operation
-	 */
-	ii = NCiocount(ncp, varp, edges, &iocount);
-
-	if(ii == -1)
-	{
-		return( putNCv_uchar(ncp, varp, start, iocount, value) );
-	}
-
-	assert(ii >= 0);
-
-
-	{ /* inline */
-	ALLOC_ONSTACK(coord, size_t, varp->ndims);
-	ALLOC_ONSTACK(upper, size_t, varp->ndims);
-	const size_t index = ii;
-
-	/* copy in starting indices */
-	(void) memcpy(coord, start, varp->ndims * sizeof(size_t));
-
-	/* set up in maximum indices */
-	set_upper(upper, start, edges, &upper[varp->ndims]);
-
-	/* ripple counter */
-	while(*coord < *upper)
-	{
-		const int lstatus = putNCv_uchar(ncp, varp, coord, iocount,
-				 value);
-		if(lstatus != NC_NOERR)
-		{
-			if(lstatus != NC_ERANGE)
-			{
-				status = lstatus;
-				/* fatal for the loop */
-				break;
-			}
-			/* else NC_ERANGE, not fatal for the loop */
-			if(status == NC_NOERR)
-				status = lstatus;
-		}
-		value += iocount;
-		odo1(start, upper, coord, &upper[index], &coord[index]);
-	}
-
-	FREE_ONSTACK(upper);
-	FREE_ONSTACK(coord);
-	} /* end inline */
-
-	return status;
-}
-
-int
-nc_put_vara_schar(int ncid, int varid,
-	 const size_t *start, const size_t *edges, const schar *value)
-{
-	int status = NC_NOERR;
-	NC *ncp;
-	const NC_var *varp;
-	int ii;
-	size_t iocount;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, start);
-	if(status != NC_NOERR)
-		return status;
-	status = NCedgeck(ncp, varp, start, edges);
-	if(status != NC_NOERR)
-		return status;
-
-	if(varp->ndims == 0) /* scalar variable */
-	{
-		return( putNCv_schar(ncp, varp, start, 1, value) );
-	}
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *start + *edges);
-		if(status != NC_NOERR)
-			return status;
-
-		if(varp->ndims == 1
-			&& ncp->recsize <= varp->len)
-		{
-			/* one dimensional && the only record variable  */
-			return( putNCv_schar(ncp, varp, start, *edges, value) );
-		}
-	}
-
-	/*
-	 * find max contiguous
-	 *   and accumulate max count for a single io operation
-	 */
-	ii = NCiocount(ncp, varp, edges, &iocount);
-
-	if(ii == -1)
-	{
-		return( putNCv_schar(ncp, varp, start, iocount, value) );
-	}
-
-	assert(ii >= 0);
-
-
-	{ /* inline */
-	ALLOC_ONSTACK(coord, size_t, varp->ndims);
-	ALLOC_ONSTACK(upper, size_t, varp->ndims);
-	const size_t index = ii;
-
-	/* copy in starting indices */
-	(void) memcpy(coord, start, varp->ndims * sizeof(size_t));
-
-	/* set up in maximum indices */
-	set_upper(upper, start, edges, &upper[varp->ndims]);
-
-	/* ripple counter */
-	while(*coord < *upper)
-	{
-		const int lstatus = putNCv_schar(ncp, varp, coord, iocount,
-				 value);
-		if(lstatus != NC_NOERR)
-		{
-			if(lstatus != NC_ERANGE)
-			{
-				status = lstatus;
-				/* fatal for the loop */
-				break;
-			}
-			/* else NC_ERANGE, not fatal for the loop */
-			if(status == NC_NOERR)
-				status = lstatus;
-		}
-		value += iocount;
-		odo1(start, upper, coord, &upper[index], &coord[index]);
-	}
-
-	FREE_ONSTACK(upper);
-	FREE_ONSTACK(coord);
-	} /* end inline */
-
-	return status;
-}
-
-int
-nc_put_vara_short(int ncid, int varid,
-	 const size_t *start, const size_t *edges, const short *value)
-{
-	int status = NC_NOERR;
-	NC *ncp;
-	const NC_var *varp;
-	int ii;
-	size_t iocount;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, start);
-	if(status != NC_NOERR)
-		return status;
-	status = NCedgeck(ncp, varp, start, edges);
-	if(status != NC_NOERR)
-		return status;
-
-	if(varp->ndims == 0) /* scalar variable */
-	{
-		return( putNCv_short(ncp, varp, start, 1, value) );
-	}
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *start + *edges);
-		if(status != NC_NOERR)
-			return status;
-
-		if(varp->ndims == 1
-			&& ncp->recsize <= varp->len)
-		{
-			/* one dimensional && the only record variable  */
-			return( putNCv_short(ncp, varp, start, *edges, value) );
-		}
-	}
-
-	/*
-	 * find max contiguous
-	 *   and accumulate max count for a single io operation
-	 */
-	ii = NCiocount(ncp, varp, edges, &iocount);
-
-	if(ii == -1)
-	{
-		return( putNCv_short(ncp, varp, start, iocount, value) );
-	}
-
-	assert(ii >= 0);
-
-
-	{ /* inline */
-	ALLOC_ONSTACK(coord, size_t, varp->ndims);
-	ALLOC_ONSTACK(upper, size_t, varp->ndims);
-	const size_t index = ii;
-
-	/* copy in starting indices */
-	(void) memcpy(coord, start, varp->ndims * sizeof(size_t));
-
-	/* set up in maximum indices */
-	set_upper(upper, start, edges, &upper[varp->ndims]);
-
-	/* ripple counter */
-	while(*coord < *upper)
-	{
-		const int lstatus = putNCv_short(ncp, varp, coord, iocount,
-				 value);
-		if(lstatus != NC_NOERR)
-		{
-			if(lstatus != NC_ERANGE)
-			{
-				status = lstatus;
-				/* fatal for the loop */
-				break;
-			}
-			/* else NC_ERANGE, not fatal for the loop */
-			if(status == NC_NOERR)
-				status = lstatus;
-		}
-		value += iocount;
-		odo1(start, upper, coord, &upper[index], &coord[index]);
-	}
-
-	FREE_ONSTACK(upper);
-	FREE_ONSTACK(coord);
-	} /* end inline */
-
-	return status;
-}
-
-int
-nc_put_vara_int(int ncid, int varid,
-	 const size_t *start, const size_t *edges, const int *value)
-{
-	int status = NC_NOERR;
-	NC *ncp;
-	const NC_var *varp;
-	int ii;
-	size_t iocount;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, start);
-	if(status != NC_NOERR)
-		return status;
-	status = NCedgeck(ncp, varp, start, edges);
-	if(status != NC_NOERR)
-		return status;
-
-	if(varp->ndims == 0) /* scalar variable */
-	{
-		return( putNCv_int(ncp, varp, start, 1, value) );
-	}
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *start + *edges);
-		if(status != NC_NOERR)
-			return status;
-
-		if(varp->ndims == 1
-			&& ncp->recsize <= varp->len)
-		{
-			/* one dimensional && the only record variable  */
-			return( putNCv_int(ncp, varp, start, *edges, value) );
-		}
-	}
-
-	/*
-	 * find max contiguous
-	 *   and accumulate max count for a single io operation
-	 */
-	ii = NCiocount(ncp, varp, edges, &iocount);
-
-	if(ii == -1)
-	{
-		return( putNCv_int(ncp, varp, start, iocount, value) );
-	}
-
-	assert(ii >= 0);
-
-
-	{ /* inline */
-	ALLOC_ONSTACK(coord, size_t, varp->ndims);
-	ALLOC_ONSTACK(upper, size_t, varp->ndims);
-	const size_t index = ii;
-
-	/* copy in starting indices */
-	(void) memcpy(coord, start, varp->ndims * sizeof(size_t));
-
-	/* set up in maximum indices */
-	set_upper(upper, start, edges, &upper[varp->ndims]);
-
-	/* ripple counter */
-	while(*coord < *upper)
-	{
-		const int lstatus = putNCv_int(ncp, varp, coord, iocount,
-				 value);
-		if(lstatus != NC_NOERR)
-		{
-			if(lstatus != NC_ERANGE)
-			{
-				status = lstatus;
-				/* fatal for the loop */
-				break;
-			}
-			/* else NC_ERANGE, not fatal for the loop */
-			if(status == NC_NOERR)
-				status = lstatus;
-		}
-		value += iocount;
-		odo1(start, upper, coord, &upper[index], &coord[index]);
-	}
-
-	FREE_ONSTACK(upper);
-	FREE_ONSTACK(coord);
-	} /* end inline */
-
-	return status;
-}
-
-int
-nc_put_vara_long(int ncid, int varid,
-	 const size_t *start, const size_t *edges, const long *value)
-{
-	int status = NC_NOERR;
-	NC *ncp;
-	const NC_var *varp;
-	int ii;
-	size_t iocount;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, start);
-	if(status != NC_NOERR)
-		return status;
-	status = NCedgeck(ncp, varp, start, edges);
-	if(status != NC_NOERR)
-		return status;
-
-	if(varp->ndims == 0) /* scalar variable */
-	{
-		return( putNCv_long(ncp, varp, start, 1, value) );
-	}
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *start + *edges);
-		if(status != NC_NOERR)
-			return status;
-
-		if(varp->ndims == 1
-			&& ncp->recsize <= varp->len)
-		{
-			/* one dimensional && the only record variable  */
-			return( putNCv_long(ncp, varp, start, *edges, value) );
-		}
-	}
-
-	/*
-	 * find max contiguous
-	 *   and accumulate max count for a single io operation
-	 */
-	ii = NCiocount(ncp, varp, edges, &iocount);
-
-	if(ii == -1)
-	{
-		return( putNCv_long(ncp, varp, start, iocount, value) );
-	}
-
-	assert(ii >= 0);
-
-
-	{ /* inline */
-	ALLOC_ONSTACK(coord, size_t, varp->ndims);
-	ALLOC_ONSTACK(upper, size_t, varp->ndims);
-	const size_t index = ii;
-
-	/* copy in starting indices */
-	(void) memcpy(coord, start, varp->ndims * sizeof(size_t));
-
-	/* set up in maximum indices */
-	set_upper(upper, start, edges, &upper[varp->ndims]);
-
-	/* ripple counter */
-	while(*coord < *upper)
-	{
-		const int lstatus = putNCv_long(ncp, varp, coord, iocount,
-				 value);
-		if(lstatus != NC_NOERR)
-		{
-			if(lstatus != NC_ERANGE)
-			{
-				status = lstatus;
-				/* fatal for the loop */
-				break;
-			}
-			/* else NC_ERANGE, not fatal for the loop */
-			if(status == NC_NOERR)
-				status = lstatus;
-		}
-		value += iocount;
-		odo1(start, upper, coord, &upper[index], &coord[index]);
-	}
-
-	FREE_ONSTACK(upper);
-	FREE_ONSTACK(coord);
-	} /* end inline */
-
-	return status;
-}
-
-int
-nc_put_vara_float(int ncid, int varid,
-	 const size_t *start, const size_t *edges, const float *value)
-{
-	int status = NC_NOERR;
-	NC *ncp;
-	const NC_var *varp;
-	int ii;
-	size_t iocount;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, start);
-	if(status != NC_NOERR)
-		return status;
-	status = NCedgeck(ncp, varp, start, edges);
-	if(status != NC_NOERR)
-		return status;
-
-	if(varp->ndims == 0) /* scalar variable */
-	{
-		return( putNCv_float(ncp, varp, start, 1, value) );
-	}
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *start + *edges);
-		if(status != NC_NOERR)
-			return status;
-
-		if(varp->ndims == 1
-			&& ncp->recsize <= varp->len)
-		{
-			/* one dimensional && the only record variable  */
-			return( putNCv_float(ncp, varp, start, *edges, value) );
-		}
-	}
-
-	/*
-	 * find max contiguous
-	 *   and accumulate max count for a single io operation
-	 */
-	ii = NCiocount(ncp, varp, edges, &iocount);
-
-	if(ii == -1)
-	{
-		return( putNCv_float(ncp, varp, start, iocount, value) );
-	}
-
-	assert(ii >= 0);
-
-
-	{ /* inline */
-	ALLOC_ONSTACK(coord, size_t, varp->ndims);
-	ALLOC_ONSTACK(upper, size_t, varp->ndims);
-	const size_t index = ii;
-
-	/* copy in starting indices */
-	(void) memcpy(coord, start, varp->ndims * sizeof(size_t));
-
-	/* set up in maximum indices */
-	set_upper(upper, start, edges, &upper[varp->ndims]);
-
-	/* ripple counter */
-	while(*coord < *upper)
-	{
-		const int lstatus = putNCv_float(ncp, varp, coord, iocount,
-				 value);
-		if(lstatus != NC_NOERR)
-		{
-			if(lstatus != NC_ERANGE)
-			{
-				status = lstatus;
-				/* fatal for the loop */
-				break;
-			}
-			/* else NC_ERANGE, not fatal for the loop */
-			if(status == NC_NOERR)
-				status = lstatus;
-		}
-		value += iocount;
-		odo1(start, upper, coord, &upper[index], &coord[index]);
-	}
-
-	FREE_ONSTACK(upper);
-	FREE_ONSTACK(coord);
-	} /* end inline */
-
-	return status;
-}
-
-int
-nc_put_vara_double(int ncid, int varid,
-	 const size_t *start, const size_t *edges, const double *value)
-{
-	int status = NC_NOERR;
-	NC *ncp;
-	const NC_var *varp;
-	int ii;
-	size_t iocount;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_readonly(ncp))
-		return NC_EPERM;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, start);
-	if(status != NC_NOERR)
-		return status;
-	status = NCedgeck(ncp, varp, start, edges);
-	if(status != NC_NOERR)
-		return status;
-
-	if(varp->ndims == 0) /* scalar variable */
-	{
-		return( putNCv_double(ncp, varp, start, 1, value) );
-	}
-
-	if(IS_RECVAR(varp))
-	{
-		status = NCvnrecs(ncp, *start + *edges);
-		if(status != NC_NOERR)
-			return status;
-
-		if(varp->ndims == 1
-			&& ncp->recsize <= varp->len)
-		{
-			/* one dimensional && the only record variable  */
-			return( putNCv_double(ncp, varp, start, *edges, value) );
-		}
-	}
-
-	/*
-	 * find max contiguous
-	 *   and accumulate max count for a single io operation
-	 */
-	ii = NCiocount(ncp, varp, edges, &iocount);
-
-	if(ii == -1)
-	{
-		return( putNCv_double(ncp, varp, start, iocount, value) );
-	}
-
-	assert(ii >= 0);
-
-
-	{ /* inline */
-	ALLOC_ONSTACK(coord, size_t, varp->ndims);
-	ALLOC_ONSTACK(upper, size_t, varp->ndims);
-	const size_t index = ii;
-
-	/* copy in starting indices */
-	(void) memcpy(coord, start, varp->ndims * sizeof(size_t));
-
-	/* set up in maximum indices */
-	set_upper(upper, start, edges, &upper[varp->ndims]);
-
-	/* ripple counter */
-	while(*coord < *upper)
-	{
-		const int lstatus = putNCv_double(ncp, varp, coord, iocount,
-				 value);
-		if(lstatus != NC_NOERR)
-		{
-			if(lstatus != NC_ERANGE)
-			{
-				status = lstatus;
-				/* fatal for the loop */
-				break;
-			}
-			/* else NC_ERANGE, not fatal for the loop */
-			if(status == NC_NOERR)
-				status = lstatus;
-		}
-		value += iocount;
-		odo1(start, upper, coord, &upper[index], &coord[index]);
-	}
-
-	FREE_ONSTACK(upper);
-	FREE_ONSTACK(coord);
-	} /* end inline */
-
-	return status;
-}
-
-
-
-
-int
-nc_get_vara_text(int ncid, int varid,
-	 const size_t *start, const size_t *edges, char *value)
-{
-	int status = NC_NOERR;
-	NC *ncp;
-	const NC_var *varp;
-	int ii;
-	size_t iocount;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type != NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, start);
-	if(status != NC_NOERR)
-		return status;
-	status = NCedgeck(ncp, varp, start, edges);
-	if(status != NC_NOERR)
-		return status;
-
-	if(varp->ndims == 0) /* scalar variable */
-	{
-		return( getNCv_text(ncp, varp, start, 1, value) );
-	}
-
-	if(IS_RECVAR(varp))
-	{
-		if(*start + *edges > NC_get_numrecs(ncp))
-			return NC_EEDGE;
-		if(varp->ndims == 1 && ncp->recsize <= varp->len)
-		{
-			/* one dimensional && the only record variable  */
-			return( getNCv_text(ncp, varp, start, *edges, value) );
-		}
-	}
-
-	/*
-	 * find max contiguous
-	 *   and accumulate max count for a single io operation
-	 */
-	ii = NCiocount(ncp, varp, edges, &iocount);
-
-	if(ii == -1)
-	{
-		return( getNCv_text(ncp, varp, start, iocount, value) );
-	}
-
-	assert(ii >= 0);
-
-
-	{ /* inline */
-	ALLOC_ONSTACK(coord, size_t, varp->ndims);
-	ALLOC_ONSTACK(upper, size_t, varp->ndims);
-	const size_t index = ii;
-
-	/* copy in starting indices */
-	(void) memcpy(coord, start, varp->ndims * sizeof(size_t));
-
-	/* set up in maximum indices */
-	set_upper(upper, start, edges, &upper[varp->ndims]);
-
-	/* ripple counter */
-	while(*coord < *upper)
-	{
-		const int lstatus = getNCv_text(ncp, varp, coord, iocount,
-				value);
-		if(lstatus != NC_NOERR)
-		{
-			if(lstatus != NC_ERANGE)
-			{
-				status = lstatus;
-				/* fatal for the loop */
-				break;
-			}
-			/* else NC_ERANGE, not fatal for the loop */
-			if(status == NC_NOERR)
-				status = lstatus;
-		}
-		value += iocount;
-		odo1(start, upper, coord, &upper[index], &coord[index]);
-	}
-
-	FREE_ONSTACK(upper);
-	FREE_ONSTACK(coord);
-	} /* end inline */
-
-	return status;
-}
-
-
-int
-nc_get_vara_uchar(int ncid, int varid,
-	 const size_t *start, const size_t *edges, uchar *value)
-{
-	int status = NC_NOERR;
-	NC *ncp;
-	const NC_var *varp;
-	int ii;
-	size_t iocount;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, start);
-	if(status != NC_NOERR)
-		return status;
-	status = NCedgeck(ncp, varp, start, edges);
-	if(status != NC_NOERR)
-		return status;
-
-	if(varp->ndims == 0) /* scalar variable */
-	{
-		return( getNCv_uchar(ncp, varp, start, 1, value) );
-	}
-
-	if(IS_RECVAR(varp))
-	{
-		if(*start + *edges > NC_get_numrecs(ncp))
-			return NC_EEDGE;
-		if(varp->ndims == 1 && ncp->recsize <= varp->len)
-		{
-			/* one dimensional && the only record variable  */
-			return( getNCv_uchar(ncp, varp, start, *edges, value) );
-		}
-	}
-
-	/*
-	 * find max contiguous
-	 *   and accumulate max count for a single io operation
-	 */
-	ii = NCiocount(ncp, varp, edges, &iocount);
-
-	if(ii == -1)
-	{
-		return( getNCv_uchar(ncp, varp, start, iocount, value) );
-	}
-
-	assert(ii >= 0);
-
-
-	{ /* inline */
-	ALLOC_ONSTACK(coord, size_t, varp->ndims);
-	ALLOC_ONSTACK(upper, size_t, varp->ndims);
-	const size_t index = ii;
-
-	/* copy in starting indices */
-	(void) memcpy(coord, start, varp->ndims * sizeof(size_t));
-
-	/* set up in maximum indices */
-	set_upper(upper, start, edges, &upper[varp->ndims]);
-
-	/* ripple counter */
-	while(*coord < *upper)
-	{
-		const int lstatus = getNCv_uchar(ncp, varp, coord, iocount,
-				value);
-		if(lstatus != NC_NOERR)
-		{
-			if(lstatus != NC_ERANGE)
-			{
-				status = lstatus;
-				/* fatal for the loop */
-				break;
-			}
-			/* else NC_ERANGE, not fatal for the loop */
-			if(status == NC_NOERR)
-				status = lstatus;
-		}
-		value += iocount;
-		odo1(start, upper, coord, &upper[index], &coord[index]);
-	}
-
-	FREE_ONSTACK(upper);
-	FREE_ONSTACK(coord);
-	} /* end inline */
-
-	return status;
-}
-
-int
-nc_get_vara_schar(int ncid, int varid,
-	 const size_t *start, const size_t *edges, schar *value)
-{
-	int status = NC_NOERR;
-	NC *ncp;
-	const NC_var *varp;
-	int ii;
-	size_t iocount;
-
-	status = NC_check_id(ncid, &ncp); 
-	if(status != NC_NOERR)
-		return status;
-
-	if(NC_indef(ncp))
-		return NC_EINDEFINE;
-
-	varp = NC_lookupvar(ncp, varid);
-	if(varp == NULL)
-		return NC_ENOTVAR; /* TODO: lost NC_EGLOBAL */
-
-	if(varp->type == NC_CHAR)
-		return NC_ECHAR;
-
-	status = NCcoordck(ncp, varp, start);
-	if(status != NC_NOERR)
-		return st