70 long int ipHi , ipLo , i;
75 " HydroLevel%3ld finds arrays, with optical depths defined? %li induced 2ph=%12.3e\n",
80 fprintf(
ioQQQ,
"up:%2ld", ipHi );
81 fprintf(
ioQQQ,
"lo" );
82 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
84 fprintf(
ioQQQ,
"%9ld", ipLo );
86 fprintf(
ioQQQ,
"\n" );
88 fprintf(
ioQQQ,
"%3ld", ipHi );
89 fprintf(
ioQQQ,
" A*esc" );
90 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
95 fprintf(
ioQQQ,
"\n" );
97 fprintf(
ioQQQ,
"%3ld", ipHi );
98 fprintf(
ioQQQ,
" A*ees" );
99 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
104 fprintf(
ioQQQ,
"\n" );
106 fprintf(
ioQQQ,
"%3ld", ipHi );
107 fprintf(
ioQQQ,
" tauin" );
108 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
112 fprintf(
ioQQQ,
"\n" );
114 fprintf(
ioQQQ,
"%3ld", ipHi );
115 fprintf(
ioQQQ,
" t tot" );
116 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
120 fprintf(
ioQQQ,
"\n" );
122 fprintf(
ioQQQ,
"%3ld", ipHi );
123 fprintf(
ioQQQ,
" Esc " );
124 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
128 fprintf(
ioQQQ,
"\n" );
130 fprintf(
ioQQQ,
"%3ld", ipHi );
131 fprintf(
ioQQQ,
" Eesc " );
132 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
136 fprintf(
ioQQQ,
"\n" );
138 fprintf(
ioQQQ,
"%3ld", ipHi );
139 fprintf(
ioQQQ,
" Dest " );
140 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
144 fprintf(
ioQQQ,
"\n" );
146 fprintf(
ioQQQ,
"%3ld", ipHi );
147 fprintf(
ioQQQ,
" A*dst" );
148 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
153 fprintf(
ioQQQ,
"\n" );
155 fprintf(
ioQQQ,
"%3ld", ipHi );
156 fprintf(
ioQQQ,
" StrkE" );
157 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
161 fprintf(
ioQQQ,
"\n" );
163 fprintf(
ioQQQ,
"%3ld", ipHi );
164 fprintf(
ioQQQ,
" B(ul)" );
165 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
170 fprintf(
ioQQQ,
"\n" );
172 fprintf(
ioQQQ,
"%3ld", ipHi );
173 fprintf(
ioQQQ,
" tcont" );
174 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
178 fprintf(
ioQQQ,
"\n" );
180 fprintf(
ioQQQ,
"%3ld", ipHi );
181 fprintf(
ioQQQ,
" C(ul)" );
182 for( ipLo=
ipH1s; ipLo < ipHi; ipLo++ )
186 fprintf(
ioQQQ,
"\n" );
190 fprintf(
ioQQQ,
" FeIIo");
193 fprintf(
ioQQQ,
"\n");
197 fprintf(
ioQQQ,
" " );
201 fprintf(
ioQQQ,
"%9ld", i );
203 fprintf(
ioQQQ,
"\n" );
221 if( (
trace.lgTrace &&
trace.lgIsoTraceFull[ipISO]) && (nelem ==
trace.ipIsoTrace[ipISO]) )
228 if( (
trace.lgIsoTraceFull[ipISO] &&
trace.lgTrace) && (nelem ==
trace.ipIsoTrace[ipISO]) )
230 fprintf(
ioQQQ,
" HLEV HGAMNC" );
233 for( i=
ipH2s; i <
iso_sp[ipISO][nelem].numLevels_local; i++ )
237 fprintf(
ioQQQ,
"\n" );
239 fprintf(
ioQQQ,
" HLEV TOTCAP" );
241 for( i=1; i <
iso_sp[ipISO][nelem].numLevels_local; i++ )
245 fprintf(
ioQQQ,
" tot");
247 fprintf(
ioQQQ,
"\n" );
249 fprintf(
ioQQQ,
" HLEV IND Rc" );
251 for( i=
ipH1s; i <
iso_sp[ipISO][nelem].numLevels_local; i++ )
255 fprintf(
ioQQQ,
"\n" );
258 fprintf(
ioQQQ,
" IND Rc LTE " );
260 for( i=
ipH1s; i <
iso_sp[ipISO][nelem].numLevels_local; i++ )
263 iso_sp[ipISO][nelem].fb[i].gamnc*
iso_sp[ipISO][nelem].fb[i].PopLTE ));
265 fprintf(
ioQQQ,
"\n" );
268 fprintf(
ioQQQ,
" HLEV HLTE" );
270 for( i=
ipH1s; i <
iso_sp[ipISO][nelem].numLevels_local; i++ )
274 fprintf(
ioQQQ,
"\n" );
277 fprintf(
ioQQQ,
" HLEVfr cion" );
279 for( i=
ipH1s; i <
iso_sp[ipISO][nelem].numLevels_local; i++ )
282 iso_sp[ipISO][nelem].fb[i].ColIoniz*
dense.EdenHCorr/
MAX2(1e-30,
iso_sp[ipISO][nelem].fb[i].RateLevel2Cont) ) );
284 fprintf(
ioQQQ,
"\n" );
287 if(
ionbal.RateRecomTot[nelem][nelem]> 0. )
289 fprintf(
ioQQQ,
" HLEVfrPhIon" );
291 for( i=
ipH1s; i <
iso_sp[ipISO][nelem].numLevels_local; i++ )
294 iso_sp[ipISO][nelem].fb[i].gamnc/
MAX2(1e-30,
iso_sp[ipISO][nelem].fb[i].RateLevel2Cont) ) );
296 fprintf(
ioQQQ,
"\n" );
299 fprintf(
ioQQQ,
" HLEV HN" );
301 for( i=
ipH1s; i <
iso_sp[ipISO][nelem].numLevels_local; i++ )
305 fprintf(
ioQQQ,
"\n" );
307 fprintf(
ioQQQ,
" HLEV b(n)" );
309 for( i=
ipH1s; i <
iso_sp[ipISO][nelem].numLevels_local; i++ )
313 fprintf(
ioQQQ,
"\n" );
315 fprintf(
ioQQQ,
" HLEV X12tot");
317 fprintf(
ioQQQ,
" Grn dest:");
319 ionbal.RateIoniz[nelem][nelem][nelem+1] ));
320 fprintf(
ioQQQ,
"\n");
328 fprintf(
ioQQQ,
" HydroLevel Z:%2ld return %s te=",
330 iso_sp[ipISO][nelem].chTypeAtomUsed );
332 fprintf(
ioQQQ,
" density=%.4e",
dense.xIonDense[nelem][nelem-ipISO] );
334 fprintf(
ioQQQ,
" simple=%.4e",
iso_sp[ipISO][nelem].xIonSimple);
338 fprintf(
ioQQQ,
" ion rate=%.4e",
ionbal.RateIonizTot(nelem,nelem-ipISO) );
340 fprintf(
ioQQQ,
" TotRec=%.4e",
ionbal.RateRecomTot[nelem][nelem-ipISO]);
342 fprintf(
ioQQQ,
" RadRec=%.4e",
iso_sp[ipISO][nelem].RadRec_effec);
343 fprintf(
ioQQQ,
"\n");