PyAtomDB APEC module¶
This modules contains the APEC code. It calls many different subroutines from throughout the PyAtomDB module. Currently largely unwritten, as APEC code needs to be tidied up for transfer.
The apec module contains routines crucial for the APEC code. This also includes some interfaces to external C libraries (or will, eventually).
Version 0.1 - initial release Adam Foster September 16th 2015
-
pyatomdb.apec.
solve_ionbal
(ionrate, recrate, init_pop=False, tau=False)¶ solve_ionbal: given a set of ionization and recombination rates, find the equilibrium ionization balance. If init_pop and tau are set, do an non-equilibrium calculation starting from init_pop and evolving for n_e * t = tau (cm^-3 s)
Parameters: ionrate : float array
the ionization rates, starting with neutral ionizing to +1
recrate : float array
the recombination rates, starting with singly ionized recombining to neutral
init_pop : float array
initial population of ions for non-equlibrium calculations. Will be renormalised to 1.
tau : float
N_e * t for the non-equilibrium ioniziation
Returns: final_pop : float array
final populations.