Merge branch 'oksanaguba/eamxx/fixer-clean' (PR #7553)

Adds (dycore) energy fixer. The infrastructure can be used for another process (and will be used for IEFLX),
but the current impl is for the dycore.
It is identical to the EAM impl, where pressure adjustment and dycore energy leaks are fixed as one
global change to temperature for each column/level as the same temperature tendency.

There is a check if fixer debug output is on ( enable_energy_fixer_debug_info ) ,
which computes global energy after the fixer to confirm
that it does match. This debug output is used in a test that has postproc *py script to
confirm that the energy is fixed up to a tolerance.

For full discussion see #7276 .

Climo results ( model vs model, with fixer and without from the same branch)
https://web.lcrc.anl.gov/public/e3sm/diagnostic_output/ac.onguba/theta/
fixerJuly24-2025/latlon-pdf-vs-default/viewer/lat_lon/index.html .

[nonBFB] for any cime runs with eamxx and dycore. Potentially climate-changing, too.

Author: oksanaguba

components/eamxx/src/share/property_checks/mass_and_energy_conservation_check.cpp

@@ -304,6 +304,7 @@ PropertyCheck::ResultAndMsg MassAndEnergyConservationCheck::check() const
 
 void MassAndEnergyConservationCheck::global_fixer(const bool & print_debug_info)
 {
+  using TPF = ekat::TeamPolicyFactory<DefaultDevice::execution_space>;
   const auto ncols = m_num_cols;
   const auto nlevs = m_num_levs;
 
@@ -328,7 +329,7 @@ void MassAndEnergyConservationCheck::global_fixer(const bool & print_debug_info)
   const auto ice_flux   = m_fields.at("ice_flux"  ).get_view<const Real*>();
   const auto heat_flux  = m_fields.at("heat_flux" ).get_view<const Real*>();
 
-  const auto policy = ExeSpaceUtils::get_default_team_policy(ncols, nlevs);
+  const auto policy = TPF::get_default_team_policy(ncols, nlevs);
 
   using namespace ekat::units;
   using namespace ShortFieldTagsNames;
@@ -485,10 +486,11 @@ compute_gas_mass_on_column (const KT::MemberType&       team,
                             const int                   nlevs,
                             const uview_1d<const Real>& pseudo_density)
 {
+  using RU = ekat::ReductionUtils<DefaultDevice::execution_space>;
   using PC = scream::physics::Constants<Real>;
   const Real gravit = PC::gravit;
 
-  return ExeSpaceUtils::parallel_reduce<Real>(team, 0, nlevs,
+  return RU::parallel_reduce<Real>(team, 0, nlevs,
                                               [&] (const int lev, Real& local_mass) {
     local_mass += pseudo_density(lev)/gravit;
   });

components/eamxx/src/share/property_checks/mass_and_energy_conservation_check.hpp

@@ -6,8 +6,8 @@
 #include "share/field/field.hpp"
 #include "share/field/field_utils.hpp"
 
-#include "ekat/kokkos/ekat_kokkos_utils.hpp"
-#include "ekat/mpi/ekat_comm.hpp"
+#include <ekat_team_policy_utils.hpp>
+#include "ekat_comm.hpp"
 
 namespace scream {