51 for( nelem=0; nelem<
LIMELM; ++nelem )
53 for( ion=0; ion<nelem+1; ++ion )
55 ionbal.CompRecoilIonRate[nelem][ion] = 0.;
56 ionbal.CompRecoilHeatRate[nelem][ion] = 0.;
57 if(
dense.xIonDense[nelem][ion] > 0. )
63 for( i=
ionbal.ipCompRecoil[nelem][ion]; i <
rfield.nflux; ++i)
68 ionbal.nCompRecoilElec[nelem-ion];
73 ionbal.CompRecoilIonRate[nelem][ion] += crsphi;
79 recoil_energy =
rfield.anu2[i] /
81 Heavy.Valence_IP_Ryd[nelem][ion];
84 ionbal.CompRecoilHeatRate[nelem][ion] += crsphi*recoil_energy;
97 for( nelem=0; nelem<
LIMELM; ++nelem )
99 for( ion=0; ion<nelem+1; ++ion )
101 ionbal.CompRecoilIonRate[nelem][ion] = 0.;
102 ionbal.CompRecoilHeatRate[nelem][ion] = 0.;
124 if(
hextra.lg_CR_B_equipartition )
129 CosRayDen =
hextra.background_density *
143 hextra.cr_energydensity = CosRayDen/
hextra.background_density *
147 sqthot = sqrt(
hextra.crtemp);
153 ionbal.CosRayHeatThermalElectrons = 5.5e-14/sqthot*CosRayDen;
156 ionbal.CosRayIonRate = 1.22e-4/sqthot*
157 log(2.72*pow(
hextra.crtemp/1e8,0.097))*CosRayDen;
164 ionbal.CosRayIonRate *= 3.;
170 ionbal.CosRayIonRate *= 10.;
184 fprintf(
ioQQQ,
" highen: cosmic ray density;%10.2e CRion rate;%10.2e CR heat rate=;%10.2e CRtemp;%10.2e\n",
190 ionbal.CosRayIonRate = 0.;
191 ionbal.CosRayHeatNeutralParticles = 0.;
216 if(
hextra.lgHextraDepth )
223 if(
hextra.turback != 0. )
229 else if(
hextra.lgHextraDensity )
235 else if(
hextra.lgHextraSS )
240 pow((
double)
hextra.HextraSSradius,-3.0),0.5);
243 else if(
DoppVel.DispScale > 0. )
248 ionbal.ExtraHeatRate = 3.45e-28 / 2.82824 * turb * turb * turb *
260 ionbal.ExtraHeatRate = 0.;
268 if(
hextra.lgNeutrnHeatOn )
278 ionbal.xNeutronHeatRate = 0.;
303 for( i=0; i <
rfield.nflux; i++ )
339 fprintf(
ioQQQ,
" HIGHEN: COOL num=%8.2e HEAT num=%8.2e\n",
351 " HIGHEN finds heating fracs= frac(compt)=%10.2e "
352 " f(pair)%10.2e totHeat=%10.2e\n",