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 #include "phycon.h"
00012
00013 void CoolVana(void)
00014 {
00015 realnum a21,
00016 a31,
00017 a32,
00018 g1,
00019 g2,
00020 g3,
00021 p2,
00022 p3,
00023 cs;
00024
00025 DEBUG_ENTRY( "CoolVana()" );
00026
00027
00028
00029
00030 a21 = 0.05f;
00031 a31 = 0.10f;
00032 a32 = 0.00;
00033 g1 = 28.;
00034 g2 = 12.;
00035 g3 = 18.;
00036
00037
00038 p3 = (realnum)atom_pop3(g1,g2,g3,g1,g2,g3,a21,a31,a32,16303.,606.,&p2,
00039 dense.xIonDense[ipVANADIUM][2], 0.,0.,0.);
00040
00041 CoolHeavy.V38830 = p2*a21*2.25e-12;
00042 CoolHeavy.V38507 = p3*a31*2.34e-12;
00043 CoolAdd("V 3",8823,CoolHeavy.V38830);
00044 CoolAdd("V 3",8507,CoolHeavy.V38507);
00045
00046
00047 a21 = 0.054f;
00048 a31 = 0.039f;
00049 a32 = 0.007f;
00050 g1 = 21.;
00051 g2 = 5.;
00052 g3 = 9.;
00053
00054
00055 p3 = (realnum)atom_pop3(g1,g2,g3,g1,g2,g3,a21,a31,a32,15159.,3437.,&p2,
00056 dense.xIonDense[ipVANADIUM][3], 0.,0.,0.);
00057
00058 CoolHeavy.V47741 = p3*a31*2.57e-12;
00059
00060 CoolHeavy.V49496 = p2*a21*2.09e-12;
00061
00062 CoolHeavy.V44p2m = p3*a32*4.74e-13;
00063 CoolAdd("V 4",7735,CoolHeavy.V47741);
00064 CoolAdd("V 4",9489,CoolHeavy.V49496);
00065 CoolAdd("V 4",42,CoolHeavy.V44p2m);
00066
00067
00068
00069
00070 PutCS(2.39,&TauLines[ipVa07130]);
00071 atom_level2(&TauLines[ipVa07130]);
00072
00073
00074
00075
00076
00077 if(phycon.te < 3.566E6)
00078 {
00079 cs = (realnum)(0.0149*phycon.te10*phycon.te05*phycon.te004*phycon.te0003);
00080 }
00081 else
00082 {
00083 cs = (realnum)(47.350653/((phycon.te40/phycon.te02)*phycon.te0002));
00084 }
00085
00086
00087 PutCS(cs,&TauLines[ipVa15172]);
00088 atom_level2(&TauLines[ipVa15172]);
00089 return;
00090 }