johnson-c / ColRadPy

Collisional-radiative python solver
Other
27 stars 10 forks source link

Use cubic interpolation of adf11 rates for solving ionization balances #25

Closed dmkriete closed 5 months ago

dmkriete commented 5 months ago

This switches the interpolation of adf11 rates from linear to cubic (the interpolation is still done in logarithmic space).

For adf11 rates this appears to give better results. See comparison below; the linear interpolation exhibits unphysical "bumps" in the curves that disappear when cubic interpolation is used. (technically the below comparison is technically energy balance, but ionization balance is a component of energy balance, so comparison is applicable to both).

Linear interpolation (old) c2+ thermalization with protium for Ti=1 0Te_old

Cubic interpolation (new) c2+ thermalization with protium for Ti=1 0Te_new