23 static double aNI[
N1_SIZE] = {2.755e-5,4.123e-5,7.536e-6,1.486e-5,4.516e-5,-2.935e-6,4.000e-6,3.751e-6,-2.176e-6,1.024e-5};
24 static double bNI[
N1_SIZE] = {-8.150e-8,-1.220e-7,-2.226e-8,-4.390e-8,-1.130e-7,8.000e-9,-1.1447e-8,-1.061e-8,5.610e-9,-3.227e-8};
25 static double cNI[
N1_SIZE] = {2.140e-4,3.272e-4,-4.944e-6,3.473e-6,-8.772e-4,1.654e-3,1.675e-3,1.123e-3,3.867e-3,3.376e-4};
27 static double NI[6][6][3];
38 double temp_max = 50e3;
59 NI[i][j][0] =
aNI[index];
60 NI[i][j][1] =
bNI[index];
61 NI[i][j][2] =
cNI[index];
72 printf(
"NI %i%i a:%e ", i, j,
NI[i][j][0]);
73 printf(
"%i%i b:%e\n", i, j,
NI[i][j][1]);
87 else if(init < 1 || init > 5 || final < 1 || final > 5)
94 temp =
MAX2(temp, temp_min);
95 temp =
MIN2(temp, temp_max);
142 static double gAr4[5]={4.,6.,4.,2.,4.};
143 static double exAr4[4]={19224.464,8.713,9605.74,0.386};
145 static long int *ipN1Pump=NULL,
146 nN1Pump=0 , lgFirst_N1=
true;
148 static bool lgFirst_N3=
true;
149 static long int *ipN3Pump=NULL,
156 double pop1, pop2,pop3,pop4,pop5;
157 ofp = fopen(
"c:\\projects\\cloudy\\trunk\\mydata\\pop.out",
"w");
160 printf(
"Can't open the file\n");
207 else if( nN1Pump > 0 )
209 ipN1Pump = (
long *)
MALLOC((
unsigned)(nN1Pump)*
sizeof(
long) );
219 ipN1Pump[nN1Pump] = i;
232 for( i=0; i<nN1Pump; ++i )
236 fprintf(
ioQQQ,
"DEBUG N %li %.3e %.3e\n",
310 atom_pop5(gAr4,exAr4,cs12,cs13,cs14,cs15,cs23,cs24,cs25,cs34,cs35,
311 cs45,a12,a13,a14,a15,a23,a24,a25,a34,a35,a45,pop,
dense.
xIonDense[
ipNITROGEN][0]);
319 fprintf(ofp,
"%lf %lf %lf %lf %lf\n", pop1, pop2, pop3, pop4, pop5);
332 p3 =
atom_pop3(4.,10.,6.,cs21,cs31,cs32,a21,a31,a32,2.769e4,1.38e4,
375 # define USE_LENNON_BURKE 1
376 # if( USE_LENNON_BURKE )
389 p3 =
atom_pop3( 9.,5.,1., cs21, cs31, cs32, a21, a31, a32 ,
407 # if( USE_LENNON_BURKE )
430 # if( USE_LENNON_BURKE )
469 else if( nN3Pump > 0 )
471 ipN3Pump = (
long *)
MALLOC((
unsigned)(nN3Pump)*
sizeof(
long) );
481 ipN3Pump[nN3Pump] = i;
492 for( i=0; i<nN3Pump; ++i )
496 fprintf(
ioQQQ,
"DEBUG C %li %.3e %.3e\n",
511 0.201 , 1.088 , 0.668 , 2.044 , pump_rate,
"N 3");