Skip to content
Snippets Groups Projects
Commit f96b7f50 authored by Olivier Sauter's avatar Olivier Sauter
Browse files

repeat change repmat for dnbi

should simplify loop for units to avoid duplication of coding
parent 6201b761
No related branches found
No related tags found
1 merge request!16Imas codecamp develop
Pipeline #17920 passed
...@@ -122,7 +122,7 @@ for iunit=1:nb_units ...@@ -122,7 +122,7 @@ for iunit=1:nb_units
p_frac=gdat(shot,'\results::dnbi:fraction'); 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.data = p_frac.data(:,1:3)*0.01;
ids_nbi.unit{iunit}.beam_power_fraction.time = p_frac.t; 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); 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.data = i_frac(:,:);
ids_nbi.unit{iunit}.beam_current_fraction.time = p_frac.t; ids_nbi.unit{iunit}.beam_current_fraction.time = p_frac.t;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment