/* This is part of the netCDF package.
   Copyright 2020 University Corporation for Atmospheric Research/Unidata
   See COPYRIGHT file for conditions of use.

   Test netcdf-4 variables with unlimited dimensions.

   Ed Hartnett
*/

#include <config.h>
#include <nc_tests.h>
#include "err_macros.h"

#ifdef TESTNCZARR
#define FILE_NAME "file://tmp_unlim_vars.file#mode=nczarr,file"
#else
#define FILE_NAME "tst_unlim_vars.nc"
#endif

#define SFC_TEMP_NAME "surface_temperature"
#define LAT_NAME "lat"
#define LAT_LEN 2
#define LON_NAME "lon"
#define LON_LEN 3
#define TIME_NAME "time"
#define NUM_TIMESTEPS 4
#define NDIM3 3
#define NDIM2 2

int
main(int argc, char **argv)
{
    int stat = NC_NOERR;
    
    printf("\n*** Testing netcdf-4 variables with unlimited dimensions.\n");
    printf("*** Testing file with one var, one unlim dim...");
    {
        int ncid, sfc_tempid;
        float data_out[NUM_TIMESTEPS][LAT_LEN][LON_LEN], data_in[NUM_TIMESTEPS][LAT_LEN][LON_LEN];
        int lat, lon, time;
        int dimids[NDIM3];
        nc_type xtype_in;
        int ndims_in, dimids_in[10], natts_in, nvars_in, unlimdimid_in;
        size_t len_in;
        char name_in[NC_MAX_NAME+1];
        size_t start[NDIM3], count[NDIM3];
        int d;


        /* Set up phony data. */
        for (time = 0; time < NUM_TIMESTEPS; time++)
            for (lat = 0; lat < LAT_LEN; lat++)
                for (lon = 0; lon < LON_LEN; lon++)
                    data_out[time][lat][lon] = 25.5f + (float)(lat + lon + time);

        /* Create a file with a 3D surface temp variable. */
        if (nc_create(FILE_NAME, NC_NETCDF4, &ncid)) ERR;

        /* Create three dims, one unlimited. */
        if (nc_def_dim(ncid, TIME_NAME, NC_UNLIMITED, &dimids[0])) ERR;
        if (nc_def_dim(ncid, LAT_NAME, LAT_LEN, &dimids[1])) ERR;
        if (nc_def_dim(ncid, LON_NAME, LON_LEN, &dimids[2])) ERR;

        /* Define a var. */
        for (d = 0; d < NDIM3; d++)
            dimids[d] = d;
        if (nc_def_var(ncid, SFC_TEMP_NAME, NC_FLOAT, NDIM3, dimids, &sfc_tempid)) ERR;

        /* Check some metadata. */
        if (nc_inq(ncid, &ndims_in, &nvars_in, &natts_in, &unlimdimid_in)) ERR;
        if (ndims_in != 3 || nvars_in != 1 || natts_in != 0 || unlimdimid_in != 0) ERR;
        if (nc_inq_var(ncid, 0, name_in, &xtype_in, &ndims_in, dimids_in,
                       &natts_in)) ERR;
        if (strcmp(name_in, SFC_TEMP_NAME) || xtype_in != NC_FLOAT ||
            ndims_in != 3 || natts_in != 0) ERR;
        for (d = 0; d < NDIM3; d++)
            if (dimids_in[d] != dimids[d]) ERR;
        if (nc_inq_dim(ncid, 0, name_in, &len_in)) ERR;
        if (len_in != 0 || strcmp(name_in, TIME_NAME)) ERR;
        if (nc_inq_dim(ncid, 1, name_in, &len_in)) ERR;
        if (len_in != LAT_LEN || strcmp(name_in, LAT_NAME)) ERR;
        if (nc_inq_dim(ncid, 2, name_in, &len_in)) ERR;
        if (len_in != LON_LEN || strcmp(name_in, LON_NAME)) ERR;

	{
	    size_t chunksizes[3];
	    int i,storage;
	    fprintf(stderr,">>> chunks =");
            if (nc_inq_var_chunking(ncid, sfc_tempid, &storage, chunksizes)) ERR;
	    for(i=0;i<3;i++) fprintf(stderr," %llu", (unsigned long long)chunksizes[i]);
	    fprintf(stderr," ; storage=%d\n", storage);
	}

        if (nc_close(ncid)) ERR;

        if ((stat=nc_open(FILE_NAME, 0, &ncid)))
	    ERR;

        /* Check metadata. */
        if (nc_inq(ncid, &ndims_in, &nvars_in, &natts_in, &unlimdimid_in)) ERR;
        if (ndims_in != 3 || nvars_in != 1 || natts_in != 0 || unlimdimid_in != 0) ERR;
        if (nc_inq_var(ncid, 0, name_in, &xtype_in, &ndims_in, dimids_in,
                       &natts_in)) ERR;
        if (strcmp(name_in, SFC_TEMP_NAME) || xtype_in != NC_FLOAT ||
            ndims_in != 3 || natts_in != 0) ERR;
        for (d = 0; d < NDIM3; d++)
            if (dimids_in[d] != dimids[d]) ERR;
        if (nc_inq_dim(ncid, 0, name_in, &len_in)) ERR;
        if (len_in != 0 || strcmp(name_in, TIME_NAME)) ERR;
        if (nc_inq_dim(ncid, 1, name_in, &len_in)) ERR;
        if (len_in != LAT_LEN || strcmp(name_in, LAT_NAME)) ERR;
        if (nc_inq_dim(ncid, 2, name_in, &len_in)) ERR;
        if (len_in != LON_LEN || strcmp(name_in, LON_NAME)) ERR;

        if (nc_close(ncid)) ERR;

        if (nc_open(FILE_NAME, NC_WRITE, &ncid)) ERR;

        /* Check metadata. */
        if (nc_inq(ncid, &ndims_in, &nvars_in, &natts_in, &unlimdimid_in)) ERR;
        if (ndims_in != 3 || nvars_in != 1 || natts_in != 0 || unlimdimid_in != 0) ERR;
        if (nc_inq_var(ncid, 0, name_in, &xtype_in, &ndims_in, dimids_in,
                       &natts_in)) ERR;
        if (strcmp(name_in, SFC_TEMP_NAME) || xtype_in != NC_FLOAT ||
            ndims_in != 3 || natts_in != 0) ERR;
        for (d = 0; d < NDIM3; d++)
            if (dimids_in[d] != dimids[d]) ERR;
        if (nc_inq_dim(ncid, 0, name_in, &len_in)) ERR;
        if (len_in != 0 || strcmp(name_in, TIME_NAME)) ERR;
        if (nc_inq_dim(ncid, 1, name_in, &len_in)) ERR;
        if (len_in != LAT_LEN || strcmp(name_in, LAT_NAME)) ERR;
        if (nc_inq_dim(ncid, 2, name_in, &len_in)) ERR;
        if (len_in != LON_LEN || strcmp(name_in, LON_NAME)) ERR;
        if (nc_close(ncid)) ERR;

        /* Write some data to the var.*/
        if (nc_open(FILE_NAME, NC_WRITE, &ncid)) ERR;
        for (d = 0; d < NDIM3; d++)
            start[d] = 0;
        count[0] = NUM_TIMESTEPS;
        count[1] = LAT_LEN;
        count[2] = LON_LEN;
        if (nc_put_vara_float(ncid, 0, start, count, (float *)data_out)) ERR;

        /* Check metadata. */
        if (nc_inq(ncid, &ndims_in, &nvars_in, &natts_in, &unlimdimid_in)) ERR;
        if (ndims_in != 3 || nvars_in != 1 || natts_in != 0 || unlimdimid_in != 0) ERR;
        if (nc_inq_var(ncid, 0, name_in, &xtype_in, &ndims_in, dimids_in,
                       &natts_in)) ERR;
        if (strcmp(name_in, SFC_TEMP_NAME) || xtype_in != NC_FLOAT ||
            ndims_in != 3 || natts_in != 0) ERR;
        for (d = 0; d < NDIM3; d++)
            if (dimids_in[d] != dimids[d]) ERR;
        if (nc_inq_dim(ncid, 0, name_in, &len_in)) ERR;
        if (len_in != NUM_TIMESTEPS || strcmp(name_in, TIME_NAME)) ERR;
        if (nc_inq_dim(ncid, 1, name_in, &len_in)) ERR;
        if (len_in != LAT_LEN || strcmp(name_in, LAT_NAME)) ERR;
        if (nc_inq_dim(ncid, 2, name_in, &len_in)) ERR;
        if (len_in != LON_LEN || strcmp(name_in, LON_NAME)) ERR;

        if (nc_close(ncid)) ERR;

        if (nc_open(FILE_NAME, NC_WRITE, &ncid)) ERR;

        /* Check data. */
        if (nc_get_vara_float(ncid, 0, start, count, (float *)data_in)) ERR;
        for (time = 0; time < NUM_TIMESTEPS; time++)
            for (lat = 0; lat < LAT_LEN; lat++)
                for (lon = 0; lon < LON_LEN; lon++)
                    if (data_in[time][lat][lon] != data_out[time][lat][lon]) ERR;

        if (nc_close(ncid)) ERR;
    }
    SUMMARIZE_ERR;
    printf("*** Testing netcdf-4 variable 2 unlimited dimensions...");
    {
#define DIM0_NAME "Speed"
#define DIM1_NAME "Height"
#define VAR_NAME "Superman"
#define VAR1_NAME "Black_Widow"
        int ncid, varid, varid1, dimid[NDIM2];
        int data = TEST_VAL_42;
        int data_in;
        size_t index[NDIM2] = {1, 1};
        char name_in[NC_MAX_NAME + 1];
        size_t len_in;

        /* Create a file with 2 unlimited dims. */
        if (nc_create(FILE_NAME, NC_NETCDF4, &ncid)) ERR;
        if (nc_def_dim(ncid, DIM0_NAME, NC_UNLIMITED, &dimid[0])) ERR;
        if (nc_def_dim(ncid, DIM1_NAME, NC_UNLIMITED, &dimid[1])) ERR;
        if (nc_def_var(ncid, VAR_NAME, NC_INT, NDIM2, dimid, &varid)) ERR;
        if (nc_def_var(ncid, VAR1_NAME, NC_INT, NDIM2, dimid, &varid1)) ERR;
        if (nc_close(ncid)) ERR;

        /* Check the file. */
        if (nc_open(FILE_NAME, NC_NOWRITE, &ncid)) ERR;
        if (nc_inq_dim(ncid, 0, name_in, &len_in)) ERR;
        if (strcmp(name_in, DIM0_NAME) || len_in != 0) ERR;
        if (nc_inq_dim(ncid, 1, name_in, &len_in)) ERR;
        if (strcmp(name_in, DIM1_NAME) || len_in != 0) ERR;
        if (nc_close(ncid)) ERR;

        /* Reopen the file and add data. */
        if (nc_open(FILE_NAME, NC_WRITE, &ncid)) ERR;
        if (nc_put_var1_int(ncid, 0, index, &data)) ERR;
        if (nc_inq_dim(ncid, 0, NULL, &len_in)) ERR;
        if (len_in != 2) ERR;
        if (nc_inq_dim(ncid, 1, NULL, &len_in)) ERR;
        if (len_in != 2) ERR;
        if (nc_get_var1_int(ncid, 0, index, &data_in)) ERR;
        if (data_in != data) ERR;
        if (nc_get_var1_int(ncid, 1, index, &data_in)) ERR;
        if (data_in != NC_FILL_INT) ERR;
        if (nc_close(ncid)) ERR;

    }
    SUMMARIZE_ERR;
    printf("*** Testing netcdf-4 variable 3 unlimited dimensions...");
    {
#define D3_DIM0_NAME "Gadgets"
#define D3_DIM1_NAME "Brains"
#define D3_DIM2_NAME "Money"
#define D3_VAR_NAME "Batman"
#define D3_VAR1_NAME "Aquaman"
        int ncid, varid, dimid[NDIM3];
        int data = TEST_VAL_42;
        int data_in;
        size_t start[NDIM3] = {1, 1, 1}, count[NDIM3] = {1, 1, 1};
        size_t index[NDIM3] = {0, 0, 1};
        char name_in[NC_MAX_NAME + 1];
        size_t len_in;

        /* Create a file with 2 unlimited dims. */
        if (nc_create(FILE_NAME, NC_NETCDF4, &ncid)) ERR;
        if (nc_def_dim(ncid, D3_DIM0_NAME, NC_UNLIMITED, &dimid[0])) ERR;
        if (nc_def_dim(ncid, D3_DIM1_NAME, NC_UNLIMITED, &dimid[1])) ERR;
        if (nc_def_dim(ncid, D3_DIM2_NAME, NC_UNLIMITED, &dimid[2])) ERR;
        if (nc_def_var(ncid, D3_VAR_NAME, NC_INT, NDIM3, dimid, &varid)) ERR;
        if (nc_def_var(ncid, D3_VAR1_NAME, NC_INT, NDIM3, dimid, &varid)) ERR;
        if (nc_close(ncid)) ERR;

        /* Check the file. */
        if (nc_open(FILE_NAME, NC_NOWRITE, &ncid)) ERR;
        if (nc_inq_dim(ncid, 0, name_in, &len_in)) ERR;
        if (strcmp(name_in, D3_DIM0_NAME) || len_in != 0) ERR;
        if (nc_inq_dim(ncid, 1, name_in, &len_in)) ERR;
        if (strcmp(name_in, D3_DIM1_NAME) || len_in != 0) ERR;
        if (nc_inq_dim(ncid, 2, name_in, &len_in)) ERR;
        if (strcmp(name_in, D3_DIM2_NAME) || len_in != 0) ERR;
        if (nc_close(ncid)) ERR;

        /* Reopen the file and add data. */
        if (nc_open(FILE_NAME, NC_WRITE, &ncid)) ERR;
        if (nc_put_vara_int(ncid, 0, start, count, &data)) ERR;
        if (nc_put_var1_int(ncid, 1, index, &data)) ERR;
        if (nc_inq_dim(ncid, 0, NULL, &len_in)) ERR;
        if (len_in != 2) ERR;
        if (nc_inq_dim(ncid, 1, NULL, &len_in)) ERR;
        if (len_in != 2) ERR;
        if (nc_inq_dim(ncid, 2, NULL, &len_in)) ERR;
        if (len_in != 2) ERR;
        if (nc_get_vara_int(ncid, 0, start, count, &data_in)) ERR;
        if (data_in != data) ERR;
        if (nc_get_vara_int(ncid, 1, start, count, &data_in)) ERR;
        if (data_in != NC_FILL_INT) ERR;
        if (nc_close(ncid)) ERR;

        /* Reopen file and check again. */
        if (nc_open(FILE_NAME, NC_NOWRITE, &ncid)) ERR;
        if (nc_inq_dim(ncid, 0, NULL, &len_in)) ERR;
        if (len_in != 2) ERR;
        if (nc_inq_dim(ncid, 1, NULL, &len_in)) ERR;
        if (len_in != 2) ERR;
        if (nc_inq_dim(ncid, 2, NULL, &len_in)) ERR;
        if (len_in != 2) ERR;
        if (nc_get_vara_int(ncid, 0, start, count, &data_in)) ERR;
        if (data_in != data) ERR;
        if (nc_get_vara_int(ncid, 1, start, count, &data_in)) ERR;
        if (data_in != NC_FILL_INT) ERR;
        if (nc_get_var1_int(ncid, 1, index, &data_in)) ERR;
        if (data_in != data) ERR;
        if (nc_close(ncid)) ERR;
    }
    SUMMARIZE_ERR;
    FINAL_RESULTS;
}
