Adding set_default_fillpatch_bc() for intermittency, and AMREX_ALWAYS…

…_ASSERT() for num_ghost and fld.num_grow() in TiogaInterface::register_solution()
Exawind · Sep 19, 2024 · 3b250e4 · 3b250e4
1 parent 74bd5ac
commit 3b250e4
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 3 deletions.
diff --git a/amr-wind/overset/TiogaInterface.cpp b/amr-wind/overset/TiogaInterface.cpp
@@ -199,6 +199,7 @@ void TiogaInterface::register_solution(
             auto& fld = repo.get_field(cvar);
             const int ncomp = fld.num_comp();
             fld.fillpatch(m_sim.time().new_time());
+            AMREX_ALWAYS_ASSERT(amrex::IntVect(num_ghost) == fld.num_grow());
             field_ops::copy(*m_qcell, fld, 0, icomp, ncomp, num_ghost);
             icomp += ncomp;
         }

diff --git a/amr-wind/physics/Intermittency.cpp b/amr-wind/physics/Intermittency.cpp
@@ -6,17 +6,19 @@
 namespace amr_wind {
 
 Intermittency::Intermittency(const CFDSim& sim)
-    : m_intermittency(sim.repo().declare_field("intermittency",1,0,1))
+    : m_intermittency(sim.repo().declare_field(
+          "intermittency", 1, sim.pde_manager().num_ghost_state()))
 {
     amrex::ParmParse pp("incflo");
     pp.query("gamma_intermittency", m_gamma);
+    auto& gamma_int_fld = sim.repo().get_field("intermittency");
+    gamma_int_fld.set_default_fillpatch_bc(sim.time());
 }
 
 /** Initialize the intermittency field at the beginning of the
  *  simulation.
  */
-void Intermittency::initialize_fields(
-    int level, const amrex::Geometry& geom)
+void Intermittency::initialize_fields(int level, const amrex::Geometry& geom)
 {
     using namespace utils;