37 static double one=1.0;
39 struc.nzonePreviousIteration = -1;
53 for(
long ion=0; ion<nelem+2; ++ion )
54 mole.set_location( nelem, ion, &(
dense.xIonDense[nelem][ion]) );
74 (*diatom)->H2_Reset();
78 for(
long ion=0; ion<=nelem; ++ion )
79 Heavy.xLyaHeavy[nelem][ion] = 0;
82 if(
radius.StopRadius[0] > 0. )
84 for(
long j=0; j <
iterations.iter_malloc; j++ )
90 fprintf(
ioQQQ,
" PROBLEM stopping radius is <= inner radius. Bailing out.\n" );
108 rfield.ipEnergyBremsThin = 0;
109 rfield.EnergyBremsThin = 0.;
113 rfield.extin_mag_B_point = 0.;
114 rfield.extin_mag_V_point = 0.;
115 rfield.extin_mag_B_extended = 0.;
116 rfield.extin_mag_V_extended = 0.;
117 rfield.EnerGammaRay = 7676.;
123 tr->Emis().TauTrack().clear();
132 fprintf(
ioQQQ,
" InitSimPostparse sets AccelGravity %.3e lgDisk?%c\n",
145 for(
long nelem=ipISO; nelem <
LIMELM; ++nelem )
147 iso_sp[ipISO][nelem].Reset();
151 iso_sp[ipISO][nelem].numLevels_max = 0;
152 iso_sp[ipISO][nelem].nCollapsed_max = 0;
153 iso_sp[ipISO][nelem].n_HighestResolved_max = 0;
155 iso_sp[ipISO][nelem].numLevels_local = 0;
156 iso_sp[ipISO][nelem].nCollapsed_local = 0;
157 iso_sp[ipISO][nelem].n_HighestResolved_local = 0;
168 if(
iso_ctrl.lgPrintNumberOfLevels )
170 fprintf(
ioQQQ,
"\n\n Number of levels in ions treated by iso sequences.\n");
171 fprintf(
ioQQQ,
" ISO Element hi-n(l-resolved) #(l-resolved) n(collapsed)\n");
175 for(
long nelem=ipISO; nelem<
LIMELM; ++nelem )
178 fprintf(
ioQQQ,
" %s %s %4li %4li %4li \n",
181 iso_sp[ipISO][nelem].n_HighestResolved_max,
182 iso_sp[ipISO][nelem].numLevels_max-
iso_sp[ipISO][nelem].nCollapsed_max,
183 iso_sp[ipISO][nelem].nCollapsed_max);
void InitSimPostparse(void)