58 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
60 if( (
dense.IonHigh[nelem] >= nelem + 1-ipISO) )
62 for( ipHi=1; ipHi <
iso_sp[ipISO][nelem].numLevels_local; ipHi++ )
65 for( ipLo=0; ipLo < ipHi - 1; ipLo++ )
68 if(
iso_sp[ipISO][nelem].trans(ipHi,ipLo).ipCont() > 0 )
70 OneLine =
iso_sp[ipISO][nelem].trans(ipHi,ipLo).Emis().pump()*
71 iso_sp[ipISO][nelem].trans(ipHi,ipLo).EnergyErg()*
72 iso_sp[ipISO][nelem].trans(ipHi,ipLo).Emis().PopOpc();
74 accel_iso[ipISO] += OneLine;
82 for( ipLo=0; ipLo <
iso_sp[ipISO][nelem].numLevels_local; ipLo++ )
86 if((*tr).ipCont() > 0 )
88 OneLine = (*tr).Emis().pump()*
90 (*tr).Emis().PopOpc();
92 accel_iso[ipISO] += OneLine;
97 for( ipHi=
iso_sp[ipISO][nelem].st[
iso_sp[ipISO][nelem].numLevels_local-1].n()+1; ipHi <
iso_ctrl.nLyman[ipISO]; ipHi++ )
101 if( (*tr).ipCont() > 0 )
103 OneLine = (*tr).Emis().pump()*
105 (*tr).Emis().PopOpc();
107 accel_iso[ipISO] += OneLine;
138 for( i=0; i <
nUTA; i++ )
155 for(
long ipSpecies=0; ipSpecies<
nSpecies; ipSpecies++ )
162 int ipHi = (*tr).ipHi();
163 if (ipHi >=
dBaseSpecies[ipSpecies].numLevels_local || (*tr).ipCont() <= 0)
165 OneLine = (*tr).EnergyErg()*(*tr).Emis().pump()*(*tr).Emis().PopOpc();
174 h2drive += (*diatom)->H2_Accel();
181 fe2drive + h2drive + AllRest;