netcdf_handler.h

#pragma once
 
#include "netcdf.h"
 
#include <cmath>
#include <fstream>
#include <iostream>
#include <string>
#include <vector>
 
class NetcdfHandler
{
   const int U_V_FILE          = 0; 
   const int Z_X_VELOCITY_FILE = 1; 
   const int Z_PROFILE_FILE    = 2; 
   const int U_V_W_FILE        = 3; 
 
   int ncid;        
   int type = U_V_FILE; 
 
   size_t dim1 = -1; 
   size_t dim2 = -1; 
   size_t dim3 = -1; 
 
   int xdim = -1; 
   int ydim = -1; 
   int zdim = -1; 
 
   float* geo1; 
   float* geo2; 
   float* geo3; 
 
   float low1, low2, low3;   
   float diff1, diff2, diff3; 
 
   float minVal;  
   float meanVal; 
   float maxVal;  
 
   bool hasCachedValues = false; 
 
   std::vector<std::vector<std::vector<float>>> currentFieldU; 
   std::vector<std::vector<std::vector<float>>> currentFieldV; 
   std::vector<std::vector<std::vector<float>>> currentFieldW; 
 
public:
   void readNCFile(const char* filename);
 
   int cacheCurrentField();
 
   void closeNCFile();
 
   void getDomainBounds(float* bounds);
 
   int getVar(char* name);
 
   int getCurrentSpeedAt(float x, float y, float z, float* value);
 
   int calculateDescriptorValues(float* values);
 
   bool isProfileFile();
};