species_mass problematic

Nicholaswogan · Oct 2, 2024 · 94c1ba3 · 94c1ba3
1 parent f177ea3
commit 94c1ba3
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/equilibrate.f90 b/src/equilibrate.f90
@@ -111,7 +111,10 @@ function create_ChemEquiAnalysis(thermopath, atoms, species, err) result(cea)
 
     cea%atoms_names = cea%dat%names_atoms(1:size(cea%dat%names_atoms)-1) ! skip E
     cea%species_names = cea%dat%names_reactants_orig
-    cea%species_mass = cea%dat%mol_weight
+    allocate(cea%species_mass(size(cea%dat%mol_weight)))
+    do i = 1,size(cea%species_mass)
+      cea%species_mass(i) = cea%dat%mol_weight(cea%dat%id_reactants(i,1))
+    enddo
 
     allocate(cea%condensate_names(cea%dat%N_cond))
     allocate(cea%gas_names(cea%dat%N_gas))