Skip to content

Commit

Permalink
scaling comment
Browse files Browse the repository at this point in the history
  • Loading branch information
jvshields committed Dec 11, 2024
1 parent d4d96d9 commit a733cbc
Showing 1 changed file with 5 additions and 4 deletions.
9 changes: 5 additions & 4 deletions stardis/radiation_field/opacities/opacities_solvers/base.py
Original file line number Diff line number Diff line change
Expand Up @@ -557,7 +557,8 @@ def calc_alan_entries(
# We want to consider grid points within a certain range of the line_nu
line_broadening = (
((line_gamma + doppler_width) * alpha) / d_nu * 20
) # Scale by alpha of the line
) # Scale by alpha of the line, though this would be good to revisit and come up with a better scaling.
# This produces the same answer as before, but really just makes large lines search very very far away.
line_broadening_range = max(10.0, line_broadening) # Force a minimum range

lower_freq_index = max(
Expand Down Expand Up @@ -699,9 +700,9 @@ def calc_alphas(
stellar_radiation_field.frequencies,
opacity_config.line,
)
stellar_radiation_field.opacities.opacities_dict["alpha_line_at_nu"] = (
alpha_line_at_nu
)
stellar_radiation_field.opacities.opacities_dict[
"alpha_line_at_nu"
] = alpha_line_at_nu
stellar_radiation_field.opacities.opacities_dict["alpha_line_at_nu_gammas"] = gammas
stellar_radiation_field.opacities.opacities_dict[
"alpha_line_at_nu_doppler_widths"
Expand Down

0 comments on commit a733cbc

Please sign in to comment.