From f96b7f508c8ead68f5f42438319f1afc767f673f Mon Sep 17 00:00:00 2001
From: Olivier Sauter <Olivier.Sauter@epfl.ch>
Date: Wed, 18 Sep 2019 20:41:54 +0200
Subject: [PATCH] repeat change repmat for dnbi
should simplify loop for units to avoid duplication of coding
---
matlab/TCV_IMAS/tcv_get_ids_nbi.m | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/matlab/TCV_IMAS/tcv_get_ids_nbi.m b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
index c0ce5639..6e900f9a 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;
--
GitLab