minor bug fix to add scope for geometry

E3SM-Project · Mar 26, 2024 · 1e61c90 · 1e61c90
1 parent 9d692dc
commit 1e61c90
Showing 1 changed file with 11 additions and 5 deletions.
diff --git a/dynamics/spam/src/models/extrudedmodel.h b/dynamics/spam/src/models/extrudedmodel.h
@@ -248,6 +248,8 @@ class ModelTendencies : public ExtrudedTendencies {
                                 FieldSet<nprognostic> &x) {
     real max_w = params.max_w;
 
+    YAKL_SCOPE(primal_geometry, this->primal_geometry);
+
     const auto &primal_topology = primal_geometry.topology;
 
     const auto wvar = x.fields_arr[WVAR].data;
@@ -271,9 +273,12 @@ class ModelTendencies : public ExtrudedTendencies {
 
   void adjust_crm_per_phys_using_vert_cfl(ModelParameters &params,
                                           FieldSet<nprognostic> &x) {
+
     // cfl = c * dt / dx
     // dt = cfl * dx / c
 
+    YAKL_SCOPE(primal_geometry, this->primal_geometry);
+
     const auto &primal_topology = primal_geometry.topology;
     const auto wvar = x.fields_arr[WVAR].data;
     const int pis = primal_topology.is;
@@ -302,11 +307,12 @@ class ModelTendencies : public ExtrudedTendencies {
     params.crm_per_phys = (int)std::ceil(params.dt_crm_phys / adj_min_dt);
     params.dtcrm = params.dt_crm_phys / params.crm_per_phys;
 
-    //  std::cout << "min_dt:         " << min_dt << "\n";
-    //  std::cout << "adj_min_dt:         " << adj_min_dt << "\n";
-    //  std::cout << "dtcrm:         " << params.dtcrm << "\n";
-    //  std::cout << "dt_crm_phys:         " << params.dt_crm_phys << "\n";
-    //  std::cout << "crm per phys:     " << params.crm_per_phys << "\n";
+    // std::cout << "adjust_crm_per_phys_using_vert_cfl: min_dt:         " << min_dt << "\n";
+    // std::cout << "adjust_crm_per_phys_using_vert_cfl: adj_min_dt:         " << adj_min_dt << "\n";
+    // std::cout << "adjust_crm_per_phys_using_vert_cfl: dtcrm:         " << params.dtcrm << "\n";
+    // std::cout << "adjust_crm_per_phys_using_vert_cfl: dt_crm_phys:         " << params.dt_crm_phys << "\n";
+    // std::cout << "adjust_crm_per_phys_using_vert_cfl: crm per phys:     " << params.crm_per_phys << "\n";
+    // yakl::fence();
   }
 
   void convert_dynamics_to_coupler_state(PamCoupler &coupler,