only compute lya damping above zmin.

bd-j · Oct 8, 2024 · 2f0588b · 2f0588b
1 parent e6af39a
commit 2f0588b
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/prospect/models/sedmodel.py b/prospect/models/sedmodel.py
@@ -797,9 +797,10 @@ def add_dla(self, wave_rest, spec):
         return spec
 
     def add_damping_wing(self, wave_rest, spec):
-        if np.any(self.params.get("igm_damping", False)):
+        zmin = 5.0
+        if (self._zred > zmin) & np.any(self.params.get("igm_damping", False)):
             x_HI = self.params.get("igm_factor", 1.0)
-            tau = tau_damping(wave_rest, self._zred, x_HI, cosmo=cosmo)
+            tau = tau_damping(wave_rest, self._zred, x_HI, zmin=zmin, cosmo=cosmo)
             spec *= np.exp(-tau)
         return spec