35 limit = nelem + 1 -
NISO;
38 limit =
MAX2( 0 , limit );
42 for( ion=0; ion < limit; ion++ )
44 if(
dense.xIonDense[nelem][ion] > 0. )
50 for( nshell=0; nshell <
Heavy.nsShells[nelem][ion]; nshell++ )
53 if( nshell==
Heavy.nsShells[nelem][ion]-1 )
56 low =
opac.ipElement[nelem][ion][nshell][0];
57 ipHi =
opac.ipElement[nelem][ion][nshell][1];
58 ipop =
opac.ipElement[nelem][ion][nshell][2];
68 for( ion=limit; ion<nelem+1; ++ion )
71 long int ipISO = nelem-ion;
78 if(
dense.xIonDense[nelem][ion] > 0. )
85 abundance =
iso_sp[ipISO][nelem].st[n].Pop();
92 abundance =
dense.xIonDense[nelem][ion];
98 iso_sp[ipISO][nelem].fb[n].ipOpac,
99 iso_sp[ipISO][nelem].fb[n].ipIsoLevNIonCon,
106 iso_sp[ipISO][nelem].fb[n].DepartCoef ,
'v' );
110 if(
iso_sp[ipISO][nelem].st[3].Pop() > 0. )
115 for(
long level =1; level <
iso_sp[ipISO][nelem].numLevels_local; level++ )
117 if( level==
iso_sp[ipISO][nelem].numLevels_max-1 )
120 else if(
iso_sp[ipISO][nelem].st[level].n() >= 5 )
135 iso_sp[ipISO][nelem].fb[level].ipOpac,
136 iso_sp[ipISO][nelem].fb[level].ipIsoLevNIonCon,
139 iso_sp[ipISO][nelem].fb[0].ipIsoLevNIonCon,
140 iso_sp[ipISO][nelem].st[level].Pop(),
142 iso_sp[ipISO][nelem].fb[level].DepartCoef , chType );