00001
00002
00003
00004 #include "cddefines.h"
00005 #include "taulines.h"
00006 #include "coolheavy.h"
00007 #include "dense.h"
00008 #include "lines_service.h"
00009 #include "atoms.h"
00010 #include "cooling.h"
00011
00012 void CoolScan(void)
00013 {
00014 double a21,
00015 a31,
00016 a32,
00017 g1,
00018 g2,
00019 g3,
00020 p3;
00021 realnum p2;
00022
00023 DEBUG_ENTRY( "CoolScan()" );
00024
00025
00026
00027
00028
00029
00030 a21 = 0.004;
00031 a31 = 6.95e-3;
00032 a32 = 5.0e-4;
00033 g1 = 15.;
00034 g2 = 5.;
00035 g3 = 21.;
00036
00037
00038 p3 = atom_pop3(g1,g2,g3,g1,g2,g3,a21,a31,a32,3504.,3407.,&p2,
00039 dense.xIonDense[ipSCANDIUM][1], 0.,0.,0.);
00040
00041
00042 CoolHeavy.Sc22p08m = p3*a31*9.56e-13;
00043
00044
00045 CoolHeavy.Sc24p1m = p2*a21*4.85e-13;
00046
00047
00048 CoolHeavy.Sc24p2m = p3*a32*4.71e-13;
00049
00050 CoolAdd("Sc 2",41,CoolHeavy.Sc24p1m);
00051 CoolAdd("Sc 2",42,CoolHeavy.Sc24p2m);
00052 CoolAdd("Sc 2",21,CoolHeavy.Sc22p08m);
00053
00054
00055
00056
00057 g1 = 10.;
00058 g2 = 2.;
00059 CoolHeavy.Sc33936 = atom_pop2(g1,g1,g2,0.03,3.66e4,
00060 dense.xIonDense[ipSCANDIUM][2])* 5.06e-12;
00061 CoolAdd("Sc 3",3933,CoolHeavy.Sc33936);
00062
00063
00064
00065
00066 PutCS(6.00,&TauLines[ipSc05231]);
00067 atom_level2(&TauLines[ipSc05231]);
00068
00069
00070 a21 = 4.94;
00071 a31 = 49.24;
00072 a32 = 4.3;
00073 g1 = 9.;
00074 g2 = 5.;
00075 g3 = 1.;
00076
00077
00078 p3 = atom_pop3(g1,g2,g3,g1,g2,g3,a21,a31,a32,28464.,40045.,&p2,
00079 dense.xIonDense[ipSCANDIUM][5], 0.,0.,0.);
00080
00081 CoolHeavy.Sc42100 = p3*a31*9.47e-12;
00082
00083 CoolHeavy.Sc45058 = p2*a21*3.93e-12;
00084
00085 CoolHeavy.Sc43595 = p3*a32*5.53e-12;
00086 CoolAdd("Sc 6",2100,CoolHeavy.Sc42100);
00087 CoolAdd("Sc 6",5054,CoolHeavy.Sc45058);
00088 CoolAdd("Sc 6",3592,CoolHeavy.Sc43595);
00089
00090
00091
00092 PutCS(0.182,&TauLines[ipSc13264]);
00093 atom_level2(&TauLines[ipSc13264]);
00094
00095 return;
00096 }