function out = model
%
% errorModel.m
%
% Model exported on Jun 13 2013, 07:46 by COMSOL 4.3.0.151.

import com.comsol.model.*
import com.comsol.model.util.*

model = ModelUtil.create('Model');

model.modelPath('/home/bauerle/Diplomarbeit/Implementierung/data');

model.name('PIPS-modified_MeasurementModel_Template_simple1.mph');

model.comments('Measurement model for generating synthetic measurement data.');

model.modelNode.create('mod1');

model.file.create('res5');
model.file.create('res6');

model.func.create('pw1', 'Piecewise');
model.func.create('im1', 'Image');
model.func('pw1').model('mod1');
model.func('pw1').set('pieces', {'11000' '15000' '1'; '15000' '19000' '-1'});
model.func('im1').model('mod1');
model.func('im1').set('sourcetype', 'model');
model.func('im1').set('modelres', 'res5');
model.func('im1').set('importedname', 'berge02.jpg');
model.func('im1').set('importedwidth', '2400 px');
model.func('im1').set('importedheight', '1600 px');
model.func('im1').set('xmax', '20000');
model.func('im1').set('ymax', '30000');

model.file('res5').resource('/tmp/cs096807/tmp7763294645619765531_copy');

model.geom.create('geom1', 3);
model.geom('geom1').feature.create('blk3', 'Block');
model.geom('geom1').feature.create('ps1', 'ParametricSurface');
model.geom('geom1').feature.create('blk4', 'Block');
model.geom('geom1').feature.create('blk5', 'Block');
model.geom('geom1').feature.create('blk6', 'Block');
model.geom('geom1').feature.create('blk7', 'Block');
model.geom('geom1').feature.create('blk8', 'Block');
model.geom('geom1').feature.create('pyr1', 'Pyramid');
model.geom('geom1').feature.create('tet1', 'Tetrahedron');
model.geom('geom1').feature.create('measurement_pt1', 'Point');
model.geom('geom1').feature.create('measurement_pt2', 'Point');
model.geom('geom1').feature.create('measurement_pt3', 'Point');
model.geom('geom1').feature.create('measurement_pt4', 'Point');
model.geom('geom1').feature.create('measurement_pt5', 'Point');
model.geom('geom1').feature.create('measurement_pt6', 'Point');
model.geom('geom1').feature.create('measurement_pt7', 'Point');
model.geom('geom1').feature.create('measurement_pt8', 'Point');
model.geom('geom1').feature.create('measurement_pt9', 'Point');
model.geom('geom1').feature.create('measurement_pt10', 'Point');
model.geom('geom1').feature.create('measurement_pt11', 'Point');
model.geom('geom1').feature.create('measurement_pt12', 'Point');
model.geom('geom1').feature.create('measurement_pt13', 'Point');
model.geom('geom1').feature.create('measurement_pt14', 'Point');
model.geom('geom1').feature.create('measurement_pt15', 'Point');
model.geom('geom1').feature.create('measurement_pt16', 'Point');
model.geom('geom1').feature.create('measurement_pt17', 'Point');
model.geom('geom1').feature.create('measurement_pt18', 'Point');
model.geom('geom1').feature.create('measurement_pt19', 'Point');
model.geom('geom1').feature.create('measurement_pt20', 'Point');
model.geom('geom1').feature.create('measurement_pt21', 'Point');
model.geom('geom1').feature.create('measurement_pt22', 'Point');
model.geom('geom1').feature.create('measurement_pt23', 'Point');
model.geom('geom1').feature.create('measurement_pt24', 'Point');
model.geom('geom1').feature.create('measurement_pt25', 'Point');
model.geom('geom1').feature.create('measurement_pt26', 'Point');
model.geom('geom1').feature.create('measurement_pt27', 'Point');
model.geom('geom1').feature.create('measurement_pt28', 'Point');
model.geom('geom1').feature.create('measurement_pt29', 'Point');
model.geom('geom1').feature.create('measurement_pt30', 'Point');
model.geom('geom1').feature.create('measurement_pt31', 'Point');
model.geom('geom1').feature.create('measurement_pt32', 'Point');
model.geom('geom1').feature.create('measurement_pt33', 'Point');
model.geom('geom1').feature.create('measurement_pt34', 'Point');
model.geom('geom1').feature.create('measurement_pt35', 'Point');
model.geom('geom1').feature.create('measurement_pt36', 'Point');
model.geom('geom1').feature.create('measurement_pt37', 'Point');
model.geom('geom1').feature.create('measurement_pt38', 'Point');
model.geom('geom1').feature.create('measurement_pt39', 'Point');
model.geom('geom1').feature.create('measurement_pt40', 'Point');
model.geom('geom1').feature.create('measurement_pt41', 'Point');
model.geom('geom1').feature.create('measurement_pt42', 'Point');
model.geom('geom1').feature.create('measurement_pt43', 'Point');
model.geom('geom1').feature.create('measurement_pt44', 'Point');
model.geom('geom1').feature.create('measurement_pt45', 'Point');
model.geom('geom1').feature.create('measurement_pt46', 'Point');
model.geom('geom1').feature.create('measurement_pt47', 'Point');
model.geom('geom1').feature.create('measurement_pt48', 'Point');
model.geom('geom1').feature.create('measurement_pt49', 'Point');
model.geom('geom1').feature.create('measurement_pt50', 'Point');
model.geom('geom1').feature.create('measurement_pt51', 'Point');
model.geom('geom1').feature.create('measurement_pt52', 'Point');
model.geom('geom1').feature.create('measurement_pt53', 'Point');
model.geom('geom1').feature.create('measurement_pt54', 'Point');
model.geom('geom1').feature.create('measurement_pt55', 'Point');
model.geom('geom1').feature.create('measurement_pt56', 'Point');
model.geom('geom1').feature.create('measurement_pt57', 'Point');
model.geom('geom1').feature.create('measurement_pt58', 'Point');
model.geom('geom1').feature.create('measurement_pt59', 'Point');
model.geom('geom1').feature.create('measurement_pt60', 'Point');
model.geom('geom1').feature.create('measurement_pt61', 'Point');
model.geom('geom1').feature.create('measurement_pt62', 'Point');
model.geom('geom1').feature.create('measurement_pt63', 'Point');
model.geom('geom1').feature.create('measurement_pt64', 'Point');
model.geom('geom1').feature.create('measurement_pt65', 'Point');
model.geom('geom1').feature('blk3').set('pos', {'10000' '15000' '-3500'});
model.geom('geom1').feature('blk3').set('base', 'center');
model.geom('geom1').feature('blk3').set('size', {'4000' '8000' '1000'});
model.geom('geom1').feature('ps1').active(false);
model.geom('geom1').feature('ps1').set('pos', {'0' '0' '-500'});
model.geom('geom1').feature('ps1').set('parmax1', '20000');
model.geom('geom1').feature('ps1').set('parmax2', '30000');
model.geom('geom1').feature('ps1').set('coord', {'s1' 's2' '1800 * mod1.im1(s1,s2)'});
model.geom('geom1').feature('blk4').active(false);
model.geom('geom1').feature('blk4').set('pos', {'0' '0' '-6000'});
model.geom('geom1').feature('blk4').set('size', {'20000' '30000' '8000'});
model.geom('geom1').feature('blk4').feature.comments('');
model.geom('geom1').feature('blk5').set('pos', {'0' '0' '-6000'});
model.geom('geom1').feature('blk5').set('size', {'20000' '30000' '6000'});
model.geom('geom1').feature('blk5').feature.comments('');
model.geom('geom1').feature('blk6').set('pos', {'0' '0' '0'});
model.geom('geom1').feature('blk6').set('size', {'20000' '30000' '2000'});
model.geom('geom1').feature('blk6').feature.comments('');
model.geom('geom1').feature('blk7').active(false);
model.geom('geom1').feature('blk7').set('pos', {'3000' '15000' '-2800'});
model.geom('geom1').feature('blk7').set('size', {'8000' '10000' '500'});
model.geom('geom1').feature('blk8').set('pos', {'4000' '1000' '-2600'});
model.geom('geom1').feature('blk8').set('size', {'5000' '9000' '2000'});
model.geom('geom1').feature('pyr1').active(false);
model.geom('geom1').feature('pyr1').set('pos', {'13000' '22000' '-2300'});
model.geom('geom1').feature('pyr1').set('a', '4000');
model.geom('geom1').feature('pyr1').set('b', '7000');
model.geom('geom1').feature('pyr1').set('h', '1500');
model.geom('geom1').feature('pyr1').set('rat', '0');
model.geom('geom1').feature('tet1').active(false);
model.geom('geom1').feature('tet1').set('p', {'7000' '14000' '7000' '7000'; '17000' '21000' '21000' '21000'; '-2000' '-1400' '-1400' '-600'});
model.geom('geom1').feature('measurement_pt1').name('Measurement Point 1');
model.geom('geom1').feature('measurement_pt1').set('p', {'5000'; '4000'; '600'});
model.geom('geom1').feature('measurement_pt2').name('Measurement Point 2');
model.geom('geom1').feature('measurement_pt2').set('p', {'5000'; '5000'; '600'});
model.geom('geom1').feature('measurement_pt3').name('Measurement Point 3');
model.geom('geom1').feature('measurement_pt3').set('p', {'5000'; '7000'; '600'});
model.geom('geom1').feature('measurement_pt4').name('Measurement Point 4');
model.geom('geom1').feature('measurement_pt4').set('p', {'5000'; '9000'; '600'});
model.geom('geom1').feature('measurement_pt5').name('Measurement Point 5');
model.geom('geom1').feature('measurement_pt5').set('p', {'5000'; '11000'; '600'});
model.geom('geom1').feature('measurement_pt6').name('Measurement Point 6');
model.geom('geom1').feature('measurement_pt6').set('p', {'5000'; '13000'; '600'});
model.geom('geom1').feature('measurement_pt7').name('Measurement Point 7');
model.geom('geom1').feature('measurement_pt7').set('p', {'5000'; '15000'; '600'});
model.geom('geom1').feature('measurement_pt8').name('Measurement Point 8');
model.geom('geom1').feature('measurement_pt8').set('p', {'5000'; '17000'; '600'});
model.geom('geom1').feature('measurement_pt9').name('Measurement Point 9');
model.geom('geom1').feature('measurement_pt9').set('p', {'5000'; '19000'; '600'});
model.geom('geom1').feature('measurement_pt10').name('Measurement Point 10');
model.geom('geom1').feature('measurement_pt10').set('p', {'5000'; '21000'; '600'});
model.geom('geom1').feature('measurement_pt11').name('Measurement Point 11');
model.geom('geom1').feature('measurement_pt11').set('p', {'5000'; '23000'; '600'});
model.geom('geom1').feature('measurement_pt12').name('Measurement Point 12');
model.geom('geom1').feature('measurement_pt12').set('p', {'5000'; '25000'; '600'});
model.geom('geom1').feature('measurement_pt13').name('Measurement Point 13');
model.geom('geom1').feature('measurement_pt13').set('p', {'5833'; '26000'; '600'});
model.geom('geom1').feature('measurement_pt14').name('Measurement Point 14');
model.geom('geom1').feature('measurement_pt14').set('p', {'6667'; '26000'; '600'});
model.geom('geom1').feature('measurement_pt15').name('Measurement Point 15');
model.geom('geom1').feature('measurement_pt15').set('p', {'7500'; '25000'; '600'});
model.geom('geom1').feature('measurement_pt16').name('Measurement Point 16');
model.geom('geom1').feature('measurement_pt16').set('p', {'7500'; '23000'; '600'});
model.geom('geom1').feature('measurement_pt17').name('Measurement Point 17');
model.geom('geom1').feature('measurement_pt17').set('p', {'7500'; '21000'; '600'});
model.geom('geom1').feature('measurement_pt18').name('Measurement Point 18');
model.geom('geom1').feature('measurement_pt18').set('p', {'7500'; '19000'; '600'});
model.geom('geom1').feature('measurement_pt19').name('Measurement Point 19');
model.geom('geom1').feature('measurement_pt19').set('p', {'7500'; '17000'; '600'});
model.geom('geom1').feature('measurement_pt20').name('Measurement Point 20');
model.geom('geom1').feature('measurement_pt20').set('p', {'7500'; '15000'; '600'});
model.geom('geom1').feature('measurement_pt21').name('Measurement Point 21');
model.geom('geom1').feature('measurement_pt21').set('p', {'7500'; '13000'; '600'});
model.geom('geom1').feature('measurement_pt22').name('Measurement Point 22');
model.geom('geom1').feature('measurement_pt22').set('p', {'7500'; '11000'; '600'});
model.geom('geom1').feature('measurement_pt23').name('Measurement Point 23');
model.geom('geom1').feature('measurement_pt23').set('p', {'7500'; '9000'; '600'});
model.geom('geom1').feature('measurement_pt24').name('Measurement Point 24');
model.geom('geom1').feature('measurement_pt24').set('p', {'7500'; '7000'; '600'});
model.geom('geom1').feature('measurement_pt25').name('Measurement Point 25');
model.geom('geom1').feature('measurement_pt25').set('p', {'7500'; '5000'; '600'});
model.geom('geom1').feature('measurement_pt26').name('Measurement Point 26');
model.geom('geom1').feature('measurement_pt26').set('p', {'8333'; '4000'; '600'});
model.geom('geom1').feature('measurement_pt27').name('Measurement Point 27');
model.geom('geom1').feature('measurement_pt27').set('p', {'9167'; '4000'; '600'});
model.geom('geom1').feature('measurement_pt28').name('Measurement Point 28');
model.geom('geom1').feature('measurement_pt28').set('p', {'10000'; '5000'; '600'});
model.geom('geom1').feature('measurement_pt29').name('Measurement Point 29');
model.geom('geom1').feature('measurement_pt29').set('p', {'10000'; '7000'; '600'});
model.geom('geom1').feature('measurement_pt30').name('Measurement Point 30');
model.geom('geom1').feature('measurement_pt30').set('p', {'10000'; '9000'; '600'});
model.geom('geom1').feature('measurement_pt31').name('Measurement Point 31');
model.geom('geom1').feature('measurement_pt31').set('p', {'10000'; '11000'; '600'});
model.geom('geom1').feature('measurement_pt32').name('Measurement Point 32');
model.geom('geom1').feature('measurement_pt32').set('p', {'10000'; '13000'; '600'});
model.geom('geom1').feature('measurement_pt33').name('Measurement Point 33');
model.geom('geom1').feature('measurement_pt33').set('p', {'10000'; '15000'; '600'});
model.geom('geom1').feature('measurement_pt34').name('Measurement Point 34');
model.geom('geom1').feature('measurement_pt34').set('p', {'10000'; '17000'; '600'});
model.geom('geom1').feature('measurement_pt35').name('Measurement Point 35');
model.geom('geom1').feature('measurement_pt35').set('p', {'10000'; '19000'; '600'});
model.geom('geom1').feature('measurement_pt36').name('Measurement Point 36');
model.geom('geom1').feature('measurement_pt36').set('p', {'10000'; '21000'; '600'});
model.geom('geom1').feature('measurement_pt37').name('Measurement Point 37');
model.geom('geom1').feature('measurement_pt37').set('p', {'10000'; '23000'; '600'});
model.geom('geom1').feature('measurement_pt38').name('Measurement Point 38');
model.geom('geom1').feature('measurement_pt38').set('p', {'10000'; '25000'; '600'});
model.geom('geom1').feature('measurement_pt39').name('Measurement Point 39');
model.geom('geom1').feature('measurement_pt39').set('p', {'10833'; '26000'; '600'});
model.geom('geom1').feature('measurement_pt40').name('Measurement Point 40');
model.geom('geom1').feature('measurement_pt40').set('p', {'11667'; '26000'; '600'});
model.geom('geom1').feature('measurement_pt41').name('Measurement Point 41');
model.geom('geom1').feature('measurement_pt41').set('p', {'12500'; '25000'; '600'});
model.geom('geom1').feature('measurement_pt42').name('Measurement Point 42');
model.geom('geom1').feature('measurement_pt42').set('p', {'12500'; '23000'; '600'});
model.geom('geom1').feature('measurement_pt43').name('Measurement Point 43');
model.geom('geom1').feature('measurement_pt43').set('p', {'12500'; '21000'; '600'});
model.geom('geom1').feature('measurement_pt44').name('Measurement Point 44');
model.geom('geom1').feature('measurement_pt44').set('p', {'12500'; '19000'; '600'});
model.geom('geom1').feature('measurement_pt45').name('Measurement Point 45');
model.geom('geom1').feature('measurement_pt45').set('p', {'12500'; '17000'; '600'});
model.geom('geom1').feature('measurement_pt46').name('Measurement Point 46');
model.geom('geom1').feature('measurement_pt46').set('p', {'12500'; '15000'; '600'});
model.geom('geom1').feature('measurement_pt47').name('Measurement Point 47');
model.geom('geom1').feature('measurement_pt47').set('p', {'12500'; '13000'; '600'});
model.geom('geom1').feature('measurement_pt48').name('Measurement Point 48');
model.geom('geom1').feature('measurement_pt48').set('p', {'12500'; '11000'; '600'});
model.geom('geom1').feature('measurement_pt49').name('Measurement Point 49');
model.geom('geom1').feature('measurement_pt49').set('p', {'12500'; '9000'; '600'});
model.geom('geom1').feature('measurement_pt50').name('Measurement Point 50');
model.geom('geom1').feature('measurement_pt50').set('p', {'12500'; '7000'; '600'});
model.geom('geom1').feature('measurement_pt51').name('Measurement Point 51');
model.geom('geom1').feature('measurement_pt51').set('p', {'12500'; '5000'; '600'});
model.geom('geom1').feature('measurement_pt52').name('Measurement Point 52');
model.geom('geom1').feature('measurement_pt52').set('p', {'13333'; '4000'; '600'});
model.geom('geom1').feature('measurement_pt53').name('Measurement Point 53');
model.geom('geom1').feature('measurement_pt53').set('p', {'14167'; '4000'; '600'});
model.geom('geom1').feature('measurement_pt54').name('Measurement Point 54');
model.geom('geom1').feature('measurement_pt54').set('p', {'15000'; '5000'; '600'});
model.geom('geom1').feature('measurement_pt55').name('Measurement Point 55');
model.geom('geom1').feature('measurement_pt55').set('p', {'15000'; '7000'; '600'});
model.geom('geom1').feature('measurement_pt56').name('Measurement Point 56');
model.geom('geom1').feature('measurement_pt56').set('p', {'15000'; '9000'; '600'});
model.geom('geom1').feature('measurement_pt57').name('Measurement Point 57');
model.geom('geom1').feature('measurement_pt57').set('p', {'15000'; '11000'; '600'});
model.geom('geom1').feature('measurement_pt58').name('Measurement Point 58');
model.geom('geom1').feature('measurement_pt58').set('p', {'15000'; '13000'; '600'});
model.geom('geom1').feature('measurement_pt59').name('Measurement Point 59');
model.geom('geom1').feature('measurement_pt59').set('p', {'15000'; '15000'; '600'});
model.geom('geom1').feature('measurement_pt60').name('Measurement Point 60');
model.geom('geom1').feature('measurement_pt60').set('p', {'15000'; '17000'; '600'});
model.geom('geom1').feature('measurement_pt61').name('Measurement Point 61');
model.geom('geom1').feature('measurement_pt61').set('p', {'15000'; '19000'; '600'});
model.geom('geom1').feature('measurement_pt62').name('Measurement Point 62');
model.geom('geom1').feature('measurement_pt62').set('p', {'15000'; '21000'; '600'});
model.geom('geom1').feature('measurement_pt63').name('Measurement Point 63');
model.geom('geom1').feature('measurement_pt63').set('p', {'15000'; '23000'; '600'});
model.geom('geom1').feature('measurement_pt64').name('Measurement Point 64');
model.geom('geom1').feature('measurement_pt64').set('p', {'15000'; '25000'; '600'});
model.geom('geom1').feature('measurement_pt65').name('Measurement Point 65');
model.geom('geom1').feature('measurement_pt65').set('p', {'15000'; '26000'; '600'});
model.geom('geom1').run;

model.view('view1').hideObjects.create('hide1');
model.view.create('view6', 'geom1');
model.view.create('view7', 3);
model.view.create('view8', 3);
model.view.create('view9', 3);
model.view.create('view10', 3);
model.view.create('view11', 3);

model.material.create('mat2');
model.material('mat2').selection.set([2]);
model.material.create('mat3');
model.material('mat3').selection.set([1 4]);
model.material.create('mat1');
model.material.create('mat4');
model.material.create('mat5');
model.material('mat5').selection.set([3]);

model.physics.create('ec', 'ConductiveMedia', 'geom1');
model.physics('ec').field('electricpotential').field('u');
model.physics('ec').feature.create('ncd1', 'NormalCurrentDensity', 2);
model.physics('ec').feature('ncd1').selection.set([1 2 3 8 22]);
model.physics('ec').feature.create('ncd2', 'NormalCurrentDensity', 2);
model.physics('ec').feature('ncd2').selection.set([4 5 7 9 23]);
model.physics('ec').feature.create('cs1', 'CurrentSource', 3);
model.physics('ec').feature('cs1').selection.set([4]);
model.physics('ec').feature.create('cs2', 'CurrentSource', 3);
model.physics('ec').feature('cs2').selection.all;

model.mesh.create('mesh1', 'geom1');
model.mesh('mesh1').feature.create('ftet1', 'FreeTet');
model.mesh('mesh1').feature('ftet1').selection.geom('geom1', 3);
model.mesh('mesh1').feature('ftet1').selection.set([2]);
model.mesh('mesh1').feature('ftet1').feature.create('size1', 'Size');
model.mesh('mesh1').feature.create('ftet2', 'FreeTet');
model.mesh('mesh1').feature('ftet2').selection.geom('geom1', 3);
model.mesh('mesh1').feature('ftet2').selection.set([1 3 4]);
model.mesh('mesh1').feature('ftet2').feature.create('size1', 'Size');

model.view('view1').set('renderwireframe', true);
model.view('view1').set('scenelight', 'off');
model.view('view1').hideObjects('hide1').set({});
model.view('view6').set('renderwireframe', true);
model.view('view6').set('transparency', 'on');

model.material('mat2').name('Air');
model.material('mat2').propertyGroup('def').func.name('Functions');
model.material('mat2').propertyGroup('def').set('relpermittivity', {'10' '0' '0' '0' '10' '0' '0' '0' '10'});
model.material('mat2').propertyGroup('def').set('electricconductivity', {'0.000001' '0' '0' '0' '0.000001' '0' '0' '0' '0.000001'});
model.material('mat3').name('Sand');
model.material('mat3').propertyGroup('def').func.name('Functions');
model.material('mat3').propertyGroup('def').set('electricconductivity', {'0.007' '0' '0' '0' '0.007' '0' '0' '0' '0.007'});
model.material('mat3').propertyGroup('def').set('relpermittivity', {'20' '0' '0' '0' '20' '0' '0' '0' '20'});
model.material('mat1').name('Sandstone');
model.material('mat1').propertyGroup('def').func.name('Functions');
model.material('mat1').propertyGroup('def').set('relpermittivity', {'7' '0' '0' '0' '7' '0' '0' '0' '7'});
model.material('mat1').propertyGroup('def').set('electricconductivity', {'0.04' '0' '0' '0' '0.04' '0' '0' '0' '0.04'});
model.material('mat4').name('Limestone');
model.material('mat4').propertyGroup('def').func.name('Functions');
model.material('mat4').propertyGroup('def').set('electricconductivity', {'0.025' '0' '0' '0' '0.025' '0' '0' '0' '0.025'});
model.material('mat4').propertyGroup('def').set('relpermittivity', {'8' '0' '0' '0' '8' '0' '0' '0' '8'});
model.material('mat5').name('Water');
model.material('mat5').propertyGroup('def').func.name('Functions');
model.material('mat5').propertyGroup('def').set('electricconductivity', {'4.0' '0' '0' '0' '4.0' '0' '0' '0' '4.0'});
model.material('mat5').propertyGroup('def').set('relpermittivity', {'80' '0' '0' '0' '80' '0' '0' '0' '80'});

model.physics('ec').prop('ShapeProperty').set('order_electricpotential', '1');
model.physics('ec').prop('EquationForm').set('form', 'Stationary');
model.physics('ec').feature('cucn1').set('rho0', '0');
model.physics('ec').feature('cucn1').set('Tref', '293.15[K]');
model.physics('ec').feature('cucn1').set('minput_electricfield_src', 'root.mod1.ec.Ex');
model.physics('ec').feature('ncd1').set('nJ', '-(1/sqrt((x-10000)^2+(y-15000)^2+(z+3500)^2))*0.007*u');
model.physics('ec').feature('ncd1').set('J0', {'sqrt((x-4000)^2+'; '0'; '0'});
model.physics('ec').feature('ncd1').name('Robin Ground');
model.physics('ec').feature('ncd2').set('nJ', '-(1/sqrt((x-10000)^2+(y-15000)^2+(z+3500)^2))*0.000001*u');
model.physics('ec').feature('ncd2').set('J0', {'sqrt((x-4000)^2+'; '0'; '0'});
model.physics('ec').feature('ncd2').name('Robin Air');
model.physics('ec').feature('cs1').selection.active(false);
model.physics('ec').feature('cs1').set('Qj', '2*mod1.pw1(y)');
model.physics('ec').feature('cs1').name('Excitation');
model.physics('ec').feature('cs2').selection.active(false);
model.physics('ec').feature('cs2').set('Qj', '1');
model.physics('ec').feature('cs2').active(false);
model.physics('ec').feature('cs2').name('Scale');

model.mesh('mesh1').feature('size').set('hauto', 3);
model.mesh('mesh1').feature('ftet1').feature('size1').set('hauto', 3);
model.mesh('mesh1').feature('ftet2').feature('size1').set('hmax', '900');
model.mesh('mesh1').feature('ftet2').feature('size1').set('hmin', '45');
model.mesh('mesh1').feature('ftet2').feature('size1').set('hminactive', false);
model.mesh('mesh1').feature('ftet2').feature('size1').set('hcurve', '0.3');
model.mesh('mesh1').feature('ftet2').feature('size1').set('hcurveactive', false);
model.mesh('mesh1').feature('ftet2').feature('size1').set('hnarrow', '0.85');
model.mesh('mesh1').feature('ftet2').feature('size1').set('hnarrowactive', false);
model.mesh('mesh1').feature('ftet2').feature('size1').set('hgrad', '1.35');
model.mesh('mesh1').feature('ftet2').feature('size1').set('hgradactive', false);
model.mesh('mesh1').feature('ftet2').feature('size1').set('hauto', '2');
model.mesh('mesh1').feature('ftet2').feature('size1').set('hmax', '900');
model.mesh('mesh1').feature('ftet2').feature('size1').set('hminactive', false);
model.mesh('mesh1').feature('ftet2').feature('size1').set('hcurveactive', false);
model.mesh('mesh1').feature('ftet2').feature('size1').set('hnarrowactive', false);
model.mesh('mesh1').feature('ftet2').feature('size1').set('hgradactive', false);
model.mesh('mesh1').run;

model.frame('material1').sorder(1);

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

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.create('i1', 'Iterative');
model.sol('sol1').feature('s1').feature('i1').feature.create('mg1', 'Multigrid');
model.sol('sol1').feature('s1').feature.remove('fcDef');

model.result.dataset.create('cpl1', 'CutPlane');
model.result.dataset.create('surf1', 'Surface');
model.result.dataset('surf1').selection.set([15 16 17 18 19]);
model.result.dataset.create('surf2', 'Surface');
model.result.dataset('surf2').selection.set([6]);
model.result.dataset.create('surf3', 'Surface');
model.result.dataset('surf3').selection.set([1 2 3 6 8]);
model.result.dataset.create('measurement_pc1', 'ParCurve3D');
model.result.dataset.create('measurement_pc2', 'ParCurve3D');
model.result.dataset.create('measurement_pc3', 'ParCurve3D');
model.result.dataset.create('measurement_pc4', 'ParCurve3D');
model.result.dataset.create('measurement_pc5', 'ParCurve3D');
model.result.dataset.create('flight_surf', 'Surface');
model.result.create('pg1', 'PlotGroup3D');
model.result('pg1').feature.create('mslc1', 'Multislice');
model.result('pg1').feature.create('surf1', 'Surface');
model.result('pg1').feature.create('surf2', 'Surface');
model.result('pg1').feature.create('surf3', 'Surface');
model.result('pg1').feature.create('surf4', 'Surface');
model.result.create('pg2', 'PlotGroup3D');
model.result('pg2').feature.create('str1', 'Streamline');
model.result('pg2').feature('str1').selection.set([1 2 3 4 5 6 7 8 9 10 11 12]);
model.result('pg2').feature('str1').feature.create('col1', 'Color');
model.result('pg2').feature.create('surf1', 'Surface');
model.result.create('pg_flight_path', 'PlotGroup3D');
model.result('pg_flight_path').feature.create('measurement_line1', 'Line');
model.result('pg_flight_path').feature.create('measurement_line2', 'Line');
model.result('pg_flight_path').feature.create('measurement_line3', 'Line');
model.result('pg_flight_path').feature.create('measurement_line4', 'Line');
model.result('pg_flight_path').feature.create('measurement_line5', 'Line');
model.result('pg_flight_path').feature.create('surf1', 'Surface');
model.result('pg_flight_path').feature.create('surf2', 'Surface');

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').feature('s1').feature('i1').set('linsolver', 'cg');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').set('prefun', 'amg');
model.sol('sol1').runAll;

model.result.dataset('cpl1').set('quickplane', 'xy');
model.result.dataset('cpl1').set('quickz', '600');
model.result.dataset('surf1').name('Oil Reservoir');
model.result.dataset('surf2').name('Ground Surface');
model.result.dataset('surf3').name('Ground');
model.result.dataset('measurement_pc1').name('Flight Lane 1');
model.result.dataset('measurement_pc1').set('exprx', '5000');
model.result.dataset('measurement_pc1').set('expry', '5000 + 20000*s');
model.result.dataset('measurement_pc1').set('exprz', '600');
model.result.dataset('measurement_pc2').name('Flight Lane 2');
model.result.dataset('measurement_pc2').set('exprx', '7500');
model.result.dataset('measurement_pc2').set('expry', '5000 + 20000*s');
model.result.dataset('measurement_pc2').set('exprz', '600');
model.result.dataset('measurement_pc3').name('Flight Lane 3');
model.result.dataset('measurement_pc3').set('exprx', '10000');
model.result.dataset('measurement_pc3').set('expry', '5000 + 20000*s');
model.result.dataset('measurement_pc3').set('exprz', '600');
model.result.dataset('measurement_pc4').name('Flight Lane 4');
model.result.dataset('measurement_pc4').set('exprx', '12500');
model.result.dataset('measurement_pc4').set('expry', '5000 + 20000*s');
model.result.dataset('measurement_pc4').set('exprz', '600');
model.result.dataset('measurement_pc5').name('Flight Lane 5');
model.result.dataset('measurement_pc5').set('exprx', '15000');
model.result.dataset('measurement_pc5').set('expry', '5000 + 20000*s');
model.result.dataset('measurement_pc5').set('exprz', '600');
model.result('pg1').name('Electric Potential');
model.result('pg1').feature('mslc1').set('titletype', 'none');
model.result('pg1').feature('mslc1').set('multiplaneymethod', 'coord');
model.result('pg1').feature('mslc1').set('ycoord', '26000');
model.result('pg1').feature('mslc1').set('multiplanezmethod', 'coord');
model.result('pg1').feature('mslc1').set('zcoord', '-6000');
model.result('pg1').feature('surf1').active(false);
model.result('pg1').feature('surf1').set('data', 'cpl1');
model.result('pg1').feature('surf2').set('data', 'surf1');
model.result('pg1').feature('surf2').feature.comments('');
model.result('pg1').feature('surf3').active(false);
model.result('pg1').feature('surf3').set('data', 'surf3');
model.result('pg1').feature('surf3').set('coloring', 'uniform');
model.result('pg1').feature('surf3').set('color', 'black');
model.result('pg1').feature('surf3').set('wireframe', true);
model.result('pg1').feature('surf4').active(false);
model.result('pg1').feature('surf4').set('data', 'surf3');
model.result('pg2').name('Electric Field');
model.result('pg2').set('titletype', 'none');
model.result('pg2').feature('str1').set('data', 'dset1');
model.result('pg2').feature('str1').set('expr', {'ec.Ex' 'ec.Ey' 'ec.Ez'});
model.result('pg2').feature('str1').set('descr', 'Electric field');
model.result('pg2').feature('str1').set('titletype', 'none');
model.result('pg2').feature('str1').set('posmethod', 'magnitude');
model.result('pg2').feature('str1').set('linetype', 'tube');
model.result('pg2').feature('str1').set('radiusexpr', '3');
model.result('pg2').feature('str1').set('tuberadiusscale', '40.00000000000001');
model.result('pg2').feature('str1').set('tuberadiusscaleactive', false);
model.result('pg2').feature('str1').feature.comments('');
model.result('pg2').feature('str1').feature('col1').set('colorlegend', false);
model.result('pg2').feature('surf1').set('data', 'surf1');
model.result('pg_flight_path').name('Flight Path Measurements');
model.result('pg_flight_path').set('titletype', 'manual');
model.result('pg_flight_path').set('title', 'Measurement data on flight path');
model.result('pg_flight_path').feature('measurement_line1').set('data', 'measurement_pc1');
model.result('pg_flight_path').feature('measurement_line1').set('linetype', 'tube');
model.result('pg_flight_path').feature('measurement_line1').set('radiusexpr', '10');
model.result('pg_flight_path').feature('measurement_line1').set('tuberadiusscale', '10');
model.result('pg_flight_path').feature('measurement_line1').set('tuberadiusscaleactive', true);
model.result('pg_flight_path').feature('measurement_line1').set('colorlegend', false);
model.result('pg_flight_path').feature('measurement_line1').set('colortablesym', true);
model.result('pg_flight_path').feature('measurement_line2').set('data', 'measurement_pc2');
model.result('pg_flight_path').feature('measurement_line2').set('linetype', 'tube');
model.result('pg_flight_path').feature('measurement_line2').set('radiusexpr', '10');
model.result('pg_flight_path').feature('measurement_line2').set('tuberadiusscale', '10');
model.result('pg_flight_path').feature('measurement_line2').set('tuberadiusscaleactive', true);
model.result('pg_flight_path').feature('measurement_line2').set('colorlegend', false);
model.result('pg_flight_path').feature('measurement_line2').set('colortablesym', true);
model.result('pg_flight_path').feature('measurement_line3').set('data', 'measurement_pc3');
model.result('pg_flight_path').feature('measurement_line3').set('linetype', 'tube');
model.result('pg_flight_path').feature('measurement_line3').set('radiusexpr', '10');
model.result('pg_flight_path').feature('measurement_line3').set('tuberadiusscale', '10');
model.result('pg_flight_path').feature('measurement_line3').set('tuberadiusscaleactive', true);
model.result('pg_flight_path').feature('measurement_line3').set('colorlegend', false);
model.result('pg_flight_path').feature('measurement_line3').set('colortablesym', true);
model.result('pg_flight_path').feature('measurement_line4').set('data', 'measurement_pc4');
model.result('pg_flight_path').feature('measurement_line4').set('linetype', 'tube');
model.result('pg_flight_path').feature('measurement_line4').set('radiusexpr', '10');
model.result('pg_flight_path').feature('measurement_line4').set('tuberadiusscale', '10');
model.result('pg_flight_path').feature('measurement_line4').set('tuberadiusscaleactive', true);
model.result('pg_flight_path').feature('measurement_line4').set('colorlegend', false);
model.result('pg_flight_path').feature('measurement_line4').set('colortablesym', true);
model.result('pg_flight_path').feature('measurement_line5').set('data', 'measurement_pc5');
model.result('pg_flight_path').feature('measurement_line5').set('linetype', 'tube');
model.result('pg_flight_path').feature('measurement_line5').set('radiusexpr', '10');
model.result('pg_flight_path').feature('measurement_line5').set('tuberadiusscale', '10');
model.result('pg_flight_path').feature('measurement_line5').set('tuberadiusscaleactive', true);
model.result('pg_flight_path').feature('measurement_line5').set('colorlegend', false);
model.result('pg_flight_path').feature('measurement_line5').set('colortablesym', true);
model.result('pg_flight_path').feature('surf1').set('data', 'flight_surf');
model.result('pg_flight_path').feature('surf1').set('coloring', 'uniform');
model.result('pg_flight_path').feature('surf1').set('color', 'black');
model.result('pg_flight_path').feature('surf1').set('wireframe', true);
model.result('pg_flight_path').feature('surf2').set('data', 'flight_surf');
model.result('pg_flight_path').feature('surf2').set('coloring', 'uniform');
model.result('pg_flight_path').feature('surf2').set('color', 'gray');

out = model;