Newer
Older
clear all
close all
clc
% This script will complete the IDS of CORSICA with the missing fields recomputing the missing quantities from available information
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';
data_path_out = '/NoTivoli/carpanes/LIU_RAP_ITER/CORSICA_ids/CORSICA_130510_extended.mat';
% Load the IDS file
IDS = load(data_path_in);
G = G_ITER_generator;
% 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
G = liug(G,P);
L = liuc(P,G);
% Extract information from CORSICA IDS
[LC,LXC,LYC] = LLXLY_IDS_corsica(L, IDS);
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
%% -------------------- ADD the missing information on the existing IDS CORSICA -----------------------
%% Consider each circuit to be composed by only one coil
% TODO NEED TO ADD THE VERTICAL STABILIZATION COIL STILL
tmp = data_coils();
for ii = 1:numel(tmp.names)
for jj=1:numel(IDS.pf_active.coil)
% Remove strange character in coil names
IDS.pf_active.coil{jj}.name = regexprep(IDS.pf_active.coil{jj}.name,'[\n\r]+','');
if strcmp(tmp.names{ii},IDS.pf_active.coil{jj}.name)
IDS.pf_active.coil{jj}.element{1}.geometry.geometry_type = 2; % Rectangle description
IDS.pf_active.coil{jj}.element{1}.geometry.rectangle.r = tmp.R(ii);
IDS.pf_active.coil{jj}.element{1}.geometry.rectangle.z = tmp.Z(ii);
IDS.pf_active.coil{jj}.element{1}.geometry.rectangle.width = tmp.dR(ii);
IDS.pf_active.coil{jj}.element{1}.geometry.rectangle.height = tmp.dZ(ii);
IDS.pf_active.coil{jj}.element{1}.turns_with_sign = tmp.N(ii);
end
end
end
%% Add the data to the circuit
Ncoils = numel(IDS.pf_active.coil);
for ii=1:Ncoils-2
IDS.pf_active.circuit{ii}.name = IDS.pf_active.coil{ii}.name;
IDS.pf_active.supply{ii}.name = IDS.pf_active.coil{ii}.name;
IDS.pf_active.circuit{ii}.current.data = IDS.pf_active.coil{ii}.current.data/IDS.pf_active.coil{ii}.element{1}.turns_with_sign;
IDS.pf_active.circuit{ii}.current.time = IDS.pf_active.coil{ii}.current.time;
IDS.pf_active.circuit{ii}.connections = zeros(2, 4*Ncoils);
IDS.pf_active.circuit{ii}.connections(1, 2*ii) =1;
IDS.pf_active.circuit{ii}.connections(2, 2*ii-1) = 1;
IDS.pf_active.circuit{ii}.connections(1, 2*Ncoils + 2*ii-1) =1;
IDS.pf_active.circuit{ii}.connections(2, 2*Ncoils + 2*ii) =1;
end
% Vertical stabilization coil needs to be added still
for ii=(Ncoils-1):Ncoils
IDS.pf_active.circuit{ii}.name = IDS.pf_active.coil{ii}.name;
IDS.pf_active.supply{ii}.name = IDS.pf_active.coil{ii}.name;
IDS.pf_active.circuit{ii}.name = IDS.pf_active.coil{ii}.name;
IDS.pf_active.circuit{ii}.current.data = zeros(LXC.nt, 1);
IDS.pf_active.circuit{ii}.current.time = LXC.t;
IDS.pf_active.circuit{ii}.connections = zeros(2, 4*Ncoils);
IDS.pf_active.circuit{ii}.connections(1, 2*ii) =1;
IDS.pf_active.circuit{ii}.connections(2, 2*ii-1) = 1;
IDS.pf_active.circuit{ii}.connections(1, 2*Ncoils + 2*ii-1) =1;
IDS.pf_active.circuit{ii}.connections(2, 2*Ncoils + 2*ii) =1;
end
%% Limiter description
tmp = data_limiter();
IDS.wall.description_2d{1}.limiter.unit{1}.outline.r = tmp.r;
IDS.wall.description_2d{1}.limiter.unit{1}.outline.z = tmp.z;
%% Vessel description
% Understand what I need to do for the double layer vessel
%% -------------- Synthetic diagnostics------------ Need to be recomputed from CORSICA flux map
IDS.magnetics.method{1}.ip.time = LXC.t;
IDS.magnetics.method{1}.ip.data = LXC.Ip;
%% Ff
tmp = data_Ff();
for ii=1:numel(tmp.r)
IDS.magnetics.flux_loop{ii}.position{1}.r = tmp.r(ii);
IDS.magnetics.flux_loop{ii}.position{1}.z = tmp.z(ii);
IDS.magnetics.flux_loop{ii}.name = tmp.name{ii};
end
for ii=1:numel(tmp.r)
IDS.magnetics.flux_loop{ii}.flux.data = -LXC.Ff(ii,:)';
IDS.magnetics.flux_loop{ii}.flux.time = LXC.t;
end
%% Bm
tmp = data_Bm();
for ii=1:numel(tmp.name)
IDS.magnetics.bpol_probe{ii}.position.r = tmp.r(ii);
IDS.magnetics.bpol_probe{ii}.position.z = tmp.z(ii);
IDS.magnetics.bpol_probe{ii}.poloidal_angle = tmp.am(ii);
IDS.magnetics.bpol_probe{ii}.name = tmp.name{ii};
end
for ii=1:numel(tmp.name)
IDS.magnetics.bpol_probe{ii}.field.data = LXC.Bm(ii,:)';
IDS.magnetics.bpol_probe{ii}.field.time = LXC.t;
end
%% Ft
IDS.magnetics.method{1}.diamagnetic_flux.data = -LXC.Ft;
IDS.magnetics.method{1}.diamagnetic_flux.time = LXC.t;
%% rBt
IDS.tf.time = LXC.t;
IDS.tf.b_field_tor_vacuum_r.time = LXC.t;
IDS.tf.b_field_tor_vacuum_r.data = LXC.rBt;
%% Store the resulting data
save(data_path_out, '-struct', 'IDS')
fprintf('\n wrote file %s \n', data_path_out);