Fixing initial state values in MPhys

adflow/mphys/mphys_adflow.py

@@ -439,6 +439,23 @@
     def _set_states(self, outputs):
         self.solver.setStates(outputs["adflow_states"])
 
+    def _setup_vectors(self, root_vectors):
+        super()._setup_vectors(root_vectors)
+
+        # Turn off extra printouts
+        self.solver.setOption("printIterations", False)
+        self.solver.adflow.inputiteration.printbcwarnings = False
+
+        # Set the AP
+        self.solver.setAeroProblem(self.ap)
+
+        # Turn printing back on
+        self.solver.setOption("printIterations", True)
+        self.solver.adflow.inputiteration.printbcwarnings = True
+
+        # 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
@@ -450,6 +467,7 @@
     def solve_nonlinear(self, inputs, outputs):
         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?

adflow/pyADflow.py

@@ -93,6 +93,12 @@ def __init__(self, comm=None, options=None, debug=False, dtype="d"):
             curDir = os.path.basename(os.path.dirname(os.path.realpath(__file__)))
             self.adflow = MExt.MExt("libadflow", curDir, debug=debug)._module
 
+        # This is a flag not exposed to the user as an option that controls
+        # printing the boundary condition warnings.  This should only
+        # be used in the MPhys layer when we make multiple calls to
+        # updateBCAllLevels.
+        self.adflow.inputiteration.printbcwarnings = True
+
         libLoadTime = time.time()
 
         # Information for base class: