27 const int nLine = 100;
31 long int iperm[
nLine];
44 for(i=0; i<
nLine; ++i)
50 if(
pressure.pres_radiation_lines_curr > 1e-30 )
55 for(
long ipISO = 0; ipISO<
NISO; ipISO++ )
57 for( nelem=ipISO; nelem <
LIMELM; nelem++ )
60 if(
dense.IonHigh[nelem] >= nelem + 1 - ipISO )
64 for( ipHi=1; ipHi <
iso_sp[ipISO][nelem].numLevels_local -
iso_sp[ipISO][nelem].nCollapsed_local; ipHi++ )
66 for( ipLo=0; ipLo < ipHi; ipLo++ )
68 if(
iso_sp[ipISO][nelem].trans(ipHi,ipLo).ipCont() <= 0 )
74 if(
iso_sp[ipISO][nelem].trans(ipHi,ipLo).Emis().PopOpc() > smallfloat &&
76 ( (
iso_sp[ipISO][nelem].trans(ipHi,ipLo).Emis().TauTot() -
iso_sp[ipISO][nelem].trans(ipHi,ipLo).Emis().TauIn()) > smallfloat ) )
82 wl[ip] =
iso_sp[ipISO][nelem].trans(ipHi,ipLo).WLAng();
100 if( (*
TauLines[i].Hi()).Pop() > 1e-30 )
125 if( (*
TauLine2[i].Hi()).Pop() > 1e-30 )
149 if( (*
HFLines[i].Hi()).Pop() > 1e-30 )
171 for(
long ipSpecies=0; ipSpecies<
nSpecies; ipSpecies++ )
179 int ipHi = (*tr).ipHi();
180 if (ipHi >=
dBaseSpecies[ipSpecies].numLevels_local || (*tr).ipCont() <= 0)
182 if( (*(*tr).Hi()).Pop() > 1e-30 )
193 wl[ip] = (*tr).WLAng();
208 RadPres1 = (*diatom)->H2_RadPress();
214 strcpy(chLab[ip],
"H2 ");
224 fprintf( ioPRESSURE,
"\n" );
243 fprintf( ioPRESSURE,
" P(Lines):" );
244 for( i=0; i <
MIN2(10,ip); i++ )
246 int ipline = iperm[i];
247 fprintf( ioPRESSURE,
"(%4.4s ", chLab[ipline]);
248 prt_wl(ioPRESSURE, wl[ipline] );
249 fprintf(ioPRESSURE,
" %.2f) ",frac[ipline]);
253 fprintf( ioPRESSURE,
"\n" );