00001
00002
00003
00004 #include "cddefines.h"
00005 #include "embesq.h"
00006 #include "taulines.h"
00007 #include "phycon.h"
00008 #include "coolheavy.h"
00009 #include "dense.h"
00010 #include "ligbar.h"
00011 #include "lines_service.h"
00012 #include "atoms.h"
00013 #include "cooling.h"
00014
00015 void CoolAlum(void)
00016 {
00017 double cs ,
00018 cs2s2p,
00019 cs2s3p;
00020 realnum p2;
00021
00022 DEBUG_ENTRY( "CoolAlum()" );
00023
00024
00025 MakeCS( TauLines[ipAlI3957] );
00026 atom_level2( TauLines[ipAlI3957]);
00027
00028
00029 MakeCS(TauLines[ipAlI3090]);
00030 atom_level2(TauLines[ipAlI3090]);
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045 PutCS(3.56,TauLines[ipT2670]);
00046
00047
00048 AtomSeqBeryllium(1.67,2.00,6.54,TauLines[ipT2670],3.67e-3);
00049 embesq.em2669 = (realnum)(atoms.PopLevels[3]*3.67e-3*7.45e-12);
00050
00051
00052
00053
00054
00055
00056 cs = 4.407*phycon.te10*phycon.te03*phycon.te01;
00057 cs = MIN2(25.0,cs);
00058 PutCS(cs*0.667,TauLines[ipT1855]);
00059 PutCS(cs*0.333,TauLines[ipT1863]);
00060 PutCS(1.0,*TauDummy);
00061 atom_level3(TauLines[ipT1863],*TauDummy,TauLines[ipT1855]);
00062
00063
00064
00065
00066
00067
00068 if(phycon.te < 1.58E5)
00069 {
00070 cs = (realnum)(0.893/(phycon.te05*phycon.te005*phycon.te001*phycon.te0002));
00071 }
00072 else
00073 {
00074 cs = (realnum)(3.1991/((phycon.te20/phycon.te04)*(phycon.te003/phycon.te0002)));
00075 }
00076
00077 PutCS(cs,TauLines[ipAl529]);
00078
00079 atom_level2(TauLines[ipAl529]);
00080
00081
00082 cs = 639.1/(phycon.sqrte*pow(phycon.te03,phycon.te003)*phycon.te001);
00083 cs = MIN2(5.5 , cs);
00084 PutCS(cs,TauLines[ipAl6366]);
00085
00086 cs = MIN2(1.10,49.37/(phycon.sqrte/phycon.te10*phycon.te02/
00087 phycon.te001));
00088 PutCS(cs,TauLines[ipAl6912]);
00089
00090 cs = MIN2(2.0,319.11/(phycon.sqrte*phycon.te10/phycon.te02/
00091 phycon.te001));
00092 PutCS(cs,*TauDummy);
00093
00094 atom_level3(TauLines[ipAl6366],TauLines[ipAl6912],*TauDummy);
00095
00096
00097
00098
00099
00100
00101 CoolHeavy.c1170 = atom_pop3(9.,5.,1.,1.044,0.145,0.463,6.63,72.9,7.79,
00102 5.92e4,6.767e4,&p2,dense.xIonDense[ipALUMINIUM][5],0.,0.,0.)*72.9*1.70e-11;
00103
00104 CoolHeavy.c2428 = p2*6.63*8.19e-12;
00105 CoolHeavy.c2125 = CoolHeavy.c1170*(7.79/72.9)*(1169.5/2124.9);
00106 CoolAdd("Al 6",1170,CoolHeavy.c1170);
00107 CoolAdd("Al 6",2428,CoolHeavy.c2428);
00108 CoolAdd("Al 6",2125,CoolHeavy.c2125);
00109
00110
00111
00112
00113 cs = MIN2(0.39,0.0459*phycon.te20/phycon.te003/phycon.te003);
00114 PutCS(cs,TauLines[ipAl8575]);
00115 cs = MIN2(1.062,0.0407*phycon.te30/phycon.te003/phycon.te003);
00116 PutCS(cs,TauLines[ipAl8370]);
00117 cs = MIN2(0.27,2.694e-3*phycon.te20*phycon.te20*phycon.te01*
00118 phycon.te003);
00119 PutCS(cs,*TauDummy);
00120 atom_level3(TauLines[ipAl8575],TauLines[ipAl8370],*TauDummy);
00121
00122
00123 PutCS(1.,TauLines[ipAl09204]);
00124 atom_level2(TauLines[ipAl09204]);
00125
00126
00127
00128
00129 cs = 0.73492 - 0.16964*phycon.alogte + 0.0096631*POW2(phycon.alogte);
00130 cs = MAX2(0.01,cs);
00131 PutCS(cs,TauLines[ipT639]);
00132 atom_level2(TauLines[ipT639]);
00133
00134
00135
00136 ligbar(13,TauLines[ipTAl550],TauLines[ipTAl48],&cs2s2p,&cs2s3p);
00137 PutCS(cs2s2p,TauLines[ipTAl550]);
00138 PutCS(cs2s2p*0.5,TauLines[ipTAl568]);
00139 PutCS(1.0,*TauDummy);
00140 atom_level3(TauLines[ipTAl568],*TauDummy,TauLines[ipTAl550]);
00141
00142 PutCS(cs2s3p,TauLines[ipTAl48]);
00143 atom_level2(TauLines[ipTAl48]);
00144 return;
00145 }