38 static bool lgFinitePop[
LIMELM];
39 static bool lgMustInit[
NISO]={
true,
true};
42 double solomon_assumed;
43 double *PumpSave=NULL;
47 if( lgMustInit[ipISO] )
50 lgFinitePop[nelem] =
true;
87 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
154 lgFinitePop[nelem] =
true;
157 else if( lgFinitePop[nelem] )
161 lgFinitePop[nelem] =
false;
171 for( ipLo=0; ipLo < ipHi; ipLo++ )
177 Transitions[ipISO][nelem][ipHi][ipLo].Emis->PopOpc = 0.;
200 lgMustInit[ipISO] =
false;
221 lgH_chem_conv =
false;
223 while( loop_H_chem < 5 && !lgH_chem_conv )
235 lgH_chem_conv =
true;
243 lgH_chem_conv =
false;
254 enum {DEBUG_LOC=
false};
258 fprintf(
ioQQQ,
"DEBUG \t%.2f\t%.3e\t%.3e\t%.3e\t%.3e\t%.3e\t%.3e\t%.3e\n",
313 enum {DEBUG_LOC=
false};
315 if(DEBUG_LOC && (fabs(error) > 1e-4) )
400 fprintf(
ioQQQ,
" Hydrogenic return %.2f ",
fnzone);
406 fprintf(
ioQQQ,
" REC, COL, GAMT= ");
409 fprintf(
ioQQQ,
"%.2e ", coltot);
410 fprintf(
ioQQQ,
"%.2e ", gamtot);
411 fprintf(
ioQQQ,
" CSUP=");
413 fprintf(
ioQQQ,
"\n");
435 double sum_atom_iso , renorm;
479 for( ipLo=ipH1s; ipLo < ipHi; ipLo++ )
495 long int in, il, is, i, ipLo, nResolved, ipFirstCollapsed=LONG_MIN;
496 char chPrtType[2][12]={
"populations",
"departure"};
498 char chSpin[3][9]= {
"singlets",
"doublets",
"triplets"};
500 #define ITEM_TO_PRINT(A_) ( lgPrtDeparCoef ? iso.DepartCoef[ipISO][nelem][A_] : StatesElem[ipISO][nelem][A_].Pop )
506 for( is = 1; is<=3; ++is)
510 else if( ipISO ==
ipHE_LIKE && is != 1 && is != 3 )
514 nResolved =
StatesElem[ipISO][nelem][ipFirstCollapsed-1].n;
519 fprintf(
ioQQQ,
" %s %s %s %s\n",
523 chPrtType[lgPrtDeparCoef]);
526 fprintf(
ioQQQ,
" n\\l=> ");
527 for( i =0; i < nResolved; ++i)
529 fprintf(
ioQQQ,
"%2ld ",i);
534 for( in = 1; in <= nResolved; ++in)
539 fprintf(
ioQQQ,
" %2ld ",in);
541 for( il = 0; il < in; ++il)
543 if( ipISO==
ipHE_LIKE && (in==2) && (il==1) && (is==3) )
563 fprintf(
ioQQQ,
" %2ld ",in);
578 double TeList[NTE] = {6000.,10000.,15000.,20000.,25000.};
584 fprintf(ioPun,
"Te\t2 3s 33s\n");
589 for( i=0; i<NTE; ++i )
593 fprintf(ioPun ,
"%.0f\t",
595 fprintf(ioPun ,
"%.2f\t",
597 fprintf(ioPun ,
"%.2f\t",
599 fprintf(ioPun ,
"%.2f\t",
601 fprintf(ioPun ,
"%.3f\t",
604 fprintf(ioPun ,
"%.1f\n",