Better output for why Newton Iters may crash. (erf-model#1938)

ewquon · Nov 7, 2024 · 583c93b · 583c93b
1 parent e237f0e
commit 583c93b
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/Source/Utils/ERF_HSE_utils.H b/Source/Utils/ERF_HSE_utils.H
@@ -54,6 +54,7 @@ namespace HSEutils
             Real lo = getRhogivenThetaPress(T, P - eps, RdOCp, qv);
             Real dFdp = 1.0 + 0.25*ieps*(hi - lo)*g*dz;
             P -= F/dFdp;
+            if (P < 1.0e3) amrex::Warning("P < 1000 [Pa]; Domain height may be too large...");
             AMREX_ALWAYS_ASSERT(P > 0.0);
 
             // Diagnose density and residual