added density floor for levden array in macro_pops; see #425 and #438

agnwinds · Nov 22, 2018 · 0a05836 · 0a05836
1 parent 2fb29e6
commit 0a05836
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/source/macro_gov.c b/source/macro_gov.c
@@ -662,8 +662,8 @@ macro_pops (xplasma, xne)
                  index_lvl++)
             {
               /* JM Nov 18 -- if statement to prevent nan in fractional populations */
-              if (this_ion_density == 0.0)
-                xplasma->levden[config[index_lvl].nden] = 0.0;
+              if (this_ion_density <= DENSITY_MIN || populations[conf_to_matrix[index_lvl]] <= DENSITY_MIN)
+                xplasma->levden[config[index_lvl].nden] = DENSITY_MIN;
               else
                 xplasma->levden[config[index_lvl].nden] = populations[conf_to_matrix[index_lvl]] / this_ion_density;
             }