cloudy  trunk
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
ion_phosi.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 /*IonPhosi derive ionization balance for phosphorus */
4 #include "cddefines.h"
5 #include "dense.h"
6 #include "ionbal.h"
7 
8 void IonPhosi(void)
9 {
10  const int NDIM = ipPHOSPHORUS+1;
11 
12  static const double dicoef[2][NDIM] = {
13  {1.62e-3,1.09e-2,3.35e-2,3.14e-2,1.27e-2,1.47e-2,1.34e-2,2.38e-2,
14  3.19e-2,7.13e-2,.08,.0796,.0134,.402,0.},
15  {0.,1.20e-2,6.59e-2,6.89e-2,.187,.129,1.04,1.12,1.40,1.00,.555,
16  1.63,.304,.298,0.}
17  };
18  static const double dite[2][NDIM] = {
19  {1.25e5,1.92e5,1.89e5,1.68e5,1.38e5,1.80e6,6.90e5,5.84e5,5.17e5,
20  6.66e5,6.00e5,5.09e5,2.91e5,2.41e7,0.},
21  {0.,1.80e4,1.59e5,8.04e4,1.71e5,1.75e6,2.15e6,2.59e6,2.91e6,2.32e6,
22  2.41e6,6.37e6,1.04e6,4.67e6,0.}
23  };
24  static const double aa[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
25  static const double bb[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
26  static const double cc[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
27  static const double dd[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
28  static const double ff[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
29  static const double ditcrt[NDIM] = {2.2e4,1.2e4,1.4e4,1.5e4,1.4e4,2.9e5,
30  1.3e5,1.1e5,9.0e4,9.0e4,9.0e4,8.3e4,6.0e4,5.0e6,1e20};
31 
32  DEBUG_ENTRY( "IonPhosi()" );
33 
34  /* phosphorus nelem=15
35  * data from sulphus
36  *
37  * rec from +13 14 15 from Arnauld et al 85 */
38  /* Pequignot and Aldrovandi Ast Ap 161, 169. */
39 
41  {
42  return;
43  }
44 
46 
47  ion_photo(ipPHOSPHORUS,false);
48 
49  /* find collisional ionization rates */
51 
52  /* get recombination coefficients */
53  ion_recomb(false,(const double*)dicoef,(const double*)dite,ditcrt,aa,bb,cc,dd,ff,ipPHOSPHORUS);
54 
55  /* solve for ionization balance */
56  ion_solver(ipPHOSPHORUS,false);
57  return;
58 }

Generated for cloudy by doxygen 1.8.3.1