function out = model
%
% toroidalgeom_userdefeq.m
%
% Model exported on Jun 21 2011, 15:26 by COMSOL 4.2.0.150.

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

model = ModelUtil.create('Model');

model.modelPath('/Users/sebastianboie/Desktop');

model.modelNode.create('mod1');

model.geom.create('geom1', 3);

model.mesh.create('mesh1', 'geom1');

model.physics.create('ht', 'BioHeat', 'geom1', {'I1' 'I2' 'I3' 'I4' 'I5' 'I6' 'I7' 'I8' 'I9' 'I10' 'I11' 'I12' 'I13' 'I14' 'I15' 'I16' 'I17' 'I18' 'I19' 'I20' 'I21' 'I22' 'I23' 'I24'});

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

model.geom('geom1').feature.create('tor1', 'Torus');
model.geom('geom1').feature('tor1').set('rmaj', '8');
model.geom('geom1').feature('tor1').set('rmin', '4');
model.geom('geom1').feature('tor1').set('axistype', 'spherical');
model.geom('geom1').runAll;

model.sol.create('sol1');
model.sol('sol1').study('std1');
model.sol('sol1').feature.create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
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('fc1').set('dtech', 'auto');
model.sol('sol1').feature('s1').feature('fc1').set('initstep', 0.01);
model.sol('sol1').feature('s1').feature('fc1').set('minstep', 1.0E-6);
model.sol('sol1').feature('s1').feature('fc1').set('maxiter', 50);
model.sol('sol1').feature('s1').feature.create('i1', 'Iterative');
model.sol('sol1').feature('s1').feature('i1').set('linsolver', 'gmres');
model.sol('sol1').feature('s1').feature('i1').set('prefuntype', 'left');
model.sol('sol1').feature('s1').feature('i1').set('rhob', 20);
model.sol('sol1').feature('s1').feature('i1').set('itrestart', 50);
model.sol('sol1').feature('s1').feature('fc1').set('linsolver', 'i1');
model.sol('sol1').feature('s1').feature('i1').feature.create('mg1', 'Multigrid');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').set('prefun', 'gmg');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').set('mcasegen', 'any');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature.create('sl1', 'SORLine');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('iter', 2);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('linerelax', 0.4);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('seconditer', 2);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('relax', 0.7);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('sl1', 'SORLine');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('iter', 2);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('linerelax', 0.4);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('seconditer', 2);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('relax', 0.7);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('cs').feature.create('d1', 'Direct');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('cs').feature('d1').set('linsolver', 'pardiso');
model.sol('sol1').feature('s1').feature('fc1').set('dtech', 'auto');
model.sol('sol1').feature('s1').feature('fc1').set('initstep', 0.01);
model.sol('sol1').feature('s1').feature('fc1').set('minstep', 1.0E-6);
model.sol('sol1').feature('s1').feature('fc1').set('maxiter', 50);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.result.create('pg1', 3);
model.result('pg1').set('data', 'dset1');
model.result('pg1').feature.create('surf1', 'Surface');
model.result('pg1').feature('surf1').set('expr', {'T'});
model.result('pg1').name('Temperature (ht)');
model.result('pg1').feature('surf1').set('colortable', 'ThermalLight');
model.result.create('pg2', 3);
model.result('pg2').set('data', 'dset1');
model.result('pg2').feature.create('iso1', 'Isosurface');
model.result('pg2').feature('iso1').set('expr', {'T'});
model.result('pg2').name('Isothermal Contours (ht)');
model.result('pg2').feature('iso1').set('colortable', 'ThermalLight');
model.result('pg2').feature('iso1').set('number', '10');
model.result('pg2').feature.create('arwv1', 'ArrowVolume');
model.result('pg2').feature('arwv1').set('expr', {'ht.gradTx' 'ht.gradTy' 'ht.gradTz'});
model.result('pg2').feature('arwv1').set('arrowlength', 'logarithmic');
model.result('pg2').feature('arwv1').set('color', 'gray');
model.result('pg1').run;

model.geom('geom1').runAll;

model.material.create('mat1');
model.material('mat1').selection.set([1]);
model.material('mat1').selection.all;
model.material('mat1').propertyGroup('def').set('thermalconductivity', {'10'});
model.material('mat1').propertyGroup('def').set('density', {'1'});
model.material('mat1').propertyGroup('def').set('heatcapacity', {'2'});

model.sol('sol1').study('std1');
model.sol('sol1').feature.remove('s1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').feature.create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
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('fc1').set('dtech', 'auto');
model.sol('sol1').feature('s1').feature('fc1').set('initstep', 0.01);
model.sol('sol1').feature('s1').feature('fc1').set('minstep', 1.0E-6);
model.sol('sol1').feature('s1').feature('fc1').set('maxiter', 50);
model.sol('sol1').feature('s1').feature.create('i1', 'Iterative');
model.sol('sol1').feature('s1').feature('i1').set('linsolver', 'gmres');
model.sol('sol1').feature('s1').feature('i1').set('prefuntype', 'left');
model.sol('sol1').feature('s1').feature('i1').set('rhob', 20);
model.sol('sol1').feature('s1').feature('i1').set('itrestart', 50);
model.sol('sol1').feature('s1').feature('fc1').set('linsolver', 'i1');
model.sol('sol1').feature('s1').feature('i1').feature.create('mg1', 'Multigrid');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').set('prefun', 'gmg');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').set('mcasegen', 'any');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature.create('sl1', 'SORLine');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('iter', 2);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('linerelax', 0.4);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('seconditer', 2);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('relax', 0.7);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('sl1', 'SORLine');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('iter', 2);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('linerelax', 0.4);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('seconditer', 2);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('relax', 0.7);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('cs').feature.create('d1', 'Direct');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('cs').feature('d1').set('linsolver', 'pardiso');
model.sol('sol1').feature('s1').feature('fc1').set('dtech', 'auto');
model.sol('sol1').feature('s1').feature('fc1').set('initstep', 0.01);
model.sol('sol1').feature('s1').feature('fc1').set('minstep', 1.0E-6);
model.sol('sol1').feature('s1').feature('fc1').set('maxiter', 50);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');
model.sol('sol1').runAll;

model.result('pg1').run;

model.physics.remove('ht');

model.material.remove('mat1');

model.view('view1').set('showlabels', 'off');
model.view('view1').set('showDirections', 'on');
model.view('view1').set('renderwireframe', 'on');
model.view('view1').set('showlabels', 'on');

model.mesh('mesh1').autoMeshSize(2);
model.mesh('mesh1').run;

model.modelNode('mod1').identifier('SDTorus');

model.physics.create('w', 'WeakFormPDE', 'geom1', {'u'});

model.study.create('std2');
model.study('std2').feature.create('time', 'Transient');

model.physics('w').prop('EquationForm').set('form', 1, 'Transient');
model.physics('w').feature.create('wfeq2', 'WeakFormPDE', 3);
model.physics('w').feature.remove('wfeq2');
model.physics('w').feature.create('cfeq1', 'CoefficientFormPDE', 3);
model.physics('w').feature.create('gfeq1', 'GeneralFormPDE', 3);
model.physics('w').feature('cfeq1').set('f', 1, 'u');
model.physics('w').feature.create('ge1', 'GlobalEquations', -1);
model.physics('w').feature('ge1').set('initialValueU', 1, 1, '1');
model.physics('w').feature('ge1').set('description', 1, 1, 'Aktivator');
model.physics('w').feature('ge1').set('name', 1, 1, 'akt');
model.physics('w').feature('ge1').set('name', 2, 1, 'inh');
model.physics('w').feature('ge1').set('equation', 2, 1, '');
model.physics('w').feature('ge1').set('initialValueU', 2, 1, '0');
model.physics('w').feature('ge1').set('initialValueUt', 2, 1, '0');
model.physics('w').feature('ge1').set('description', 2, 1, '');
model.physics('w').feature('ge1').set('description', 2, 1, 'Inhibitor');
model.physics('w').feature('ge1').set('initialValueU', 2, 1, '1');
model.physics('w').feature('ge1').set('equation', 1, 1, '-u_t+u-1/3*u*u*u-v+D*div*grad(u)');
model.physics('w').feature('ge1').set('equation', 2, 1, '-v_t+eps*(u+beta-gamma*v)');
model.physics('w').feature.remove('gfeq1');
model.physics('w').feature.remove('cfeq1');

model.variable.create('var1');
model.variable('var1').set('u', 'u');
model.variable('var1').remove('u');

model.param.set('eps', '0');
model.param.descr('eps', 'Zeitskalentrennung');
model.param.set('beta', '1.35');

model.variable('var1').set('act', 'u');
model.variable('var1').remove('act');
model.variable('var1').set('u', '1');
model.variable('var1').set('v', '1');
model.variable('var1').set('u', '0');
model.variable('var1').set('v', '0');

model.param.set('gamma', '0');
model.param.set('eps', '0.04');

model.physics('w').feature('wfeq1').featureInfo('info').lock('root.SDTorus.w.wfeq1.weak$1', {''});
model.physics('w').feature('wfeq1').set('weak', 1, '0');

model.sol.create('sol2');
model.sol('sol2').study('std2');
model.sol('sol2').feature.create('st1', 'StudyStep');
model.sol('sol2').feature('st1').set('study', 'std2');
model.sol('sol2').feature('st1').set('studystep', 'time');
model.sol('sol2').feature.create('v1', 'Variables');
model.sol('sol2').feature.create('t1', 'Time');
model.sol('sol2').feature('t1').set('tlist', 'range(0,0.1,1)');
model.sol('sol2').feature('t1').set('plot', 'off');
model.sol('sol2').feature('t1').set('plotgroup', 'pg1');
model.sol('sol2').feature('t1').set('plotfreq', 'tout');
model.sol('sol2').feature('t1').set('probesel', 'all');
model.sol('sol2').feature('t1').set('probes', {});
model.sol('sol2').feature('t1').set('probefreq', 'tsteps');
model.sol('sol2').feature('t1').set('control', 'time');
model.sol('sol2').feature('t1').feature.create('fc1', 'FullyCoupled');
model.sol('sol2').feature('t1').feature.remove('fcDef');
model.sol('sol2').attach('std2');

model.result.numerical.create('gev1', 'EvalGlobal');
model.result.numerical('gev1').set('data', 'dset2');
model.result.numerical('gev1').set('expr', 'akt');
model.result.numerical('gev1').name('Aktivator');
model.result.numerical.create('gev2', 'EvalGlobal');
model.result.numerical('gev2').set('data', 'dset2');
model.result.numerical('gev2').set('expr', 'inh');
model.result.numerical('gev2').name('Inhibitor');
model.result.create('pg3', 1);
model.result('pg3').set('data', 'dset2');
model.result('pg3').feature.create('glob1', 'Global');
model.result('pg3').feature('glob1').set('expr', {'akt' 'inh'});
model.result.create('pg4', 3);
model.result('pg4').set('data', 'dset2');
model.result('pg4').feature.create('slc1', 'Slice');
model.result('pg4').feature('slc1').set('expr', 'u');
model.result('pg3').run;

model.study.remove('std1');

model.physics('w').feature('ge1').set('equation', 1, 1, '-u_t+u-1/3*u*u*u-v+D*div*gradient(u)');

model.sol('sol2').study('std2');
model.sol('sol2').feature.remove('t1');
model.sol('sol2').feature.remove('v1');
model.sol('sol2').feature.remove('st1');
model.sol('sol2').feature.create('st1', 'StudyStep');
model.sol('sol2').feature('st1').set('study', 'std2');
model.sol('sol2').feature('st1').set('studystep', 'time');
model.sol('sol2').feature.create('v1', 'Variables');
model.sol('sol2').feature.create('t1', 'Time');
model.sol('sol2').feature('t1').set('tlist', 'range(0,0.1,1)');
model.sol('sol2').feature('t1').set('plot', 'off');
model.sol('sol2').feature('t1').set('plotgroup', 'pg3');
model.sol('sol2').feature('t1').set('plotfreq', 'tout');
model.sol('sol2').feature('t1').set('probesel', 'all');
model.sol('sol2').feature('t1').set('probes', {});
model.sol('sol2').feature('t1').set('probefreq', 'tsteps');
model.sol('sol2').feature('t1').set('control', 'time');
model.sol('sol2').feature('t1').feature.create('fc1', 'FullyCoupled');
model.sol('sol2').feature('t1').feature.remove('fcDef');
model.sol('sol2').attach('std2');

model.result('pg3').run;

model.physics('w').feature.remove('ge1');
model.physics.remove('w');
model.physics.create('c', 'CoefficientFormPDE', 'geom1', {'u' 'v'});

model.study.create('std3');
model.study('std3').feature.create('time', 'Transient');
model.study.remove('std2');

model.physics('c').feature('cfeq1').set('f', 1, 'u*u*u');
model.physics('c').feature('cfeq1').set('f', 2, '0');
model.physics('c').feature('cfeq1').set('c', 4, {'0' '0' '0' '0' '0' '0' '0' '0' '0'});
model.physics('c').feature('cfeq1').set('c', 1, {'1' '0' '0' '0' '1' '0' '0' '0' '1'});
model.physics('c').feature('cfeq1').set('c', 2, {'0' '0' '0' '0' '0' '0' '0' '0' '0'});
model.physics('c').feature('cfeq1').set('c', 1, {'1' '0' '0' '0' '1' '0' '0' '0' '1'});
model.physics('c').feature('cfeq1').set('da', 1, '-1');
model.physics('c').feature('cfeq1').set('da', 4, '-1');
model.physics('c').feature('cfeq1').set('a', 1, '1');
model.physics('c').feature('cfeq1').set('c', 2, {'0' '0' '0' '0' '0' '0' '0' '0' '0'});
model.physics('c').feature('cfeq1').set('da', 1, '-3');
model.physics('c').feature('cfeq1').set('a', 1, '3');
model.physics('c').feature('cfeq1').set('a', 3, '-3');
model.physics('c').feature('cfeq1').set('a', 2, '0');
model.physics('c').feature('cfeq1').set('a', 4, '-1');
model.physics('c').feature('cfeq1').set('c', 1, {'-0.12' '0' '0' '0' '-0.12' '0' '0' '0' '-0.12'});

model.mesh('mesh1').autoMeshSize(1);
model.mesh('mesh1').run;

model.sol.create('sol1');
model.sol('sol1').study('std3');
model.sol('sol1').feature.create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std3');
model.sol('sol1').feature('st1').set('studystep', 'time');
model.sol('sol1').feature.create('v1', 'Variables');
model.sol('sol1').feature.create('t1', 'Time');
model.sol('sol1').feature('t1').set('tlist', 'range(0,0.1,1)');
model.sol('sol1').feature('t1').set('plot', 'off');
model.sol('sol1').feature('t1').set('plotfreq', 'tout');
model.sol('sol1').feature('t1').set('probesel', 'all');
model.sol('sol1').feature('t1').set('probes', {});
model.sol('sol1').feature('t1').set('probefreq', 'tsteps');
model.sol('sol1').feature('t1').set('control', 'time');
model.sol('sol1').feature('t1').feature.create('fc1', 'FullyCoupled');
model.sol('sol1').feature('t1').feature.remove('fcDef');
model.sol('sol1').attach('std3');

model.result.create('pg1', 3);
model.result('pg1').set('data', 'dset1');
model.result('pg1').feature.create('slc1', 'Slice');
model.result('pg1').feature('slc1').set('expr', 'u');

model.mesh('mesh1').autoMeshSize(4);
model.mesh('mesh1').run;
model.mesh('mesh1').autoMeshSize(3);
model.mesh('mesh1').run;
model.mesh('mesh1').run;

model.sol('sol1').study('std3');
model.sol('sol1').feature.remove('t1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').feature.create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std3');
model.sol('sol1').feature('st1').set('studystep', 'time');
model.sol('sol1').feature.create('v1', 'Variables');
model.sol('sol1').feature.create('t1', 'Time');
model.sol('sol1').feature('t1').set('tlist', 'range(0,0.1,1)');
model.sol('sol1').feature('t1').set('plot', 'off');
model.sol('sol1').feature('t1').set('plotgroup', 'pg1');
model.sol('sol1').feature('t1').set('plotfreq', 'tout');
model.sol('sol1').feature('t1').set('probesel', 'all');
model.sol('sol1').feature('t1').set('probes', {});
model.sol('sol1').feature('t1').set('probefreq', 'tsteps');
model.sol('sol1').feature('t1').set('control', 'time');
model.sol('sol1').feature('t1').feature.create('fc1', 'FullyCoupled');
model.sol('sol1').feature('t1').feature.remove('fcDef');
model.sol('sol1').attach('std3');
model.sol('sol1').runAll;

model.result('pg1').run;
model.result('pg1').set('data', 'none');
model.result('pg1').run;
model.result('pg1').set('data', 'dset1');
model.result('pg1').run;
model.result('pg1').set('view', 'view1');

model.view('view1').set('showgrid', 'on');
model.view('view1').set('transparency', 'on');

model.physics('c').feature('init1').set('u', 1, '1');
model.physics('c').feature('init1').set('v', 1, '1');
model.physics('c').feature('init1').set('u', 1, '2');

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');

model.sol('sol1').study('std3');
model.sol('sol1').feature.remove('t1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').feature.create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std3');
model.sol('sol1').feature('st1').set('studystep', 'time');
model.sol('sol1').feature.create('v1', 'Variables');
model.sol('sol1').feature.create('t1', 'Time');
model.sol('sol1').feature('t1').set('tlist', 'range(0,0.1,1)');
model.sol('sol1').feature('t1').set('plot', 'off');
model.sol('sol1').feature('t1').set('plotgroup', 'pg1');
model.sol('sol1').feature('t1').set('plotfreq', 'tout');
model.sol('sol1').feature('t1').set('probesel', 'all');
model.sol('sol1').feature('t1').set('probes', {});
model.sol('sol1').feature('t1').set('probefreq', 'tsteps');
model.sol('sol1').feature('t1').set('control', 'time');
model.sol('sol1').feature('t1').feature.create('fc1', 'FullyCoupled');
model.sol('sol1').feature('t1').feature.remove('fcDef');
model.sol('sol1').attach('std3');

model.geom('geom1').feature('tor1').setIndex('ax2', '2*pi', 0);
model.geom('geom1').feature('tor1').setIndex('ax2', '2*pi', 1);
model.geom('geom1').feature('tor1').set('rot', '360');
model.geom('geom1').feature('tor1').setIndex('ax2', '360', 0);
model.geom('geom1').feature('tor1').setIndex('ax2', '3360', 1);
model.geom('geom1').feature('tor1').setIndex('ax2', '360', 1);
model.geom('geom1').feature('tor1').set('createselection', 'on');
model.geom('geom1').feature('tor1').set('type', 'surface');
model.geom('geom1').runAll;

model.view('view1').set('transparency', 'off');

model.sol('sol1').study('std3');
model.sol('sol1').feature.remove('t1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').feature.create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std3');
model.sol('sol1').feature('st1').set('studystep', 'time');
model.sol('sol1').feature.create('v1', 'Variables');
model.sol('sol1').feature.create('t1', 'Time');
model.sol('sol1').feature('t1').set('tlist', 'range(0,0.1,1)');
model.sol('sol1').feature('t1').set('plot', 'off');
model.sol('sol1').feature('t1').set('plotgroup', 'pg1');
model.sol('sol1').feature('t1').set('plotfreq', 'tout');
model.sol('sol1').feature('t1').set('probesel', 'all');
model.sol('sol1').feature('t1').set('probes', {});
model.sol('sol1').feature('t1').set('probefreq', 'tsteps');
model.sol('sol1').feature('t1').set('control', 'time');
model.sol('sol1').feature('t1').feature.create('fc1', 'FullyCoupled');
model.sol('sol1').feature('t1').feature.remove('fcDef');
model.sol('sol1').attach('std3');
model.sol('sol1').runAll;

model.geom('geom1').feature('tor1').set('angle', '0');
model.geom('geom1').feature('tor1').setIndex('ax2', '0', 0);
model.geom('geom1').feature('tor1').setIndex('ax2', '0', 1);
model.geom('geom1').feature('tor1').set('createselection', 'off');
model.geom('geom1').run('');
model.geom('geom1').feature('tor1').set('angle', '360');
model.geom('geom1').runAll;
model.geom('geom1').feature('tor1').setIndex('ax2', '360', 0);
model.geom('geom1').feature('tor1').setIndex('ax2', '180', 1);
model.geom('geom1').runAll;
model.geom('geom1').runAll;
model.geom('geom1').run('tor1');
model.geom('geom1').feature.create('csur1', 'ConvertToSurface');
model.geom('geom1').run('tor1');
model.geom('geom1').feature('csur1').set('createselection', 'on');
model.geom('geom1').feature('tor1').set('createselection', 'on');
model.geom('geom1').feature.remove('csur1');
model.geom('geom1').runAll;
model.geom('geom1').run('tor1');
model.geom('geom1').run('tor1');
model.geom('geom1').feature('tor1').set('type', 'solid');
model.geom('geom1').runAll;
model.geom('geom1').feature('tor1').set('type', 'surface');
model.geom('geom1').runAll;

model.sol('sol1').study('std3');
model.sol('sol1').feature.remove('t1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').feature.create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std3');
model.sol('sol1').feature('st1').set('studystep', 'time');
model.sol('sol1').feature.create('v1', 'Variables');
model.sol('sol1').feature.create('t1', 'Time');
model.sol('sol1').feature('t1').set('tlist', 'range(0,0.1,1)');
model.sol('sol1').feature('t1').set('plot', 'off');
model.sol('sol1').feature('t1').set('plotgroup', 'pg1');
model.sol('sol1').feature('t1').set('plotfreq', 'tout');
model.sol('sol1').feature('t1').set('probesel', 'all');
model.sol('sol1').feature('t1').set('probes', {});
model.sol('sol1').feature('t1').set('probefreq', 'tsteps');
model.sol('sol1').feature('t1').set('control', 'time');
model.sol('sol1').feature('t1').feature.create('fc1', 'FullyCoupled');
model.sol('sol1').feature('t1').feature.remove('fcDef');
model.sol('sol1').attach('std3');
model.sol('sol1').runAll;

model.view('view1').set('scenelight', 'on');
model.view('view1').set('transparency', 'on');
model.view('view1').set('renderwireframe', false);
model.view('view1').set('transparency', 'off');
model.view('view1').set('renderwireframe', false);
model.view('view1').set('scenelight', 'on');
model.view('view1').set('drawhidestatus', 'ignore');

model.sol('sol1').study('std3');
model.sol('sol1').feature.remove('t1');
model.sol('sol1').feature.remove('v1');
model.sol('sol1').feature.remove('st1');
model.sol('sol1').feature.create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std3');
model.sol('sol1').feature('st1').set('studystep', 'time');
model.sol('sol1').feature.create('v1', 'Variables');
model.sol('sol1').feature.create('t1', 'Time');
model.sol('sol1').feature('t1').set('tlist', 'range(0,0.1,1)');
model.sol('sol1').feature('t1').set('plot', 'off');
model.sol('sol1').feature('t1').set('plotgroup', 'pg1');
model.sol('sol1').feature('t1').set('plotfreq', 'tout');
model.sol('sol1').feature('t1').set('probesel', 'all');
model.sol('sol1').feature('t1').set('probes', {});
model.sol('sol1').feature('t1').set('probefreq', 'tsteps');
model.sol('sol1').feature('t1').set('control', 'time');
model.sol('sol1').feature('t1').feature.create('fc1', 'FullyCoupled');
model.sol('sol1').feature('t1').feature.remove('fcDef');
model.sol('sol1').attach('std3');
model.sol('sol1').runAll;
model.sol('sol1').updateSolution;

model.name('torus_start.mph');

model.result('pg1').run;
model.result('pg1').set('frametype', 'spatial');
model.result.dataset('dset1').set('frametype', 'spatial');
model.result.dataset('dset1').run;

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature('time').set('tlist', 'range(0,0.1,10)');
model.study('std3').feature('time').set('plot', 'on');

model.geom('geom1').feature('tor1').setIndex('ax2', '0', 1);
model.geom('geom1').feature('tor1').setIndex('ax2', '0', 0);
model.geom('geom1').runAll;
model.geom('geom1').feature('tor1').set('createselection', 'off');

model.sol('sol1').feature.create('st2', 'StudyStep');
model.sol('sol1').feature('st2').name('Compile Equations: Time Dependent (2)');
model.sol('sol1').feature.remove('st2');
model.sol('sol1').feature('v1').feature('SDTorus_u').set('scalemethod', 'auto');
model.sol('sol1').feature('v1').feature('SDTorus_v').set('scalemethod', 'auto');
model.sol('sol1').feature.create('s1', 'Stationary');

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');

model.sol('sol1').feature.remove('t1');
model.sol('sol1').runAll;

model.view('view1').set('transparency', 'on');
model.view('view1').set('renderwireframe', false);

model.geom('geom1').runAll;

model.sol('sol1').runAll;

model.result('pg1').feature.remove('slc1');
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').feature.create('vol1', 'Volume');
model.result('pg1').feature.remove('vol1');
model.result('pg1').run;
model.result('pg1').feature.create('vol1', 'Volume');

model.sol('sol1').runAll;

model.result('pg1').run;

model.view('view1').set('transparency', 'off');

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');

model.sol('sol1').feature.create('t1', 'Time');
model.sol('sol1').feature('s1').active(false);
model.sol('sol1').runAll;

model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').feature('vol1').set('wireframe', 'on');
model.result('pg1').feature('vol1').set('descr', 'Dependent variable u ');
model.result('pg1').feature('vol1').set('expr', 'uroot.u');
model.result('pg1').feature('vol1').set('descr', 'uroot.u ');
model.result('pg1').feature('vol1').set('expr', 'uroot.uroot.v');
model.result('pg1').feature('vol1').set('data', 'parent');
model.result('pg1').feature('vol1').set('expr', 'uroot.u');

model.sol('sol1').runAll;

model.result('pg1').run;
model.result('pg1').set('windowtitle', 'Graphics');
model.result('pg1').set('window', 'window1');
model.result('pg1').set('windowtitle', 'Plot 1');
model.result('pg1').set('window', 'graphics');
model.result('pg1').set('frametype', 'mesh');

model.sol('sol1').runAll;

model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').feature('vol1').active(false);
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').feature.create('surf1', 'Surface');
model.result('pg1').run;
model.result('pg1').feature.remove('surf1');
model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').feature('vol1').active(true);
model.result('pg1').run;
model.result('pg1').feature('vol1').set('data', 'dset1');

model.view('view1').set('scenelight', 'on');
model.view('view1').set('transparency', 'off');

model.result('pg1').set('renderdatacached', true);
model.result('pg1').set('lastinputmode', 'evaluate');
model.result('pg1').set('inputmode', 'evaluate');
model.result('pg1').set('lastinputmode', 'evaluate');
model.result('pg1').set('inputmode', 'evaluate');
model.result('pg1').set('lastinputmode', 'evaluate');
model.result('pg1').set('inputmode', 'linefirst');
model.result('pg1').set('lastinputmode', 'linefirst');
model.result('pg1').set('inputmode', 'linefirst');
model.result('pg1').set('lastinputmode', 'linefirst');
model.result('pg1').set('inputmode', 'evaluate');
model.result('pg1').feature.create('slc1', 'Slice');
model.result('pg1').feature('slc1').set('expr', 'uroot.u');
model.result('pg1').feature('slc1').set('descr', 'uroot.u');
model.result('pg1').feature('slc1').set('unit', '');
model.result('pg1').feature.create('iso1', 'Isosurface');
model.result('pg1').feature('iso1').set('expr', 'uroot.u');
model.result('pg1').feature('iso1').set('descr', 'uroot.u');
model.result('pg1').feature('iso1').set('unit', '');
model.result('pg1').feature.create('vol2', 'Volume');
model.result('pg1').feature('vol2').set('expr', 'uroot.u');
model.result('pg1').feature('vol2').set('descr', 'uroot.u');
model.result('pg1').feature('vol2').set('unit', '');
model.result('pg1').feature.create('surf1', 'Surface');
model.result('pg1').feature('surf1').set('expr', 'uroot.u');
model.result('pg1').feature('surf1').set('descr', 'uroot.u');
model.result('pg1').feature('surf1').set('unit', '');
model.result.export.create('play1', 'Player');
model.result.export('play1').set('plotgroup', 'pg1');
model.result.export('play1').set('autoplay', true);
model.result('pg1').set('window', 'graphics');
model.result('pg1').run;
model.result('pg1').set('window', 'graphics');
model.result('pg1').set('windowtitle', 'Graphics');
model.result.export('play1').set('showframe', 1);
model.result.export('play1').set('autoplay', false);
model.result.export('play1').set('showframe', 1);
model.result('pg1').feature.remove('slc1');
model.result.remove('pg1');

model.sol('sol1').runAll;

model.result('pg2').set('window', 'window2');
model.result('pg2').set('windowtitle', 'Plot 2');
model.result('pg2').run;

model.sol('sol1').feature('t1').set('control', 'time');

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature('time').set('tlist', 'range(0,0.1,10)');

model.name('torus_start.mph');

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');

model.sol('sol1').runAll;

model.result('pg2').set('window', 'window2');
model.result('pg2').set('windowtitle', 'Plot 2');
model.result('pg2').run;

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');

model.result('pg2').set('window', 'window2');
model.result('pg2').set('windowtitle', 'Plot 2');
model.result('pg2').run;
model.result.create('pg3', 'PlotGroup3D');
model.result('pg3').run;
model.result('pg3').set('data', 'dset1');

model.variable('var1').remove('u');
model.variable('var1').remove('v');

model.result('pg3').run;
model.result('pg3').set('frametype', 'spatial');
model.result('pg2').set('windowtitle', 'Plot 2');
model.result('pg3').set('windowtitle', 'Graphics');
model.result('pg3').set('window', 'window3');
model.result('pg3').set('windowtitle', 'Plot 3');
model.result('pg3').set('window', 'graphics');
model.result('pg3').run;
model.result('pg3').set('titleactive', 'off');

model.sol('sol1').runAll;

model.result('pg2').set('window', 'window2');
model.result('pg2').set('windowtitle', 'Plot 2');
model.result('pg2').run;
model.result('pg2').active(false);
model.result('pg2').set('window', 'window2');
model.result('pg2').set('windowtitle', 'Plot 2');
model.result('pg2').run;

model.sol('sol1').runAll;

model.result('pg3').run;

model.view('view1').set('transparency', 'off');

model.result('pg3').set('legendpos', 'rightdouble');
model.result('pg3').set('view', 'view1');
model.result('pg3').set('titleactive', 'on');
model.result('pg3').set('edges', 'off');
model.result('pg3').run;
model.result('pg3').set('edges', 'on');
model.result('pg3').run;
model.result('pg3').feature.create('surf1', 'Surface');
model.result('pg3').run;
model.result('pg3').feature.remove('surf1');
model.result('pg3').run;
model.result('pg3').feature.create('vol1', 'Volume');
model.result('pg3').feature('vol1').set('descractive', 'on');
model.result('pg3').feature('vol1').set('data', 'dset1');
model.result('pg3').run;
model.result('pg3').feature('vol1').set('smooth', 'everywhere');
model.result('pg3').feature('vol1').set('recover', 'ppr');
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').feature.create('mesh1', 'Mesh');
model.result('pg3').feature('mesh1').set('elemtype2', 'tri');

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');

model.sol('sol1').runAll;

model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').feature('vol1').set('expr', 'v');
model.result('pg3').feature('vol1').set('solnum', '3');
model.result('pg3').feature('vol1').set('expr', 'u');
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;

model.physics('c').feature('init1').set('u', 1, '0');

model.sol('sol1').runAll;

model.result('pg3').run;
model.result('pg3').set('view', 'auto');
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').feature('vol1').set('solnum', '21');
model.result('pg3').run;

model.physics('c').feature('init1').set('v', 1, '0');

model.sol('sol1').runAll;

model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;

model.coordSystem('sys1').set('frametype', 'material');

model.sol('sol1').runAll;

model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').set('data', 'none');
model.result('pg3').run;
model.result('pg3').set('data', 'dset1');
model.result('pg3').run;
model.result.dataset.remove('cln1');
model.result.dataset('dset1').run;
model.result('pg3').run;
model.result.dataset('dset1').run;

model.sol('sol1').runAll;

model.result('pg3').run;
model.result.dataset('dset1').set('scalefactor', '0');
model.result.dataset('dset1').run;
model.result.dataset('dset1').set('scalefactor', '1');
model.result.dataset('dset1').run;
model.result.dataset.create('fun1', 'Function2D');
model.result.dataset.create('mesh1', 'Mesh');
model.result.dataset('mesh1').run;
model.result.dataset('mesh1').run;
model.result.dataset('fun1').run;
model.result.dataset.remove('fun1');
model.result.dataset('dset1').run;
model.result.dataset('mesh1').run;
model.result.dataset.create('ps1', 'ParSurface');
model.result.dataset.remove('ps1');

model.sol('sol1').runAll;

model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').feature('vol1').set('expr', 'v');
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').feature('vol1').set('solnum', '2');
model.result('pg3').run;

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');

model.sol('sol1').runAll;

model.result('pg3').run;

model.geom('geom1').runAll;

model.result.create('pg4', 3);
model.result('pg4').set('data', 'mesh1');
model.result('pg4').feature.create('mesh1', 'Mesh');

model.mesh('mesh1').automatic(false);
model.mesh('mesh1').run;

model.physics('c').feature('init1').set('u', 1, '10');
model.physics('c').feature('init1').set('v', 1, '2');

model.sol('sol1').runAll;

model.result('pg3').run;
model.result('pg4').run;
model.result('pg4').run;
model.result('pg4').run;
model.result('pg4').run;
model.result.remove('pg4');
model.result('pg3').run;
model.result.table.create('tbl1', 'Table');
model.result.table.remove('tbl1');
model.result.dataset('dset1').set('phase', '0.5');
model.result.dataset('dset1').run;
model.result.dataset('dset1').run;
model.result.dataset('dset1').set('phase', '45');
model.result.dataset('dset1').run;
model.result.dataset('dset1').set('phase', '0.5');
model.result.dataset.create('int1', 'Integral');
model.result.dataset.remove('int1');
model.result.numerical.create('gev1', 'EvalGlobal');
model.result.table.create('tbl1', 'Table');
model.result.table('tbl1').comments('Global Evaluation 1 (u)');
model.result.numerical('gev1').set('table', 'tbl1');
model.result.numerical('gev1').set('descractive', 'on');

model.sol('sol1').runAll;

model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').feature('vol1').set('data', 'none');
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').feature('vol1').set('expr', 'uv');
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;

model.sol('sol1').runAll;

model.result('pg3').run;
model.result.export.create('data1', 'Data');
model.result.export('data1').setIndex('expr', 'u', 0);
model.result.export('data1').setIndex('expr', 'v', 1);
model.result.export('data1').set('filename', '/Users/sebastianboie/Desktop/text_torus_soln.txt');
model.result.export('data1').run;

model.geom('geom1').runAll;

model.mesh('mesh1').run;
model.mesh.remove('mesh1');

model.geom('geom1').run('tor1');
model.geom('geom1').runAll;
model.geom('geom1').run('tor1');
model.geom('geom1').run('fin');
model.geom('geom1').run('tor1');
model.geom('geom1').feature.create('csur1', 'ConvertToSurface');
model.geom('geom1').feature('csur1').selection('input').set({'tor1'});
model.geom('geom1').feature('csur1').set('keep', 'off');
model.geom('geom1').run('csur1');
model.geom('geom1').runPre('csur1');
model.geom('geom1').run('csur1');
model.geom('geom1').run('fin');
model.geom('geom1').run('csur1');

model.result('pg3').run;
model.result('pg3').active(false);
model.result('pg3').run;
model.result.create('pg4', 'PlotGroup2D');
model.result('pg4').run;
model.result('pg4').active(false);
model.result('pg4').run;
model.result('pg3').run;
model.result('pg3').active(true);
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').feature('vol1').active(false);
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').feature.create('surf1', 'Surface');
model.result('pg3').feature('surf1').set('data', 'dset1');

model.sol('sol1').runAll;

model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').feature('surf1').active(false);
model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').feature('vol1').active(true);
model.result('pg3').run;

model.geom('geom1').feature('csur1').active(false);

model.sol('sol1').feature('v1').feature('SDTorus_u').set('scalemethod', 'auto');
model.sol('sol1').run('v1');

model.result('pg3').run;

model.sol('sol1').feature('v1').set('control', 'user');
model.sol('sol1').runAll;

model.result('pg3').run;
model.result.export('data1').run;
model.result('pg3').run;
model.result('pg3').feature('vol1').set('expr', 'u');
model.result('pg3').run;
model.result('pg3').feature('vol1').set('data', 'dset1');
model.result('pg3').run;
model.result('pg3').run;
model.result.numerical('gev1').set('table', 'tbl1');

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');

model.geom('geom1').run('tor1');
model.geom('geom1').feature('tor1').set('createselection', 'on');
model.geom('geom1').run('tor1');
model.geom('geom1').run('tor1');
model.geom('geom1').runAll;

model.mesh.create('mesh1', 'geom1');
model.mesh('mesh1').automatic(false);
model.mesh('mesh1').automatic(true);
model.mesh('mesh1').autoMeshSize(3);
model.mesh('mesh1').run;

model.result('pg3').run;

model.sol('sol1').runAll;

model.result('pg3').run;

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature('time').mglevel.create('mgl1');
model.study('std3').feature('time').mglevel('mgl1').set('geomselection', 'geom1');
model.study('std3').feature('time').mglevel('mgl1').set('physselection', 'c');
model.study('std3').feature('time').mglevel('mgl1').set('usesol', 'on');
model.study('std3').feature('time').mglevel('mgl1').set('notstudy', 'std3');

model.sol('sol1').runAll;

model.result('pg3').run;

model.sol('sol1').feature('v1').set('initmethod', 'sol');
model.sol('sol1').feature('v1').set('initsol', 'sol1');
model.sol('sol1').feature('v1').set('solnum', 'auto');
model.sol('sol1').feature('v1').set('scalemethod', 'init');
model.sol('sol1').feature('v1').set('notsol', 'sol1');
model.sol('sol1').runAll;

model.result('pg3').run;

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature('time').mglevel('mgl1').set('geomselection', 'geom1');
model.study('std3').feature('time').mglevel('mgl1').set('physselection', 'c');
model.study('std3').feature('time').mglevel('mgl1').set('usesol', 'on');
model.study('std3').feature('time').mglevel('mgl1').set('notsolmethod', 'sol');

model.sol('sol1').runAll;

model.result('pg3').run;

model.study('std3').feature('time').mglevel('mgl1').set('geomselection', 'geom1');
model.study('std3').feature('time').mglevel('mgl1').set('physselection', 'c');
model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature('time').mglevel('mgl1').set('geomselection', 'geom1');
model.study('std3').feature('time').mglevel('mgl1').set('physselection', 'c');
model.study('std3').feature('time').mglevel.remove('mgl1');
model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature('time').set('tlist', 'range(0,360,10)');
model.study('std3').feature('time').set('plotfreq', 'tsteps');
model.study('std3').feature('time').set('probesel', 'all');
model.study('std3').feature('time').set('autoremesh', 'on');
model.study('std3').feature('time').set('optimization', 'on');

model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').run;

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature('time').set('rtolactive', 'on');
model.study('std3').feature('time').set('tlist', 'range(0,0.1,10) range(0,0.1,10)');

model.result('pg3').run;

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature('time').set('tlist', 'range(0,0.1,10)');

model.sol('sol1').runAll;

model.result('pg3').run;

model.sol('sol1').feature('v1').set('notsolmethod', 'sol');
model.sol('sol1').feature('t1').set('control', 'user');
model.sol('sol1').feature('t1').feature.remove('fcDef');
model.sol('sol1').feature.remove('t1');
model.sol('sol1').runAll;

model.result('pg3').run;

model.sol('sol1').feature.create('t1', 'Time');
model.sol('sol1').feature('t1').feature('fcDef').set('plot', 'on');
model.sol('sol1').runAll;

model.result('pg3').run;
model.result.export('data1').run;

model.sol('sol1').runAll;

model.result('pg3').run;

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature('time').set('plot', 'off');
model.study('std3').feature('time').set('usesol', 'on');
model.study('std3').feature('time').set('notstudy', 'std3');
model.study('std3').feature('time').set('notsolnum', 'all');

model.sol('sol1').runAll;

model.result('pg3').run;
model.result.numerical('gev1').set('table', 'tbl1');

model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature('time').set('geomselection', 'geom1');
model.study('std3').feature('time').set('physselection', 'c');
model.study('std3').feature.remove('time');

model.sol('sol1').feature('v1').set('solnum', 'auto');
model.sol('sol1').run('v1');

model.result('pg3').run;
model.result.table.remove('tbl1');

model.sol('sol1').runAll;

model.result('pg3').run;

model.variable.remove('var1');

model.view('view1').set('renderwireframe', 'off');

model.material.create('mat1');
model.material.remove('mat1');

model.physics('c').feature('init1').featureInfo('info').name('Equation View');
model.physics('c').feature('init1').set('ut', 1, '1');
model.physics('c').feature('init1').set('vt', 1, '-1');

model.sol('sol1').runAll;

model.result('pg3').run;
model.result('pg3').run;
model.result('pg3').active(false);
model.result('pg3').run;
model.result.numerical('gev1').set('expr', 'u v');
model.result.table.create('tbl1', 'Table');
model.result.table('tbl1').comments('Global Evaluation 1 (u v)');
model.result.numerical('gev1').set('table', 'tbl1');
model.result.numerical('gev1').set('expr', 'u');
model.result.table.create('tbl2', 'Table');
model.result.table('tbl2').comments('Global Evaluation 1 (u)');
model.result.numerical('gev1').set('table', 'tbl2');

model.geom('geom1').feature('tor1').active(false);

model.sol('sol1').feature('v1').feature('SDTorus_u').set('scalemethod', 'init');
model.sol('sol1').feature('v1').feature('SDTorus_v').set('scalemethod', 'init');
model.sol('sol1').feature('v1').feature('SDTorus_v').set('out', 'on');
model.sol('sol1').runAll;

model.batch.create('p1', 'Parametric');
model.batch('p1').study('std3');
model.batch.remove('p1');

model.sol('sol1').runAll;

model.geom('geom1').feature('tor1').active(true);
model.geom('geom1').feature('csur1').active(true);
model.geom('geom1').run('tor1');
model.geom('geom1').run('tor1');
model.geom('geom1').feature('csur1').active(false);

model.mesh('mesh1').run;

model.sol('sol1').run('st1');
model.sol('sol1').runAll;

out = model;