% 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');
g2=rect2(1.0E-7,1.0E-7,'base','corner','pos',[-1.0E-7,4.0E-7]);
g3=rect2(4.0E-7,1.0E-7,'base','corner','pos',[0,4.0E-7]);
g4=rect2(1.0E-7,1.0E-7,'base','corner','pos',[4.0E-7,4.0E-7]);
g5=rect2(1.0E-7,4.0E-7,'base','corner','pos',[4.0E-7,0]);
g6=rect2(1.0E-7,1.0E-7,'base','corner','pos',[4.0E-7,-1.0E-7]);
g7=rect2(4.0E-7,1.0E-7,'base','corner','pos',[0,-1.0E-7]);
g8=rect2(1.0E-7,1.0E-7,'base','corner','pos',[-1.0E-7,-1.0E-7]);
g3=move(g3,[0,5.000000000000005E-8]);
g7=move(g7,[0,5.0E-8]);
% Import CAD data
garr = geomimport('D:\Yu_che_Ho\3d\111.mphbin', ...
                  'repairtol','1.0E-5');
[g9,g10,g11]=deal(garr{:});
g12=rect2('400e-9','400e-9','base','corner','pos',{'0','0'},'rot','0');

% Analyzed geometry
clear s
s.objs={g12};
s.name={'R1'};
s.tags={'g12'};

fem.draw=struct('s',s);
fem.geom=geomcsg(fem);
g13=rect2(1.0E-7,1.0E-7,'base','corner','pos',[-1.0E-7,4.0E-7]);
g14=rect2(4.0E-7,1.0E-7,'base','corner','pos',[0,4.0E-7]);
g15=rect2(1.0E-7,1.0E-7,'base','corner','pos',[4.0E-7,4.0E-7]);
g16=rect2(1.0E-7,4.0E-7,'base','corner','pos',[4.0E-7,0]);
g17=rect2(1.0E-7,1.0E-7,'base','corner','pos',[4.0E-7,-1.0E-7]);
g18=rect2(4.0E-7,1.0E-7,'base','corner','pos',[0,-1.0E-7]);
g19=rect2(1.0E-7,1.0E-7,'base','corner','pos',[-1.0E-7,-1.0E-7]);
g20=rect2(1.0E-7,4.0E-7,'base','corner','pos',[-1.0E-7,0]);
g13=move(g13,[5.0E-8,5.000000000000005E-8]);
g13=move(g13,[-0.5e-7,-0.5e-7]);

% Analyzed geometry
clear s
s.objs={g12,g13,g14,g15,g16,g17,g18,g19,g20};
s.name={'R1','R2','R3','R4','R5','R6','R7','R8','R9'};
s.tags={'g12','g13','g14','g15','g16','g17','g18','g19','g20'};

fem.draw=struct('s',s);
fem.geom=geomcsg(fem);
parr={point2(0,0)};
g21=geomcoerce('point',parr);
parr={point2(2.0E-7,2.0E-7)};
g22=geomcoerce('point',parr);

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

s.objs={g12,g13,g14,g15,g16,g17,g18,g19,g20};
s.name={'R1','R2','R3','R4','R5','R6','R7','R8','R9'};
s.tags={'g12','g13','g14','g15','g16','g17','g18','g19','g20'};

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

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

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

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

% 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.name = 'rfweh';
appl.module = 'RF';
appl.border = 'on';
appl.assignsuffix = '_rfweh';
clear prop
prop.field='TES';
prop.inputvar='lambda';
appl.prop = prop;
clear bnd
bnd.type = {'cont','SC'};
bnd.ind = [2,2,2,1,2,1,2,1,2,1,1,1,1,2,1,2,1,1,1,1,2,2,2,2];
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;

% 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',{'scEz'}, ...
                  'outcomp',{'scEz'}, ...
                  'blocksize','auto', ...
                  'linsolver','pardiso');

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

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

% Plot solution
postplot(fem, ...
         'tridata',{'Ez','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, z component [V/m]', ...
         'axis',[-2.900521591295089E-7,6.900521780148026E-7,-1.5197633321853858E-7,5.519763521038322E-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',[-3.026896076027163E-7,7.0268962648801E-7,-1.5197633321853858E-7,5.519763521038322E-7]);

% Plot solution
postplot(fem, ...
         'tridata',{'normE_rfweh','cont','internal','unit','V/m'}, ...
         'trimap','jet(1024)', ...
         'title','Surface: Electric field, norm [V/m]', ...
         'axis',[-2.900521591295089E-7,6.900521780148026E-7,-1.5197633321853858E-7,5.519763521038322E-7]);

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

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

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

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

% Integrate
I4=postint(fem,'normE_rfweh', ...
           'unit','V/m', ...
           'recover','off', ...
           'dl',[9], ...
           'edim',0);

% Integrate
I5=postint(fem,'normE_rfweh', ...
           'unit','V/m', ...
           'recover','off', ...
           'dl',[9], ...
           'edim',0);

% Integrate
I6=postint(fem,'normE_rfweh', ...
           'unit','V/m', ...
           'recover','off', ...
           'dl',[9], ...
           'edim',0);

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

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

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