netcdf_handler.h

#ifndef NETCDF_HANDLER
#define NETCDF_HANDLER
 
#include "netcdf.h"
#include <math.h>
#include <vector>
#include <string>
#include <iostream>
#include <fstream>
 
 
class NetcdfHandler {
    const int U_V_FILE = 0,
            Z_X_VELOCITY_FILE = 1,
            Z_PROFILE_FILE = 2,
            U_V_W_FILE = 3;
        
        //NcFile *nf;
        int ncid;
        int type = U_V_FILE;
        size_t dim1 = -1, dim2 = -1, dim3 = -1;
        int xdim = -1, ydim = -1, zdim = -1;
        //std::vector<float> geo1, geo2;
        float* geo1;
        float* geo2;
        float* geo3;
        float low1, low2, low3, diff1, diff2, diff3;
        float minVal, meanVal, maxVal;
        bool hasCachedValues = false;
        std::vector<std::vector<std::vector<float>>> currentFieldU, currentFieldV, currentFieldW;
    public:
        /*
        * Open the given NetCDF file and extract basic information about its contents.
        */
        void readNCFile(const char *filename);
 
        /*
        * Read the entire current field of the opened NetCDF file and store in memory. Any lookup of
        * current values using getCurrentSpeedAt() afterwards will be based on the cached values.
        */
        int cacheCurrentField();
 
 
        /*
        * Close the opened NetCDF file. If cacheCurrentField() has been called, lookup of current speed values can be
        * done after closing the file.
        */
        void closeNCFile();
 
        /*
        * Get a 6 element matrix containing the lower and upper bounds for the x, y and z axes of the
        * opened NetCDF file.
        */
        void getDomainBounds(float* bounds);
 
        int getVar(char* name);
        // getVar(char* name, std::vector<size_t> start, std::vector<size_t> count, float value[][2]);
        int getCurrentSpeedAt(float x, float y, float z, float *value);
        int calculateDescriptorValues(float *values);
        bool isProfileFile();
 
};
#endif