grid aborts due to element non-conservation
The attached grid crashes with c10_branch@4322 and newmole@4322. The situation is better with c10_branch since there only 4 grid points crash:
IONIZATION PARAMETER= -1.000000 LOG
METALS= 3.250000 LOG GRAINS
IONIZATION PARAMETER= -0.500000 LOG
METALS= 2.750000 LOG GRAINS
IONIZATION PARAMETER= -0.500000 LOG
METALS= 3.750000 LOG GRAINS
IONIZATION PARAMETER= 0.000000 LOG
METALS= 3.750000 LOG GRAINS
On newmole the following 8 grid points crash:
IONIZATION PARAMETER= -1.000000 LOG
METALS= 3.250000 LOG GRAINS
IONIZATION PARAMETER= -1.000000 LOG
METALS= 3.750000 LOG GRAINS
IONIZATION PARAMETER= -0.500000 LOG
METALS= 2.750000 LOG GRAINS
IONIZATION PARAMETER= -0.500000 LOG
METALS= 3.250000 LOG GRAINS
IONIZATION PARAMETER= -0.500000 LOG
METALS= 3.750000 LOG GRAINS
IONIZATION PARAMETER= 0.000000 LOG
METALS= 2.250000 LOG GRAINS
IONIZATION PARAMETER= 0.000000 LOG
METALS= 3.250000 LOG GRAINS
IONIZATION PARAMETER= 0.000000 LOG
METALS= 3.750000 LOG GRAINS
In all cases the aborts were due to non-conservation of nuclei or densities exceeding MAX_DENSITY. The code was compiled with g++ 4.4.3 and MPI enabled.