Show
Ignore:
Timestamp:
05/10/08 09:03:02 (4 months ago)
Author:
peter
Message:

Merging all changes from mainline upto r2033, except r1902.

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • branches/c08_branch/source/iso_photo.cpp

    r1732 r2034  
    1111#include "gammas.h" 
    1212#include "iso.h" 
    13 /*lint -e661 Possible access of out-of-bounds pointer*/ 
    14 /*lint -e662 (Warning -- Possible creation of out-of-bounds pointer  */ 
    1513 
    1614void iso_photo( 
     
    7169 
    7270        /* option to print ground state photoionization rates */ 
    73         if( /*nelem==0*/ trace.lgTrace && trace.lgHBug && false ) 
     71        if( trace.lgTrace && trace.lgIsoTraceFull[ipISO] && (nelem == trace.ipIsoTrace[ipISO]) ) 
    7472        { 
    7573                GammaPrt(iso.ipIsoLevNIonCon[ipISO][nelem][0], 
    7674                        rfield.nflux, 
    7775                        iso.ipOpac[ipISO][nelem][0], 
    78                   ioQQQ, 
    79                   iso.gamnc[ipISO][nelem][0],iso.gamnc[ipISO][nelem][0]*0.05); 
    80         } 
    81  
    82         { 
    83                 /*@-redef@*/ 
    84                 enum {DEBUG_LOC=false}; 
    85                 /*@+redef@*/ 
    86                 if( DEBUG_LOC ) 
    87                 { 
    88                         /* option to print ground state photoionization rates */ 
    89                         if( DEBUG_LOC && nelem==ipHELIUM && ipISO==ipH_LIKE && nzone > 290) 
    90                         { 
    91                                 GammaPrt(iso.ipIsoLevNIonCon[ipISO][nelem][0], 
    92                                         rfield.nflux, 
    93                                         iso.ipOpac[ipISO][nelem][0], 
    94                                 ioQQQ, 
    95                                 iso.gamnc[ipISO][nelem][0],iso.gamnc[ipISO][nelem][0]*0.05); 
    96                         } 
    97                 } 
     76                        ioQQQ, 
     77                        iso.gamnc[ipISO][nelem][0], 
     78                        iso.gamnc[ipISO][nelem][0]*0.05); 
    9879        } 
    9980 
     
    148129                if( DEBUG_LOC ) 
    149130                { 
    150                         n = 1; 
    151                         if( ipISO==ipHE_LIKE && nelem==ipHELIUM && nzone > 200 ) 
    152                         { 
    153                                 GammaPrt(iso.ipIsoLevNIonCon[ipISO][nelem][n], 
    154                                         iso.ipIsoLevNIonCon[ipISO][nelem][ipH1s]-1, 
    155                                         iso.ipOpac[ipISO][nelem][n], 
    156                                         ioQQQ, 
    157                                         iso.gamnc[ipISO][nelem][n], 
    158                                         iso.gamnc[ipISO][nelem][n]*0.05); 
    159                         } 
    160                 } 
    161         } 
    162  
    163         { 
    164                 /*@-redef@*/ 
    165                 enum {DEBUG_LOC=false}; 
    166                 /*@+redef@*/ 
    167                 if( DEBUG_LOC ) 
    168                 { 
    169                         if( ipISO==1 && nelem==1 && nzone > 200 ) 
    170                                 fprintf(ioQQQ,"DEBUG isophoto\t%.2f\t%.4e\n", 
    171                                 fnzone,  
    172                                 iso.gamnc[ipISO][nelem][1]); 
    173                 } 
    174         } 
    175  
    176         { 
    177                 /*@-redef@*/ 
    178                 enum {DEBUG_LOC=false}; 
    179                 /*@+redef@*/ 
    180                 if( DEBUG_LOC ) 
    181                 { 
    182                         if( nelem==0 ) 
     131                        if( nelem==ipHYDROGEN ) 
    183132                        { 
    184133                                fprintf(ioQQQ," buggbugg hphotodebugg%li\t%.2e\t%.2e\t%.2e\t%.2e\t%.2e\t%.2e\n", 
     
    232181        return; 
    233182} 
    234 /*lint +e661 Possible access of out-of-bounds pointer*/ 
    235 /*lint +e662 (Warning -- Possible creation of out-of-bounds pointer  */