provide depth point indices in warning message

jvshields · Oct 14, 2024 · 39d13e8 · 39d13e8
1 parent 97acfad
commit 39d13e8
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/stardis/radiation_field/opacities/opacities_solvers/util.py b/stardis/radiation_field/opacities/opacities_solvers/util.py
@@ -58,7 +58,7 @@ def sigma_file(tracing_lambdas, temperatures, fpath, opacity_source=None):
         )  # Scaling from Stancil 1994 table
         if np.any(sigmas == 0):
             logger.warning(
-                "Outside of interpolation range for H2+ BF cross-sections in some part of the atmosphere. Assuming 0 opacity from H2+ BF for these points."
+                f"Outside of interpolation range for H2+ BF cross-sections at depth points {np.unique(np.where(sigmas == 0)[0])}. Assuming 0 opacity from H-FF for these depth points."
             )
     elif (
         opacity_source == "Hminus_ff"
@@ -87,7 +87,7 @@ def sigma_file(tracing_lambdas, temperatures, fpath, opacity_source=None):
         )
         if np.any(sigmas == 0):
             logger.warning(
-                "Outside of interpolation range for H- FF cross-sections in some part of the atmosphere. Assuming 0 opacity from H-FF for these points."
+                f"Outside of interpolation range for H- FF cross-sections at depth points {np.unique(np.where(sigmas == 0)[0])}. Assuming 0 opacity from H-FF for these depth points."
             )
 
     elif (