Show
Ignore:
Timestamp:
03/22/08 14:45:49 (10 months ago)
Author:
gary
Message:

BUGFIX - He+ + H0 -> He0 + H+ was counted as the reverse reaction

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • trunk/source/atmdat_char_tran.cpp

    r1771 r1879  
    9090                } 
    9191 
    92                 /* >>chng 04 apr 27, the above included only the radiative charge transfer from 
     92                /* The above included only the radiative charge transfer from 
    9393                 * Stancil et al 1998.  must explicitly add on the ct fitted by Kingdon & Ferland, 
    94                  * the code must have originally had the KF fit, but got hosed at some point  
    95                  * these are two different sets of physics  
    96                  * He0 + H+ -> He+ + H0 */ 
     94                 * The process H0 + He+ -> He0 + H+ */ 
    9795                if( phycon.te > 6000. ) 
    98                         atmdat.HCharExcIonOf[ipHELIUM][0] += 7.47e-15*pow(phycon.te/1e4,2.06)* 
     96                        atmdat.HCharExcRecTo[ipHELIUM][0] += 7.47e-15*pow(phycon.te/1e4,2.06)* 
    9997                        (1.+9.93*sexp(3.89e-4*phycon.te) ); 
    10098