From 583c93b438ee395e8a53eb95f72936eeeed0f922 Mon Sep 17 00:00:00 2001
From: "Aaron M. Lattanzi" <103702284+AMLattanzi@users.noreply.github.com>
Date: Thu, 7 Nov 2024 14:04:44 -0800
Subject: [PATCH] Better output for why Newton Iters may crash. (#1938)

---
 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