% 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('7.5','30','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');

% Analyzed geometry
clear s
s.objs={g1};
s.name={'CYL1'};
s.tags={'g1'};

fem.draw=struct('s',s);
fem.geom=geomcsg(fem);
g2=block3('1','7.5','15','base','corner','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g3=block2('1','7.5','15','base','corner','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g4=block3('0.1','7.5','15','base','corner','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g5=block3('0.2','7.5','15','base','corner','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g6=block3('0.1','7.5','15','base','corner','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g7=block3('0.1','7.5','15','base','corner','pos',{'0','0','0'},'axis',{'0','0','15'},'rot','0');
g8=block3('0.1','7.5','15','base','corner','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g9=block3('0.1','7.5','15','base','corner','pos',{'0','0','15'},'axis',{'0','0','1'},'rot','0');
g10=block3('0.1','7.5','15','base','corner','pos',{'0','0','15'},'axis',{'0','0','1'},'rot','180');
g11=block3('0.1','7.5','15','base','corner','pos',{'0','0','15'},'axis',{'0','0','1'},'rot','0');
g12=block3('0.1','7.5','15','base','corner','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g13=block3('0.1','7.5','15','base','corner','pos',{'0','0','15'},'axis',{'0','0','1'},'rot','0');

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

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

% (Default values are not included)

% 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';
bnd.ind = [1,1,1,1,1,1];
appl.bnd = bnd;
fem.appl{1} = appl;
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,g13};
s.name={'CYL1','BLK1'};
s.tags={'g1','g13'};

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 $)

% Geometry
g2=cylinder3('7.5','19','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g3=block3('0.1','7.5','15','base','corner','pos',{'0','0','4'},'axis',{'0','0','1'},'rot','0');

% Geometry objects
clear s
s.objs={g2,g3};
s.name={'CYL1','BLK1'};
s.tags={'g2','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 $)

% Constants
fem.const = {};

% 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;

% Load geometry and mesh from component
[garr geom mesh]=loadgeommesh('/home/netcat/comsol35/bin/probe.mph');
[g1,g4,g5,g6,g7]=deal(garr{:});

% Geometry 2
flclear fem

% Geometry objects
clear s
s.objs={g1,g4,g5,g6};
s.name={'CYL2','CYL1','EXT1','CYL3'};
s.tags={'g1','g4','g5','g6'};

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

% Geometry 3
flclear fem

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

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

% Geometry 2
g8=cylinder3('0.205','2.07','pos',{'0.635','0.635','-0.4'},'axis',{'0','0','1'},'rot','0');
g9=cylinder3('0.065','2.36','pos',{'0.635','0.635','-0.4'},'axis',{'0','0','1'},'rot','0');
g5=rotate(g5,1.5707963267948966,[1,0,0],[0,0,0]);
g6=rotate(g6,1.5707963267948966,[1,0,0],[0,0,0]);
g8=rotate(g8,1.5707963267948966,[1,0,0],[0,0,0]);
g9=rotate(g9,1.5707963267948966,[1,0,0],[0,0,0]);
[g10,g11,g12,g13]=geomcopy({g5,g6,g8,g9});
g5=rotate(g5,-1.5707963267948966,[1,0,0],[0,0,0]);
g6=rotate(g6,-1.5707963267948966,[1,0,0],[0,0,0]);
g8=rotate(g8,-1.5707963267948966,[1,0,0],[0,0,0]);
g9=rotate(g9,-1.5707963267948966,[1,0,0],[0,0,0]);
g5=rotate(g5,1.5707963267948966,[1,0,0],[0,0,0]);
g6=rotate(g6,1.5707963267948966,[1,0,0],[0,0,0]);
g8=rotate(g8,1.5707963267948966,[1,0,0],[0,0,0]);
g9=rotate(g9,1.5707963267948966,[1,0,0],[0,0,0]);
g5=rotate(g5,1.5707963267948966,[0,0,-1],[0,0,0]);
g6=rotate(g6,1.5707963267948966,[0,0,-1],[0,0,0]);
g8=rotate(g8,1.5707963267948966,[0,0,-1],[0,0,0]);
g9=rotate(g9,1.5707963267948966,[0,0,-1],[0,0,0]);
[g18,g19,g20,g21]=geomcopy({g5,g6,g8,g9});

% Geometry 1
[g30,g31,g32,g33]=geomcopy({g18,g19,g20,g21});
g30=move(g30,[7.67,0,9.5]);
g31=move(g31,[7.67,0,9.5]);
g32=move(g32,[7.67,0,9.5]);
g33=move(g33,[7.67,0,9.5]);

% Geometry 2
g5=rotate(g5,3.141592653589793,[0,0,1],[0,0,0]);
g6=rotate(g6,3.141592653589793,[0,0,1],[0,0,0]);
g8=rotate(g8,3.141592653589793,[0,0,1],[0,0,0]);
g9=rotate(g9,3.141592653589793,[0,0,1],[0,0,0]);
[g34,g35,g36,g37]=geomcopy({g5,g6,g8,g9});

% Geometry 1
[g38,g39,g40,g41]=geomcopy({g34,g35,g36,g37});
g38=move(g38,[-7.67,0,9.5]);
g39=move(g39,[-7.67,0,9.5]);
g40=move(g40,[-7.67,0,9.5]);
g41=move(g41,[-7.67,0,9.5]);
xfem.fem{3}=fem;
fem=xfem.fem{1};

% Geometry objects
clear s
s.objs={g2,g3,g30,g31,g32,g33,g38,g39,g40,g41};
s.name={'CYL1','BLK1','EXT1','CYL2','CYL3','CYL4','EXT2','CYL5', ...
  'CYL6','CYL7'};
s.tags={'g2','g3','g30','g31','g32','g33','g38','g39','g40','g41'};

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

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

% Geometry objects
clear s
s.objs={g5,g6,g8,g9};
s.name={'EXT1','CYL3','CYL1','CYL2'};
s.tags={'g5','g6','g8','g9'};

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

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

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

fem.draw=struct('s',s);
xfem.fem{3}=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';
bnd.ind = [1,1,1,1,1,1];
appl.bnd = bnd;
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};

% 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{2} = fem;

fem=xfem.fem{3};
fem.sdim = {'x','y'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{3} = 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
g30=move(g30,[0,0.635,0]);
g31=move(g31,[0,0.635,0]);
g32=move(g32,[0,0.635,0]);
g33=move(g33,[0,0.635,0]);
g38=move(g38,[0,-0.635,0]);
g39=move(g39,[0,-0.635,0]);
g40=move(g40,[0,-0.635,0]);
g41=move(g41,[0,-0.635,0]);
fem=xfem.fem{1};

% Analyzed geometry
clear s
s.objs={g2,g3,g30,g31,g32,g33,g38,g39,g40,g41};
s.name={'CYL1','BLK1','EXT1','CYL2','CYL3','CYL4','EXT2','CYL5', ...
  'CYL6','CYL7'};
s.tags={'g2','g3','g30','g31','g32','g33','g38','g39','g40','g41'};

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={g5,g6,g8,g9};
s.name={'EXT1','CYL3','CYL1','CYL2'};
s.tags={'g5','g6','g8','g9'};

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

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

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

fem.draw=struct('s',s);
xfem.fem{3}=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.inport = {0,0,1};
bnd.type = {'E0','cont','port'};
bnd.modespec = {'userdef','userdef','coaxial'};
bnd.ind = [1,1,1,1,2,2,1,2,2,1,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2, ...
  2,2,2,2,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,1,2,2,2,2, ...
  2,2,2,2,1,2,2,1,2,2,2,2,2,2,2,1,2,2,2,1,1,2,1,1,1,2,1,2,2,2,2,2,2,2,2, ...
  2,2,2,2,2,2,2,1,1,2,2,1,1,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2, ...
  1,2,2,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,2,2,2,2,2,2,1,2,2,2,2,1,1,3,1];
appl.bnd = bnd;
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};

% 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{2} = fem;

fem=xfem.fem{3};
fem.sdim = {'x','y'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{3} = 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={g5,g6,g8,g9};
s.name={'EXT1','CYL3','CYL1','CYL2'};
s.tags={'g5','g6','g8','g9'};

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

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

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

fem.draw=struct('s',s);
xfem.fem{3}=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};
bnd.type = {'E0','cont','port'};
bnd.modespec = {'userdef','userdef','coaxial'};
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,1,1,1,1,1,1,2, ...
  2,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,1,1,2,2,1,2,2,2,2, ...
  1,1,1,1,1,2,2,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,2,2,2, ...
  2,1,1,2,1,1,2,1,1,2,2,1,1,2,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,1,1,1,1,1,1,1,1,1,1,1,2,1,1,2,1,1,2,1,1,1,1,1,1,1,1,1,3,1];
appl.bnd = bnd;
clear equ
equ.sigma = {0,'10^-25',4.8};
equ.epsilonr = {1,2.1,80.10};
equ.mur = {1,'1.2567*10^-6',0.999992};
equ.usage = {0,1,1};
equ.ind = [1,2,1,1,2,1,3,2,1,1,1,1,1,2,1,2,1,1,1,2,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};

% 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{2} = fem;

fem=xfem.fem{3};
fem.sdim = {'x','y'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{3} = 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={g5,g6,g8,g9};
s.name={'EXT1','CYL3','CYL1','CYL2'};
s.tags={'g5','g6','g8','g9'};

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

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

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

fem.draw=struct('s',s);
xfem.fem{3}=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};
bnd.type = {'E0','cont','port'};
bnd.modespec = {'userdef','userdef','coaxial'};
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,1,1,1,1,1,1,2, ...
  2,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,1,1,2,2,1,2,2,2,2, ...
  1,1,1,1,1,2,2,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,2,2,2, ...
  2,1,1,2,1,1,2,1,1,2,2,1,1,2,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,1,1,1,1,1,1,1,1,1,1,1,2,1,1,2,1,1,2,1,1,1,1,1,1,1,1,1,3,1];
appl.bnd = bnd;
clear equ
equ.sigma = {0,'10^-25',4.8};
equ.epsilonr = {1,2.1,80.1};
equ.mur = {1,'1.2567*10^-6',0.999992};
equ.usage = {0,1,1};
equ.ind = [1,2,1,1,2,1,3,2,1,1,1,1,1,2,1,2,1,1,1,2,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};

% 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{2} = fem;

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

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

% Multiphysics
xfem=multiphysics(xfem);