NetCdfSeafloorBathymetry.h

#pragma once
 
#include "util/NetCDFParser.h"
 
#include <marenv/seafloor/BathymetryField.h>
#include <cstdint>
#include <memory>
#include <string>
#include <vector>
 
namespace environmentmodels
{
 
class NetCdfSeafloorBathymetry : public marenv::seafloor::BathymetryField
{
public:
   explicit NetCdfSeafloorBathymetry(const std::string& filename);
 
   [[nodiscard]] std::unique_ptr<NetCdfSeafloorBathymetry> MakeVisualTerrain(std::uint32_t seed) const;
 
   [[nodiscard]] marenv::Status GetDepth(const double pos[2], double& depthOut) override;
 
   [[nodiscard]] marenv::Status ComputeSurfaceUnitNormal(const double pos[2], double unitNormalOut[3]) override;
 
   [[nodiscard]] marenv::Status GetSoilState(const double pos[3], marenv::seafloor::SoilState& soilOut) override;
 
   const netcdf_sinmod::Grid& GetGrid() const { return m_grid; }
 
   const double* GetDepthMap() const { return m_depthMap.data(); }
 
   double GetLengthX() const { return m_length_x; }
 
   double GetLengthY() const { return m_length_y; }
 
private:
   NetCdfSeafloorBathymetry(const netcdf_sinmod::Grid& grid, std::vector<double> depthMap);
 
   double GetDepthInternal(const double pos[2]) const;
 
   double DepthSlope(const double pos[2], int axis) const;
 
   netcdf_sinmod::Grid m_grid {};
   std::vector<double> m_depthMap; 
   double              m_length_x = 0.0;
   double              m_length_y = 0.0;
};
 
} // namespace environmentmodels