% COMSOL Multiphysics Model M-file
% Generated by COMSOL 3.5a (COMSOL 3.5.0.603, $Date: 2008/12/03 17:02:19 $)
% Some geometry objects are stored in a separate file.
% The name of this file is given by the variable 'flbinaryfile'.

clear
flclear fem

% COMSOL version
clear vrsn
vrsn.name = 'COMSOL 3.5';
vrsn.ext = 'a';
vrsn.major = 0;
vrsn.build = 603;
vrsn.rcs = '$Name:  $';
vrsn.date = '$Date: 2008/12/03 17:02:19 $';
fem.version = vrsn;

flbinaryfile='circ3_1inc.mphm';

% Geometry
clear draw
g10=flbinary('g10','draw',flbinaryfile);
g11=flbinary('g11','draw',flbinaryfile);
g12=flbinary('g12','draw',flbinaryfile);
g13=flbinary('g13','draw',flbinaryfile);
g14=flbinary('g14','draw',flbinaryfile);
g15=flbinary('g15','draw',flbinaryfile);
g4=flbinary('g4','draw',flbinaryfile);
g5=flbinary('g5','draw',flbinaryfile);
g6=flbinary('g6','draw',flbinaryfile);
g7=flbinary('g7','draw',flbinaryfile);
g8=flbinary('g8','draw',flbinaryfile);
g9=flbinary('g9','draw',flbinaryfile);
draw.p.objs = {g10,g11,g12,g13,g14,g15,g4,g5,g6,g7,g8,g9};
draw.p.name = {'PT1','PT10','PT11','PT12','PT2','PT3','PT4','PT5','PT6','PT7', ...
  'PT8','PT9'};
draw.p.tags = {'g10','g11','g12','g13','g14','g15','g4','g5','g6','g7','g8', ...
  'g9'};
g1=flbinary('g1','draw',flbinaryfile);
g16=flbinary('g16','draw',flbinaryfile);
g3=flbinary('g3','draw',flbinaryfile);
draw.s.objs = {g1,g16,g3};
draw.s.name = {'C1','C2','CO1'};
draw.s.tags = {'g1','g16','g3'};
fem.draw = draw;
fem.geom = geomcsg(fem);

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

% (Default values are not included)

% Application mode 1
clear appl
appl.mode.class = 'ConductiveMediaDC';
appl.assignsuffix = '_dc';
clear prop
prop.elemdefault='Lag1';
appl.prop = prop;
clear pnt
pnt.name = {'','sink','source'};
pnt.Qj0 = {0,-1.0,1.0};
pnt.ind = [3,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];
appl.pnt = pnt;
clear bnd
bnd.type = {'nJ0','cont'};
bnd.ind = [1,1,1,1,2,2,2,2,1,1,2,2,2,2,1,1,1,2,2,2,2,1,1,1,1,1,1,1];
appl.bnd = bnd;
clear equ
equ.sigma = {1.0,10.0,1.0,10.0};
equ.ind = [1,2,3,4];
appl.equ = equ;
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);

% KJE added 2009Jul23
fem.xmesh = meshextend(fem);
for isig = 1:5
    sig2 = isig * 10;
    equ.sigma{2} = sig2;
    equ.sigma
    fem.sol = femstatic(fem);
    save(isig,1:539) = fem.sol.u;
end