00001
00002
00003
00004 #include "cddefines.h"
00005 #include "coolheavy.h"
00006 #include "taulines.h"
00007 #include "dense.h"
00008 #include "phycon.h"
00009 #include "lines_service.h"
00010 #include "atoms.h"
00011 #include "cooling.h"
00012
00013 void CoolSodi(void)
00014 {
00015 double a12,
00016 a13,
00017 a23,
00018 cs,
00019 cs13,
00020 cs23,
00021 p3;
00022 realnum p2;
00023
00024 DEBUG_ENTRY( "CoolSodi()" );
00025
00026
00027
00028
00029 cs = 2.12e-2*(phycon.te70/phycon.te02);
00030 PutCS(cs,TauLines[ipT5895]);
00031 atom_level2(TauLines[ipT5895]);
00032
00033
00034
00035
00036 cs = MIN2(0.40,0.198*phycon.te05*phycon.te01);
00037 cs = MAX2(0.35,cs);
00038 PutCS(cs,TauLines[ipfsNa373]);
00039 atom_level2(TauLines[ipfsNa373]);
00040
00041
00042
00043
00044 PutCS(0.802,TauLines[ipfsNa490]);
00045
00046 PutCS(0.273,TauLines[ipfsNa421]);
00047
00048 PutCS(.111,*TauDummy);
00049 atom_level3(TauLines[ipfsNa490],TauLines[ipfsNa421],*TauDummy);
00050
00051
00052
00053
00054
00055
00056
00057
00058 cs = 0.919;
00059 cs13 = 0.359;
00060 cs23 = 1.915;
00061 a12 = 8.16e-3;
00062 a13 = 3.81;
00063 a23 = 2.68;
00064
00065 p3 = atom_pop3(4.,10.,6.,cs,cs13,cs23,a12,a13,a23,6.96e4,3.58e4,&p2,
00066 dense.xIonDense[ipSODIUM][4],0.,0.,0.);
00067 CoolHeavy.c1365 = p3*a13*1.46e-11;
00068 CoolHeavy.c4017 = p3*a23*4.95e-12;
00069 CoolHeavy.c2067 = p2*a12*9.63e-12;
00070 CoolAdd("Na 5",1365,CoolHeavy.c1365);
00071 CoolAdd("Na 5",4017,CoolHeavy.c4017);
00072 CoolAdd("Na 5",2067,CoolHeavy.c2067);
00073
00074
00075 cs = MIN2(0.77,2.346/(phycon.te10*phycon.te02*phycon.te001));
00076 PutCS(cs,TauLines[ipxNa6143]);
00077 cs = MIN2(2.15,6.934/(phycon.te10*phycon.te03/phycon.te001/
00078 phycon.te001));
00079 PutCS(cs,TauLines[ipxNa6862]);
00080 cs = MIN2(0.53,1.518/(phycon.te10*phycon.te01*phycon.te003*
00081 phycon.te003));
00082 PutCS(cs,*TauDummy);
00083
00084 atom_level3(TauLines[ipxNa6143],TauLines[ipxNa6862],*TauDummy);
00085
00086
00087
00088 cs = MIN2(0.2876,2.603e-3/(phycon.sqrte/phycon.te10));
00089 cs = MAX2(0.1,cs);
00090
00091 CoolHeavy.c2569 = atom_pop3(9.,5.,1.,1.38,0.173,cs,1.68,16.9,5.27,5.01e4,
00092 5.60e4,&p2,dense.xIonDense[ipSODIUM][5],0.,0.,0.)*5.27*6.70e-12;
00093
00094 CoolHeavy.c1357 = CoolHeavy.c2569*(16.9/5.27)*(2972./1357.);
00095 CoolHeavy.c2972 = p2*1.68*6.70e-12;
00096
00097 CoolAdd("Na 6",2569,CoolHeavy.c2569);
00098 CoolAdd("Na 6",1357,CoolHeavy.c1357);
00099 CoolAdd("Na 6",2972,CoolHeavy.c2972);
00100
00101
00102
00103
00104 PutCS(1.5,TauLines[ipxNa0746]);
00105 atom_level2(TauLines[ipxNa0746]);
00106 return;
00107 }