CableSurface.h

#pragma once
 
#include <fhsim_environment/BodySubmergence.h>
#include <marenv/EnvironmentFacade.h>
 
#include <algorithm>
#include <cmath>
#include <memory>
 
namespace cable_surface
{
 
template<class Element>
void SetSurface(marenv::EnvironmentFacade& environment, double maxWaveElevation, double T, const double posA[3], const double posB[3], Element& element)
{
   const double halfSize = element.SurfaceHalfSize();
   const double top      = std::min(posA[2], posB[2]) - halfSize; // shallowest point, z down
   const double bottom   = std::max(posA[2], posB[2]) + halfSize;
   if (top > maxWaveElevation || bottom < -maxWaveElevation) {
      element.SetSurfaceElevation(0.0, 0.0);
      return;
   }
   double etaA = 0.0, etaB = 0.0;
   if (environment.GetSurfaceElevation(T, posA, etaA) != marenv::Status::OK)
      etaA = 0.0;
   if (environment.GetSurfaceElevation(T, posB, etaB) != marenv::Status::OK)
      etaB = 0.0;
   element.SetSurfaceElevation(etaA, etaB);
}
 
// The gated surface queries (R66, R67; P-off, MARE-0343) are fhsim_environment's (ENV-0090, MARE-0344).
using environment::body_submergence::CylinderWetFraction;
using environment::body_submergence::MeanWaveAcceleration;
using environment::body_submergence::Submergence;
using environment::body_submergence::WaveAcceleration;
using environment::body_submergence::WavesWithElevation;
 
template<class Element>
void SetWaveAcceleration(marenv::wave::WaveField* waves, double T, const double posA[3], const double posB[3], Element& element)
{
   double mean[3];
   MeanWaveAcceleration(waves, T, posA, posB, mean);
   element.SetWaveAcceleration(mean);
}
 
} // namespace cable_surface