Changeset 1200

Show
Ignore:
Timestamp:
06/18/07 20:08:39 (17 months ago)
Author:
gary
Message:

BUGFIX - in calculation of 21 cm equilibrium timescale, H0 collision rate coefficient is multiplied by xIonDense[ipHELIUM][ipHYDROGEN]. this should only affect 21 cm timescale, not spin temperature.

Files:
1 modified

Legend:

Unmodified
Added
Removed

  • trunk/source/radius_increment.cpp

    r1159 r1200  
    231231        } 
    232232 
    233         /* H 21 cm equilibrium timescale, H21cm returns H (not e) collisional deexcitation rate (not cs) */ 
    234         t = H21cm_H_atom( phycon.te )* dense.xIonDense[ipHELIUM][ipHYDROGEN] + 
     233        /* H 21 cm equilibrium timescale, H21cm returns H (not e) collisional  
     234         * deexcitation rate (not cs) */ 
     235        t = H21cm_H_atom( phycon.te )* dense.xIonDense[ipHYDROGEN][0] + 
    235236                /* >>chng 02 feb 14, add electron term as per discussion in */ 
    236237                /* >>refer      H1      21cm    Liszt, H., 2001, A&A, 371, 698 */