00001
00002
00003
00004 #include "cddefines.h"
00005 #include "trace.h"
00006 #include "dense.h"
00007 #include "ionbal.h"
00008 #include "mole.h"
00009 #include "oxy.h"
00010 #include "atmdat.h"
00011
00012 void IonNelem(
00013
00014 bool lgPrintIt,
00015
00016 long int nelem)
00017 {
00018 DEBUG_ENTRY( "IonNelem()" );
00019
00020 if( dense.lgElmtOn[nelem] )
00021 {
00022 ion_zero(nelem);
00023
00024 ion_photo(nelem,lgPrintIt);
00025
00026
00027 ion_collis(nelem);
00028
00029
00030 ion_CX( nelem );
00031
00032
00033 ion_recomb(lgPrintIt,nelem);
00034
00035
00036 ion_solver(nelem,lgPrintIt);
00037
00038 if( trace.lgTrace && trace.lgHeavyBug )
00039 {
00040 fprintf( ioQQQ, " IonNelem nelem\t%li\tfnzone\t%6.2f\tfrac\t",
00041 nelem, fnzone );
00042 for( int i=0; i < nelem+2; i++ )
00043 {
00044 fprintf( ioQQQ, "\t%10.3e", dense.xIonDense[nelem][i]/
00045 dense.gas_phase[nelem] );
00046 }
00047 fprintf( ioQQQ, "\n" );
00048 }
00049 }
00050
00051 return;
00052 }