00001
00002
00003
00004 #include "cddefines.h"
00005 #include "dense.h"
00006 #include "ionbal.h"
00007
00008 void IonLithi(void)
00009 {
00010 const int NDIM = ipLITHIUM+1;
00011
00012 static const double dicoef[2][NDIM] = { {2.54e-3,6.15e-3,0.}, {4.42e-2,5.88e-2,0.} };
00013 static const double dite[2][NDIM] = { {1.57e5,1.41e5,0.}, {3.74e5,1.41e5,0.} };
00014 static const double ditcrt[NDIM] = {1.2e4,1.2e4,1e20};
00015 static const double aa[NDIM] = {0.,0.,0.};
00016 static const double bb[NDIM] = {0.,0.,0.};
00017 static const double cc[NDIM] = {0.,0.,0.};
00018 static const double dd[NDIM] = {0.,0.,0.};
00019 static const double ff[NDIM] = {0.1,0.1,0.};
00020
00021 DEBUG_ENTRY( "IonLithi()" );
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037 if( !dense.lgElmtOn[ipLITHIUM] )
00038 {
00039 return;
00040 }
00041
00042
00043 ion_zero(ipLITHIUM);
00044
00045 ion_photo(ipLITHIUM,false);
00046
00047
00048 ion_collis(ipLITHIUM);
00049
00050
00051 ion_recomb(false,(const double*)dicoef,(const double*)dite,ditcrt,aa,bb,cc,dd,ff,ipLITHIUM);
00052
00053
00054 ion_solver(ipLITHIUM,false);
00055 return;
00056 }