Handle ionpot_ev > Emin

lukeshingles · Nov 6, 2024 · d347457 · d347457
1 parent f121a5a
commit d347457
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/pynonthermal/spencerfano.py b/pynonthermal/spencerfano.py
@@ -360,8 +360,8 @@ def add_ionisation(self, Z, ion_stage, n_ion):
         dfcollion_thision = self.dfcollion.query("Z == @Z and ion_stage == @ion_stage", inplace=False)
 
         for index, shell in dfcollion_thision.iterrows():
-            assert shell.ionpot_ev >= self.engrid[0]
-            self._add_ionisation_shell(n_ion, shell)
+            if shell.ionpot_ev >= self.engrid[0]:
+                self._add_ionisation_shell(n_ion, shell)
 
     def calculate_n_e(self):
         # number density of free electrons [cm-^3]