4198 raise exceptions on negative fuel ion charge density and temperature profile elements

process/models/physics/bootstrap_current.py

@@ -1470,6 +1470,10 @@ def bootstrap_fraction_sauter(self, plasma_profile: float) -> float:
 
         # Calculate electron and ion temperature profiles
         tempe = plasma_profile.teprofile.profile_y
+        # Check for any negative temperature in profile: always fatal eventually,
+        # report explicitly at source
+        if (tempe < 0).any():
+            raise ValueError("Negative temperature in plasma profile")
         tempi = (
             physics_variables.temp_plasma_ion_vol_avg_kev
             / physics_variables.temp_plasma_electron_vol_avg_kev

process/models/physics/physics.py

@@ -1161,6 +1161,12 @@ def plasma_composition():
             - znimp
         )
 
+        # Negative znfuel can occur during solution and is always fatal,
+        # frequently resulting in a confusing bootstrap current error later.
+        # Catch early and explicitly instead
+        if znfuel < 0.0:
+            raise ValueError(f"znfuel is negative: {znfuel}")
+
         # ======================================================================
 
         # Fuel ion density, nd_plasma_fuel_ions_vol_avg