28 bool lgUpdateFineOpac )
45 fprintf(
ioQQQ,
" RT_line_all called\n" );
79 ipISO = 0; nelem = 0;ipLo = 0;
81 dTau =
Transitions[ipISO][nelem][ipHi][ipLo].Emis->PopOpc *
82 Transitions[ipISO][nelem][ipHi][ipLo].Emis->opacity /
85 aa = log(1. + dTau ) /
SDIV(dTau);
86 fprintf(
ioQQQ,
"DEBUG dTau\t%.2f\t%.5e\t%.5e\t%.5e\n",
fnzone,dTau,
100 for( nelem=ipISO; nelem<
LIMELM; ++nelem )
138 for( ipLo=0; ipLo < ipHi; ++ipLo )
142 if(
Transitions[ipISO][nelem][ipHi][ipLo].ipCont < 0 )
145 double SavePopOpc =
Transitions[ipISO][nelem][ipHi][ipLo].Emis->PopOpc;
147 Transitions[ipISO][nelem][ipHi][ipLo].Emis->PopOpc *= factor;
162 enum {DEBUG_LOC=
false};
163 if( DEBUG_LOC && nelem==1&& ipLo==0 )
165 fprintf(
ioQQQ,
"DEBUG pdest\t%3li\t%.2f\t%.3e\n",
168 Transitions[ipISO][nelem][ipHi][ipLo].Emis->Pdest);
173 Transitions[ipISO][nelem][ipHi][ipLo].Emis->PopOpc = SavePopOpc;
177 if( ipISO > 0 && lgDoEsc )
196 if(
Transitions[ipISO][nelem][ipHi][ipLo].ipCont <= 0 )
234 for( ipISO=0; ipISO<
NISO; ipISO++ )
236 for( nelem=ipISO; nelem<
LIMELM; nelem++ )
249 for( ipLo=0; ipLo < ipHi; ipLo++ )
251 if(
Transitions[ipISO][nelem][ipHi][ipLo].ipCont <= 0 )
256 if(
Transitions[ipISO][nelem][ipHi][ipLo].Emis->Pesc<1. )
290 fprintf(
ioQQQ,
" RT_line_all calls P8446 who found pmph31=%10.2e\n",
314 for( ipISO=ipH_LIKE; ipISO<
NISO; ++ipISO )
317 for( nelem=ipISO; nelem<
LIMELM; ++nelem )
332 ((
realnum)SaveLyaPdest[ipISO][nelem] +
335 ((
realnum)SaveLyaPesc[ipISO][nelem] +
359 Transitions[ipISO][nelem][ipHi][ipLo].Emis->pump = 0.;
366 enum {DEBUG_LOC=
false};
368 if( DEBUG_LOC &&
nzone>433 )
398 for( i=0; i <
nUTA; i++ )
406 UTALines[i].Lo->Pop =
dense.
xIonDense[UTALines[i].Hi->nelem-1][UTALines[i].Hi->IonStg-1];
407 UTALines[i].Hi->Pop = 0.;
408 RT_line_one(&UTALines[i] , lgDoEsc , lgUpdateFineOpac,
true);
417 for( nelem=ipISO; nelem<
LIMELM; ++nelem )