65 limit_hi =
MAX2( 1 , limit_hi );
69 if( !
conv.nPres2Ioniz &&
gv.lgDustOn() )
75 limit_lo =
dense.IonLow[nelem];
81 for( ion=limit_lo; ion < limit_hi; ion++ )
84 for( ns=0; ns <
Heavy.nsShells[nelem][ion]; ns++ )
87 if( (ns==(
Heavy.nsShells[nelem][ion]-1) ||
opac.lgRedoStatic) )
90 iplow =
opac.ipElement[nelem][ion][ns][0];
91 iphi =
opac.ipElement[nelem][ion][ns][1];
92 ipop =
opac.ipElement[nelem][ion][ns][2];
98 ionbal.PhotoRate_Shell[nelem][ion][ns][0] =
101 &photoHeat )*
ionbal.lgPhotoIoniz_On;
115 ns = (
Heavy.nsShells[nelem][ion]-1);
117 ionbal.PhotoRate_Shell[nelem][ion][ns][0] +=
ionbal.CompRecoilIonRate[nelem][ion];
119 ionbal.PhotoRate_Shell[nelem][ion][ns][2] +=
ionbal.CompRecoilHeatRate[nelem][ion];
129 opac.ipElement[nelem][ion][ns][0],
130 opac.ipElement[nelem][ion][ns][1],
131 opac.ipElement[nelem][ion][ns][2],
133 ionbal.PhotoRate_Shell[nelem][ion][ns][0],
137 for( ns=0; ns <
Heavy.nsShells[nelem][0]; ns++ )
142 for( ion=0; ion <
dense.IonHigh[nelem]; ion++ )
144 if(
Heavy.nsShells[nelem][ion] > ns )
146 fprintf(
ioQQQ,
" %8.1e",
ionbal.PhotoRate_Shell[nelem][ion][ns][0] );
161 long ns = 6, ion = 1;
162 ionbal.PhotoRate_Shell[nelem][ion][ns][0] +=
ca.dstCala;
171 ionbal.PhotoRate_Shell[nelem][ion][ns][0] =
172 (
HMRATE((1e-10)*3.0,0,0)*(
hmi.UV_Cont_rel2_Habing_TH85_face*
173 exp(-(3.0*
rfield.extin_mag_V_point))/1.66));
175 ionbal.PhotoRate_Shell[nelem][ion][ns][1] = 0.;
176 ionbal.PhotoRate_Shell[nelem][ion][ns][2] = 0.;
255 opac.ipo3exc[2] , 1., &dummy ));
259 opac.ipo3exc3[2] , 1., &dummy ));
263 opac.ipo1exc[2] , 1., &dummy ));
266 aeff = 0.0263 +
oxy.d5007r;
271 enum {DEBUG_LOC=
false};
274 fprintf(
ioQQQ,
"pop rel %.1e rate %.1e grnd rate %.1e\n",
284 (1. -
oxy.poiexc) +
oxy.d6300*
oxy.poiexc;
288 (1. -
oxy.poiii2 -
oxy.poiii3) +
oxy.d5007r*
oxy.poiii2 +
359 const int NDIM =
ipIRON+1;
361 static const double fyield[NDIM+1] = {.34,.34,.35,.35,.36,.37,.37,.38,.39,.40,
362 .41,.42,.43,.44,.45,.46,.47,.47,.48,.48,.49,.49,.11,.75,0.,0.,0.};
364 long int i, limit, limit2;
385 for( i=limit; i < limit2; i++ )