cloudy  trunk
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
abund.cpp
Go to the documentation of this file.
1 /* This file is part of Cloudy and is copyright (C)1978-2022 by Gary J. Ferland and
2  * others. For conditions of distribution and use see copyright notice in license.txt */
3 #include "cddefines.h"
4 #include "abund.h"
6 
8 {
9  DEBUG_ENTRY( "t_abund::zero()" );
10  for( long nelem=0; nelem < LIMELM; nelem++ )
11  {
12  /* these are depletion scale factors */
13  depset[nelem] = 1.;
14  /*begin sanity check */
15  if( depset[nelem] == 0. )
16  {
17  fprintf( ioQQQ, " ZERO finds insane abundance or depletion.\n" );
18  fprintf( ioQQQ, " atomic number=%6ld abundance=%10.2e depletion=%10.2e\n",
19  nelem, solar[nelem], depset[nelem] );
20  ShowMe();
22  }
23  /*end sanity check */
24  }
25 
26  /* typical ISM depletion factors, subjective mean of Cowie and Songaila
27  * and Jenkins
28  * */
29  Depletion[0] = 1.;
30  Depletion[1] = 1.;
31  Depletion[2] = .16f;
32  Depletion[3] = .6f;
33  Depletion[4] = .13f;
34  Depletion[5] = 0.4f;
35  Depletion[6] = 1.0f;
36  Depletion[7] = 0.6f;
37  Depletion[8] = .3f;
38  Depletion[9] = 1.f;
39  Depletion[10] = 0.2f;
40  Depletion[11] = 0.2f;
41  Depletion[12] = 0.01f;
42  Depletion[13] = 0.03f;
43  Depletion[14] = .25f;
44  Depletion[15] = 1.0f;
45  Depletion[16] = 0.4f;
46  Depletion[17] = 1.0f;
47  Depletion[18] = .3f;
48  Depletion[19] = 1e-4f;
49  Depletion[20] = 5e-3f;
50  Depletion[21] = 8e-3f;
51  Depletion[22] = 6e-3f;
52  Depletion[23] = 6e-3f;
53  Depletion[24] = 5e-2f;
54  Depletion[25] = 0.01f;
55  Depletion[26] = 0.01f;
56  Depletion[27] = 0.01f;
57  Depletion[28] = .1f;
58  Depletion[29] = .25f;
59 
60  lgDepln = false;
61  ScaleMetals = 1.;
62 
63 }
realnum depset[LIMELM]
Definition: abund.h:245
realnum Depletion[LIMELM]
Definition: abund.h:245
FILE * ioQQQ
Definition: cddefines.cpp:7
t_abund abund
Definition: abund.cpp:5
realnum ScaleMetals
Definition: abund.h:254
#define EXIT_FAILURE
Definition: cddefines.h:168
#define cdEXIT(FAIL)
Definition: cddefines.h:482
void zero()
Definition: abund.cpp:7
const int LIMELM
Definition: cddefines.h:308
#define DEBUG_ENTRY(funcname)
Definition: cddefines.h:723
bool lgDepln
Definition: abund.h:251
int fprintf(const Output &stream, const char *format,...)
Definition: service.cpp:1121
Definition: abund.h:189
void ShowMe(void)
Definition: service.cpp:205
realnum solar[LIMELM]
Definition: abund.h:210