cloudy  trunk
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
ion_lithi.cpp
Go to the documentation of this file.
1 /* This file is part of Cloudy and is copyright (C)1978-2008 by Gary J. Ferland and
2  * others. For conditions of distribution and use see copyright notice in license.txt */
3 /*IonLithi compute ionization balance for lithium */
4 #include "cddefines.h"
5 #include "dense.h"
6 #include "ionbal.h"
7 
8 void IonLithi(void)
9 {
10  const int NDIM = ipLITHIUM+1;
11 
12  static const double dicoef[2][NDIM] = { {2.54e-3,6.15e-3,0.}, {4.42e-2,5.88e-2,0.} };
13  static const double dite[2][NDIM] = { {1.57e5,1.41e5,0.}, {3.74e5,1.41e5,0.} };
14  static const double ditcrt[NDIM] = {1.2e4,1.2e4,1e20};
15  static const double aa[NDIM] = {0.,0.,0.};
16  static const double bb[NDIM] = {0.,0.,0.};
17  static const double cc[NDIM] = {0.,0.,0.};
18  static const double dd[NDIM] = {0.,0.,0.};
19  static const double ff[NDIM] = {0.1,0.1,0.};
20 
21  DEBUG_ENTRY( "IonLithi()" );
22 
23  /* lithium nelem=3
24  * data are for carbon
25  *
26  * fluorescent fields, first dim is stage of ionization, sec is shell
27  *
28  * rates from Shull and van Steenberg, Ap.J. Sup 48, 95. */
29  /* DATA GRDEFF/0.10,0.10,0.10/
30  * GRDEFF is fraction of recombinations to ground state, used for
31  * outward diffuse fields
32  *
33  * rec from +3, +4 from Arnaud et al Ast Ap Sup 60 425. (1985)
34  * rec from fully ionized uses Seaton '79 in ionrat */
35  /* Pequignot and Aldrovandi Ast Ap 161, 169. */
36 
37  if( !dense.lgElmtOn[ipLITHIUM] )
38  {
39  return;
40  }
41 
42  /* zero out ionization balance arrays */
44 
45  ion_photo(ipLITHIUM,false);
46 
47  /* find collisional ionization rates */
49 
50  /* get recombination coefficients */
51  ion_recomb(false,(const double*)dicoef,(const double*)dite,ditcrt,aa,bb,cc,dd,ff,ipLITHIUM);
52 
53  /* solve for ionization balance */
54  ion_solver(ipLITHIUM,false);
55  return;
56 }

Generated for cloudy by doxygen 1.8.1.1