/* ------- file: -------------------------- atom.h ------------------------- */


#ifndef __ATOM_H__
#define __ATOM_H__

/* Defines atomic data for H and He
 *
 * Sara Seager
 * Last modified: Feb. 27, 1999
 */

#define Lam2s1sH   8.22458  /* H 2s-1s two photon rate in [s^-1] */
#define Lam2s1sHe  51.3  /* HeI 2s-1s two photon rate in [s^-1] */

#define LyalphaH   1.215670e-07  /* H Lyman alpha wavelength in [m] */
#define LyalphaHe   5.843344e-08 /* HeI 2^1p - 1^1s wavelength in [m] */

#define EionH2s  5.446605e-19 /* H 2s ionization energy in [J] */
#define EionHe2s 6.363254e-19 /* HeI 2s ionization energy in [J] */

#define E2s1sH  1.63403067e-18 /* H 2s energy from the ground state [J]     */
#define E2s1sHe 3.30301387e-18 /* HeI 2s energy from the ground state in [J] */

#define E2p2sHe 9.64908313e-20 /* HeI 2p - 2s energy difference   in  [J] */

/* --- Associated function prototypes -------------------------------------- */

double SahaBoltz(double gi, double gj, double ne, double E_ion, double T);
double Boltzmann(double gi, double gj, double E, double T);
void evaluate(int choice, double z, double y[], double fvec[], double dfdz[], double **fjac, int n);

#endif /* !__ATOM_H__ */

/* ------- end ---------------------------- atom.h ------------------------- */