model = ModelUtil.create('Model');

model.modelPath('C:\Users\Administrator\Desktop\heterogen');

model.name('heterogeneity1.mph');

model.param.set('mug', '1e-5[Pa*s]');
model.param.set('k', '1e-14[m^2]');

model.modelNode.create('mod1');

model.geom.create('geom1', 2);
model.geom('geom1').feature.create('sq1', 'Square');
model.geom('geom1').feature('sq1').set('size', '10');
model.geom('geom1').run;

model.mesh.create('mesh1', 'geom1');
model.mesh('mesh1').feature.create('map1', 'Map');
model.mesh('mesh1').feature('map1').feature.create('dis1', 'Distribution');
model.mesh('mesh1').feature('map1').feature('dis1').selection.all;


model.mesh('mesh1').feature('map1').feature('dis1').set('numelem', '10');
model.mesh('mesh1').run;

model.view.create('view2', 3);
model.view.create('view3', 2);
model.view.create('view4', 3);
model.view.create('view5', 2);

model.physics.create('c', 'CoefficientFormPDE', 'geom1');
model.physics('c').field('dimensionless').field('p');
model.physics('c').field('dimensionless').component({'p'});
model.physics('c').feature.create('dir1', 'DirichletBoundary', 1);
model.physics('c').feature('dir1').selection.set([1]);
model.physics('c').feature.create('dir2', 'DirichletBoundary', 1);
model.physics('c').feature('dir2').selection.set([4]);

% model.mesh.create('mesh1', 'geom1');
% model.mesh('mesh1').feature.create('map1', 'Map');
% model.mesh('mesh1').feature('map1').feature.create('dis1', 'Distribution');
% model.mesh('mesh1').feature('map1').feature('dis1').selection.all;

model.result.table.create('evl2', 'Table');

model.view('view1').axis.set('xmin', '-0.5');
model.view('view1').axis.set('xmax', '10.5');
model.view('view1').axis.set('ymin', '-1.6550002098083496');
model.view('view1').axis.set('ymax', '11.655000686645508');
model.view('view3').axis.set('xmin', '0.13068175315856934');
model.view('view3').axis.set('xmax', '10.869318008422852');
model.view('view3').axis.set('ymin', '0.5499999523162842');
model.view('view3').axis.set('ymax', '10.449999809265137');
model.view('view5').axis.set('xmin', '-0.5125000476837158');
model.view('view5').axis.set('xmax', '10.226136207580566');
model.view('view5').axis.set('ymin', '0.17272724211215973');
model.view('view5').axis.set('ymax', '10.07272720336914');

model.physics('c').prop('Units').set('DependentVariableQuantity', 'elasticity');
model.physics('c').prop('Units').set('CustomSourceTermUnit', 'Pa/s');
model.physics('c').feature('cfeq1').set('c', {'k*p/mug' '0' '0' 'k*p/mug'});
model.physics('c').feature('cfeq1').set('f', '0');
model.physics('c').feature('init1').set('p', '10e6');
model.physics('c').feature('dir1').set('r', '10e6');
model.physics('c').feature('dir2').set('r', '1e6');

% model.mesh('mesh1').feature('map1').feature('dis1').set('numelem', '10');
% model.mesh('mesh1').run;

model.result.table('evl2').name('Evaluation 2D');
model.result.table('evl2').comments('Interactive 2D values');

model.study.create('std1');
model.study('std1').feature.create('stat', 'Stationary');
ModelUtil.showProgress(true);

model.sol.create('sol1');
model.sol('sol1').study('std1');
model.sol('sol1').attach('std1');
model.sol('sol1').feature.create('st1', 'StudyStep');
model.sol('sol1').feature.create('v1', 'Variables');
model.sol('sol1').feature.create('s1', 'Stationary');
model.sol('sol1').feature('s1').feature.create('fc1', 'FullyCoupled');
model.sol('sol1').feature('s1').feature.remove('fcDef');

model.result.create('pg1', 'PlotGroup2D');
model.result('pg1').feature.create('surf1', 'Surface');
model.result('pg1').feature.create('mesh1', 'Mesh');

model.sol('sol1').attach('std1');
model.sol('sol1').feature('st1').name('Compile Equations: Stationary');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').feature('v1').set('control', 'stat');
model.sol('sol1').feature('s1').set('control', 'stat');
model.sol('sol1').runAll;

model.result('pg1').feature('mesh1').set('titletype', 'none');
model.result('pg1').feature('mesh1').set('elemcolor', 'none');

mphplot(model,'pg1','rangenum',1);