diff --git a/amr-wind/wind_energy/ABLStats.cpp b/amr-wind/wind_energy/ABLStats.cpp
index e9e3a1149b..26cf92cc81 100644
--- a/amr-wind/wind_energy/ABLStats.cpp
+++ b/amr-wind/wind_energy/ABLStats.cpp
@@ -55,11 +55,15 @@ void ABLStats::initialize()
         pp.query("normal_direction", m_normal_dir);
         AMREX_ASSERT((0 <= m_normal_dir) && (m_normal_dir < AMREX_SPACEDIM));
         pp.query("kappa", m_kappa);
+        pp.get("reference_temperature", m_ref_theta);
+        pp.query("stats_do_energy_budget", m_do_energy_budget);
+    }
+
+    {
+        amrex::ParmParse pp("incflo");
         amrex::Vector<amrex::Real> gravity{0.0, 0.0, -9.81};
         pp.queryarr("gravity", gravity);
         m_gravity = utils::vec_mag(gravity.data());
-        pp.get("reference_temperature", m_ref_theta);
-        pp.query("stats_do_energy_budget", m_do_energy_budget);
     }
 
     // Get normal direction and associated stuff
diff --git a/amr-wind/wind_energy/ABLWallFunction.cpp b/amr-wind/wind_energy/ABLWallFunction.cpp
index 9675145cae..5e7a123c6b 100644
--- a/amr-wind/wind_energy/ABLWallFunction.cpp
+++ b/amr-wind/wind_energy/ABLWallFunction.cpp
@@ -16,6 +16,11 @@ namespace amr_wind {
 ABLWallFunction::ABLWallFunction(const CFDSim& sim)
     : m_sim(sim), m_mesh(sim.mesh())
 {
+    {
+        amrex::ParmParse pp("incflo");
+        pp.queryarr("gravity", m_gravity);
+    }
+
     amrex::ParmParse pp("ABL");
 
     pp.query("kappa", m_mo.kappa);
@@ -25,7 +30,6 @@ ABLWallFunction::ABLWallFunction(const CFDSim& sim)
     pp.query("mo_beta_h", m_mo.beta_h);
     pp.query("surface_roughness_z0", m_mo.z0);
     pp.query("normal_direction", m_direction);
-    pp.queryarr("gravity", m_gravity);
     AMREX_ASSERT((0 <= m_direction) && (m_direction < AMREX_SPACEDIM));
 
     if (pp.contains("log_law_height")) {