/* readvis5d.c */


/*
 *  This is an example  program to show how to read the header and 
 *  data arrays from vis5d data files.  The structure of this program 
 *  comes from v5dstats.c 
 */


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

/*  need to have the path set to the v5d.h file in the vis5d directory tree */

#include "../vis5d-4.1/src/v5d.h"



int main( int argc, char *argv[] )
{
   v5dstruct v;
   int time, var;

  /*
   *  Check the inputs
   */

   if (argc!=2) {
      printf("Usage:\n");
      printf("   v5dstats file\n");
      exit(0);
   }

  /*
   *  Open the input vis5d file and read in the vis5d header.  The header
   *  variables are stored in the structure v.
   */

   if (!v5dOpenFile( argv[1], &v )) {
      printf("Error: couldn't open %s for reading\n", argv[1] );
      exit(0);
   }
   
  /*
   *  Now the data arrays can be read in.  The data arrays are stored for
   *  each time and each variable.  You send the time and the variable of
   *  the data array you want to read to v5dReadGrid and it returns the
   *  needed data array.  Space is opened up for the data in the malloc
   *  command.  If you want more than one data array available at once,
   *  you have to create a new variable to store the data.
   */

   for (time=0; time<v.NumTimes; time++) {

      for (var=0; var<v.NumVars; var++) {
	 int i, nrncnl;
	 float *data;
	 float min, max, sum, sumsum;
	 int missing, good;

        /*
         *  nrncnl --  Total number of data points.  (Number of rows *
         *             Number of columns * Number of levels)
         */

	 nrncnl = v.Nr * v.Nc * v.Nl[var];

        /*
         *  Open up the space for the data.
         */

	 data = (float *) malloc( nrncnl * sizeof(float) );

        /*
         *  Read in the data, given the time and the variable wanted.
         */

	 if (!v5dReadGrid( &v, time, var, data )) {
	    printf("Error while reading grid (time=%d,var=%s)\n",
		   time+1, v.VarName[var] );
	    exit(0);
	 }

        /*
         *  The (0,0) point is the Northern most and Western most point.
         *  The level increase from bottom to top.
         *  The Volume Imaging Lidar is located at the center of the
         *  image.  The center point is at ((nr + 1) / 2, (nc + 1) / 2).
         *  The data points are determined by the following formula:
         *  row+nr*(col+nc*lev).  Therefore, the point where the VIL is 
         *  located is (nr + 1) / 2 + nr * (((nc + 1) / 2) + nc * lev),
         *  where lev is the points above the VIL.
         */

        /*  
         *  This is an example section which counts the number of missing
         *  data points in the data array.  Missing data points are given
         *  the value of MISSING which is defined in the file v5d.h.

	 min = MISSING;
	 max = -MISSING;
	 missing = 0;
	 good = 0;
	 sum = 0.0;
	 sumsum = 0.0;

	 for (i=0;i<nrncnl;i++) {
	    if ( IS_MISSING(data[i]) ) {
	       missing++;
	    }
	    else {
	       good++;
	       if (data[i]<min) {
		  min = data[i];
	       }
	       if (data[i]>max) {
		  max = data[i];
	       }
	       sum += data[i];
	       sumsum += data[i]*data[i];
	    }
	 }

         */

	 free( data );
      }
      printf("\n");
   }

   v5dCloseFile( &v );
   return 0;
}