54 for( i=0; i < 6; i++ )
74 " P8446 called for first time, finds A*escape prob 3-1 =%10.2e\n",
90 fprintf(
ioQQQ,
" P8446 tin, tout=%10.2e%10.2e\n", tin, tout );
115 " P8446 finds Lbeta, OI widths=%10.2e%10.2e and esc prob=%10.2e%10.2e esAB=%10.2e\n",
124 xoi = opacoi/(opacoi + opaclb);
125 xlb = opaclb/(opacoi + opaclb);
135 " P8446 opac Lb, OI=%10.2e%10.2e X Lb, OI=%10.2e%10.2e FLb, OI=%10.2e%10.2e\n",
136 opaclb, opacoi, xlb, xoi, flb, foi );
172 fprintf(
ioQQQ,
" P8446 PMPH31=%10.2e EscP3-1%10.2e ESCH31=%10.2e\n",
282 static realnum g[6]={9.,3.,9.,3.,15.,9};
329 c12 = c21*g[1]/g[0]*e12;
330 c13 = c31*g[2]/g[0]*e12*e23;
331 c14 = c41*g[3]/g[0]*e12*e23*e34;
332 c15 = c51*g[4]/g[0]*e12*e23*e34*e45;
333 c16 = c61*g[5]/g[0]*e12*e23*e34*e45*e56;
340 c23 = c32*g[2]/g[1]*e23;
341 c24 = c42*g[3]/g[1]*e23*e34;
342 c25 = c52*g[4]/g[1]*e23*e34*e45;
343 c26 = c62*g[5]/g[1]*e23*e34*e45*e56;
349 c34 = c43*g[3]/g[2]*e34;
350 c35 = c53*g[4]/g[2]*e34*e45;
351 c36 = c63*g[5]/g[2]*e34*e45*e56;
356 c45 = c54*g[4]/g[3]*e45;
357 c46 = c64*g[5]/g[3]*e45*e56;
360 c56 = c65*g[5]/g[4]*e56;
365 simple = (c16 +
atoms.
pmpo15)/(c61 + c62 + c64 + a65 + a64 + a62);
369 for( i=1; i < 6; i++ )
379 for( i=0; i < 6; i++ )
390 zz[1][1] = c21 + c23 + c24 + c25 + c26 + a21;
391 zz[2][1] = -c32 - a32;
394 zz[5][1] = -c62 - a62;
399 zz[2][2] = c31 + c32 + c34 + c35 + c36 + a32;
400 zz[3][2] = -c43 - a43;
401 zz[4][2] = -c53 - a53;
408 zz[3][3] = c41 + c42 + c43 + c45 + c46 + a41 + a43;
410 zz[5][3] = -c64 - a64;
417 zz[4][4] = c51 + c52 + c53 + c54 + c56 + a51 + a53;
418 zz[5][4] = -c65 - a65;
426 zz[5][5] = c61 + c62 + c63 + c64 + c65 + a65 + a64 + a62;
429 for( j=0; j < 6; j++ )
431 for( i=0; i < 6; i++ )
433 amat[i][j] = zz[i][j];
441 getrs_wrapper(
'N', 6, 1, (
double*)amat, 6, ipiv, bvec, 6, &ner);
447 fprintf(
ioQQQ,
" oi_level_pops: dgetrs finds singular or ill-conditioned matrix\n" );
453 for( i=0; i < 6; i++ )
459 for( i=0; i < 6; i++ )
480 e12*e23*e34*e45*e56);
482 fprintf(
ioQQQ,
" oilevl finds levl pop" );
485 fprintf(
ioQQQ,
"\n" );
487 fprintf(
ioQQQ,
" oilevl finds dep coef" );
489 fprintf(
ioQQQ,
"%11.3e", deptoi[i] );
490 fprintf(
ioQQQ,
"\n" );
498 fprintf(
ioQQQ,
" OILEVL finds negative population" );
501 fprintf(
ioQQQ,
"\n" );
503 fprintf(
ioQQQ,
" simple 5 =%10.2e\n", simple );