Newer
Older
function ids_struct_out = tcv_get_ids_flux_loop(shot, ids_structures)
Francesco Carpanese
committed
Francesco Carpanese
committed
mdsopen(shot);
tmp = gdat(shot, 'tcv_idealloop("FLUX")');
names = tmp.dim{2};
time = tmp.dim{1};
data = tmp.data;
% Preallocate output structure
Nprobes = size(tmp.data,2);
Francesco Carpanese
committed
ids_struct_out(1:Nprobes) = ids_structures;
Francesco Carpanese
committed
Francesco Carpanese
committed
for ii=1:Nprobes
ids_struct_out{ii}.name = [names{ii}];
ids_struct_out{ii}.position{1}.r = mdsvalue('STATIC("R_F" )[$1]',ids_struct_out{ii}.name);
ids_struct_out{ii}.position{1}.z = mdsvalue('STATIC("Z_F" )[$1]',ids_struct_out{ii}.name);
ids_struct_out{ii}.flux.data = data(:,ii);
ids_struct_out{ii}.flux.time = time;
ids_struct_out{ii}.flux.data_error_upper = 0.0012;
ids_struct_out{ii}.flux.data_error_lower = 0.0012;
%ids_struct_out{ii}.position.phi = 0; % TO BE FOUND