Change deep water check

Author: tridelat

src/wave_types.cpp

@@ -11,6 +11,14 @@
 #include <iomanip>
 
 
+bool is_in_deep_water(double wavenumber, double water_depth) {
+    if (wavenumber * water_depth > 89.4) {
+        return true;
+    } else {
+        return false;
+    }
+}
+
 double GetEta(const Eigen::Vector3d& position,
               double time,
               double omega,
@@ -73,7 +81,7 @@ Eigen::Vector3d GetWaterVelocity(const Eigen::Vector3d& position,
 
     // get water velocity
     auto water_velocity = Eigen::Vector3d(0.0, 0.0, 0.0);
-    if (2 * M_PI / wavenumber > water_depth || wavenumber * water_depth > 500.0) {
+    if (is_in_deep_water(wavenumber, water_depth)) {
         // deep water
         water_velocity[0] =
             omega * amplitude * std::exp(wavenumber * z_pos) * cos(wavenumber * x_pos - omega * time + phase);
@@ -105,7 +113,7 @@ Eigen::Vector3d GetWaterAcceleration(const Eigen::Vector3d& position,
 
     // get water velocity
     auto water_acceleration = Eigen::Vector3d(0.0, 0.0, 0.0);
-    if (2 * M_PI / wavenumber > water_depth || wavenumber * water_depth > 500.0) {
+    if (is_in_deep_water(wavenumber, water_depth)) {
         // deep water
         water_acceleration[0] =
             omega * omega * amplitude * std::exp(wavenumber * z_pos) * sin(wavenumber * x_pos - omega * time + phase);