Skip to content
Snippets Groups Projects

Add corsica liuqe complete ids from meq

Merged Olivier Sauter requested to merge add_corsica_liuqe_complete_ids_from_meq into master
12 files
+ 2778
3
Compare changes
  • Side-by-side
  • Inline
Files
12
+ 61
0
 
function G = G_ITER_generator()
 
G = struct();
 
 
tmp = data_Bm();
 
 
G.rm = tmp.r;
 
G.zm = tmp.z;
 
G.am = tmp.am;
 
G.nm = numel(G.zm);
 
 
tmp = data_Ff();
 
G.rf = tmp.r;
 
G.zf = tmp.z;
 
G.nf = numel(G.zf);
 
 
tmp = data_limiter();
 
G.rl = tmp.r;
 
G.zl = tmp.z;
 
 
% TODO vessel is ready but will be added in a second moment
 
G.rv = [];
 
G.zv = [];
 
 
% Point for the gaps
 
tmp = data_rW();
 
G.rW = tmp.r;
 
G.zW = tmp.z;
 
G.oW = tmp.oW;
 
G.aW = tmp.aW;
 
 
%% Coils
 
tmp = data_coils();
 
R = tmp.R;
 
Z = tmp.Z;
 
dR = tmp.dR;
 
dZ = tmp.dZ;
 
N = tmp.N;
 
dima= tmp.names;
 
 
%% Sort by CORSICA ordering for the structure liux
 
% Not all che coils are used in this
 
% Matrix to be given to eval_trap
 
vvtmp = [R';Z';dR';dZ'];
 
% Sort for Corsica convention
 
vv = zeros(size(vvtmp));
 
sort_list = [3,4,2,5,1,6,7,8,9,10,11,12, 13, 14];
 
for ii=1:numel(sort_list)
 
vv(:,ii) = vvtmp(:,sort_list(ii));
 
end
 
N = N(sort_list);
 
% Compute the windings and connection matrix
 
[rw,zw,Twa] = gen_filament(vv,N,2);
 
G.dima = dima(sort_list);
 
 
G.rw = rw;
 
G.zw = zw;
 
G.Twa = Twa;
 
G.na = numel(R);
 
G.Tvs = [];
 
 
end
Loading