sirocco-rt / sirocco

This is the repository for Sirocco, the radiative transfer code used to model winds in AGN and other systems
GNU General Public License v3.0
30 stars 24 forks source link

LTE_Tr (and LTE_Te) does not calcuate ionization for macro atoms in LTE #1103

Open kslong opened 1 month ago

kslong commented 1 month ago

LTE is used to initialize a Python cycle, but in LTE mode even in subsequent cycles macro atoms are ionized based on estimators due to photon passages. As a result, macro atoms are not held in LTE for the LTE_Tr and LTE_te modes. This should not be the case. The LTE modes should force all ions, macro or otherwise, to LTE based on Tr or Te.

(Note that it is unclear to me how some of the other "simple" modes operate)>