69 sprintf( chField7,
"####%3ld",
nzone );
73 sprintf( chField7,
"###%4ld",
nzone );
76 fprintf(
ioQQQ,
" %7.7s %cTe:",chField7, chLet);
78 fprintf(
ioQQQ,
" Hden:");
80 fprintf(
ioQQQ,
" Ne:");
84 fprintf(
ioQQQ,
" R-R0:");
86 fprintf(
ioQQQ,
" dR:");
88 fprintf(
ioQQQ,
" NTR:%3ld Htot:",
conv.nPres2Ioniz);
90 fprintf(
ioQQQ,
" T912:");
92 fprintf(
ioQQQ,
"###\n");
96 fprintf(
ioQQQ,
" H:%.2e %.2e 2H2/H: %.2e He: %.2e %.2e %.2e\n",
117 fprintf(
ioQQQ,
" Abun:" );
118 for( i=0; i <
LIMELM; i++ )
122 fprintf(
ioQQQ,
"\n" );
127 if( !
wind.lgStatic() )
131 if(
wind.AccelTotalOutward == 0. )
134 fac =
wind.AccelTotalOutward;
136 " Dynamics wind V:%.3e km/s a(grav):%.2e a(tot):%.2e Fr(cont):%6.3f "
140 wind.AccelTotalOutward,
142 wind.AccelLine/fac );
160 hatmic +=
mole.species[mol].den*
mole_global.list[mol]->nAtom[elHydrogen];
165 fprintf(
ioQQQ,
" Hydrogen ");
168 fprintf(
ioQQQ,
" H+o/Hden");
171 fprintf(
ioQQQ,
" H- H2");
175 fprintf(
ioQQQ,
" H2+ HeH+");
177 fprintf(
ioQQQ,
" Ho+ ColD");
180 fprintf(
ioQQQ,
"\n");
185 fprintf(
ioQQQ,
" Hydrogen " );
188 fprintf(
ioQQQ,
" H+o/Hden");
191 fprintf(
ioQQQ,
" H- H2");
193 fprintf(
ioQQQ,
" H2+");
195 fprintf(
ioQQQ,
" H3+");
197 fprintf(
ioQQQ,
"\n" );
205 (*diatom)->H2_PrtDepartCoef();
208 if(
prt.lgPrintHeating )
210 fprintf(
ioQQQ,
" ");
216 fprintf(
ioQQQ,
"\n");
226 double Pcon_nT = (Tincid + Tdiff) /
SPEEDLIGHT;
229 if(
prt.lgPrintHeating )
231 fprintf(
ioQQQ,
" ");
233 fprintf(
ioQQQ,
" ");
235 fprintf(
ioQQQ,
" BoundCom");
237 fprintf(
ioQQQ,
" Extra:");
239 fprintf(
ioQQQ,
" Pairs:");
241 fprintf(
ioQQQ,
" H-lines\n");
247 fprintf(
ioQQQ,
" Helium " );
248 for( i=0; i < 3; i++ )
253 fprintf(
ioQQQ,
" HeI 2s3S");
256 fprintf(
ioQQQ,
" Comp H,C");
259 fprintf(
ioQQQ ,
" Fill Fac");
261 fprintf(
ioQQQ ,
" Gam1/tot");
263 fprintf(
ioQQQ,
"\n");
268 fprintf(
ioQQQ,
" Helium " );
273 fprintf(
ioQQQ,
" Comp H,C");
276 fprintf(
ioQQQ ,
" Fill Fac");
278 fprintf(
ioQQQ ,
" Gam1/tot");
280 fprintf(
ioQQQ,
"\n");
286 if(
prt.lgPrintHeating )
293 fprintf(
ioQQQ,
" ");
296 fprintf(
ioQQQ,
" Lines:");
299 fprintf(
ioQQQ,
" Compton:");
301 fprintf(
ioQQQ,
" FFHeatig");
303 fprintf(
ioQQQ,
"\n");
310 fprintf(
ioQQQ,
" He singlet n " );
320 fprintf(
ioQQQ,
" He tripl" );
331 fprintf(
ioQQQ,
"\n" );
339 for( nelem=ipISO; nelem<
LIMELM; ++nelem )
341 if(
dense.lgElmtOn[nelem] )
343 if(
iso_sp[ipISO][nelem].lgPrtLevelPops )
347 if(
iso_sp[ipISO][nelem].lgPrtDepartCoef )
359 fprintf(
ioQQQ,
" Pressure NgasTgas");
361 fprintf(
ioQQQ,
" P(total)");
363 fprintf(
ioQQQ,
" P( gas )");
365 fprintf(
ioQQQ,
" P(Radtn)");
367 fprintf(
ioQQQ,
" Rad accl");
369 fprintf(
ioQQQ,
" ForceMul");
371 fprintf(
ioQQQ,
"\n" );
373 fprintf(
ioQQQ ,
" Texc(La) ");
376 fprintf(
ioQQQ ,
" T(I con)");
378 fprintf(
ioQQQ ,
" T(D con)");
380 fprintf(
ioQQQ ,
" T(U tot)");
383 fprintf(
ioQQQ ,
" nT (c+d)");
386 fprintf(
ioQQQ ,
" Prad/Gas");
389 fprintf(
ioQQQ ,
" Pmag/Gas");
391 fprintf(
ioQQQ,
"\n" );
393 if(
gv.lgGrainPhysicsOn )
395 for(
size_t nd=0; nd <
gv.bin.size(); nd++ )
400 chQHMark = (char)((
gv.bin[nd]->lgQHeat &&
gv.bin[nd]->lgUseQHeat ) ?
'*' :
' ');
401 fprintf(
ioQQQ,
"%-12.12s%c DustTemp",
gv.bin[nd]->chDstLab, chQHMark);
403 fprintf(
ioQQQ,
" Pot Volt");
405 fprintf(
ioQQQ,
" Chrg (e)");
407 fprintf(
ioQQQ,
" drf cm/s");
409 fprintf(
ioQQQ,
" Heating:");
411 fprintf(
ioQQQ,
" Frac tot");
413 fprintf(
ioQQQ,
"\n" );
421 fprintf(
ioQQQ,
" Molecules CH/Ctot:");
423 fprintf(
ioQQQ,
" CH+/Ctot");
425 fprintf(
ioQQQ,
" CO/Ctot:");
427 fprintf(
ioQQQ,
" CO+/Ctot");
429 fprintf(
ioQQQ,
" H2O/Otot");
431 fprintf(
ioQQQ,
" OH/Ototl");
433 fprintf(
ioQQQ,
"\n");
438 (*diatom)->H2_Prt_Zone();
444 fprintf(
ioQQQ,
" Lithium " );
445 for( i=0; i < 4; i++ )
449 fprintf(
ioQQQ,
" Berylliu" );
450 for( i=0; i < 5; i++ )
456 fprintf(
ioQQQ,
" sec ion:" );
458 fprintf(
ioQQQ,
"\n" );
461 if(
prt.lgPrintHeating )
463 fprintf(
ioQQQ,
" " );
464 for( i=0; i < 3; i++ )
468 fprintf(
ioQQQ,
" " );
470 for( i=0; i < 4; i++ )
474 fprintf(
ioQQQ,
"\n" );
481 fprintf(
ioQQQ,
" Boron " );
482 for( i=0; i < 6; i++ )
486 fprintf(
ioQQQ,
"\n" );
489 if(
prt.lgPrintHeating )
491 fprintf(
ioQQQ,
" " );
492 for( i=0; i < 5; i++ )
496 fprintf(
ioQQQ,
"\n" );
501 fprintf(
ioQQQ,
" Carbon " );
502 for( i=0; i < 7; i++ )
507 fprintf(
ioQQQ,
" H2O+/O " );
509 fprintf(
ioQQQ,
" OH+/Otot" );
512 fprintf(
ioQQQ,
" Hex(tot)" );
514 fprintf(
ioQQQ,
"\n" );
517 if(
prt.lgPrintHeating )
519 fprintf(
ioQQQ,
" " );
524 fprintf(
ioQQQ,
"\n" );
528 fprintf(
ioQQQ,
" Nitrogen " );
529 for( i=1; i <= 8; i++ )
533 fprintf(
ioQQQ,
" O2/Ototl" );
535 fprintf(
ioQQQ,
" O2+/Otot" );
537 fprintf(
ioQQQ,
"\n" );
540 if(
prt.lgPrintHeating )
542 fprintf(
ioQQQ,
" " );
547 fprintf(
ioQQQ,
"\n" );
552 fprintf(
ioQQQ,
" Oxygen " );
553 for( i=1; i <= 9; i++ )
557 fprintf(
ioQQQ,
"\n" );
560 if(
prt.lgPrintHeating )
562 fprintf(
ioQQQ,
" " );
567 fprintf(
ioQQQ,
"\n" );
575 if(
dense.lgElmtOn[nelem] )
580 for( i=0; i < nelem+2; i++ )
584 fprintf(
ioQQQ,
"\n" );
587 if(
prt.lgPrintHeating )
589 fprintf(
ioQQQ,
" " );
590 for( i=0; i < nelem+1; i++ )
594 fprintf(
ioQQQ,
"\n" );
602 if(
dense.lgElmtOn[nelem] )
608 ishift =
MAX2(0,
dense.IonHigh[nelem]-LINE+1);
613 for( i=0; i < LINE; i++ )
617 fprintf(
ioQQQ,
"\n" );
620 if(
prt.lgPrintHeating )
622 fprintf(
ioQQQ,
" " );
623 for( i=0; i < LINE; i++ )
628 fprintf(
ioQQQ,
"\n" );