Show
Ignore:
Timestamp:
03/09/08 17:20:59 (9 months ago)
Author:
rjrw
Message:

Remove some more likely counterproductive tests on iteration.

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • branches/newmole/source/mole_reactions.cpp

    r1830 r1839  
    870870 
    871871        binding_energy = rate->b; 
    872         /*>>chng 06 nov 28 only include source from molecules if we have an estimated first  
    873          * solution - first test is that we have called mole at least twice, 
    874          * second test is that we are on a later iteration */ 
    875         if( conv.nTotalIoniz > 1 || iteration > 1 ) 
    876         { 
    877                 for( nd=0; nd < gv.nBin; nd++ ) 
    878                 { 
    879                         exponent += exp(-binding_energy/gv.bin[nd]->tedust); 
    880                 }        
    881         } 
     872        for( nd=0; nd < gv.nBin; nd++ ) 
     873        { 
     874                exponent += exp(-binding_energy/gv.bin[nd]->tedust); 
     875        }        
    882876        number_of_sites = 1.5e15; 
    883877 
     
    11371131STATIC double hopexch(struct mole_reaction_s *rate) 
    11381132{ 
    1139         if ( iteration==1 && conv.lgSearch ) 
    1140                 return 0.; 
    1141         else 
    1142                 return atmdat.HCharExcRecTo[ipOXYGEN][0]; 
     1133        return atmdat.HCharExcRecTo[ipOXYGEN][0]; 
    11431134} 
    11441135STATIC double hpoexch(struct mole_reaction_s *rate) 
    11451136{ 
    1146         if ( iteration==1 && conv.lgSearch ) 
    1147                 return 0.; 
    1148         else 
    1149                 return atmdat.HCharExcIonOf[ipOXYGEN][0]; 
     1137        return atmdat.HCharExcIonOf[ipOXYGEN][0]; 
    11501138} 
    11511139STATIC double hmattach(struct mole_reaction_s *rate)