Small fixes

species/plot/plot_mcmc.py

@@ -404,7 +404,7 @@ def plot_posterior(
             o_h_ratio = np.zeros(samples.shape[0])
             c_o_ratio = np.zeros(samples.shape[0])
 
-            for sample_item in samples:
+            for i,sample_item in enumerate(samples):
                 abund_dict = {}
                 for line_item in line_species:
                     abund_dict[line_item] = sample_item[abund_index[line_item]]

species/util/radtrans_util.py

@@ -195,7 +195,11 @@ def retrieval_spectrum(
             model_param["albedo"] = sample[indices["albedo"]]
 
     elif len(cloud_species) > 0:
-        model_param["fsed"] = sample[indices["fsed"]]
+        if "fsed" in indices:
+            model_param["fsed"] = sample[indices["fsed"]]
+        else:
+            for item in cloud_species:
+                model_param[f"fsed_{item}"] = sample[indices[f"fsed_{item}"]]
         model_param["sigma_lnorm"] = sample[indices["sigma_lnorm"]]
 
         if "kzz" in indices:

species/util/retrieval_util.py

@@ -900,6 +900,8 @@ def calc_spectrum_clear(
             pressure,
             Pquench_carbon=p_quench,
         )
+        # Extract the mean molecular weight
+        mmw = abund_in["MMW"]
 
     elif chemistry == "free":
         # Free abundances