% 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=rect2('400e-9','400e-9','base','corner','pos',{'0','0'},'rot','0');
parr={point2(0,0)};
g2=geomcoerce('point',parr);
parr={point2(2.0E-7,2.0E-7)};
g3=geomcoerce('point',parr);
g4=rect2(1.0E-7,1.0E-7,'base','corner','pos',[-1.0E-7,4.0E-7]);
g5=rect2(4.0E-7,1.0E-7,'base','corner','pos',[0,4.0E-7]);
g6=rect2(1.0E-7,1.0E-7,'base','corner','pos',[4.0E-7,4.0E-7]);
g7=rect2(1.0E-7,4.0E-7,'base','corner','pos',[4.0E-7,0]);
g8=rect2(1.0E-7,1.0E-7,'base','corner','pos',[4.0E-7,-1.0E-7]);
g9=rect2(4.0E-7,1.0E-7,'base','corner','pos',[0,-1.0E-7]);
g10=rect2(1.0E-7,1.0E-7,'base','corner','pos',[-1.0E-7,-1.0E-7]);
g11=rect2(1.0E-7,4.0E-7,'base','corner','pos',[-1.0E-7,0]);

% Analyzed geometry
clear p s
p.objs={g3};
p.name={'PT1'};
p.tags={'g3'};

s.objs={g1,g4,g5,g6,g7,g8,g9,g10,g11};
s.name={'R1','R2','R3','R4','R5','R6','R7','R8','R9'};
s.tags={'g1','g4','g5','g6','g7','g8','g9','g10','g11'};

fem.draw=struct('p',p,'s',s);
fem.geom=geomcsg(fem);
g12=geomcomp({g1,g5,g6,g7,g8,g9},'ns',{'g1','g5','g6','g7','g8','g9'},'sf','g1+g5+g6+g7+g8+g9','edge','none');
g4=move(g4,[-5.0E-8,5.0E-8]);
g10=move(g10,[-5.0E-8,5.0E-8]);
g11=move(g11,[-5.0E-8,5.0E-8]);
g4=move(g4,[0,-5.0E-8]);
g10=move(g10,[0,-5.0E-8]);
g11=move(g11,[0,-5.0E-8]);
g4=move(g4,[0.5e-7,0]);
g10=move(g10,[0.5e-7,0]);
g11=move(g11,[0.5e-7,0]);
g13=geomcomp({g4,g10,g11},'ns',{'g4','g10','g11'},'sf','g4+g10+g11','edge','none');

% Create assembly and imprints
clear p s
p.objs={g3};
p.name={'PT1'};
p.tags={'g3'};

s.objs={g12,g13};
s.name={'CO1','CO2'};
s.tags={'g12','g13'};

fem.draw=struct('p',p,'s',s);
fem.geom=geomgroup(fem,'imprint','on',...
  'paircand',{'CO1','CO2'});

% Create assembly and imprints
clear p s
p.objs={g3};
p.name={'PT1'};
p.tags={'g3'};

s.objs={g12,g13};
s.name={'CO1','CO2'};
s.tags={'g12','g13'};

fem.draw=struct('p',p,'s',s);
fem.geom=geomgroup(fem,'imprint','on',...
  'paircand',{'CO1','CO2'});

% Analyzed geometry
clear p s
p.objs={g3};
p.name={'PT1'};
p.tags={'g3'};

s.objs={g12,g13};
s.name={'CO1','CO2'};
s.tags={'g12','g13'};

fem.draw=struct('p',p,'s',s);
fem.geom=geomgroup(fem,'imprint','off','paircand','none');

% Initialize mesh
fem.mesh=meshinit(fem, ...
                  'hauto',2);

% Refine mesh
fem.mesh=meshrefine(fem, ...
                    'mcase',0, ...
                    'rmethod','regular');

% Refine mesh
fem.mesh=meshrefine(fem, ...
                    'mcase',0, ...
                    'rmethod','regular');

% (Default values are not included)

% Application mode 1
clear appl
appl.mode.class = 'InPlaneWaves';
appl.module = 'RF';
appl.assignsuffix = '_rfwe';
clear prop
prop.inputvar='lambda';
appl.prop = prop;
clear bnd
bnd.type = {'SC','cont'};
bnd.ind = [1,1,1,2,1,2,1,1,1,1,1,1,1,2,1,2,1,2,1,2,2,2,2,1,1,1,1];
clear pair
pair.type = 'E';
pair.E0 = {{0;0;1}};
pair.pair = 'Pair 1';
bnd.pair = pair;
appl.bnd = bnd;
clear equ
equ.Sd = {{'lambdaS_rfwe';'lambdaS_rfwe'},{'lambdaS_rfwe';'Sdy_guess_rfwe'}, ...
  {'Sdx_guess_rfwe';'lambdaS_rfwe'},{'Sdx_guess_rfwe';'Sdy_guess_rfwe'}};
equ.coordOn = {{1;1},{1;0},{0;1},{0;0}};
equ.Stype = {'coord','coord','coord','none'};
equ.ind = [1,2,1,3,4,3,1,2,1];
appl.equ = equ;
fem.appl{1} = appl;
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;

% Boundary settings
clear bnd
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];
bnd.dim = {'Ez'};

% Boundary pairs
clear pair
pair{1}.type= 'identity';
pair{1}.name= 'Pair 1';
pair{1}.src.dl = [11,13,15];
pair{1}.src.operator = 'src2dst_ip2';
pair{1}.dst.dl = [8,9,10];
pair{1}.dst.operator = 'dst2src_ip2';
bnd.pair = pair;

fem.bnd = bnd;

% Point settings
clear pnt
pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1];

% Point pairs
clear pair
pair{1}.type= 'identity';
pair{1}.name= 'Pair 1';
pair{1}.src.dl = [5,6,7,8];
pair{1}.src.operator = 'src2dst_ip1';
pair{1}.dst.dl = [9,10,11,12];
pair{1}.dst.operator = 'dst2src_ip1';
pnt.pair = pair;

fem.pnt = pnt;

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

% Multiphysics
fem=multiphysics(fem);

% Extend mesh
fem.xmesh=meshextend(fem);

% Solve problem
fem.sol=femstatic(fem, ...
                  'complexfun','on', ...
                  'solcomp',{'Ez'}, ...
                  'outcomp',{'Ez'}, ...
                  'blocksize','auto', ...
                  'linsolver','pardiso');

% Save current fem structure for restart purposes
fem0=fem;

% Plot solution
postplot(fem, ...
         'tridata',{'Ez','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, z component [V/m]', ...
         'axis',[-4.5340287361663264E-7,8.534028782910715E-7,-2.1590613029363458E-7,6.640938799901311E-7]);

% (Default values are not included)

% Application mode 1
clear appl
appl.mode.class = 'InPlaneWaves';
appl.module = 'RF';
appl.assignsuffix = '_rfwe';
clear prop
prop.inputvar='lambda';
appl.prop = prop;
clear bnd
bnd.type = {'SC','cont'};
bnd.ind = [1,1,1,2,1,2,1,1,1,1,1,1,1,2,1,2,1,2,1,2,2,2,2,1,1,1,1];
clear pair
pair.type = 'E';
pair.E0 = {{0;0;1}};
pair.pair = 'Pair 1';
bnd.pair = pair;
appl.bnd = bnd;
clear equ
equ.Sd = {{'lambdaS_rfwe';'lambdaS_rfwe'},{'lambdaS_rfwe';'Sdy_guess_rfwe'}, ...
  {'Sdx_guess_rfwe';'lambdaS_rfwe'},{'Sdx_guess_rfwe';'Sdy_guess_rfwe'}};
equ.coordOn = {{1;1},{1;0},{0;1},{0;0}};
equ.Stype = {'coord','coord','coord','none'};
equ.ind = [1,2,1,3,4,3,1,2,1];
appl.equ = equ;
appl.var = {'lambda0','400e-9'};
fem.appl{1} = appl;
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;

% Boundary settings
clear bnd
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];
bnd.dim = {'Ez'};

% Boundary pairs
clear pair
pair{1}.type= 'identity';
pair{1}.name= 'Pair 1';
pair{1}.src.dl = [11,13,15];
pair{1}.src.operator = 'src2dst_ip2';
pair{1}.dst.dl = [8,9,10];
pair{1}.dst.operator = 'dst2src_ip2';
bnd.pair = pair;

fem.bnd = bnd;

% Point settings
clear pnt
pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1];

% Point pairs
clear pair
pair{1}.type= 'identity';
pair{1}.name= 'Pair 1';
pair{1}.src.dl = [5,6,7,8];
pair{1}.src.operator = 'src2dst_ip1';
pair{1}.dst.dl = [9,10,11,12];
pair{1}.dst.operator = 'dst2src_ip1';
pnt.pair = pair;

fem.pnt = pnt;

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

% Multiphysics
fem=multiphysics(fem);

% Extend mesh
fem.xmesh=meshextend(fem);

% Solve problem
fem.sol=femstatic(fem, ...
                  'complexfun','on', ...
                  'solcomp',{'Ez'}, ...
                  'outcomp',{'Ez'}, ...
                  'blocksize','auto', ...
                  'linsolver','pardiso');

% Save current fem structure for restart purposes
fem0=fem;

% Plot solution
postplot(fem, ...
         'tridata',{'Ez','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, z component [V/m]', ...
         'axis',[-4.6643617566358444E-7,8.664361803380232E-7,-2.1590613029363458E-7,6.640938799901311E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'Hx_rfwe','cont','internal','unit','A/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Magnetic field, x component [A/m]', ...
         'axis',[-4.6643617566358444E-7,8.664361803380232E-7,-2.1590613029363458E-7,6.640938799901311E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'Hx_rfwe','cont','internal','unit','A/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Magnetic field, x component [A/m]', ...
         'axis',[-4.534028736166325E-7,8.534028782910713E-7,-2.4253382003997273E-7,6.907215697364692E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'Hx_rfwe','cont','internal','unit','A/m'}, ...
         'trimap','jet(1024)', ...
         'lindata',{'Ez','cont','internal','unit','V/m'}, ...
         'linmap','jet(1024)', ...
         'title','Surface: Magnetic field, x component [A/m]   Boundary: Electric field, z component [V/m]', ...
         'axis',[-4.6643617566358444E-7,8.664361803380232E-7,-2.1590613029363458E-7,6.640938799901311E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'Hx_rfwe','cont','internal','unit','A/m'}, ...
         'trimap','jet(1024)', ...
         'lindata',{'Ez','cont','internal','unit','V/m'}, ...
         'linmap','jet(1024)', ...
         'title','Surface: Magnetic field, x component [A/m]   Boundary: Electric field, z component [V/m]', ...
         'axis',[-4.534028736166325E-7,8.534028782910713E-7,-2.7627166763372394E-7,7.244594173302204E-7]);

% Plot solution
postplot(fem, ...
         'lindata',{'Ez','cont','internal','unit','V/m'}, ...
         'linmap','jet(1024)', ...
         'title','Boundary: Electric field, z component [V/m]', ...
         'axis',[-4.534028736166325E-7,8.534028782910713E-7,-2.7627166763372394E-7,7.244594173302204E-7]);

% Plot solution
postplot(fem, ...
         'lindata',{'Ez','cont','internal','unit','V/m'}, ...
         'linmap','jet(1024)', ...
         'title','Boundary: Electric field, z component [V/m]', ...
         'axis',[-4.534028736166325E-7,8.534028782910713E-7,-2.3598006028030677E-7,6.841678099768032E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'Hx_rfwe','cont','internal','unit','A/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Magnetic field, x component [A/m]', ...
         'axis',[-4.534028736166325E-7,8.534028782910713E-7,-2.3598006028030677E-7,6.841678099768032E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'normH_rfwe','cont','internal','unit','A/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Magnetic field, norm [A/m]', ...
         'axis',[-4.6643617566358444E-7,8.664361803380232E-7,-2.1590613029363458E-7,6.640938799901311E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'normH_rfwe','cont','internal','unit','A/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Magnetic field, norm [A/m]', ...
         'axis',[-4.534028736166325E-7,8.534028782910713E-7,-2.4164921777478174E-7,6.898369674712782E-7]);

% Integrate
I1=postint(fem,'Ez', ...
           'unit','V/m', ...
           'recover','off', ...
           'dl',[15], ...
           'edim',0);

% Geometry
g14=geomcomp({g12},'ns',{'g12'},'sf','g12','edge','none');
g15=geomcomp({g14},'ns',{'g14'},'sf','g14','edge','none');
g16=geomcomp({g13},'ns',{'g13'},'sf','g13','edge','none');

% Create assembly and imprints
clear p s
p.objs={g3};
p.name={'PT1'};
p.tags={'g3'};

s.objs={g15,g16};
s.name={'CO1','CO3'};
s.tags={'g15','g16'};

fem.draw=struct('p',p,'s',s);
fem.geom=geomgroup(fem,'imprint','on',...
  'paircand',{'CO1','CO3'});

% Create assembly and imprints
clear p s
p.objs={g3};
p.name={'PT1'};
p.tags={'g3'};

s.objs={g15,g16};
s.name={'CO1','CO3'};
s.tags={'g15','g16'};

fem.draw=struct('p',p,'s',s);
fem.geom=geomgroup(fem,'imprint','on',...
  'paircand',{'CO1','CO3'});

% Analyzed geometry
clear p s
p.objs={g3};
p.name={'PT1'};
p.tags={'g3'};

s.objs={g15,g16};
s.name={'CO1','CO3'};
s.tags={'g15','g16'};

fem.draw=struct('p',p,'s',s);
fem.geom=geomgroup(fem,'imprint','off','paircand','none');

% Initialize mesh
fem.mesh=meshinit(fem, ...
                  'hauto',2);

% Refine mesh
fem.mesh=meshrefine(fem, ...
                    'mcase',0, ...
                    'rmethod','regular');

% (Default values are not included)

% Application mode 1
clear appl
appl.mode.class = 'InPlaneWaves';
appl.module = 'RF';
appl.assignsuffix = '_rfwe';
clear prop
prop.inputvar='lambda';
appl.prop = prop;
clear bnd
bnd.type = {'SC','cont'};
bnd.ind = [1,1,1,2,1,2,1,1,1,1,1,1,1,2,1,2,1,2,1,2,2,2,2,1,1,1,1];
clear pair
pair.type = 'E';
pair.E0 = {{0;0;1}};
pair.pair = 'Pair 1';
bnd.pair = pair;
appl.bnd = bnd;
clear equ
equ.Sd = {{'lambdaS_rfwe';'lambdaS_rfwe'},{'lambdaS_rfwe';'Sdy_guess_rfwe'}, ...
  {'Sdx_guess_rfwe';'lambdaS_rfwe'},{'Sdx_guess_rfwe';'Sdy_guess_rfwe'}};
equ.coordOn = {{1;1},{1;0},{0;1},{0;0}};
equ.Stype = {'coord','coord','coord','none'};
equ.ind = [1,2,1,3,4,3,1,2,1];
appl.equ = equ;
appl.var = {'lambda0','400e-9'};
fem.appl{1} = appl;
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;

% Boundary settings
clear bnd
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];
bnd.dim = {'Ez'};

% Boundary pairs
clear pair
pair{1}.type= 'identity';
pair{1}.name= 'Pair 1';
pair{1}.src.dl = [11,13,15];
pair{1}.src.operator = 'src2dst_ip2';
pair{1}.dst.dl = [8,9,10];
pair{1}.dst.operator = 'dst2src_ip2';
bnd.pair = pair;

fem.bnd = bnd;

% Point settings
clear pnt
pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1];

% Point pairs
clear pair
pair{1}.type= 'identity';
pair{1}.name= 'Pair 1';
pair{1}.src.dl = [5,6,7,8];
pair{1}.src.operator = 'src2dst_ip1';
pair{1}.dst.dl = [9,10,11,12];
pair{1}.dst.operator = 'dst2src_ip1';
pnt.pair = pair;

fem.pnt = pnt;

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

% Multiphysics
fem=multiphysics(fem);

% Extend mesh
fem.xmesh=meshextend(fem);

% Solve problem
fem.sol=femstatic(fem, ...
                  'complexfun','on', ...
                  'solcomp',{'Ez'}, ...
                  'outcomp',{'Ez'}, ...
                  'blocksize','auto', ...
                  'linsolver','pardiso');

% Save current fem structure for restart purposes
fem0=fem;

% Plot solution
postplot(fem, ...
         'tridata',{'normH_rfwe','cont','internal','unit','A/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Magnetic field, norm [A/m]', ...
         'axis',[-4.5340287361663264E-7,8.534028782910715E-7,-2.1590613029363458E-7,6.640938799901311E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'normE_rfwe','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, norm [V/m]', ...
         'axis',[-4.534028736166325E-7,8.534028782910713E-7,-2.4164921777478174E-7,6.898369674712782E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'Ez','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, z component [V/m]', ...
         'axis',[-4.534028736166325E-7,8.534028782910713E-7,-2.3684323839499443E-7,6.850309880914909E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'Ez','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, z component [V/m]', ...
         'axis',[-4.534028736166325E-7,8.534028782910713E-7,-2.3340981483266122E-7,6.815975645291577E-7]);

% (Default values are not included)

% Application mode 1
clear appl
appl.mode.class = 'InPlaneWaves';
appl.module = 'RF';
appl.assignsuffix = '_rfwe';
clear prop
prop.inputvar='lambda';
appl.prop = prop;
clear bnd
bnd.type = {'SC','cont'};
bnd.ind = [1,1,1,2,1,2,1,1,1,1,1,1,1,2,1,2,1,2,1,2,2,2,2,1,1,1,1];
clear pair
pair.type = 'E';
pair.E0 = {{0;0;1}};
pair.pair = 'Pair 1';
bnd.pair = pair;
appl.bnd = bnd;
clear equ
equ.Sd = {{'lambdaS_rfwe';'lambdaS_rfwe'},{'lambdaS_rfwe';'Sdy_guess_rfwe'}, ...
  {'Sdx_guess_rfwe';'lambdaS_rfwe'},{'Sdx_guess_rfwe';'Sdy_guess_rfwe'}};
equ.coordOn = {{1;1},{1;0},{0;1},{0;0}};
equ.Stype = {'coord','coord','coord','none'};
equ.ind = [1,2,1,3,4,3,1,2,1];
appl.equ = equ;
appl.var = {'lambda0','400e-9'};
fem.appl{1} = appl;
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;

% Boundary settings
clear bnd
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];
bnd.dim = {'Ez'};

% Boundary pairs
clear pair
pair{1}.type= 'identity';
pair{1}.name= 'Pair 1';
pair{1}.src.dl = [11,13,15];
pair{1}.src.operator = 'src2dst_ip2';
pair{1}.dst.dl = [8,9,10];
pair{1}.dst.operator = 'dst2src_ip2';
bnd.pair = pair;

fem.bnd = bnd;

% Point settings
clear pnt
pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1];

% Point pairs
clear pair
pair{1}.type= 'identity';
pair{1}.name= 'Pair 1';
pair{1}.src.dl = [5,6,7,8];
pair{1}.src.operator = 'src2dst_ip1';
pair{1}.dst.dl = [9,10,11,12];
pair{1}.dst.operator = 'dst2src_ip1';
pnt.pair = pair;

fem.pnt = pnt;

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

% Multiphysics
fem=multiphysics(fem);

% Extend mesh
fem.xmesh=meshextend(fem);

% Solve problem
fem.sol=femstatic(fem, ...
                  'complexfun','on', ...
                  'solcomp',{'Ez'}, ...
                  'outcomp',{'Ez'}, ...
                  'blocksize','auto', ...
                  'linsolver','gmres');

% Save current fem structure for restart purposes
fem0=fem;

% Plot solution
postplot(fem, ...
         'tridata',{'Ez','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, z component [V/m]', ...
         'axis',[1.1832464284298788E-7,2.816753618314509E-7,1.6909387420551286E-7,2.7909387549098357E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'normE_rfwe','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, norm [V/m]', ...
         'axis',[-1.4647386903682136E-7,5.464738737112583E-7,-4.6579699922066994E-9,4.5284571968870287E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'normE_rfwe','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, norm [V/m]', ...
         'axis',[-1.396979784623616E-7,5.396979831367986E-7,-1.5542974646390053E-8,4.637307243428862E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'Ez','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, z component [V/m]', ...
         'axis',[-1.3969797846236168E-7,5.396979831367987E-7,-1.554297464639004E-8,4.637307243428862E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'Ez','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, z component [V/m]', ...
         'axis',[-1.4224730664847483E-7,5.422473113229118E-7,-1.554297464639004E-8,4.637307243428862E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'scEz','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Scattered electric field, z component [V/m]', ...
         'axis',[-1.4224730664847483E-7,5.422473113229118E-7,-1.554297464639004E-8,4.637307243428862E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'scEz','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Scattered electric field, z component [V/m]', ...
         'axis',[-1.403353105088899E-7,5.403353151833269E-7,-1.554297464639004E-8,4.637307243428862E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'scEz','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Scattered electric field, z component [V/m]', ...
         'axis',[-1.403353105088899E-7,5.403353151833269E-7,-1.554297464639004E-8,4.637307243428862E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'scEz','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Scattered electric field, z component [V/m]', ...
         'axis',[-1.403353105088899E-7,5.403353151833269E-7,-1.554297464639004E-8,4.637307243428862E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'Ez','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, z component [V/m]', ...
         'axis',[-1.403353105088899E-7,5.403353151833269E-7,-1.554297464639004E-8,4.637307243428862E-7]);

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

% (Default values are not included)

% Application mode 1
clear appl
appl.mode.class = 'InPlaneWaves';
appl.module = 'RF';
appl.sshape = 2;
appl.assignsuffix = '_rfwe';
clear prop
prop.inputvar='lambda';
appl.prop = prop;
clear bnd
bnd.type = {'SC','cont'};
bnd.ind = [1,1,1,2,1,2,1,1,1,1,1,1,1,2,1,2,1,2,1,2,2,2,2,1,1,1,1];
clear pair
pair.type = 'E';
pair.E0 = {{0;0;1}};
pair.pair = 'Pair 1';
bnd.pair = pair;
appl.bnd = bnd;
clear equ
equ.coordOn = {{1;1},{1;0},{0;1},{0;0}};
equ.Stype = {'coord','coord','coord','none'};
equ.ind = [1,2,1,3,4,3,1,2,1];
appl.equ = equ;
appl.var = {'lambda0','400e-9'};
fem.appl{1} = appl;
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;

% Boundary settings
clear bnd
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];
bnd.dim = {'Ez'};

% Boundary pairs
clear pair
pair{1}.type= 'identity';
pair{1}.name= 'Pair 1';
pair{1}.src.dl = [11,13,15];
pair{1}.src.operator = 'src2dst_ip2';
pair{1}.dst.dl = [8,9,10];
pair{1}.dst.operator = 'dst2src_ip2';
bnd.pair = pair;

fem.bnd = bnd;

% Point settings
clear pnt
pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1];

% Point pairs
clear pair
pair{1}.type= 'identity';
pair{1}.name= 'Pair 1';
pair{1}.src.dl = [5,6,7,8];
pair{1}.src.operator = 'src2dst_ip1';
pair{1}.dst.dl = [9,10,11,12];
pair{1}.dst.operator = 'dst2src_ip1';
pnt.pair = pair;

fem.pnt = pnt;

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

% Multiphysics
fem=multiphysics(fem);