diff --git a/Source/reactions/Castro_react.cpp b/Source/reactions/Castro_react.cpp
index 7d005e726c..8d39095ff3 100644
--- a/Source/reactions/Castro_react.cpp
+++ b/Source/reactions/Castro_react.cpp
@@ -208,9 +208,8 @@ Castro::react_state(MultiFab& s, MultiFab& r, Real time, Real dt, const int stra
 
         const auto dx = geom.CellSizeArray();
 #ifdef MODEL_PARSER
-        const auto problo = geom.ProbLoArray();
-#endif
         const auto geomdata = geom.data();
+#endif
 
 #if defined(AMREX_USE_GPU)
         ParallelFor(bx, [=] AMREX_GPU_DEVICE (int i, int j, int k)
diff --git a/Source/scf/scf_relax.cpp b/Source/scf/scf_relax.cpp
index 7b1899ee13..674a509f89 100644
--- a/Source/scf/scf_relax.cpp
+++ b/Source/scf/scf_relax.cpp
@@ -243,11 +243,10 @@ Castro::do_hscf_solve()
                 [=] AMREX_GPU_DEVICE (int i, int j, int k) -> ReduceTuple
                 {
                     const auto *dx = geomdata.CellSize();
-                    const auto *problo = geomdata.ProbLo();
 
                     // The below assumes we are rotating on the z-axis.
 
-                    GpuArray<Real, 3>r;
+                    GpuArray<Real, 3>r = {0.0};
                     position(i, j, k, geomdata, r);
 
                     // Do a trilinear interpolation to find the contribution from
@@ -369,7 +368,6 @@ Castro::do_hscf_solve()
                 [=] AMREX_GPU_DEVICE (int i, int j, int k) -> ReduceTuple
                 {
                     const auto *dx = geomdata.CellSize();
-                    const auto *problo = geomdata.ProbLo();
 
                     // The below assumes we are rotating on the z-axis.
 
@@ -441,9 +439,6 @@ Castro::do_hscf_solve()
                     // enthalpy + gravitational potential + rotational potential = const
                     // We already have the constant, so our goal is to construct the enthalpy field.
 
-                    const auto *dx = geomdata.CellSize();
-                    const auto *problo = geomdata.ProbLo();
-
                     GpuArray<Real, 3> r = {0.0};
                     position(i, j, k, geomdata, r);
 
@@ -618,8 +613,6 @@ Castro::do_hscf_solve()
                 {
                     Real dM = 0.0, dK = 0.0, dU = 0.0, dE = 0.0;
 
-                    const auto* problo = geomdata.ProbLo();
-
                     GpuArray<Real, 3> r = {0.0};
                     position(i, j, k, geomdata, r);