tcv_ids_supply.m

function [ids_struct_out,ids_struct_out_description] =  tcv_ids_supply(shot, ids_structures, gdat_params, varargin)
%
% [ids_struct_out,ids_struct_out_description] =  tcv_get_ids_supply(shot, ids_structures, gdat_params, varargin);
%
% Get power supply/coils names for each circuit
%
% error_bar: from gdat_params.error_bar
%            'delta' (default): error_bar to be added inserted in "upper" only as mentioned in description
%            'delta_with_lower' : error_bar (abs) inserted in both lower and upper
%            'added': value already added to data: upper/lower = data +/- error_bar
%

error_bar = 'delta';
if exist('gdat_params') && isfield(gdat_params,'error_bar') && ~isempty(gdat_params.error_bar)
  error_bar = gdat_params.error_bar;
end

%% Get power supply/coils names for each circuit.
[tcv_circuit_info] =  tcv_ids_pf_active_definition();

% Preallocate memory and get data
ids_struct_out(1:tcv_circuit_info.ntotpowersupplies) = ids_structures;
for ii=1:tcv_circuit_info.ntotpowersupplies
  if shot ==-1
    warning('TCV:IDS:NoTimeData','no time data loaded for shot %d',shot);
    tmpdata.dim{1} = [];
    tmpdata.data = [];
    % empty dummy
    ids_struct_out_description{ii}.current =  sprintf('Not loaded for shot %d',shot);

  else
    tmpdata = gdat_tcv(shot,['' tcv_circuit_info.mds_paths{ii} '']); % Get current
    ids_struct_out_description{ii}.current = ['abs value from ' tmpdata.data_fullpath];
  end
  
  %ids_struct_out{ii}.current.data = abs(tmpdata.data); %%%%%%% FF:WHY ABS?? %%%%%%%
  ids_struct_out{ii}.current.data = tmpdata.data; %% FF removed abs()
  ids_struct_out{ii}.current.time = tmpdata.dim{1};
  ids_struct_out{ii}.name = tcv_circuit_info.power_supply_names{ii}{1};
  ids_struct_out_description{ii}.name = ...
      ['aa=tcv_ids_pf_active_definition, then from aa.power_supply_names{' num2str(ii) '}{1}'];
end

fixed_error_fraction = 0.03;
switch error_bar
 case 'delta'
  for ii=1:tcv_circuit_info.ntotpowersupplies
    ids_struct_out{ii}.current.data_error_upper = fixed_error_fraction.*max(abs(ids_struct_out{ii}.current.data)) ...
        .* ones(size(ids_struct_out{ii}.current.data));
    ids_struct_out_description{ii}.current_data_error_upper = ['from fixed_fraction (' num2str(fixed_error_fraction) ...
                    ') *max() in case: ' error_bar];
    ids_struct_out_description{ii}.current_data_error_lower = ['not provided since symmetric'];
  end
 case 'delta_with_lower'
  for ii=1:tcv_circuit_info.ntotcircuits
    ids_struct_out{ii}.current.data_error_upper = fixed_error_fraction.*max(abs(ids_struct_out{ii}.current.data)) ...
        .* ones(size(ids_struct_out{ii}.current.data));
    ids_struct_out{ii}.current.data_error_lower = ids_struct_out{ii}.current.data_error_upper;
    ids_struct_out_description{ii}.current_data_error_upper = ['from fixed_fraction (' num2str(fixed_error_fraction) ...
                    ') *max() in case: ' error_bar];
    ids_struct_out_description{ii}.current_data_error_lower = ids_struct_out_description{ii}.current_data_error_upper;
  end
 case 'added'
  for ii=1:tcv_circuit_info.ntotcircuits
    ids_struct_out{ii}.current.data_error_upper = ids_struct_out{ii}.current.data + ...
          fixed_error_fraction.*max(abs(ids_struct_out{ii}.current.data)).*ones(size(ids_struct_out{ii}.current.data));
    ids_struct_out{ii}.current.data_error_lower = ids_struct_out{ii}.current.data - ...
          fixed_error_fraction.*max(abs(ids_struct_out{ii}.current.data)).*ones(size(ids_struct_out{ii}.current.data));
    ids_struct_out_description{ii}.current_data_error_upper = ['from data + fixed_fraction (' num2str(fixed_error_fraction) ...
                    ') *max() in case: ' error_bar];
    ids_struct_out_description{ii}.current_data_error_lower = ['from data - fixed_fraction (' num2str(fixed_error_fraction) ...
                    ') *max() in case: ' error_bar];
  end
 otherwise
  error(['tcv_ids_bpol_loop: error_bar option not known: ' error_bar])
end