Remove 90 deg Latitude Constraint from Hurricane Forcing (#872)

* Removing the constraint that latitude must be 90 degrees for the hurricane forcing. * Small update for formatting. * Ran clang-format... * Putting back what I deleted with clang-format... Nice job, Matt... --------- Co-authored-by: Matt Churchfield <[email protected]>
Exawind · Jul 12, 2023 · 05367b3 · 05367b3
1 parent ef0a198
commit 05367b3
Showing 1 changed file with 4 additions and 6 deletions.
diff --git a/amr-wind/equation_systems/icns/source_terms/HurricaneForcing.cpp b/amr-wind/equation_systems/icns/source_terms/HurricaneForcing.cpp
@@ -21,14 +21,12 @@ HurricaneForcing::HurricaneForcing(const CFDSim& sim) : m_mesh(sim.mesh())
         amrex::ParmParse pp("CoriolisForcing");
         amrex::Real rot_time_period = 86400.0;
         pp.query("rotational_time_period", rot_time_period);
-        m_coriolis_factor = 2.0 * utils::two_pi() / rot_time_period;
-        amrex::Print() << "Geostrophic forcing: Coriolis factor = "
-                       << m_coriolis_factor << std::endl;
         amrex::Real latitude = 90.0;
         pp.query("latitude", latitude);
-        AMREX_ALWAYS_ASSERT(
-            amrex::Math::abs(latitude - 90.0) <
-            static_cast<amrex::Real>(vs::DTraits<float>::eps()));
+        m_coriolis_factor = (2.0 * utils::two_pi() / rot_time_period) *
+                            std::sin(utils::radians(latitude));
+        amrex::Print() << "Geostrophic forcing: Coriolis factor = "
+                       << m_coriolis_factor << std::endl;
     }
 
     {