%% Create the geometry 1 for the T-REX experiment 
% uses the upper ellipse and constant radius center column.

C_column_length=400;%mm
C_column_radius=10; %mm
O_electrode_length=400;%mm

% Outer vacuum vessel 
splid=1;
geomcells(splid).Dval=0; % V
geomcells(splid).order=3;
geomcells(splid).dim=2;
geomcells(splid).name='vacuum vessel';
geomcells(splid).type=0;
geomcells(splid).periodic=0;
rmax=50; %mm
lw=linspace(40,rmax,8);
rw=linspace(250,480,20);
uw=flip(linspace(0,rmax,30));
geomcells(splid).points=flip([250*ones(size(lw')) lw';
                   rw' rmax*ones(size(rw'));
                   480*ones(size(uw')) uw']/1e3);
    
 
% Central electrode
splid=2;
geomcells(splid).Dval=-20000; % V
geomcells(splid).order=3;
geomcells(splid).dim=2;
%splines(1).epsge=1e-6;
%splines(1).epsce=1e-6;
geomcells(splid).name='vacuum vessel';
geomcells(splid).type=0;
geomcells(splid).periodic=0;
lw=linspace(-.5,0,10)+C_column_radius;
rw=linspace(270.5,C_column_length-0.5,60);
uw=flip(linspace(0,1,10))*C_column_radius;
geomcells(splid).points=[ 270*ones(size(lw')) lw';
                   rw' C_column_radius*ones(size(rw'));
                   C_column_length*ones(size(uw')) uw']/1e3;

         
% Dielectric ring
splid=4;
geomcells(splid).Dval=0; % V
geomcells(splid).order=3;
geomcells(splid).dim=2;
%splines(1).epsge=1e-6;
%splines(1).epsce=1e-6;
geomcells(splid).name='vacuum vessel';
geomcells(splid).type=2;
geomcells(splid).periodic=0;
lw=flip(linspace(C_column_radius+.1,17.5,10));
rw=linspace(250,C_column_length-0.5,60);
uw=flip(linspace(0,C_column_radius,10));
geomcells(splid).points=[ 
                   250*ones(size(lw')) lw';
                   rw' C_column_radius*ones(size(rw'));
                   C_column_length*ones(size(uw')) uw']/1e3;
               
% Bottom plate
splid=3;
geomcells(splid).Dval=0; % V
geomcells(splid).order=3;
geomcells(splid).dim=2;
%splines(1).epsge=1e-6;
%splines(1).epsce=1e-6;
geomcells(splid).name='vacuum vessel';
geomcells(splid).type=0;
geomcells(splid).periodic=0;
rw=linspace(240,250,10);
lw=linspace(C_column_radius-.5,17.5,40);
uw=flip(linspace(240,250,10));
geomcells(splid).points=flip([
                         250*ones(size(lw')) lw';]/1e3);
               
               
               


               
% Outer electrode
splid=5;
geomcells(splid).Dval=0; % V
geomcells(splid).order=3;
geomcells(splid).dim=2;
%splines(1).epsge=1e-6;
%splines(1).epsce=1e-6;
geomcells(splid).name='vacuum vessel';
geomcells(splid).type=0;
geomcells(splid).periodic=0;
base=linspace(250,276,80);

%arc
radius=700; % mm
r_c=716.87399;
z_c=290.26217;

arcz=linspace(277,O_electrode_length);
arcr=r_c-sqrt(radius^2-(arcz-z_c).^2);

%tiltedwall

axialend=linspace(arcr(end)+0.5,39.5,15);
upperend=flip(linspace(250,arcz(end),60));

geomcells(splid).points=flip([base' 17*ones(size(base))';
                        arcz' arcr';
                        arcz(end)*ones(size(axialend')) axialend';
                        upperend' 40*ones(size(upperend'))]/1e3);
               
               
               
               %% Plots
f=figure;
for k=1:length(geomcells)
    plothandle=plot(geomcells(k).points(:,1), geomcells(k).points(:,2),'k-x','linewidth',1.5);
    hold on
    %geomcells(k).points=[geomcells(k).Z; geomcells(k).R];
    order=geomcells(k).order;
    knots=linspace(0,1,size(geomcells(k).points,1)-(order-2));
    knots=augknt(knots, order);
    sizec=size(geomcells(k).points,1);
    order=length(knots)-sizec(end);
    coeffs=geomcells(k).points';
    pp=spmak(knots,coeffs);
    s=linspace(0,1,1000);
    fittedpos=fnval(pp,s);
    plot(fittedpos(1,:),fittedpos(2,:),'x-')
end


legend(plothandle,{'Gun geometry'},'location','southwest')
f.PaperUnits='centimeters';
f.PaperSize=[12,8];
xlabel('z [m]')
ylabel('r [m]')

axis equal

savegeomtoh5('exp_arc_with_vessel_geom.h5',geomcells,1e-2,true);
% print(f,name,'-dpdf','-fillpage')
% savefig(f,name)
% set(f, 'Color', 'w');
% export_fig(f,name,'-eps')
hold off