Skip to content

Commit

Permalink
Update fiasco/ions.py
Browse files Browse the repository at this point in the history
Co-authored-by: Will Barnes <will.t.barnes@gmail.com>
  • Loading branch information
jwreep and wtbarnes authored Aug 9, 2024
1 parent 68ddf50 commit ed6252a
Showing 1 changed file with 0 additions and 3 deletions.
3 changes: 0 additions & 3 deletions fiasco/ions.py
Original file line number Diff line number Diff line change
Expand Up @@ -1584,13 +1584,10 @@ def free_bound_radiative_loss(self) -> u.erg * u.cm**3 / u.s:
E_obs = recombined._fblvl['E_obs']*const.h*const.c
E_th = recombined._fblvl['E_th']*const.h*const.c
E_fb = np.where(E_obs==0*u.erg, E_th, E_obs)

wvl_n0 = const.h * const.c / (recombined.ip - E_fb[0])
wvl_n1 = (recombined._fblvl['n'][0] + 1)**2 /(const.Ryd * z**2)

g_fb0 = self._gaunt_factor_free_bound_total(ground_state=True)
g_fb1 = self._gaunt_factor_free_bound_total(ground_state=False)

term1 = g_fb0 * np.exp(-const.h*const.c/(const.k_B * self.temperature * wvl_n0))
term2 = g_fb1 * np.exp(-const.h*const.c/(const.k_B * self.temperature * wvl_n1))

Expand Down

0 comments on commit ed6252a

Please sign in to comment.