Made reduceIO write states at writeTime only. (#579)

mdolab · Jan 30, 2024 · 69e5bc8 · 69e5bc8
1 parent 0f42fe5
commit 69e5bc8
Showing 1 changed file with 34 additions and 30 deletions.
diff --git a/src/adjoint/DASolver/DASolver.C b/src/adjoint/DASolver/DASolver.C
@@ -9118,44 +9118,48 @@ void DASolver::writeAdjStates(const label writeMesh)
         Write only the adjoint states
     */
 
-    // volVector states
-    forAll(stateInfo_["volVectorStates"], idxI)
+    if (runTimePtr_->writeTime())
     {
-        const word stateName = stateInfo_["volVectorStates"][idxI];
-        volVectorField& state =
-            const_cast<volVectorField&>(meshPtr_->thisDb().lookupObject<volVectorField>(stateName));
 
-        state.write();
-    }
+        // volVector states
+        forAll(stateInfo_["volVectorStates"], idxI)
+        {
+            const word stateName = stateInfo_["volVectorStates"][idxI];
+            volVectorField& state =
+                const_cast<volVectorField&>(meshPtr_->thisDb().lookupObject<volVectorField>(stateName));
 
-    // volScalar states
-    forAll(stateInfo_["volScalarStates"], idxI)
-    {
-        const word stateName = stateInfo_["volScalarStates"][idxI];
-        volScalarField& state =
-            const_cast<volScalarField&>(meshPtr_->thisDb().lookupObject<volScalarField>(stateName));
+            state.write();
+        }
 
-        state.write();
-    }
+        // volScalar states
+        forAll(stateInfo_["volScalarStates"], idxI)
+        {
+            const word stateName = stateInfo_["volScalarStates"][idxI];
+            volScalarField& state =
+                const_cast<volScalarField&>(meshPtr_->thisDb().lookupObject<volScalarField>(stateName));
 
-    // model states
-    forAll(stateInfo_["modelStates"], idxI)
-    {
-        const word stateName = stateInfo_["modelStates"][idxI];
-        volScalarField& state =
-            const_cast<volScalarField&>(meshPtr_->thisDb().lookupObject<volScalarField>(stateName));
+            state.write();
+        }
 
-        state.write();
-    }
+        // model states
+        forAll(stateInfo_["modelStates"], idxI)
+        {
+            const word stateName = stateInfo_["modelStates"][idxI];
+            volScalarField& state =
+                const_cast<volScalarField&>(meshPtr_->thisDb().lookupObject<volScalarField>(stateName));
 
-    // surfaceScalar states
-    forAll(stateInfo_["surfaceScalarStates"], idxI)
-    {
-        const word stateName = stateInfo_["surfaceScalarStates"][idxI];
-        surfaceScalarField& state =
-            const_cast<surfaceScalarField&>(meshPtr_->thisDb().lookupObject<surfaceScalarField>(stateName));
+            state.write();
+        }
+
+        // surfaceScalar states
+        forAll(stateInfo_["surfaceScalarStates"], idxI)
+        {
+            const word stateName = stateInfo_["surfaceScalarStates"][idxI];
+            surfaceScalarField& state =
+                const_cast<surfaceScalarField&>(meshPtr_->thisDb().lookupObject<surfaceScalarField>(stateName));
 
-        state.write();
+            state.write();
+        }
     }
 
     scalar endTime = runTimePtr_->endTime().value();