machine limit the output

AMReX-Astro · Dec 13, 2024 · 83c6103 · 83c6103
1 parent 92aa762
commit 83c6103
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/screening/screen.H b/screening/screen.H
@@ -173,7 +173,11 @@ number_t debye_huckel (const plasma_state_t<number_t>& state,
 
     // eq. A1
     number_t h_DH = z1z2 * gcem::sqrt(3 * amrex::Math::powi<3>(Gamma_e) * z_ratio);
-    return admath::exp(h_DH);
+
+    // machine limit the output
+    constexpr amrex::Real h_max = 300.e0_rt;
+    h = admath::min(h_DH, h_max);
+    return admath::exp(h);
 }
 
 #if SCREEN_METHOD == SCREEN_METHOD_screen5