38 CollisIonizatCoolingTotal,
39 dCollisIonizatCoolingTotalDT,
43 CollisIonizatCoolingDT,
50 valarray<double> CollisIonizatCoolingArr,
51 CollisIonizatCoolingDTArr,
56 long int nlo_heat_max , nhi_heat_max;
125 CollisIonizatCoolingTotal = 0.;
126 dCollisIonizatCoolingTotalDT = 0.;
139 CollisIonizatCooling =
142 CollisIonizatCoolingTotal += CollisIonizatCooling;
146 CollisIonizatCoolingDT = CollisIonizatCooling*
149 dCollisIonizatCoolingTotalDT += CollisIonizatCoolingDT;
153 CollisIonizatCoolingArr[n] = CollisIonizatCooling;
154 CollisIonizatCoolingDTArr[n] = CollisIonizatCoolingDT;
168 double dcool = dCollisIonizatCoolingTotalDT *
EN1RYD *
192 fprintf(
ioQQQ,
"DEBUG coll ioniz cool contributors:");
195 if( CollisIonizatCoolingArr[n] /
SDIV( CollisIonizatCoolingTotal ) > 0.1 )
196 fprintf(
ioQQQ,
" %2li %.1e",
198 CollisIonizatCoolingArr[n]/
SDIV( CollisIonizatCoolingTotal ) );
201 fprintf(
ioQQQ,
"DEBUG coll ioniz derivcontributors:");
204 if( CollisIonizatCoolingDTArr[n] /
SDIV( dCollisIonizatCoolingTotalDT ) > 0.1 )
205 fprintf(
ioQQQ,
" %2li %.1e",
207 CollisIonizatCoolingDTArr[n]/
SDIV( dCollisIonizatCoolingTotalDT ) );
227 edenHCorr_IonAbund = collider*
dense.
xIonDense[nelem][nelem+1-ipISO];
256 iso.
RadRecCool[ipISO][nelem] += SaveRadRecCool[0] * edenIonAbund;
276 ThinCoolingSum += thin;
280 enum {DEBUG_LOC=
false};
288 fprintf(
ioQQQ,
"\t%.2f",SaveRadRecCool[n]/ThinCoolingSum);
305 ThinCoolingCaseB = (-25.859117 +
317 ThinCoolingCaseB = (-25.859117 +
328 ThinCoolingCaseB =
POW3(1.+nelem-ipISO)*pow(10.,ThinCoolingCaseB)/(
phycon.
te/
POW2(1.+nelem-ipISO) );
331 RecCoolExtra = ThinCoolingCaseB - ThinCoolingSum;
335 ThinCoolingCaseB = 0.;
341 RecCoolExtra =
MAX2(0., RecCoolExtra );
360 HeatExcited += SavePhotoHeat[n];
361 if( SavePhotoHeat[n] > biggest )
363 biggest = SavePhotoHeat[n];
369 enum {DEBUG_LOC=
false};
370 if( DEBUG_LOC && ipISO==0 && nelem==0 &&
nzone > 700)
375 fprintf(
ioQQQ,
"ipISO = %li nelem=%li ipbig=%li biggest=%g HeatExcited=%.2e ctot=%.2e\n",
384 fprintf(
ioQQQ,
"DEBUG phot heat%2li\t%.2e\t%.2e\t%.2e\t%.2e\t%.2e\n",
386 SavePhotoHeat[n]/HeatExcited,
436 enum {DEBUG_LOC=
false};
437 if( DEBUG_LOC && ipISO==0 && nelem==5 )
439 fprintf(
ioQQQ,
" ipISO=%li nelem=%li ctot = %.2e\n",
443 fprintf(
ioQQQ,
"sum\t%.2e\t%.2e\t%.2e\n",
447 fprintf(
ioQQQ,
"sum\tp ht\tr cl\ti cl\n");
452 fprintf(
ioQQQ,
"%li\t%.2e\t%.2e\t%.2e\t%.2e\t%.2e\t%.2e \n",
461 fprintf(
ioQQQ,
" \n");
489 StatesElem[ipISO][nelem][ipHi].Pop)*edenHCorr_IonAbund*
508 if( hlone < heat_max )
520 enum {DEBUG_LOC=
false};
524 fprintf(
ioQQQ,
"%li %li %.2e\n", nlo_heat_max, nhi_heat_max, heat_max );
539 StatesElem[ipISO][nelem][ipHi].Pop)* edenHCorr_IonAbund*
566 StatesElem[ipISO][nelem][ipHi].Pop)*edenHCorr_IonAbund*
575 StatesElem[ipISO][nelem][ipHi].Pop)*edenHCorr_IonAbund*
594 StatesElem[ipISO][nelem][ipHi].Pop)*edenHCorr_IonAbund*
619 iso.
dLTot[ipISO][nelem] = -dCdT_all;
624 enum {DEBUG_LOC=
false};
629 fprintf(
ioQQQ,
"%.2e la %.2f restly %.2f barest %.2f hrest %.2f\n",
641 enum {DEBUG_LOC=
false};
644 if( DEBUG_LOC && (nelem==1 || nelem==0) )
652 fprintf(
ioQQQ,
"%li\t%li\t%g\t%g\t%g\t%g\tT=\t%g\t%g\t%g\t%g\n", nelem,n,
677 fprintf(
ioQQQ,
"%li\t%li\t%g\t%g\t%g\t%g\tT=\t%g\t%g\t%g\t%g\n", nelem,n,
686 SaveInducCool[n]+SaveRadRecCool[n]*edenIonAbund ,
690 SaveRadRecCool[n]*edenIonAbund );
697 #if defined(__HP_aCC)