Merge pull request #207 from natashabatalha/dev

Critical hot fix for correct He mass
natashabatalha · Jul 8, 2024 · bbc2b67 · bbc2b67
2 parents b2b4521 + a82531e
commit bbc2b67
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/picaso/atmsetup.py b/picaso/atmsetup.py
@@ -330,7 +330,9 @@ def get_weights(self, molecule):
                     el, num = sep
                 #default isotope
                 if iso_num=='main':
-                    iso_num = list(ele[el].isotopes.keys())[0] 
+                    #select the main isotope off according to that with the highest relative abundance
+                    main_iso = np.argmax([ele[el].isotopes[i].abundance for i in ele[el].isotopes.keys()])
+                    iso_num = list(ele[el].isotopes.keys())[main_iso] 
                 totmass += ele[el].isotopes[iso_num].mass*float(num)
 
             weights[i]=totmass