22 if(
dense.EdenSet > 0.f )
26 dense.eden_from_metals = 1.;
29 fprintf(
ioQQQ,
" eden_sum zn: %.2f eden set to: %.4e\n",
fnzone,
dense.EdenSet );
31 else if(
dense.EdenFraction > 0.f )
35 dense.eden_from_metals = 1.;
39 fprintf(
ioQQQ,
" eden_sum zn: %.2f eden ratio set to: %.4e, eden is: %.4e\n",
50 double sum_all_ions = 0.;
51 double sum_metals = 0.;
54 eden_ions[nelem] = 0.;
55 for(
long ion=1; ion <= nelem+1; ion++ )
56 eden_ions[nelem] += ion*
dense.xIonDense[nelem][ion];
58 sum_all_ions += eden_ions[nelem];
60 sum_metals += eden_ions[nelem];
62 sum_all_ions +=
deut.xIonDense[1];
63 dense.EdenTrue += sum_all_ions;
79 dense.EdenTrue +=
gv.TotalEden*
gv.lgGrainElectrons;
87 " eden_sum zn: %.2f current: %.4e new true: %.4e ions: %.4e mole: %.4e"
88 " grain: %.4e extra: %.4e LaOTS: %.4e\n",
94 gv.TotalEden*
gv.lgGrainElectrons,
103 fprintf(
ioQQQ,
" eden_sum H -Ne:" );
104 else if( nelem == 10 )
105 fprintf(
ioQQQ,
" eden_sum Na-Ca:" );
106 else if( nelem == 20 )
107 fprintf(
ioQQQ,
" eden_sum Sc-Zn:" );
108 fprintf(
ioQQQ,
" %.4e", eden_ions[nelem] );
110 fprintf(
ioQQQ,
"\n" );