diff --git a/ForMoSA/nested_sampling/nested_sampling.py b/ForMoSA/nested_sampling/nested_sampling.py
index 72f3f32..f216a5e 100755
--- a/ForMoSA/nested_sampling/nested_sampling.py
+++ b/ForMoSA/nested_sampling/nested_sampling.py
@@ -213,9 +213,9 @@ def loglike(theta, theta_index, global_params, main_file, for_plot='no'):
                 if len(transm_obs_ns_u) != 0: # Merge the transmissions (if necessary)
                     transm_obs_ns_u = np.concatenate((transm_obs_ns_u, transm_obs[ind_spectro]))
                 if len(star_flx_obs_ns_u) != 0: # Merge star fluxes (if necessary)
-                    star_flx_obs_ns_u = np.concatenate((star_flx_obs_ns_u, star_flx_obs[ind_grid_spectro_sel]),axis=0)
+                    star_flx_obs_ns_u = np.concatenate((star_flx_obs_ns_u, star_flx_obs[ind_spectro]),axis=0)
                 if len(system_obs) != 0: # Merge systematics model (if necessary)
-                    system_obs_ns_u = np.concatenate((system_obs_ns_u, system_obs[ind_grid_spectro_sel]), axis=0)
+                    system_obs_ns_u = np.concatenate((system_obs_ns_u, system_obs[ind_spectro]), axis=0)
                 wav_obs_photo_ns_u = np.concatenate((wav_obs_photo_ns_u, wav_obs_photo[ind_photo]))
                 flx_obs_photo_ns_u = np.concatenate((flx_obs_photo_ns_u, flx_obs_photo[ind_photo]))
                 err_obs_photo_ns_u = np.concatenate((err_obs_photo_ns_u, err_obs_photo[ind_photo]))