Fixes for ALM restart

Exawind · Dec 10, 2024 · b2efa43 · b2efa43
1 parent 6acdd39
commit b2efa43
Show file tree

Hide file tree

Showing 3 changed files with 22 additions and 0 deletions.
diff --git a/amr-wind/wind_energy/actuator/Actuator.cpp b/amr-wind/wind_energy/actuator/Actuator.cpp
@@ -121,6 +121,16 @@ void Actuator::post_regrid_actions()
 {
     BL_PROFILE("amr-wind::actuator::Actuator::post_regrid_actions");
 
+    amrex::Vector<int> act_proc_count(amrex::ParallelDescriptor::NProcs(), 0);
+    for (auto& act : m_actuators) {
+        auto& info = act->info();
+        if (info.root_proc == -1) {
+            info.procs =
+                utils::determine_influenced_procs(m_sim.mesh(), info.bound_box);
+            utils::determine_root_proc(act->info(), act_proc_count);
+        }
+    }
+
     for (auto& act : m_actuators) {
         act->determine_influenced_procs();
     }

diff --git a/amr-wind/wind_energy/actuator/ActuatorModel.H b/amr-wind/wind_energy/actuator/ActuatorModel.H
@@ -39,6 +39,8 @@ public:
 
     virtual const ActInfo& info() const = 0;
 
+    virtual ActInfo& info() = 0;
+
     virtual void read_inputs(const utils::ActParser&) = 0;
 
     virtual void determine_influenced_procs() = 0;
@@ -116,6 +118,9 @@ public:
     //! Return the core info object for this actuator instance
     const ActInfo& info() const override { return m_data.info(); }
 
+    //! Return the core info object for this actuator instance
+    ActInfo& info() override { return m_data.info(); }
+
     //! Return the meta info object for this actuator instance
     const auto& meta() const { return m_data.meta(); }
 

diff --git a/amr-wind/wind_energy/actuator/turbine/fast/FastIface.cpp b/amr-wind/wind_energy/actuator/turbine/fast/FastIface.cpp
@@ -202,6 +202,7 @@ void FastIface::init_turbine(const int local_id)
 
     case SimMode::restart: {
         fast_restart_turbine(fi);
+        prepare_netcdf_file(fi);
         break;
     }
     }
@@ -367,6 +368,12 @@ void FastIface::prepare_netcdf_file(FastTurbine& fi)
     }
 
     const std::string fname = m_output_dir + "/" + fi.tlabel + ".nc";
+
+    // Don't overwrite existing
+    if (amrex::FileSystem::Exists(fname)) {
+        return;
+    }
+
     auto ncf = ncutils::NCFile::create(fname, NC_CLOBBER | NC_NETCDF4);
     const std::string nt_name = "num_time_steps";
     const std::string np_name = "num_vel_points";