2 test on c checks whether collisions are possible at this temperature, should add photo excitation
2 these pops ARE NOT defined below
1 update ct to Kimura et al. (1996)
1 above rate not intended for very low temperatures - find ref for low-T rate, probably is 1e-9 like above
0 these should be values at 1e5 K
2 not currently used - include as deexcitation process
2 this is redundant with contents of oxygen line arrays use them instead when removing this, make sure all line intensity predictions also go into oi line arrays
2 what are we trying to print here?
2 add term for protons from Rouef, E., & Le Bourlot, J. 1990, A&A, 236, 515
1 add neutral helium Staemmler, V., & Flower, D. R. 1991, J. Phys. B, 24, 2343
3 change to atom_level3
1 add this as a Lya excitation process
2 put supra thermal excitation into excitation of electronic bands
2 - insert reference
2 - this check may also too liberal.
2 find a transition probability for this 2^3P0 - 2^3P1 transition. It will require a bit of trickery to insert into the rate matrix, because of the fact that the lower level has a higher index. See discussion "Energy order within 2 3P" near the top of helike.c
2 find a transition probability for this 2^3P1 - 2^3P2 transition. It will require a bit of trickery to insert into the rate matrix, because of the fact that the lower level has a higher index. See discussion "Energy order within 2 3P" near the top of helike.c
2 this algorithm might be more efficient with Brent
2 should collisional heating/cooling be included here?
2 use something like log(ThermRatio) + log(delta) ????
2 this algorithm might be more efficient with Brent
1 add part of hard heat to secondaries
2 find correct high-energy limit for these
2 initialization already done at this point
1 streamline all of this, using transition::Zero and then setting dangerously large negative numbers.
2 following always true, why? either remove test or use it - it is here to save time - this step routine is called very often
2 this should be done with new populations after converged soln
2 add charged heavy elements
1 stupid - why this test on opacity_abs ? - we only get here if we already passed above test on this very thing
2 Need to include divergence terms in pressure balance if flux index is != 0
2 Need to include divergence terms in pressure balance if flux index is != 0
2 this is the total binding energy of the molecules, and is negative, the energy need to get back to free atoms never set and only appears in print statements
0 add this press term due to cosmic rays - hextra.cr_energydensity
0 code has variable magnetic energydensity and pressure, which are equal, as they must be - del one or the other
3 consider making the fine opacity array a double. with a float, the opacity itself often becomes a denormalized number, it then becomes significant when multiplied by dr - can cause numerical noise. this is why the coarse opacity array is a double
1 set temp of 1e4K will be too coarse a line for PDRs where H2 line overlap is very important
2 should use ConEmitLocal for all continua, but not followed by rfield.DiffuseEscape - put that at the end. Once continua all bundled this way, it will be easy to save them as a function of depth and then do exact rt
2 Replace this constant with the appropriate macro, if any
2 add fegrain to outward beams, but within main formalism by including grains in all x-ray processes
1 dynamics; this test assumes that neg vel are subsonic, so that sobolev length would overestimate the optical depth, since ion is at most present over computed slab, and possibly more.
1 rewrite so that this checks on size not sign of windv
- check rebinning of Tlusty models
- why was it necessary to change stars_tlusty.in? (change from r43 to r50?)
- check all interpolation modes of CoStar
- compare models with original code, dump atmospheres!
- check all Edges arrays...
- update Doxygen documentation
2 - should this include mass in grains?
2 - should this include mass in grain mantle ice deposits?
xray - secondaries from incident electrons still need to be added in
xray - primary, secondary, auger electrons need to be added into suprathermals
2 change this to CS_new
2 file opacities are disabled for now - reinstate this when arrays settle down