From 44e7625f54f5c196fac8a9baddb18893d7d12aed Mon Sep 17 00:00:00 2001
From: Eliot Quon <eliot.quon@nrel.gov>
Date: Fri, 8 Sep 2023 21:53:08 -0600
Subject: [PATCH] Add initHSE call for init_type==ideal

---
 Source/ERF.cpp | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/Source/ERF.cpp b/Source/ERF.cpp
index 3dbf25359..d49de1435 100644
--- a/Source/ERF.cpp
+++ b/Source/ERF.cpp
@@ -793,6 +793,7 @@ ERF::init_only (int lev, Real time)
         // The base state is initialized from WRF wrfinput data, output by
         // ideal.exe or real.exe
         init_from_wrfinput(lev);
+        if (init_type == "ideal") initHSE();
 
     } else if (init_type == "metgrid") {
         // The base state is initialized from data output by WPS metgrid;
@@ -802,13 +803,13 @@ ERF::init_only (int lev, Real time)
     } else if (init_type == "uniform") {
         // Initialize a uniform background field and base state based on the
         // problem-specified reference density and temperature
-        initHSE(lev);
         init_uniform(lev);
+        initHSE(lev);
     } else {
         // No background flow initialization specified, initialize the
         // background field to be equal to the base state, calculated from the
         // problem-specific erf_init_dens_hse
-        initHSE(lev);
+        initHSE(lev); // need to call this first
         init_from_hse(lev);
     }