36 fprintf(
ioQQQ,
" RT_tau_inc called.\n" );
56 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
60 int ion = nelem+1-ipISO;
69 ipISO , nelem , -1 , ipHi );
71 for( ipLo=0; ipLo < ipHi; ipLo++ )
76 if(
Transitions[ipISO][nelem][ipHi][ipLo].ipCont <= 0 )
79 save =
Transitions[ipISO][nelem][ipHi][ipLo].Emis->PopOpc;
80 Transitions[ipISO][nelem][ipHi][ipLo].Emis->PopOpc *= factor;
83 ipISO , nelem , ipHi , ipLo );
85 Transitions[ipISO][nelem][ipHi][ipLo].Emis->PopOpc = save;
98 -1 ,ipISO , nelem , ipHi );
130 for( i=0; i <
nUTA; i++ )
136 UTALines[i].Lo->Pop =
dense.
xIonDense[UTALines[i].Hi->nelem-1][UTALines[i].Hi->IonStg-1];
137 UTALines[i].Hi->Pop = 0.;
149 if( save<=0. )
continue;
157 dense.
xIonDense[HFLines[i].Hi->nelem-1][HFLines[i].Hi->IonStg-1] = save;
179 -10 , -10 , -10 , i );
191 fprintf(
ioQQQ,
" RT_tau_inc updated optical depths:\n" );
196 fprintf(
ioQQQ,
" RT_tau_inc returns.\n" );