fix for sirocco-rt#498, don't overwrite macro_pops values

source/matrix_ion.c

@@ -307,16 +307,20 @@ matrix_ion_populations (xplasma, mode)
       {
         newden[nn] = xplasma->density[nn] / elem_dens[ion[nn].z];
       }
-
-      for (mm = 0; mm < nrows; mm++)    // inner loop over the elements of the population array
+
+      /* if the ion is "simple" then find it's calculated ionization state in populations array */
+      else 
       {
-        if (xion[mm] == nn)     // if this element contains the population of the ion is question
+
+        for (mm = 0; mm < nrows; mm++)    // inner loop over the elements of the population array
         {
-          newden[nn] = populations[mm]; // get the population
+          if (xion[mm] == nn)     // if this element contains the population of the ion is question
+          {
+            newden[nn] = populations[mm]; // get the population
+          }
         }
       }
 
-
       if (newden[nn] < DENSITY_MIN)     // this wil also capture the case where population doesnt have a value for this ion
         newden[nn] = DENSITY_MIN;
     }