38 fprintf(
ioQQQ,
" RT_tau_inc called.\n" );
43 conv.lgLastSweepThisZone =
true;
53 for(
long i=0; i<
rfield.nfine; ++i )
56 rfield.fine_opt_depth[i] += tauzone;
58 rfield.trans_coef_total_stale =
true;
88 for(
long nelem=ipISO; nelem <
LIMELM; nelem++ )
92 int ion = nelem+1-ipISO;
94 if( ion <=
dense.IonHigh[nelem] &&
dense.xIonDense[nelem][ion] >
dense.density_low_limit )
99 for( ipLo=0; ipLo <
iso_sp[ipISO][nelem].numLevels_local; ipLo++ )
106 for( ipHi=1; ipHi <
iso_sp[ipISO][nelem].numLevels_local; ipHi++ )
108 for( ipLo=0; ipLo < ipHi; ipLo++ )
110 if(
iso_sp[ipISO][nelem].trans(ipHi,ipLo).ipCont() <= 0 )
120 for( ipHi=
iso_sp[ipISO][nelem].st[
iso_sp[ipISO][nelem].numLevels_local-1].n()+1; ipHi <
iso_ctrl.nLyman[ipISO]; ipHi++ )
123 (*tr).Emis().PopOpc() =
iso_sp[ipISO][nelem].st[0].Pop();
153 for( i=0; i <
nUTA; i++ )
169 if(
save<=0. )
continue;
185 (*diatom)->H2_RT_tau_inc();
188 for(
long ipSpecies=0; ipSpecies<
nSpecies; ipSpecies++ )
196 int ipHi = (*tr).ipHi();
197 if (ipHi >=
dBaseSpecies[ipSpecies].numLevels_local || (*tr).ipCont() <= 0)
199 int ipLo = (*tr).ipLo();
207 if(
wind.lgStatic() )
215 fprintf(
ioQQQ,
" RT_tau_inc updated optical depths:\n" );
220 fprintf(
ioQQQ,
" RT_tau_inc returns.\n" );
void RT_line_one_tauinc(const TransitionProxy &t, long int mas_species, long int mas_ion, long int mas_hi, long int mas_lo, realnum DopplerWidth)