| 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); |
| 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 ) |