50 if(
prt.lgPrtCitations )
54 fprintf(
ioQQQ,
"\n" );
56 fprintf(
ioQQQ,
"\n\n" );
61 fprintf(
ioQQQ,
" Lo");
63 fprintf(
ioQQQ,
"=%6.2fcm Hi-Con:",
66 fprintf(
ioQQQ,
" Ryd E(hi):");
68 fprintf(
ioQQQ,
"Ryd E(hi): %9.2f MeV\n",
rfield.egamry*0.0000136 );
107 fprintf(
ioQQQ,
" L(nu>1ryd):%9.4f Average nu:",pionl);
109 fprintf(
ioQQQ,
" L( X-ray):%9.4f L(BalC):%9.4f Q(Balmer C):%9.4f\n",
110 prt.xpow, pballog, qballog );
113 fprintf(
ioQQQ,
"\n\n WARNING - the continuum has a luminosity %.2e times greater than the sun.\n",
115 fprintf(
ioQQQ,
" WARNING - Is this correct? Check the luminosity commands.\n\n\n");
120 fprintf(
ioQQQ,
" I(nu>1ryd):%9.4f Average nu:",pionl);
122 fprintf(
ioQQQ,
" I( X-ray):%9.4f I(BalC):%9.4f Phi(BalmrC):%9.4f\n",
158 " Q(1.0-1.8):%9.4f Q(1.8-4.0):%9.4f Q(4.0-20):"
159 "%9.4f Q(20--):%9.4f Ion pht flx:",
169 " phi(1.0-1.8):%7.4f phi(1.8-4.0):%7.3f phi(4.0-20):"
170 "%7.3f phi(20--):%7.3f Ion pht flx:",
177 fprintf(
ioQQQ,
"\n");
185 if(
rfield.flux[0][ip2500-1] > 1e-28 )
197 alfox = log(ratio)/log(
rfield.anu[ip2kev-1]/
rfield.anu[ip2500-1]);
209 if(
prt.GammaLumin > 0. )
211 gpowl = log10(
prt.GammaLumin) +
radius.pirsq;
220 if(
prt.pradio > 0. )
222 radpwl = log10(
prt.pradio) +
radius.pirsq;
232 " L(gam ray):%9.4f Q(gam ray):%9.4f L(Infred):%9.4f Alf(ox):%9.4f Total lumin:%9.4f\n",
233 gpowl, qgaml, radpwl, alfox, log10(
continuum.TotalLumin) +
239 " I(gam ray):%9.4f phi(gam r):%9.4f I(Infred):%9.4f Alf(ox):%9.4f Total inten:%9.4f\n",
240 gpowl, qgaml, radpwl, alfox, log10(
continuum.TotalLumin) +
248 absbol = 4.75 - 2.5*solar;
265 bolc = absbol - absv;
267 " log L/Lsun:%9.4f Abs bol mg:%9.4f Abs V mag:%9.4f Bol cor:%9.4f nuFnu(Bbet):%9.4f\n",
278 for( i=0; i <
rfield.nflux; i++ )
309 if(
phycon.TEnerDen > (1.05*tcomp) )
319 fprintf(
ioQQQ,
" U(1.0----):");
321 fprintf(
ioQQQ,
" U(4.0----):");
323 fprintf(
ioQQQ,
" T(En-Den):");
325 fprintf(
ioQQQ,
" T(Comp):");
327 fprintf(
ioQQQ,
" nuJnu(912A):");
329 fprintf(
ioQQQ,
"\n");
332 fprintf(
ioQQQ,
" Occ(FarIR):");
334 fprintf(
ioQQQ,
" Occ(H n=6):");
345 fprintf(
ioQQQ,
" Occ(1Ryd):");
347 fprintf(
ioQQQ,
" Occ(4R):");
349 fprintf(
ioQQQ,
" Occ (Nu-hi):");
351 fprintf(
ioQQQ,
"\n");
355 fprintf(
ioQQQ,
" Tbr(FarIR):");
357 fprintf(
ioQQQ,
" Tbr(H n=6):");
367 fprintf(
ioQQQ,
" Tbr(1Ryd):");
369 fprintf(
ioQQQ,
" Tbr(4R):");
371 fprintf(
ioQQQ,
" Tbr (Nu-hi):");
373 fprintf(
ioQQQ,
"\n");
378 " >>>The radio brightness temperature is very large,%10.2eK at%10.2ecm. Is this physical???\n",
379 tradio, 9.115e-6/
rfield.anu[0] );
383 fprintf(
ioQQQ,
" \n" );