Show
Ignore:
Timestamp:
05/06/08 16:35:36 (8 months ago)
Author:
rjrw
Message:

Further refactor newton_step to more purely numerical function.

Remove more dead or leftover code and data.

Files:
1 modified

Legend:

Unmodified
Added
Removed

  • branches/newmole/source/molcol.cpp

    r1780 r2028  
    1919 
    2020        DEBUG_ENTRY( "molcol()" ); 
    21  
    22         if( strcmp(chLabel,"PRIN") == 0 ) 
    23         { 
    24                 /* total hydrogen column density, all forms */ 
    25                 fprintf( ioMEAN, "\n                                                     Log10 Column density (cm^-2)\n"); 
    26                 fprintf( ioMEAN, "   Htot  :"); 
    27                 fprintf( ioMEAN, "%7.3f",log10(MAX2(SMALLFLOAT,colden.colden[ipCOL_HTOT]))); 
    28                 fprintf( ioMEAN, "   HII   :"); 
    29                 fprintf( ioMEAN, "%7.3f",log10(MAX2(SMALLFLOAT,colden.colden[ipCOL_Hp]))); 
    30                 fprintf( ioMEAN, "   HI    :"); 
    31                 fprintf( ioMEAN, "%7.3f",log10(MAX2(SMALLFLOAT,colden.colden[ipCOL_H0]))); 
    32                 fprintf( ioMEAN, "\n"); 
    33                 fprintf( ioMEAN, "   H-    :"); 
    34                 fprintf( ioMEAN, "%7.3f",log10(MAX2(SMALLFLOAT,colden.colden[ipCOL_HMIN]))); 
    35                 fprintf( ioMEAN, "   H2g   :"); 
    36                 fprintf( ioMEAN, "%7.3f",log10(MAX2(SMALLFLOAT,colden.colden[ipCOL_H2g]))); 
    37                 fprintf( ioMEAN, "   H2*   :"); 
    38                 fprintf( ioMEAN, "%7.3f",log10(MAX2(SMALLFLOAT,colden.colden[ipCOL_H2s]))); 
    39                 fprintf( ioMEAN, "   H2+   :"); 
    40                 fprintf( ioMEAN, "%7.3f",log10(MAX2(SMALLFLOAT,colden.colden[ipCOL_H2p]))); 
    41                 fprintf( ioMEAN, "   HeH+  :"); 
    42                 fprintf( ioMEAN, "%7.3f",log10(MAX2(SMALLFLOAT,colden.colden[ipCOL_HeHp] ))); 
    43                 fprintf( ioMEAN, "\n"); 
    44                 fprintf( ioMEAN, "   H3+   :"); 
    45                 fprintf( ioMEAN, "%7.3f",log10(MAX2(SMALLFLOAT,colden.colden[ipCOL_H3p] ))); 
    46                 fprintf( ioMEAN, "\n"); 
    47         } 
    4821 
    4922        /* call large H2 and CO column density routines which will do their jobs */ 
    … …  
    7245        else if( strcmp(chLabel,"PRIN") == 0 ) 
    7346        { 
    74                 /*  print the molecular column densities 
    75                  * want to print all the molecules, not including the atoms/ions 
    76                  * that are part of the co solver.  use first to print them all */ 
    77                 /*for( i=0; i < mole.num_comole_calc; i++ )*/ 
     47                /*  print the molecular column densities */ 
    7848                int j=0; 
    7949                for( i=0; i < mole.num_calc; i++ ) 
    … …  
    8454                                fprintf( ioMEAN, "%7.3f",log10(MAX2(SMALLFLOAT,mole.list[i]->column ))); 
    8555                                j++;  
     56                                fixit();  /* Separate elemental chem more clearly -- would need to refactor mole.nElem */ 
    8657                                if( j == 8 ) 
    8758                                {