44 sp->
fb[0].gamnc =
GammaBn(sp->
fb[0].ipIsoLevNIonCon,
47 sp->
fb[0].xIsoLevNIonRyd,
48 &sp->
fb[0].RecomInducRate,
49 &sp->
fb[0].RecomInducCool_Coef,
57 ionbal.PhotoRate_Shell[nelem][nelem-ipISO][0][0] = sp->
fb[
ipH1s].gamnc;
64 ASSERT(
ionbal.CompRecoilIonRate[nelem][nelem-ipISO]>=0. &&
65 ionbal.CompRecoilHeatRate[nelem][nelem-ipISO]>= 0. );
66 sp->
fb[0].gamnc +=
ionbal.CompRecoilIonRate[nelem][nelem-ipISO];
67 sp->
fb[0].PhotoHeat +=
ionbal.CompRecoilHeatRate[nelem][nelem-ipISO];
70 ionbal.PhotoRate_Shell[nelem][nelem-ipISO][0][0] +=
ionbal.CompRecoilIonRate[nelem][nelem-ipISO];
72 ionbal.PhotoRate_Shell[nelem][nelem-ipISO][0][2] +=
ionbal.CompRecoilHeatRate[nelem][nelem-ipISO];
75 if(
trace.lgTrace &&
trace.lgIsoTraceFull[ipISO] && (nelem ==
trace.ipIsoTrace[ipISO]) )
82 sp->
fb[0].gamnc*0.05);
98 if(
hydro.lgHInducImp )
102 sp->
fb[n].ipIsoLevNIonCon,
105 sp->
fb[n].xIsoLevNIonRyd,
106 &sp->
fb[n].RecomInducRate,
107 &sp->
fb[n].RecomInducCool_Coef,
121 sp->
fb[n].RecomInducRate = 0.;
122 sp->
fb[n].RecomInducCool_Coef = 0.;
133 enum {DEBUG_LOC=
false};
139 fprintf(
ioQQQ,
" buggbugg hphotodebugg%li\t%.2e\t%.2e\t%.2e\t%.2e\t%.2e\t%.2e\n",
153 if(
opac.lgCaseB_no_photo )
157 sp->
fb[n].gamnc = 0.;
158 sp->
fb[n].RecomInducRate = 0.;
159 sp->
fb[n].RecomInducCool_Coef = 0.;
165 enum {DEBUG_LOC=
false};
167 if( DEBUG_LOC && ipISO==1 && nelem==5)
172 sp->
fb[n].RecomInducRate = 0.;
179 fprintf(
ioQQQ,
" iso_photo, ipISO%2ld nelem%2ld low, hi=",ipISO,nelem);
183 fprintf(
ioQQQ,
" total=");
185 fprintf(
ioQQQ,
"\n");
void GammaPrt(long int ipLoEnr, long int ipHiEnr, long int ipOpac, FILE *ioFILE, double total, double threshold)
double GammaBn(long int ipLoEnr, long int ipHiEnr, long int ipOpac, double thresh, double *ainduc, double *rcool, t_phoHeat *photoHeat)