% COMSOL Multiphysics Model M-file
% Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $)

flclear fem

% COMSOL version
clear vrsn
vrsn.name = 'COMSOL 3.5';
vrsn.ext = '';
vrsn.major = 0;
vrsn.build = 494;
vrsn.rcs = '$Name:  $';
vrsn.date = '$Date: 2008/09/19 16:09:48 $';
fem.version = vrsn;

% Geometry
g1=cylinder3('0.075','0.19','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g2=block3('0.001','0.075','0.15','base','corner','pos',{'0','0','0'},'axis',{'0','0','0.04'},'rot','0');
g3=block3('0.0010','0.075','0.15','base','corner','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g3=move(g3,[0,-0.075,0]);
g3=move(g3,[-0.0005,0.001,0]);
g3=move(g3,[0,0,0.04]);

% Geometry objects
clear s
s.objs={g1,g3};
s.name={'CYL1','BLK1'};
s.tags={'g1','g3'};

fem.draw=struct('s',s);

% (Default values are not included)

% Application mode 1
clear appl
appl.mode.class = 'ElectromagneticWaves';
appl.module = 'RF';
appl.assignsuffix = '_rfw';
clear edg
edg.Im0 = {};
edg.I0 = {};
edg.name = {};
edg.type = {};
edg.ind = [];
appl.edg = edg;
clear bnd
bnd.inport = {};
bnd.eta = {};
bnd.murbnd = {};
bnd.kdir = {};
bnd.hport = {};
bnd.I0 = {};
bnd.srcpnt = {};
bnd.mbsource = {};
bnd.Z_ref = {};
bnd.type = {};
bnd.Js0 = {};
bnd.nummodetype = {};
bnd.E0 = {};
bnd.applmode = {};
bnd.nucutoff = {};
bnd.A0 = {};
bnd.index = {};
bnd.kper = {};
bnd.V0 = {};
bnd.rectmodetype = {};
bnd.modenum = {};
bnd.scsource = {};
bnd.lporttype = {};
bnd.nu0 = {};
bnd.matparams = {};
bnd.modespec = {};
bnd.epsilonrbnd = {};
bnd.srcaxis = {};
bnd.circmodenum = {};
bnd.chsrcdst = {};
bnd.Pport = {};
bnd.pertype = {};
bnd.symmetryz = {};
bnd.beta = {};
bnd.wport = {};
bnd.symmetryx = {};
bnd.symmetryy = {};
bnd.Es = {};
bnd.usermodetype = {};
bnd.name = {};
bnd.curofforbeta = {};
bnd.sigmabnd = {};
bnd.wavetype = {};
bnd.sparammethod = {};
bnd.portnr = {};
bnd.nbnd = {};
bnd.Pphase = {};
bnd.symtypez = {};
bnd.symtypey = {};
bnd.H0 = {};
bnd.symtypex = {};
bnd.ah = {};
bnd.ind = [];
appl.bnd = bnd;
clear equ
equ.Sd = {};
equ.magconstrel = {};
equ.srcpnt = {};
equ.M = {};
equ.S0 = {};
equ.gporder = {};
equ.coordOn = {};
equ.P = {};
equ.sigma = {};
equ.dinit = {};
equ.epsilonr = {};
equ.name = {};
equ.dr = {};
equ.rOn = {};
equ.cporder = {};
equ.mur = {};
equ.normfH = {};
equ.Br = {};
equ.n = {};
equ.init = {};
equ.Stype = {};
equ.R0 = {};
equ.elconstrel = {};
equ.fH = {};
equ.matparams = {};
equ.Dr = {};
equ.usage = {};
equ.srcaxis = {};
equ.user = {};
equ.ind = [];
appl.equ = equ;
fem.appl{1} = appl;
fem.sdim = {'x','y','z'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;

% ODE Settings
clear ode
clear units;
units.basesystem = 'SI';
ode.units = units;
fem.ode=ode;

% Multiphysics
fem=multiphysics(fem);
% COMSOL Multiphysics Model M-file
% Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $)

% Geometry

% Geometry objects
clear s
s.objs={g1,g3};
s.name={'CYL1','BLK1'};
s.tags={'g1','g3'};

fem.draw=struct('s',s);
% COMSOL Multiphysics Model M-file
% Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $)

% Using multiple geometries
flclear xfem
% Keep global fields
fields={'functions','const','event','globalexpr','xmesh','version'};
for ii=1:length(fields)
  if isfield(fem,fields{ii})
    xfem.(fields{ii})=fem.(fields{ii});
    fem=rmfield(fem,fields{ii});
  end
end
if isfield(fem,'sol')
  fem=rmfield(fem,'sol');
end
xfem.fem{1}=fem;

% Geometry 3
carr={curve2([0,0.2],[0,0],[1,1])};
g2=geomcoerce('curve',carr);
gg=geomedit(g2);
gg{1}=beziercurve2([0,0.0127],[0,0],[1,1]);
g4=geomedit(g2,gg);
carr={curve2([0,0],[0,0.01],[1,1])};
g5=geomcoerce('curve',carr);
gg=geomedit(g5);
gg{1}=beziercurve2([0,0],[0,0.0127],[1,1]);
g6=geomedit(g5,gg);
[g7]=geomcopy({g4});
[g8]=geomcopy({g7});
g8=move(g8,[0,0.0127]);
[g9]=geomcopy({g6});
[g10]=geomcopy({g9});
g10=move(g10,[0.0127,0]);
g11=geomcoerce('solid',{g4,g6,g8,g10});

% Geometry 2
g12=cylinder3('0.00205','0.015','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g13=cylinder3('0.00065','0.018','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');

% Geometry 3
g14=ellip2(0.0050,0.0050,'base','corner','pos',[-0.03,0]);
g15=ellip2('0.00205','0.00205','base','center','pos',{'0.000635','0.000635'},'rot','0');
g15=move(g15,[-6.35E-4,-0.008585]);
g16=ellip2('0.00205','0.00205','base','center','pos',{'0.00635','-0.00795'},'rot','0');
g17=ellip2('0.00205','0.00205','base','center','pos',{'0.00635','0.00635'},'rot','0');
g19=chamfer(g11,'dist1',0.0020,'point',[1]);
g20=chamfer(g19,'dist1',0.0020,'point',[2]);
g21=chamfer(g20,'dist1',0.0020,'point',[5]);
g22=chamfer(g21,'dist1',0.0020,'point',[7]);
g23=ellip2(0.0050,0.0050,'base','corner','pos',[-0.03,0]);
g24=ellip2('0.0013','0.0013','base','corner','pos',{'-0.03','0'},'rot','0');
g24=move(g24,[0.030700000000000005,0.0114]);
g25=ellip2('0.0013','0.0013','base','corner','pos',{'0.00185','0.0114'},'rot','0');
g26=ellip2('0.0013','0.0013','base','corner','pos',{'0.0018','0.0114'},'rot','0');
g27=ellip2('0.0013','0.0013','base','corner','pos',{'0.0017','0.0114'},'rot','0');
g28=ellip2('0.0013','0.0013','base','corner','pos',{'0.0017','0.0017'},'rot','0');
g29=ellip2('0.0013','0.0013','base','corner','pos',{'0.0017','0.0015'},'rot','0');
g30=ellip2('0.0013','0.0013','base','corner','pos',{'0.0015','0.0015'},'rot','0');
g31=ellip2('0.0013','0.0013','base','corner','pos',{'0.0013','0.0013'},'rot','0');
g32=ellip2('0.0013','0.0013','base','corner','pos',{'0.00125','0.00125'},'rot','0');
g33=ellip2('0.0013','0.0013','base','corner','pos',{'0.0012','0.0012'},'rot','0');
[g34]=geomcopy({g33});
[g35]=geomcopy({g34});
g35=move(g35,[0.0012,0.0115]);
g36=ellip2('0.0013','0.0013','base','corner','pos',{'0.0012','0.0115'},'rot','0');
g37=ellip2('0.0013','0.0013','base','corner','pos',{'0.0012','0.010'},'rot','0');
g38=ellip2('0.0013','0.0013','base','corner','pos',{'0.0012','0.0098'},'rot','0');
g39=ellip2('0.0013','0.0013','base','corner','pos',{'0.0012','0.0101'},'rot','0');
g40=ellip2('0.0013','0.0013','base','corner','pos',{'0.0012','0.0089'},'rot','0');
[g41]=geomcopy({g40});
[g42]=geomcopy({g41});
g42=move(g42,[0.0086,0]);
g43=ellip2('0.0013','0.0013','base','corner','pos',{'0.0089','0.0089'},'rot','0');
[g44]=geomcopy({g33});
[g46]=geomcopy({g33});
[g47]=geomcopy({g46});
g47=move(g47,[0.0086,0]);
g48=ellip2('0.0013','0.0013','base','corner','pos',{'0.0089','0.0012'},'rot','0');
g49=geomcomp({g22,g17,g33,g40,g43,g48},'ns',{'g22','g17','g33','g40','g43','g48'},'sf','g22-g17-g33-g40-g43-g48','edge','none');
g50=extrude(g49,'distance',[0.0017],'scale',[1;1],'displ',[0;0],'twist',[0],'face','none','wrkpln',[0 1 0;0 0 1;0 0 0]);

% Geometry 1

% Geometry 3
g51=extrude(g49,'distance',[0.0017],'scale',[1;1],'displ',[0;0],'twist',[0],'face','none','wrkpln',[0 1 0;0 0 1;0 0 0]);

% Geometry 2
g12=move(g12,[0.00635,0.00635,0]);
g13=move(g13,[0.00635,0.00635,0]);
g12=move(g12,[0,0,0.0017]);
g13=move(g13,[0,0,0.0017]);

% Geometry 4
g52=ellip2(0.1,0.1,'base','corner','pos',[0,-0.2]);
g53=ellip2('0.0025','0.0025','base','center','pos',{'0','0'},'rot','0');
g54=ellip2(2.5E-4,2.5E-4,'base','corner','pos',[-0.0055,0.0015]);
g55=ellip2('0.0023','0.0023','base','center','pos',{'0','0'},'rot','0');
g56=geomcomp({g53,g55},'ns',{'g53','g55'},'sf','g53-g55','edge','none');
g57=extrude(g56,'distance',[0.04],'scale',[1;1],'displ',[0;0],'twist',[0],'face','none','wrkpln',[0 1 0;0 0 1;0 0 0]);

% Geometry 2
g57=move(g57,[0.00635,0.00635,0.04]);
g57=move(g57,[0,0,-0.08]);
[g58,g59,g60,g61]=geomcopy({g12,g13,g51,g57});

% Geometry 1
[g62,g63,g64,g65]=geomcopy({g58,g59,g60,g61});
g62=move(g62,[0,0,0]);
g63=move(g63,[0,0,0]);
g64=move(g64,[0,0,0]);
g65=move(g65,[0,0,0]);
g62=rotate(g62,1.5707963267948966,[0,1,0],[0,0,0]);
g63=rotate(g63,1.5707963267948966,[0,1,0],[0,0,0]);
g64=rotate(g64,1.5707963267948966,[0,1,0],[0,0,0]);
g65=rotate(g65,1.5707963267948966,[0,1,0],[0,0,0]);
g62=move(g62,[-0.075,-0.003,0.095]);
g63=move(g63,[-0.075,-0.003,0.095]);
g64=move(g64,[-0.075,-0.003,0.095]);
g65=move(g65,[-0.075,-0.003,0.095]);
g62=move(g62,[-0.0017,-0.001,0]);
g63=move(g63,[-0.0017,-0.001,0]);
g64=move(g64,[-0.0017,-0.001,0]);
g65=move(g65,[-0.0017,-0.001,0]);
g62=move(g62,[0,-0.001,0]);
g63=move(g63,[0,-0.001,0]);
g64=move(g64,[0,-0.001,0]);
g65=move(g65,[0,-0.001,0]);
g62=move(g62,[0,-0.001,0]);
g63=move(g63,[0,-0.001,0]);
g64=move(g64,[0,-0.001,0]);
g65=move(g65,[0,-0.001,0]);
g62=move(g62,[0,-0.0005,0]);
g63=move(g63,[0,-0.0005,0]);
g64=move(g64,[0,-0.0005,0]);
g65=move(g65,[0,-0.0005,0]);
g62=move(g62,[0,-0.0002,0]);
g63=move(g63,[0,-0.0002,0]);
g64=move(g64,[0,-0.0002,0]);
g65=move(g65,[0,-0.0002,0]);
g62=move(g62,[0,0.0003,0]);
g63=move(g63,[0,0.0003,0]);
g64=move(g64,[0,0.0003,0]);
g65=move(g65,[0,0.0003,0]);
g62=move(g62,[0,0.0001,0]);
g63=move(g63,[0,0.0001,0]);
g64=move(g64,[0,0.0001,0]);
g65=move(g65,[0,0.0001,0]);
g62=move(g62,[0,-0.00005,0]);
g63=move(g63,[0,-0.00005,0]);
g64=move(g64,[0,-0.00005,0]);
g65=move(g65,[0,-0.00005,0]);
g62=move(g62,[0,0,0.00635]);
g63=move(g63,[0,0,0.00635]);
g64=move(g64,[0,0,0.00635]);
g65=move(g65,[0,0,0.00635]);
[g66,g67,g68,g69]=geomcopy({g58,g59,g60,g61});
g66=move(g66,[0,0,0]);
g67=move(g67,[0,0,0]);
g68=move(g68,[0,0,0]);
g69=move(g69,[0,0,0]);
g66=rotate(g66,-1.5707963267948966,[0,1,0],[0,0,0]);
g67=rotate(g67,-1.5707963267948966,[0,1,0],[0,0,0]);
g68=rotate(g68,-1.5707963267948966,[0,1,0],[0,0,0]);
g69=rotate(g69,-1.5707963267948966,[0,1,0],[0,0,0]);
g66=move(g66,[0.0767,0,0.101]);
g67=move(g67,[0.0767,0,0.101]);
g68=move(g68,[0.0767,0,0.101]);
g69=move(g69,[0.0767,0,0.101]);
g66=move(g66,[0,-0.005,0]);
g67=move(g67,[0,-0.005,0]);
g68=move(g68,[0,-0.005,0]);
g69=move(g69,[0,-0.005,0]);
g66=move(g66,[0,-0.00125,0]);
g67=move(g67,[0,-0.00125,0]);
g68=move(g68,[0,-0.00125,0]);
g69=move(g69,[0,-0.00125,0]);
g66=move(g66,[0,-0.00005,0]);
g67=move(g67,[0,-0.00005,0]);
g68=move(g68,[0,-0.00005,0]);
g69=move(g69,[0,-0.00005,0]);
g62=move(g62,[0,0.00000000000000000000003,0]);
g63=move(g63,[0,0.00000000000000000000003,0]);
g64=move(g64,[0,0.00000000000000000000003,0]);
g65=move(g65,[0,0.00000000000000000000003,0]);
g68=move(g68,[0,0,-0.01735]);
g69=move(g69,[0,0,-0.01735]);
g66=move(g66,[0,0,-0.01735]);
g67=move(g67,[0,0,-0.01735]);
g68=move(g68,[0,0,0.005]);
g69=move(g69,[0,0,0.005]);
g66=move(g66,[0,0,0.005]);
g67=move(g67,[0,0,0.005]);
fem=xfem.fem{1};

% Geometry objects
clear s
s.objs={g1,g3,g62,g63,g64,g65,g68,g69,g66,g67};
s.name={'CYL1','BLK1','CYL2','CYL3','EXT1','EXT2','EXT3','EXT4', ...
  'CYL4','CYL5'};
s.tags={'g1','g3','g62','g63','g64','g65','g68','g69','g66','g67'};

fem.draw=struct('s',s);
xfem.fem{1}=fem;

% Geometry 2
flclear fem

% Geometry objects
clear s
s.objs={g12,g13,g51,g57};
s.name={'CYL1','CYL2','EXT1','EXT2'};
s.tags={'g12','g13','g51','g57'};

fem.draw=struct('s',s);
xfem.fem{2}=fem;

% Geometry 3
flclear fem

% Geometry objects
clear s
s.objs={g49};
s.name={'CO2'};
s.tags={'g49'};

fem.draw=struct('s',s);
xfem.fem{3}=fem;

% Geometry 4
flclear fem

% Geometry objects
clear s
s.objs={g56};
s.name={'CO1'};
s.tags={'g56'};

fem.draw=struct('s',s);
xfem.fem{4}=fem;

% (Default values are not included)

fem=xfem.fem{1};

% Application mode 1
clear appl
appl.mode.class = 'ElectromagneticWaves';
appl.module = 'RF';
appl.assignsuffix = '_rfw';
clear edg
edg.Im0 = {};
edg.I0 = {};
edg.name = {};
edg.type = {};
edg.ind = [];
appl.edg = edg;
clear bnd
bnd.inport = {};
bnd.eta = {};
bnd.murbnd = {};
bnd.kdir = {};
bnd.hport = {};
bnd.I0 = {};
bnd.srcpnt = {};
bnd.mbsource = {};
bnd.Z_ref = {};
bnd.type = {};
bnd.Js0 = {};
bnd.nummodetype = {};
bnd.E0 = {};
bnd.applmode = {};
bnd.nucutoff = {};
bnd.A0 = {};
bnd.index = {};
bnd.kper = {};
bnd.V0 = {};
bnd.rectmodetype = {};
bnd.modenum = {};
bnd.scsource = {};
bnd.lporttype = {};
bnd.nu0 = {};
bnd.matparams = {};
bnd.modespec = {};
bnd.epsilonrbnd = {};
bnd.srcaxis = {};
bnd.circmodenum = {};
bnd.chsrcdst = {};
bnd.Pport = {};
bnd.pertype = {};
bnd.symmetryz = {};
bnd.beta = {};
bnd.wport = {};
bnd.symmetryx = {};
bnd.symmetryy = {};
bnd.Es = {};
bnd.usermodetype = {};
bnd.name = {};
bnd.curofforbeta = {};
bnd.sigmabnd = {};
bnd.wavetype = {};
bnd.sparammethod = {};
bnd.portnr = {};
bnd.nbnd = {};
bnd.Pphase = {};
bnd.symtypez = {};
bnd.symtypey = {};
bnd.H0 = {};
bnd.symtypex = {};
bnd.ah = {};
bnd.ind = [];
appl.bnd = bnd;
clear equ
equ.Sd = {};
equ.magconstrel = {};
equ.srcpnt = {};
equ.M = {};
equ.S0 = {};
equ.gporder = {};
equ.coordOn = {};
equ.P = {};
equ.sigma = {};
equ.dinit = {};
equ.epsilonr = {};
equ.name = {};
equ.dr = {};
equ.rOn = {};
equ.cporder = {};
equ.mur = {};
equ.normfH = {};
equ.Br = {};
equ.n = {};
equ.init = {};
equ.Stype = {};
equ.R0 = {};
equ.elconstrel = {};
equ.fH = {};
equ.matparams = {};
equ.Dr = {};
equ.usage = {};
equ.srcaxis = {};
equ.user = {};
equ.ind = [];
appl.equ = equ;
fem.appl{1} = appl;
fem.sdim = {'x','y','z'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{1} = fem;

fem=xfem.fem{2};
fem.sdim = {'x','y','z'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{2} = fem;

fem=xfem.fem{3};
fem.sdim = {'x','y'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{3} = fem;

fem=xfem.fem{4};
fem.sdim = {'x','y'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{4} = fem;

% ODE Settings
clear ode
clear units;
units.basesystem = 'SI';
ode.units = units;
xfem.ode=ode;

% Multiphysics
xfem=multiphysics(xfem);
% COMSOL Multiphysics Model M-file
% Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $)

% Geometry 5
g2=ellip2(0.1,0.1,'base','corner','pos',[0,0]);
g4=ellip2('0.00205','0.00205','base','center','pos',{'0','0'},'rot','0');
g5=ellip2(2.5E-4,5.0E-4,'base','corner','pos',[-0.0045,0]);
g6=ellip2('0.00065','0.00065','base','center','pos',{'0','0'},'rot','0');
g7=geomcomp({g4,g6},'ns',{'g4','g6'},'sf','g4-g6','edge','none');
g8=extrude(g7,'distance',[0.015],'scale',[1;1],'displ',[0;0],'twist',[0],'face','none','wrkpln',[0 1 0;0 0 1;0 0 0]);

% Geometry 1
g8=rotate(g8,1.5707963267948966,[0,1,0],[0,0,0]);
g8=move(g8,[-0.075,0,0]);
g8=move(g8,[0,0,0.095]);
g63=move(g63,[0.0001,0,0]);
g8=move(g8,[0.0001,0,0]);
g8=move(g8,[0.0001,0,0]);
g8=move(g8,[0.0001,0,0]);
g8=move(g8,[0.0001,0,0]);
g67=move(g67,[-0.0003,0,0]);
fem=xfem.fem{1};

% Analyzed geometry
clear s
s.objs={g68,g67,g69,g64,g63,g1,g65,g3,g8};
s.name={'EXT3','CYL5','EXT4','EXT1','CYL3','CYL1','EXT2','BLK1', ...
  'EXT5'};
s.tags={'g68','g67','g69','g64','g63','g1','g65','g3','g8'};

fem.draw=struct('s',s);
fem.geom=geomcsg(fem);
g67=move(g67,[0.0001,0,0]);

% Analyzed geometry
clear s
s.objs={g68,g67,g69,g64,g63,g1,g65,g3,g8};
s.name={'EXT3','CYL5','EXT4','EXT1','CYL3','CYL1','EXT2','BLK1', ...
  'EXT5'};
s.tags={'g68','g67','g69','g64','g63','g1','g65','g3','g8'};

fem.draw=struct('s',s);
fem.geom=geomcsg(fem);
g67=move(g67,[0.0001,0,0]);

% Analyzed geometry
clear s
s.objs={g68,g67,g69,g64,g63,g1,g65,g3,g8};
s.name={'EXT3','CYL5','EXT4','EXT1','CYL3','CYL1','EXT2','BLK1', ...
  'EXT5'};
s.tags={'g68','g67','g69','g64','g63','g1','g65','g3','g8'};

fem.draw=struct('s',s);
fem.geom=geomcsg(fem);
g67=move(g67,[0,-0.00005,0]);
g68=move(g68,[0,-0.00005,0]);
g69=move(g69,[0,-0.00005,0]);

% Analyzed geometry
clear s
s.objs={g68,g67,g69,g64,g63,g1,g65,g3,g8};
s.name={'EXT3','CYL5','EXT4','EXT1','CYL3','CYL1','EXT2','BLK1', ...
  'EXT5'};
s.tags={'g68','g67','g69','g64','g63','g1','g65','g3','g8'};

fem.draw=struct('s',s);
fem.geom=geomcsg(fem);
g68=move(g68,[0,5.0E-5,0]);
g69=move(g69,[0,5.0E-5,0]);
g67=move(g67,[0,5.0E-5,0]);
g8=move(g8,[-0.0001,0,0]);

% Analyzed geometry
clear s
s.objs={g68,g67,g69,g64,g63,g1,g65,g3,g8};
s.name={'EXT3','CYL5','EXT4','EXT1','CYL3','CYL1','EXT2','BLK1', ...
  'EXT5'};
s.tags={'g68','g67','g69','g64','g63','g1','g65','g3','g8'};

fem.draw=struct('s',s);
fem.geom=geomcsg(fem);
g63=move(g63,[0.0001,0,0]);
g63=move(g63,[0.0001,0,0]);
[g9]=geomcopy({g8});
[g10]=geomcopy({g9});
g10=move(g10,[0.15,0,0]);
g10=move(g10,[-0.075,0,0]);
g10=move(g10,[0.06,0,0]);
g68=move(g68,[0,-0.00005,0]);
g67=move(g67,[0,-0.00005,0]);
g69=move(g69,[0,-0.00005,0]);
g10=move(g10,[0,-0.00005,0]);
g10=move(g10,[-0.0003,0,0]);
g10=move(g10,[-0.0003,0,0]);
g67=move(g67,[-0.0001,0,0]);
g67=move(g67,[-0.0001,0,0]);
g10=move(g10,[0,0.00005,0]);

% Analyzed geometry
clear s
s.objs={g68,g67,g69,g64,g63,g1,g65,g3,g8,g10};
s.name={'EXT3','CYL5','EXT4','EXT1','CYL3','CYL1','EXT2','BLK1', ...
  'EXT5','EXT6'};
s.tags={'g68','g67','g69','g64','g63','g1','g65','g3','g8','g10'};

fem.draw=struct('s',s);
fem.geom=geomcsg(fem);

% Removing geometries
xfem.fem={xfem.fem{1},xfem.fem{2},xfem.fem{3},xfem.fem{4}};

% Geometry 1

% Geometry objects
clear s
s.objs={g68,g67,g69,g64,g63,g1,g65,g3};
s.name={'EXT3','CYL5','EXT4','EXT1','CYL3','CYL1','EXT2','BLK1'};
s.tags={'g68','g67','g69','g64','g63','g1','g65','g3'};

fem.draw=struct('s',s);
xfem.fem{1}=fem;

% Geometry 2
fem=xfem.fem{2};

% Geometry objects
clear s
s.objs={g57,g51,g12,g13};
s.name={'EXT2','EXT1','CYL1','CYL2'};
s.tags={'g57','g51','g12','g13'};

fem.draw=struct('s',s);
xfem.fem{2}=fem;

% Geometry 3
fem=xfem.fem{3};

% Geometry objects
clear s
s.objs={g49};
s.name={'CO2'};
s.tags={'g49'};

fem.draw=struct('s',s);
xfem.fem{3}=fem;

% Geometry 4
fem=xfem.fem{4};

% Geometry objects
clear s
s.objs={g56};
s.name={'CO1'};
s.tags={'g56'};

fem.draw=struct('s',s);
xfem.fem{4}=fem;

% (Default values are not included)

fem=xfem.fem{1};

% Application mode 1
clear appl
appl.mode.class = 'ElectromagneticWaves';
appl.module = 'RF';
appl.gporder = 4;
appl.cporder = 2;
appl.sshape = 2;
appl.assignsuffix = '_rfw';
clear bnd
bnd.type = {'E0','cont'};
bnd.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1, ...
  1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,1,1,1,2,1,1,2,2,1,1,1,1,1,1,1,1, ...
  1,1,1,1,1,1,1,1,1,1,2,2,2,1,1,2,1,1,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, ...
  1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1,1];
appl.bnd = bnd;
clear equ
equ.sigma = {0,0.005,0.0000000000000000000000001};
equ.epsilonr = {1,80.1,2.1};
equ.mur = {1,0.999992,0.0000012567};
equ.usage = {0,1,1};
equ.ind = [1,1,2,3,1,1,1,3,1,1];
appl.equ = equ;
fem.appl{1} = appl;
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{1} = fem;

fem=xfem.fem{2};
fem.sdim = {'x','y','z'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{2} = fem;

fem=xfem.fem{3};
fem.sdim = {'x','y'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{3} = fem;

fem=xfem.fem{4};
fem.sdim = {'x','y'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{4} = fem;

% ODE Settings
clear ode
clear units;
units.basesystem = 'SI';
ode.units = units;
xfem.ode=ode;

% Multiphysics
xfem=multiphysics(xfem);
% COMSOL Multiphysics Model M-file
% Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $)

% Geometry 1
g2=cylinder3('0.00205','0.15','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g4=cylinder3('0.00205','0.015','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g4=rotate(g4,1.5707963267948966,[0,1,0],[0,0,0]);
g4=move(g4,[-0.0747,-3.9976955110403E-19,0.095]);
fem=xfem.fem{1};

% Analyzed geometry
clear s
s.objs={g68,g67,g69,g64,g1,g63,g65,g3,g4};
s.name={'EXT3','CYL5','EXT4','EXT1','CYL1','CYL3','EXT2','BLK1', ...
  'CYL2'};
s.tags={'g68','g67','g69','g64','g1','g63','g65','g3','g4'};

fem.draw=struct('s',s);
fem.geom=geomcsg(fem);
g5=cylinder3('0.00205','0.015','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g5=rotate(g5,1.5707963267948966,[0,1,0],[0,0,0]);
g5=move(g5,[0.0747,-5.59719371545642E-18,0.095]);
g5=move(g5,[-0.014,0,0]);
g5=move(g5,[-0.002,0,0]);
g5=move(g5,[0.001,0,0]);

% Analyzed geometry
clear s
s.objs={g68,g67,g69,g64,g1,g63,g65,g3,g4,g5};
s.name={'EXT3','CYL5','EXT4','EXT1','CYL1','CYL3','EXT2','BLK1', ...
  'CYL2','CYL4'};
s.tags={'g68','g67','g69','g64','g1','g63','g65','g3','g4','g5'};

fem.draw=struct('s',s);
fem.geom=geomcsg(fem);
xfem.fem{1}=fem;

% Geometry 2
fem=xfem.fem{2};

% Geometry objects
clear s
s.objs={g57,g12,g51,g13};
s.name={'EXT2','CYL1','EXT1','CYL2'};
s.tags={'g57','g12','g51','g13'};

fem.draw=struct('s',s);
xfem.fem{2}=fem;

% Geometry 3
fem=xfem.fem{3};

% Geometry objects
clear s
s.objs={g49};
s.name={'CO2'};
s.tags={'g49'};

fem.draw=struct('s',s);
xfem.fem{3}=fem;

% Geometry 4
fem=xfem.fem{4};

% Geometry objects
clear s
s.objs={g56};
s.name={'CO1'};
s.tags={'g56'};

fem.draw=struct('s',s);
xfem.fem{4}=fem;

% (Default values are not included)

fem=xfem.fem{1};

% Application mode 1
clear appl
appl.mode.class = 'ElectromagneticWaves';
appl.module = 'RF';
appl.gporder = 4;
appl.cporder = 2;
appl.assignsuffix = '_rfw';
clear bnd
bnd.type = {'E0','cont'};
bnd.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1, ...
  1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,1,1,1,2,1,1,2,2,2,1,1,1,1,1,1,1, ...
  1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,1,1,2,1,1,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1, ...
  1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1,1];
appl.bnd = bnd;
clear equ
equ.sigma = {0,0.005,0.0000000000000000000000001};
equ.epsilonr = {1,80.1,2.1};
equ.mur = {1,0.999992,0.0000012567};
equ.usage = {0,1,1};
equ.ind = [1,1,2,3,1,1,1,1,2,1,1,1];
appl.equ = equ;
fem.appl{1} = appl;
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{1} = fem;

fem=xfem.fem{2};
fem.sdim = {'x','y','z'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{2} = fem;

fem=xfem.fem{3};
fem.sdim = {'x','y'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{3} = fem;

fem=xfem.fem{4};
fem.sdim = {'x','y'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{4} = fem;

% ODE Settings
clear ode
clear units;
units.basesystem = 'SI';
ode.units = units;
xfem.ode=ode;

% Multiphysics
xfem=multiphysics(xfem);

% COMSOL Multiphysics Model M-file
% Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $)

% Geometry 2
fem=xfem.fem{2};

% Geometry objects
clear s
s.objs={g57,g51,g12,g13};
s.name={'EXT2','EXT1','CYL1','CYL2'};
s.tags={'g57','g51','g12','g13'};

fem.draw=struct('s',s);
xfem.fem{2}=fem;

% Geometry 3
fem=xfem.fem{3};

% Geometry objects
clear s
s.objs={g49};
s.name={'CO2'};
s.tags={'g49'};

fem.draw=struct('s',s);
xfem.fem{3}=fem;

% Geometry 4
fem=xfem.fem{4};

% Geometry objects
clear s
s.objs={g56};
s.name={'CO1'};
s.tags={'g56'};

fem.draw=struct('s',s);
xfem.fem{4}=fem;

% (Default values are not included)

fem=xfem.fem{1};

% Application mode 1
clear appl
appl.mode.class = 'ElectromagneticWaves';
appl.module = 'RF';
appl.gporder = 4;
appl.cporder = 2;
appl.sshape = 2;
appl.assignsuffix = '_rfw';
clear bnd
bnd.inport = {0,0,1,0,0,0};
bnd.type = {'E0','cont','port','port','port','port'};
bnd.usermodetype = {'TE','TE','TEM','TEM','TEM','TEM'};
bnd.portnr = {1,1,1,2,3,4};
bnd.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1, ...
  1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,3,1,1,2,1,1,2,2,2,1,1,1,1,4,1,1, ...
  1,1,1,1,1,1,5,1,1,1,1,2,2,2,2,1,1,2,1,1,2,2,6,1,1,1,1,1,1,1,1,1,1,1,1, ...
  1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1,1];
appl.bnd = bnd;
clear equ
equ.sigma = {0,0.005,0.0000000000000000000000001};
equ.epsilonr = {1,80.1,2.1};
equ.mur = {1,0.999992,0.0000012567};
equ.usage = {0,1,1};
equ.ind = [1,1,2,3,1,1,1,1,3,1,1,1];
appl.equ = equ;
fem.appl{1} = appl;
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{1} = fem;

fem=xfem.fem{2};
fem.sdim = {'x','y','z'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{2} = fem;

fem=xfem.fem{3};
fem.sdim = {'x','y'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{3} = fem;

fem=xfem.fem{4};
fem.sdim = {'x','y'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{4} = fem;

% ODE Settings
clear ode
clear units;
units.basesystem = 'SI';
ode.units = units;
xfem.ode=ode;

% Multiphysics
xfem=multiphysics(xfem);