Topics

Number of ionizing photons

vianney.lebouteiller@...
 

Hi there, 
is there a variable that tracks the number of ionizing photons in the C++ code? I can't seem to find an output quantity for it in the various save options and I would like to avoid inferring it from the spectra. 
Cheers,
Vianney

chris_morisset
 

Hi Vianney,
That's printed out in the main output file, just after the printing of the input file:

           8228CellPeak1.00E+00   Lo 3.04e-09=  29.98m   Hi-Con:7.13E+06 Ryd   E(hi):7.35E+06Ryd     E(hi):     100.00 MeV
           L(nu>1ryd):  39.5117   Average nu:1.490E+00   L( X-ray):  16.2184   L(BalC):  39.3644     Q(Balmer C):  50.2404
           Q(1.0-1.8):  49.9066   Q(1.8-4.0):  49.2849   Q(4.0-20):  46.8881   Q(20--):  25.5603     Ion pht flx:1.990E+11
           L(gam ray):   0.0000   Q(gam ray):   0.0000   L(Infred):  38.0213   Alf(ox):   0.0000     Total lumin:  39.7534
           log L/Lsun:   6.1704   Abs bol mg: -10.6835   Abs V mag:  -6.4535   Bol cor:  -4.2300     nuFnu(Bbet):  38.1230
           U(1.0----):2.654E-02   U(4.0----):2.051E-05   T(En-Den):1.493E+01   T(Comp):4.790E+04     nuJnu(912A):7.668E+00
           Occ(FarIR):8.274E-07   Occ(H n=6):8.690E-14   Occ(1Ryd):3.536E-16   Occ(4R):2.549E-20     Occ (Nu-hi):0.000E+00
           Tbr(FarIR):3.977E-10   Tbr(H n=6):3.806E-10   Tbr(1Ryd):5.582E-11   Tbr(4R):1.613E-14     Tbr (Nu-hi):0.000E+00

Just add the "Q(1.0-1.8):  49.9066   Q(1.8-4.0):  49.2849   Q(4.0-20):  46.8881   Q(20--):  25.5603" values.

In a pyCloudy object, the 4 values are in pyCloudy.CloudyModel.Q and the sum is in pyCloudy.CloudyModel.Q0

Saludos,
Christophe

vianney.lebouteiller@...
 

Hi Christophe,
thanks for your reply, it's good to know. What about as a function of depth? 
Cheers,
Vianney

vianney.lebouteiller@...
 

Quick follow-up from a previous email exchange, looks like the best solution for now is to do 
save last ionizing continuum 1 0 every "output_file.dat"
in order to save the continuum >1Ryd, then integrate and convert to U or Q. It's better than saving the full spectrum but it still creates a lot of data. It's fine I guess.
Cheers,
Vianney