/*********************************************************************
 *   Copyright 1993, UCAR/Unidata
 *   See netcdf/COPYRIGHT file for copying and redistribution conditions.
 *   $Id: varget.c 4998 2007-11-16 03:14:03Z epourmal $
 *********************************************************************/

#include <stdio.h>
#include <stdlib.h>

#include "h4config.h"
#ifdef H4_HAVE_NETCDF
#include "netcdf.h"
#else
#include "hdf4_netcdf.h"
#endif

#include "testcdf.h"		/* defines in-memory test cdf structure */
#include "error.h"
#include "tests.h"
#include "alloc.h"
#include "emalloc.h"
#ifdef HDF
#include "hdf.h"
#endif


/*
 * Test ncvarget
 *    check that proper call worked after ncvarput
 *    try with negative coords, edges, check error
 *    try with too-large coords, edges, check error
 *    try with bad variable handle, check error
 *    try in define mode, check error
 *    try with bad netCDF handle, check error
 */
void
test_ncvarget(path)
     char *path;		/* name of writable netcdf file to open */
{
    int nerrs = 0;
    static char pname[] = "test_ncvarget";
    int cdfid;			/* netcdf id */
    int iv;			/* variable id */
    struct cdfhc {		/* a hypercube with generic values */
	long cor[H4_MAX_VAR_DIMS];	/* netcdf coordinates for lower corner */
	long edg[H4_MAX_VAR_DIMS];	/* netcdf edge lengths to upper corner */
	void *vals;		/* pointer to block of values */
    } hc;			/* test hypercube */

    (void) fprintf(stderr, "*** Testing %s ...\t", &pname[5]);

    if ((cdfid = ncopen(path, NC_WRITE)) == -1) {
	error("%s: ncopen failed", pname);
	return;
    }
    /* opened in data mode, try putting and getting hypercubes of each type */
    nerrs += test_varputget (cdfid);

    /* try putting hypercube and getting various interior slabs */
    nerrs += test_slabs (cdfid);
    
    /* find a variable with at least one dimension */
    iv = 0;
    while (test.vars[iv].ndims <= 0 && iv < test.nvars)
      iv++;
    if (iv < test.nvars) {	/* iv is varid of variable with dimensions */
	int tmp;
	/* set coords */
	int id;			/* dimension id */
	for (id = 0; id < test.vars[iv].ndims; id++) {
	    hc.cor[id] = 0;
	    hc.edg[id] = 1;
	}
	/* get space for vals */
	hc.vals = emalloc(nctypelen(test.vars[iv].type) + 8);

	id = test.vars[iv].ndims - 1;
	tmp = hc.cor[id];
	hc.cor[id] = -1;	/* try negative coordinate, should fail */
	if(ncvarget (cdfid, iv, hc.cor, hc.edg, hc.vals) != -1) {
	    error("%s: ncvarget should fail for negative corner", pname);
	    ncclose(cdfid); return;
	}
	hc.cor[id] = tmp;
	tmp = hc.edg[id];
	hc.edg[id] = -1;	/* try negative edge, should fail */
	if(ncvarget (cdfid, iv, hc.cor, hc.edg, hc.vals) != -1) {
	    error("%s: ncvarget should fail for negative edge", pname);
	    ncclose(cdfid); return;
	}
	hc.edg[id] = tmp;

	{ 
		long mqv = test.vars[iv].ndims -1 ;
		int dim = test.vars[iv].dims[mqv] ;

	tmp = hc.cor[mqv];
	hc.cor[mqv] = test.dims[dim].size; /* try big coordinate, should fail */
	if(ncvarget (cdfid, iv, hc.cor, hc.edg, hc.vals) != -1) {
	    error("%s: ncvarget should fail for too-high coordinate", pname);
	    ncclose(cdfid); return;
	}
	hc.cor[mqv] = tmp;
	tmp = hc.edg[mqv];
	hc.edg[mqv] = test.dims[dim].size + 1; /* try big edge, should fail */
	if(ncvarget (cdfid, iv, hc.cor, hc.edg, hc.vals) != -1) {
	    error("%s: ncvarget should fail for too-high edge", pname);
	    ncclose(cdfid); return;
	}
	hc.edg[mqv] = tmp;
	} /* mqv block */

	if (ncredef(cdfid) == -1) {
	    error("%s: ncredef failed", pname);
	    ncclose(cdfid); return;
	}
	/* try in define mode, should fail */
	if(ncvarget (cdfid, iv, hc.cor, hc.edg, hc.vals) != -1) {
	    error("%s: ncvarget should fail in define mode",
		  pname);
	    ncclose(cdfid); return;
	}
	if (ncendef (cdfid) == -1) {
	    error("%s: ncendef failed", pname);
	    ncclose(cdfid); return;
	}
    }
    /* try with bad variable handle, should fail */
    if(ncvarget (cdfid, -1, hc.cor, hc.edg, hc.vals) != -1) {
	error("%s: ncvarget should fail for bad variable handle", pname);
	ncclose(cdfid); return;
    }
    if (ncclose (cdfid) == -1) {
	error("%s: ncclose failed", pname);
	return;
    }
    /* try with bad netCDF handle, should fail */
    if(ncvarget (cdfid, 0, hc.cor, hc.edg, hc.vals) != -1) {
	error("%s: ncvarget failed to report bad netcdf handle", pname);
	nerrs++;
    }
    Free ((char *)hc.vals);
    if (nerrs > 0)
      (void) fprintf(stderr,"FAILED! ***\n");
    else
      (void) fprintf(stderr,"ok ***\n");
}
