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.