From 4d5fa6dffa0956e49a0b5e47751660b8e4019c78 Mon Sep 17 00:00:00 2001
From: AMLattanzi <AMLattanzi@lbl.gov>
Date: Thu, 7 Nov 2024 13:41:12 -0800
Subject: [PATCH] Better output for why Newton Iters may crash.

---
 Source/Utils/ERF_HSE_utils.H | 1 +
 1 file changed, 1 insertion(+)

diff --git a/Source/Utils/ERF_HSE_utils.H b/Source/Utils/ERF_HSE_utils.H
index ac63109b3..2cab30782 100644
--- 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