00001
00002
00003
00004 #include "cddefines.h"
00005 #include "taulines.h"
00006 #include "coolheavy.h"
00007 #include "lines_service.h"
00008 #include "dense.h"
00009 #include "atoms.h"
00010 #include "cooling.h"
00011 #include "phycon.h"
00012
00013
00014 void CoolChro(void)
00015 {
00016 double a21,
00017 a31,
00018 a32;
00019 realnum p2,
00020 p3;
00021
00022 DEBUG_ENTRY( "CoolChro()" );
00023
00024
00025
00026
00027
00028 CoolHeavy.Cr3l21 = atom_pop2(25.,25.,9.,0.05,2.47e4,dense.xIonDense[ipCHROMIUM][2])*
00029 3.41e-12;
00030 CoolAdd("Cr 3",5828,CoolHeavy.Cr3l21);
00031
00032
00033
00034
00035 a21 = 0.053;
00036 a31 = 0.102;
00037 a32 = 0.00;
00038
00039
00040 p3 = (realnum)atom_pop3(28.,12.,18.,28.,12.,18.,a21,a31,a32,19795.,1356.,&p2,
00041 dense.xIonDense[ipCHROMIUM][3],0.,0.,0.);
00042
00043 CoolHeavy.Cr4l31 = p3*a31*2.92e-12;
00044
00045 CoolHeavy.Cr4l21 = p2*a21*2.74e-12;
00046 CoolAdd("Cr 4",6801,CoolHeavy.Cr4l31);
00047 CoolAdd("Cr 4",7267,CoolHeavy.Cr4l21);
00048
00049
00050
00051
00052 a21 = 0.157;
00053 a31 = 0.048;
00054 a32 = 0.016;
00055
00056
00057 p3 = (realnum)atom_pop3(21.,5.,9.,21.,5.,9.,a21,a31,a32,18028.,3842.,&p2,dense.xIonDense[ipCHROMIUM][4],
00058 0.,0.,0.);
00059
00060 CoolHeavy.Cr5l31 = p3*a31*3.02e-12;
00061
00062 CoolHeavy.Cr5l21 = p2*a21*2.49e-12;
00063
00064 CoolHeavy.Cr5l32 = p2*a32*5.31e-13;
00065 CoolAdd("Cr 5",6577,CoolHeavy.Cr5l31);
00066 CoolAdd("Cr 5",37,CoolHeavy.Cr5l32);
00067 CoolAdd("Cr 5",7979,CoolHeavy.Cr5l21);
00068
00069 return;
00070 }