From 7dbd18c25d5b8709504c7693d7a62ed5db406742 Mon Sep 17 00:00:00 2001 From: Josh Shields Date: Mon, 14 Oct 2024 11:51:14 -0400 Subject: [PATCH] fill value -> 0 for LinearNDInterpolator --- stardis/radiation_field/opacities/opacities_solvers/util.py | 2 ++ 1 file changed, 2 insertions(+) diff --git a/stardis/radiation_field/opacities/opacities_solvers/util.py b/stardis/radiation_field/opacities/opacities_solvers/util.py index 7e297e80..09ad5c1a 100644 --- a/stardis/radiation_field/opacities/opacities_solvers/util.py +++ b/stardis/radiation_field/opacities/opacities_solvers/util.py @@ -46,6 +46,7 @@ def sigma_file(tracing_lambdas, temperatures, fpath, opacity_source=None): linear_interp_h2plus_bf = LinearNDInterpolator( np.vstack([file_waves_mesh.ravel(), file_temps_mesh.ravel()]).T, file_cross_sections.flatten(), + fill_value=0, ) lambdas, temps = np.meshgrid(tracing_lambdas, temperatures) sigmas = ( @@ -68,6 +69,7 @@ def sigma_file(tracing_lambdas, temperatures, fpath, opacity_source=None): linear_interp_hminus_ff = LinearNDInterpolator( np.vstack([file_waves_mesh.ravel(), file_thetas_mesh.ravel()]).T, file_values.flatten(), + fill_value=0, ) lambdas, thetas = np.meshgrid(tracing_lambdas, 5040 / temperatures) sigmas = (