cloudy  trunk
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
ion_neon.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 /*IonNeon ionization balance for neon */
4 #include "cddefines.h"
5 #include "dense.h"
6 #include "trace.h"
7 #include "ionbal.h"
8 
9 void IonNeon(void)
10 {
11  const int NDIM = ipNEON+1;
12 
13  static const double dicoef[2][NDIM] = {
14  {9.77e-4,2.65e-3,3.69e-3,1.12e-2,2.44e-2,3.02e-2,6.10e-3,2.52e-1,1.44e-1,0.},
15  {.073,.242,1.01,.391,2.52,.445,.254,.304,.296,0.}
16  };
17  static const double dite[2][NDIM] = {
18  {3.11e5,2.84e5,2.24e5,2.7e5,3.09e5,2.83e5,1.68e5,1.4e7,1.5e7,0.},
19  {2.06e5,3.07e5,2.94e5,5.50e5,9.91e5,1.73e6,6.13e5,1.80e6,2.24e6,0.}
20  };
21  static const double ditcrt[NDIM] = {3.0e4,3.3e4,3.3e4,3.5e4,3.6e4,3.6e4,2.9e4,1.5e6,3.8e6,1e20};
22  static const double aa[NDIM] = {0.,0.0129,3.6781,-0.0254,-0.0141,19.9280,5.4751,0.,0.,0.};
23  static const double bb[NDIM] = {0.,-0.1779,14.1481,5.5365,33.8479,235.0536,203.9751,0.,0.,0.};
24  static const double cc[NDIM] = {0.,0.9353,17.1175,17.0727,43.1608,152.5096,86.9016,0.,0.,0.};
25  static const double dd[NDIM] = {0.,-0.0682,-0.5017,-0.7225,-1.6072,9.1413,-7.4568,0.,0.,0.};
26  static const double ff[NDIM] = {0.1,0.4516,0.2313,0.1702,0.1942,0.1282,2.5145,0.,0.,0.};
27 
28  DEBUG_ENTRY( "IonNeon()" );
29 
30  /* neon, atomic number 10 */
31  if( !dense.lgElmtOn[ipNEON] )
32  {
33  return;
34  }
35 
37 
38  ion_photo(ipNEON,false);
39 
40  /* find collisional ionization rates */
42 
43  /* get recombination coefficients */
44  ion_recomb(false,(const double*)dicoef,(const double*)dite,ditcrt,aa,bb,cc,dd,ff,ipNEON);
45 
46  /* solve for ionization balance */
47  ion_solver(ipNEON,false);
48 
49  if( trace.lgTrace && trace.lgHeavyBug )
50  {
51  fprintf( ioQQQ, " IonNeon returns; frac=" );
52  for( int i=0; i < 10; i++ )
53  {
54  fprintf( ioQQQ, "%10.3e", dense.xIonDense[ipNEON][i]/
56  }
57  fprintf( ioQQQ, "\n" );
58  }
59  return;
60 }

Generated for cloudy by doxygen 1.8.3.1