fix: ensure that net_field_totals reducers are created even when hydr…

…o=false
lanl · Dec 2, 2023 · f29e0d2 · f29e0d2
1 parent 5621d23
commit f29e0d2
Showing 1 changed file with 5 additions and 5 deletions.
diff --git a/src/fluid/fluid.cpp b/src/fluid/fluid.cpp
@@ -233,6 +233,11 @@ std::shared_ptr<StateDescriptor> Initialize(ParameterInput *pin) {
     physics->AddField(c::ye, mcons_scalar);
   }
 
+  AllReduce<std::vector<Real>> net_field_totals;
+  AllReduce<std::vector<Real>> net_field_totals_2;
+  physics->AddParam<>("net_field_totals", net_field_totals, true);
+  physics->AddParam<>("net_field_totals_2", net_field_totals_2, true);
+
   // If fluid is not active, still don't add reconstruction variables
   if (!active) return physics;
 
@@ -318,11 +323,6 @@ std::shared_ptr<StateDescriptor> Initialize(ParameterInput *pin) {
         HstSum, ReduceMom, "total X" + std::to_string(d + 1) + " momentum"));
   }
 
-  AllReduce<std::vector<Real>> net_field_totals;
-  AllReduce<std::vector<Real>> net_field_totals_2;
-  physics->AddParam<>("net_field_totals", net_field_totals, true);
-  physics->AddParam<>("net_field_totals_2", net_field_totals_2, true);
-
   params.Add(parthenon::hist_param_key, hst_vars);
 
   // Fill Derived and Estimate Timestep