Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
function [ids_nbi,ids_nbi_description,varargout] = tcv_get_ids_nbi(shot,ids_nbi_empty, gdat_params,varargin);
%
% [ids_nbi,ids_nbi_description,varargout] = tcv_get_ids_nbi(shot,ids_nbi_empty,varargin);
%
%
% gdat_params: gdat_data.gdat_params to get all params passed from original call, in particular error_bar options
%
[ids_nbi, params_nbi] = tcv_ids_headpart(shot,ids_nbi_empty,'nbi','homogeneous_time',0,varargin{:});
% As a general rule, for a new substructure under the main ids, construct a local structure like:
% "global_quantities" with subfields being the relevant data to get and a local structure:
% "global_quantities_desc" which contains the same subfields themselves containing the gdat string aftre shot used
%
pow=gdat_tcv(shot,'powers');
nb_units = 2; % assume 2 units: 1st NBH and DNBI
ids_nbi.unit(1:nb_units) = ids_nbi.unit(1); % copy empty structure for all units, then fill in
for iunit=1:length(nb_units)
switch iunit
case 1
ids_nbi.unit{iunit}.identifier = 'NBH1';
ids_nbi.unit{iunit}.name = '25keV 1st NBH source';
ids_nbi.unit{iunit}.power_launched.time = pow.nbi.t;
ids_nbi.unit{iunit}.power_launched.data = pow.nbi.data;
case 2
ids_nbi.unit{iunit}.identifier = 'DNBI';
ids_nbi.unit{iunit}.name = 'diagnostic NBI';
% ids_nbi.unit{iunit}.power_launched.data = pow.nbi.data; % DNBI power time values and time
% ids_nbi.unit{iunit}.power_launched.time = ;
otherwise
warning(['iunit = ' num2str(iunit) ' not yet defined']);
end
end