diff --git a/matlab/TCV_IMAS/tcv_get_ids_nbi.m b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
index c0ce563920b6ded3847d4c86b90b68aa28561115..6e900f9a6aae9be9339b67a2ed4b564e1f3aed28 100644
--- a/matlab/TCV_IMAS/tcv_get_ids_nbi.m
+++ b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
@@ -122,7 +122,7 @@ for iunit=1:nb_units
     p_frac=gdat(shot,'\results::dnbi:fraction');
     ids_nbi.unit{iunit}.beam_power_fraction.data = p_frac.data(:,1:3)*0.01;
     ids_nbi.unit{iunit}.beam_power_fraction.time = p_frac.t;
-    i_frac = p_frac.data(:,1:3).*[1 2 3];
+    i_frac = p_frac.data(:,1:3).*repmat([1 2 3],size(p_frac.data,1),1); % to be compatible with older matlab version .*[1 2 3] not ok
     i_frac = i_frac.*1./sum(i_frac);
     ids_nbi.unit{iunit}.beam_current_fraction.data = i_frac(:,:);
     ids_nbi.unit{iunit}.beam_current_fraction.time = p_frac.t;