1 get following to work
0 document this problem with molecules - time dependent
2 change the name of this routine to reflect the fact that it is entirely induced.
2 How should these affect locally emitted continuum calculated in rtdiffuse? Just subtract these rates from the emitted ones?
2 test on c checks whether collisions are possible at this temperature, should add photo excitation
2 these pops ARE NOT defined below
2 following needed to get badbugs/bug8.in to work
2 following needed to get badbugs/bug5.in to work
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 this assert is not passed if error made much smaller. This error should be related to a check on convergence of the molecular networks and their sum rules, with a criteria used here and there
0 PvH prefers test against err/10
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
2 - ground term is actually a fix level system, the vectors are created, with pointers ipFe1_54m , ipFe1_111m, must add collision date, use larger model atom
2 following to stop insane FeX strengths >>chng 96 jul 11, removed 1 / 10 factor, so using real data, 90.01 cs = cs * 0.1 >>chng 96 jun 03, transferred following >>chng 97 jan 31, I give up on this mess, use cs of unity
2 update atomic data to Chidichimo et al 1999 AASup 137, 175
2 * following not in cooling function
2 - insert reference
2 - this check may also too liberal.
2 this should be setable at run time as per Robin Williams email rshock = 4e16;
2 Switch to supersonic when bad enough failure for STRONGD -- need to improve this logic when the p(rho) data is cleaner
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
2 - put H2Lines in outward beams in RT_diffuse
2 put supra thermal excitation into excitation of electronic bands
1 add part of hard heat to secondaries
2 find correct high-energy limit for these
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 - put in actual composition dependent Tad - this is only valid for bare surfaces - not ice - for ice Tad is 555K hmi.Tad is binding energy expressed as a temperature note that hmi.Tad is set to 800. in zero tau_nu the first equation in section 2.5 equation one paragraph before equation 2 at low grain temperatures all end in para, J=0
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
2 process is net source term for H(n=3) states, must be added in
2 process is net ionization term for H(n=3) states
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
0 update photoelectric opacity for H2 to include real cross sections and energies. this is not a higher priority because when H2 is formed there can be very little ionizing radiation. this process must be trivial compared with the Solomon process following reference gives cross section for all energies >>refer H2 photo cs Yan, M., Sadeghpour, H.R., & Dalgarno, A., 1998, ApJ, 496, 1044 Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914
0 must include heating, Compton ionization
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? reverse of this reaction i not in detailed balance,why? GS
2 this rate drives numerical instability in such models as secondary1 and 2.in
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
2 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
2 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
2 must add process H2+ + H- => H2 + H, Dalgarno&Lepp 87
2 put in H2+ + gamma => H + H+
1 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
2 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
2 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
2 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
2 equivalent reaction for H2* is not included in chemistry, Big h2 does not include this reaction, what to do? GS
0 had to comment following test out - NA change to hmole caused massive prints
2 put this line back in! EFF = 43.6 / (43.6 + 0.27*COLFAC) R11 = FE(11)*FLUX(IPFE10)*3.122E-12*EFF * 1 ( PFE11A*ESCINC(T353(1),1E-4)/(1.+17.0/5.3*T353(3) ) + 2 PFE11B*ESCINC(T353(1)/3.,1E-4)/(1.+11.0/12.3*T353(3) ) ) contribution to Fe11 from continuum fluorescence call linadd( R11 , 7892 , 'Pump','i')
2 put this in
2 initialization already done at this point
1 streamline all of this, using TransitionZero and then setting dangerously large negative numbers.
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
2 - should this include mass in grains?
2 - should this include mass in grain mantle ice deposits?
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 This doesn't really seem to be the expression above!!!
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
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 should iso.ipLyaRedist[ipHE_LIKE] use ipLY_A as does H-like?
2 file opacities are disabled for now - reinstate this when arrays settle down