00001
00002
00003
00004 #include "cddefines.h"
00005 #include "taulines.h"
00006 #include "phycon.h"
00007 #include "lines_service.h"
00008 #include "dense.h"
00009 #include "atoms.h"
00010 #include "cooling.h"
00011
00012 void CoolNick(void)
00013 {
00014 double tused;
00015 realnum rate;
00016
00017 DEBUG_ENTRY( "CoolNick()" );
00018
00019
00020
00021 tused = MAX2( 30. , phycon.te );
00022 tused = MIN2( 3000. , phycon.te );
00023 tused /= 100.;
00024
00025
00026
00027 rate = (realnum)(1.2e-7 * dense.eden +
00028
00029
00030 8.0e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0] );
00031 LineConvRate2CS( TauLines[ipNi1_7m] , rate );
00032
00033
00034 rate = (realnum)(9.3e-8 * dense.eden +
00035
00036
00037 5.3e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0] );
00038 LineConvRate2CS( TauLines[ipNi1_11m] , rate );
00039
00040 rate = (realnum)(1.2e-7 * dense.eden +
00041
00042
00043 6.9e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0] );
00044 (*(*TauDummy).Hi()).g() = (*TauLines[ipNi1_11m].Hi()).g();
00045 LineConvRate2CS( *TauDummy , rate );
00046
00047 (*(*TauDummy).Hi()).g() = 0.;
00048
00049 atom_level3(TauLines[ipNi1_7m],TauLines[ipNi1_11m],*TauDummy);
00050
00051 return;
00052 }