We need to create SmnSmn when using MYNN.

erf-model · Oct 24, 2024 · 49cbb8e · 49cbb8e
1 parent 596871b
commit 49cbb8e
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/Source/ERF_make_new_arrays.cpp b/Source/ERF_make_new_arrays.cpp
@@ -362,7 +362,8 @@ ERF::update_diffusive_arrays (int lev, const BoxArray& ba, const DistributionMap
                            (solverChoice.turbChoice[lev].pbl_type        !=       PBLType::None) );
     bool l_use_kturb   = ( (solverChoice.turbChoice[lev].les_type        != LESType::None)   ||
                            (solverChoice.turbChoice[lev].pbl_type        != PBLType::None) );
-    bool l_use_ddorf   = (  solverChoice.turbChoice[lev].les_type        == LESType::Deardorff);
+    bool l_use_ke      = (  (solverChoice.turbChoice[lev].les_type       == LESType::Deardorff) ||
+                            (solverChoice.turbChoice[lev].pbl_type       == PBLType::MYNN25) );
     bool l_use_moist   = (  solverChoice.moisture_type != MoistureType::None  );
 
     BoxArray ba12 = convert(ba, IntVect(1,1,0));
@@ -429,7 +430,7 @@ ERF::update_diffusive_arrays (int lev, const BoxArray& ba, const DistributionMap
     if (l_use_kturb) {
         eddyDiffs_lev[lev] = std::make_unique<MultiFab>(ba, dm, EddyDiff::NumDiffs, 2);
         eddyDiffs_lev[lev]->setVal(0.0);
-        if(l_use_ddorf) {
+        if(l_use_ke) {
             SmnSmn_lev[lev] = std::make_unique<MultiFab>( ba, dm, 1, 0 );
         } else {
             SmnSmn_lev[lev] = nullptr;