Skip to content
Snippets Groups Projects
Commit b7b79499 authored by Francesco Carpanese's avatar Francesco Carpanese Committed by Olivier Sauter
Browse files

Ongoing conversion from COCOS 11 to COCOS 17 for CORSICA

parent 0c8e9e49
No related branches found
No related tags found
1 merge request!41Add corsica liuqe complete ids from meq
...@@ -7,11 +7,18 @@ clc ...@@ -7,11 +7,18 @@ clc
data_path_in = '/NoTivoli/carpanes/LIU_RAP_ITER/CORSICA_ids/CORSICA_130510.mat'; data_path_in = '/NoTivoli/carpanes/LIU_RAP_ITER/CORSICA_ids/CORSICA_130510.mat';
%path_IDS = '/NoTivoli/carpanes/LIU_RAP_ITER/CORSICA_ids/CORSICA_130506.mat'; %path_IDS = '/NoTivoli/carpanes/LIU_RAP_ITER/CORSICA_ids/CORSICA_130506.mat';
data_path_out = '/NoTivoli/carpanes/LIU_RAP_ITER/CORSICA_ids/CORSICA_130510_extended.mat'; data_path_out = '/NoTivoli/carpanes/LIU_RAP_ITER/CORSICA_ids/CORSICA_130510_extended_COCOS17.mat';
% Load the IDS file % Load the IDS file
IDS = load(data_path_in); IDS = load(data_path_in);
%% Correct COCOS convention on original data
% Put homogenous = 0 on magnetic but 1 on the coils.
% Sign of the angle of the magnetic probes.
%%
G = G_ITER_generator; G = G_ITER_generator;
% Get the structure from generating files % Get the structure from generating files
P = liupiter(1000,'uuu',ones(G.na,1),'vvv',ones(G.nm,1), 'www', ones(G.nf,1)); % Some default setting that could be removed probably to use the liupiter directly P = liupiter(1000,'uuu',ones(G.na,1),'vvv',ones(G.nm,1), 'www', ones(G.nf,1)); % Some default setting that could be removed probably to use the liupiter directly
...@@ -118,12 +125,13 @@ IDS.tf.time = LXC.t; ...@@ -118,12 +125,13 @@ IDS.tf.time = LXC.t;
IDS.tf.b_field_tor_vacuum_r.time = LXC.t; IDS.tf.b_field_tor_vacuum_r.time = LXC.t;
IDS.tf.b_field_tor_vacuum_r.data = LXC.rBt; IDS.tf.b_field_tor_vacuum_r.data = LXC.rBt;
%% Convert from 11 to 17
IDS = ids_generic_cocos_nodes_transformation_symbolic(IDS.magnetics, 'magnetics', 11, 17);
IDS = ids_generic_cocos_nodes_transformation_symbolic(IDS.magnetics, 'pf_active', 11, 17);
IDS = ids_generic_cocos_nodes_transformation_symbolic(IDS.magnetics, 'tf', 11, 17);
%% Store the resulting data %% Store the resulting data
save(data_path_out, '-struct', 'IDS') save(data_path_out, '-struct', 'IDS')
fprintf('\n wrote file %s \n', data_path_out); fprintf('\n wrote file %s \n', data_path_out);
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment