diff --git a/adflow/mphys/mphys_adflow.py b/adflow/mphys/mphys_adflow.py
index 6eb6a532f..2c0ee6101 100644
--- a/adflow/mphys/mphys_adflow.py
+++ b/adflow/mphys/mphys_adflow.py
@@ -559,6 +559,11 @@ def solve_nonlinear(self, inputs, outputs):
                         print("###############################################################")
                         print("# Solve Failed, not attempting a clean restart")
                         print("###############################################################")
+
+                    # write the solution so that we can diagnose
+                    solver.writeSolution(baseName="analysis_fail", number=self.solution_counter)
+                    self.solution_counter += 1
+
                     raise AnalysisError("ADFLOW Solver Fatal Fail")
 
                 else: