cloudy  trunk
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
ion_titan.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 /*IonTitan determine ionization balance for titanium */
4 #include "cddefines.h"
5 #include "dense.h"
6 #include "ionbal.h"
7 
8 void IonTitan(void)
9 {
10  const int NDIM = ipTITANIUM+1;
11 
12  static const double dicoef[2][NDIM] = {
13  {1.58e-3,8.38e-3,1.54e-2,3.75e-2,0.117,0.254,0.291,0.150,0.140,0.100,
14  0.200,0.240,0.260,0.190,0.120,0.350,0.066,0.10,0.13,0.23,0.14,0.},
15  {.456,.323,.310,.411,.359,.0975,.229,4.20,3.30,5.30,1.50,0.700,.600,
16  .5,1.,0.,7.8,6.3,5.5,3.6,4.9,0.}
17  };
18  static const double dite[2][NDIM] = {
19  {6.00e4,1.94e5,3.31e5,4.32e5,6.28e5,7.50e5,7.73e5,2.62e5,2.50e5,2.57e5,2.84e5,
20  8.69e5,4.21e5,4.57e5,2.85e5,8.18e5,1.51e6,1.30e6,1.19e6,1.09e6,9.62e5,0.},
21  {8.97e4,1.71e5,2.73e5,3.49e5,5.29e5,4.69e5,6.54e5,1.32e6,1.33e6,1.41e6,1.52e6,
22  1.51e6,1.82e6,1.84e6,2.31e6,0.,9.98e6,9.98e6,1.00e7,1.10e7,8.34e6,0.}
23  };
24  static const double ditcrt[NDIM] = {6e3,2e4,4e4,5e4,7e4,8e4,8e4,3e4,
25  3e4,3e4,3e4,9e4,4e4,5e4,3e4,9e5,2e5,2e5,2e5,2e5,1e5,1e20};
26  static const double aa[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
27  0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
28  static const double bb[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
29  0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
30  static const double cc[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
31  0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
32  static const double dd[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
33  0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
34  static const double ff[NDIM] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
35  0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
36 
37  DEBUG_ENTRY( "IonTitan()" );
38 
39  /* titanium nelem=22
40  * based on iron
41  *
42  * rates from Shull and van Steenberg, Ap.J. Sup 48, 95. */
43 
44  /* Fe rates from Woods et al. Ap.J. 249, 399.
45  * rec from +23, 24 25 from Arnauld et al 85 */
46 
47  /* Pequignot and Aldrovandi Ast Ap 161, 169. */
48 
49  if( !dense.lgElmtOn[ipTITANIUM] )
50  {
51  return;
52  }
53 
55 
56  ion_photo(ipTITANIUM,false);
57 
58  /* find collisional ionization rates */
60 
61  /* get recombination coefficients */
62  ion_recomb(false,(const double*)dicoef,(const double*)dite,ditcrt,aa,bb,cc,dd,ff,ipTITANIUM);
63 
64  /* solve for ionization balance */
65  ion_solver(ipTITANIUM,false);
66  return;
67 }

Generated for cloudy by doxygen 1.8.4