Fixing initial state values in MPhys

adflow/mphys/mphys_adflow.py

@@ -468,6 +468,28 @@
     def _set_states(self, outputs):
         self.solver.setStates(outputs["adflow_states"])
 
+    def _setup_vectors(self, root_vectors):
+        super()._setup_vectors(root_vectors)
+        solver = self.solver
+
+        # Get the default flags
+        printIterationsDefault = solver.getOption("printIterations")
+        printBCWarningsDefault = solver.adflow.inputiteration.printbcwarnings
+
+        # Turn off extra printouts
+        solver.setOption("printIterations", False)
+        solver.adflow.inputiteration.printbcwarnings = False
+
+        # Set the AP
+        self.solver.setAeroProblem(self.ap)
+
+        # Set values back to the default
+        solver.setOption("printIterations", printIterationsDefault)
+        solver.adflow.inputiteration.printbcwarnings = printBCWarningsDefault
+
+        # Set the OM states to be equal to the free stream
+        self.set_val("adflow_states", self.solver.getStates())
+
     def apply_nonlinear(self, inputs, outputs, residuals):
         solver = self.solver
         ap = self.ap
@@ -481,6 +503,7 @@
             print_func_call(self)
         solver = self.solver
         ap = self.ap
+
         if self._do_solve:
             setAeroProblem(solver, ap, self.ap_vars, inputs=inputs, outputs=outputs, print_dict=False)
             ap.solveFailed = False  # might need to clear this out?
@@ -1029,8 +1052,11 @@
         printIterationsDefault = solver.getOption("printIterations")
         printBCWarningsDefault = solver.adflow.inputiteration.printbcwarnings
 
+        # Turn off extra printouts
         solver.setOption("printIterations", False)
-        solver.adflow.inputiteration.printbcwarnings = False  # Turn off extra printouts
+        solver.adflow.inputiteration.printbcwarnings = False
+
+        # Set the aeroproblem
         solver.setAeroProblem(ap)
 
         # Reset back to true to preserve normal ADflow printout structure