From 0c9747a4e6279f312ab4dfee80c2effe6e316ca1 Mon Sep 17 00:00:00 2001 From: MatthieuR18 Date: Fri, 13 Sep 2024 14:32:47 +0200 Subject: [PATCH] MR: small index correction for the star_flx and systematics --- ForMoSA/nested_sampling/nested_sampling.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) 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]))