Changeset 2025

Show
Ignore:
Timestamp:
05/05/08 09:02:52 (2 weeks ago)
Author:
peter
Message:

source/grains.cpp:

Use safe lower bound for PAH abundance in case neutral hydrogen density gets very low or zero.

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved

  • trunk/source/grains.cpp

    r1960 r2025  
    15981598                        * PAHs are strong near the ionization front and weak in the ionized and molecular gas */ 
    15991599                        /* >>chng 04 sep 28, propto atomic fraction */ 
    1600                         GrnStdDpth_v = dense.xIonDense[ipHYDROGEN][0]/dense.gas_phase[ipHYDROGEN]; 
     1600                        GrnStdDpth_v = max(1.e-10,dense.xIonDense[ipHYDROGEN][0]/dense.gas_phase[ipHYDROGEN]); 
    16011601                } 
    16021602                else if( strcmp( gv.chPAH_abundance_fcn , "CON" )==0 )