Add KE_min param for ABL

ewquon · Oct 7, 2024 · da9afc9 · da9afc9
1 parent f8bcb64
commit da9afc9
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 4 deletions.
diff --git a/Exec/ABL/ERF_prob.H b/Exec/ABL/ERF_prob.H
@@ -13,6 +13,7 @@ struct ProbParm : ProbParmDefaults {
   amrex::Real A_0   = 1.0;
   amrex::Real KE_0  = 0.1;
 
+  amrex::Real KE_min = 1e-13;
   amrex::Real KE_decay_height = -1;
   amrex::Real KE_decay_order  =  1;
 

diff --git a/Exec/ABL/ERF_prob.cpp b/Exec/ABL/ERF_prob.cpp
@@ -17,6 +17,7 @@ Problem::Problem(const amrex::Real* problo, const amrex::Real* probhi)
   pp.query("T_0", parms.T_0);
   pp.query("A_0", parms.A_0);
   pp.query("KE_0", parms.KE_0);
+  pp.query("KE_min", parms.KE_min);
   pp.query("KE_decay_height", parms.KE_decay_height);
   pp.query("KE_decay_order", parms.KE_decay_order);
 
@@ -130,19 +131,21 @@ Problem::init_custom_pert(
         if (parms_d.KE_decay_height > 0) {
             // scale initial SGS kinetic energy with height
             state_pert(i, j, k, RhoKE_comp) *= max(
-                std::pow(1 - min(z/parms_d.KE_decay_height,1.0), parms_d.KE_decay_order),
-                1e-12);
+                std::pow(1.0 - min(z/parms_d.KE_decay_height,1.0), parms_d.KE_decay_order),
+                0.0);
         }
+        state_pert(i, j, k, RhoKE_comp) = max(state_pert(i, j, k, RhoKE_comp), r_hse(i,j,k)*parms_d.KE_min);
     }
     if (state_pert.nComp() > RhoQKE_comp) {
         // PBL
         state_pert(i, j, k, RhoQKE_comp) = r_hse(i,j,k) * 2.0 * parms_d.KE_0;
         if (parms_d.KE_decay_height > 0) {
             // scale initial SGS kinetic energy with height
             state_pert(i, j, k, RhoQKE_comp) *= max(
-                std::pow(1 - min(z/parms_d.KE_decay_height,1.0), parms_d.KE_decay_order),
-                1e-12);
+                std::pow(1.0 - min(z/parms_d.KE_decay_height,1.0), parms_d.KE_decay_order),
+                0.0);
         }
+        state_pert(i, j, k, RhoQKE_comp) = max(state_pert(i, j, k, RhoQKE_comp), 2.0*r_hse(i,j,k)*parms_d.KE_min);
     }
 
     if (use_moisture) {