function out = model
%
% model.m
%
% Model exported on Oct 29 2012, 21:54 by COMSOL 4.3.1.115.

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

model = ModelUtil.create('Model');

model.modelPath('C:\Users\Superman\AppData\Local\Opera\Opera\temporary_downloads');

model.modelNode.create('mod1');

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

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

model.physics.create('spf', 'LaminarFlow', 'geom1');

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

model.geom('geom1').run('');
model.geom('geom1').feature.create('wp1', 'WorkPlane');
model.geom('geom1').feature('wp1').geom.feature.create('r1', 'Rectangle');
model.geom('geom1').feature('wp1').geom.feature('r1').setIndex('size', '4', 0);
model.geom('geom1').feature('wp1').geom.feature('r1').setIndex('size', '3', 1);
model.geom('geom1').feature('wp1').geom.runAll;
model.geom('geom1').run('wp1');
model.geom('geom1').feature.create('ext1', 'Extrude');
model.geom('geom1').feature('ext1').setIndex('distance', '5.', 0);
model.geom('geom1').runAll;
model.geom('geom1').run;

model.material.create('mat1');
model.material('mat1').name('Air');
model.material('mat1').propertyGroup('def').set('relpermeability', '1');
model.material('mat1').propertyGroup('def').set('relpermittivity', '1');
model.material('mat1').propertyGroup('def').set('dynamicviscosity', 'eta(T[1/K])[Pa*s]');
model.material('mat1').propertyGroup('def').set('ratioofspecificheat', '1.4');
model.material('mat1').propertyGroup('def').set('electricconductivity', '0[S/m]');
model.material('mat1').propertyGroup('def').set('heatcapacity', 'Cp(T[1/K])[J/(kg*K)]');
model.material('mat1').propertyGroup('def').set('density', 'rho(pA[1/Pa],T[1/K])[kg/m^3]');
model.material('mat1').propertyGroup('def').set('thermalconductivity', 'k(T[1/K])[W/(m*K)]');
model.material('mat1').propertyGroup('def').set('soundspeed', 'cs(T[1/K])[m/s]');
model.material('mat1').propertyGroup('def').func.create('eta', 'Piecewise');
model.material('mat1').propertyGroup('def').func('eta').set('funcname', 'eta');
model.material('mat1').propertyGroup('def').func('eta').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('eta').set('extrap', 'constant');
model.material('mat1').propertyGroup('def').func('eta').set('pieces', {'200.0' '1600.0' '-8.38278E-7+8.35717342E-8*T^1-7.69429583E-11*T^2+4.6437266E-14*T^3-1.06585607E-17*T^4'});
model.material('mat1').propertyGroup('def').func.create('Cp', 'Piecewise');
model.material('mat1').propertyGroup('def').func('Cp').set('funcname', 'Cp');
model.material('mat1').propertyGroup('def').func('Cp').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('Cp').set('extrap', 'constant');
model.material('mat1').propertyGroup('def').func('Cp').set('pieces', {'200.0' '1600.0' '1047.63657-0.372589265*T^1+9.45304214E-4*T^2-6.02409443E-7*T^3+1.2858961E-10*T^4'});
model.material('mat1').propertyGroup('def').func.create('rho', 'Analytic');
model.material('mat1').propertyGroup('def').func('rho').set('funcname', 'rho');
model.material('mat1').propertyGroup('def').func('rho').set('args', {'pA' 'T'});
model.material('mat1').propertyGroup('def').func('rho').set('expr', 'pA*0.02897/8.314/T');
model.material('mat1').propertyGroup('def').func('rho').set('dermethod', 'manual');
model.material('mat1').propertyGroup('def').func('rho').set('argders', {'pA' 'd(pA*0.02897/8.314/T,pA)'; 'T' 'd(pA*0.02897/8.314/T,T)'});
model.material('mat1').propertyGroup('def').func.create('k', 'Piecewise');
model.material('mat1').propertyGroup('def').func('k').set('funcname', 'k');
model.material('mat1').propertyGroup('def').func('k').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('k').set('extrap', 'constant');
model.material('mat1').propertyGroup('def').func('k').set('pieces', {'200.0' '1600.0' '-0.00227583562+1.15480022E-4*T^1-7.90252856E-8*T^2+4.11702505E-11*T^3-7.43864331E-15*T^4'});
model.material('mat1').propertyGroup('def').func.create('cs', 'Analytic');
model.material('mat1').propertyGroup('def').func('cs').set('funcname', 'cs');
model.material('mat1').propertyGroup('def').func('cs').set('args', {'T'});
model.material('mat1').propertyGroup('def').func('cs').set('expr', 'sqrt(1.4*287*T)');
model.material('mat1').propertyGroup('def').func('cs').set('dermethod', 'manual');
model.material('mat1').propertyGroup('def').func('cs').set('argders', {'T' 'd(sqrt(1.4*287*T),T)'});
model.material('mat1').propertyGroup('def').addInput('temperature');
model.material('mat1').propertyGroup('def').addInput('pressure');

model.physics('spf').feature.create('inl1', 'Inlet', 2);
model.physics('spf').feature('inl1').selection.all;
model.physics('spf').feature('inl1').selection.set([1]);
model.physics('spf').feature('inl1').set('BoundaryCondition', 1, 'LaminarInflow');
model.physics('spf').feature('inl1').set('Uav', 1, '1');
model.physics('spf').feature.create('out1', 'Outlet', 2);
model.physics('spf').feature('out1').selection.set([4]);

model.mesh('mesh1').feature.create('ftet1', 'FreeTet');
model.mesh('mesh1').run;

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.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('itrestart', 50);
model.sol('sol1').feature('s1').feature('i1').set('rhob', 20);
model.sol('sol1').feature('s1').feature('i1').set('maxlinit', 200);
model.sol('sol1').feature('s1').feature('i1').set('nlinnormuse', 'on');
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('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
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.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', 'spf.U');
model.result('pg1').feature('surf1').set('descr', 'Velocity magnitude');
model.result.create('pg2', 3);
model.result('pg2').set('data', 'dset1');
model.result('pg2').feature.create('slc1', 'Slice');
model.result('pg2').feature('slc1').set('expr', 'spf.U');
model.result('pg2').feature('slc1').set('descr', 'Velocity magnitude');
model.result('pg1').set('windowtitle', 'Graphics');
model.result('pg2').set('windowtitle', 'Graphics');
model.result('pg1').run;

model.physics('spf').feature('out1').selection.set([6]);

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.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('itrestart', 50);
model.sol('sol1').feature('s1').feature('i1').set('rhob', 20);
model.sol('sol1').feature('s1').feature('i1').set('maxlinit', 200);
model.sol('sol1').feature('s1').feature('i1').set('nlinnormuse', 'on');
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('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
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.remove('fcDef');
model.sol('sol1').attach('std1');
model.sol('sol1').runAll;

model.result('pg1').run;

model.physics('spf').feature('inl1').set('LaminarInflowOption', 1, 'V0');
model.physics('spf').feature('inl1').set('V0', 1, '50');

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.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('itrestart', 50);
model.sol('sol1').feature('s1').feature('i1').set('rhob', 20);
model.sol('sol1').feature('s1').feature('i1').set('maxlinit', 200);
model.sol('sol1').feature('s1').feature('i1').set('nlinnormuse', 'on');
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('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
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.remove('fcDef');
model.sol('sol1').attach('std1');
model.sol('sol1').runAll;

model.result('pg1').run;

model.name('hydrogen.mph');

model.result('pg2').run;

model.geom('geom1').run('ext1');
model.geom('geom1').feature.create('blk1', 'Block');
model.geom('geom1').run('blk1');
model.geom('geom1').feature('blk1').setIndex('ax3', '1', 0);
model.geom('geom1').feature('blk1').setIndex('ax3', '1', 1);
model.geom('geom1').run('blk1');
model.geom('geom1').feature('blk1').setIndex('ax3', '0', 0);
model.geom('geom1').feature('blk1').setIndex('ax3', '0', 1);
model.geom('geom1').feature('blk1').setIndex('ax3', '0', 2);
model.geom('geom1').feature('blk1').setIndex('pos', '1', 0);
model.geom('geom1').feature('blk1').setIndex('pos', '1', 1);
model.geom('geom1').feature('blk1').setIndex('pos', '1', 2);
model.geom('geom1').run('blk1');
model.geom('geom1').run('blk1');
model.geom('geom1').feature.remove('blk1');
model.geom('geom1').run('ext1');
model.geom('geom1').feature.create('blk1', 'Block');
model.geom('geom1').feature('blk1').setIndex('pos', '1', 0);
model.geom('geom1').feature('blk1').setIndex('pos', '1', 2);
model.geom('geom1').run('blk1');
model.geom('geom1').feature('blk1').set('type', 'surface');
model.geom('geom1').run;

model.physics('spf').feature('out1').selection.set([7]);

model.mesh('mesh1').feature.create('ftet2', 'FreeTet');
model.mesh('mesh1').run;

model.geom('geom1').run('blk1');
model.geom('geom1').feature.create('spl1', 'Split');
model.geom('geom1').feature('spl1').selection('input').set({'blk1' 'ext1'});
model.geom('geom1').runAll;
model.geom('geom1').run;

model.physics('spf').feature('out1').selection.set([7]);

model.mesh('mesh1').feature.create('ftet3', 'FreeTet');
model.mesh('mesh1').run;
model.mesh('mesh1').feature.remove('ftet2');
model.mesh('mesh1').feature.remove('ftet1');

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.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('itrestart', 50);
model.sol('sol1').feature('s1').feature('i1').set('rhob', 20);
model.sol('sol1').feature('s1').feature('i1').set('maxlinit', 200);
model.sol('sol1').feature('s1').feature('i1').set('nlinnormuse', 'on');
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('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
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.remove('fcDef');
model.sol('sol1').attach('std1');

model.result('pg1').run;

model.material.remove('mat1');
model.material.create('mat1');
model.material('mat1').name('Air');
model.material('mat1').propertyGroup('def').set('relpermeability', '1');
model.material('mat1').propertyGroup('def').set('relpermittivity', '1');
model.material('mat1').propertyGroup('def').set('dynamicviscosity', 'eta(T[1/K])[Pa*s]');
model.material('mat1').propertyGroup('def').set('ratioofspecificheat', '1.4');
model.material('mat1').propertyGroup('def').set('electricconductivity', '0[S/m]');
model.material('mat1').propertyGroup('def').set('heatcapacity', 'Cp(T[1/K])[J/(kg*K)]');
model.material('mat1').propertyGroup('def').set('density', 'rho(pA[1/Pa],T[1/K])[kg/m^3]');
model.material('mat1').propertyGroup('def').set('thermalconductivity', 'k(T[1/K])[W/(m*K)]');
model.material('mat1').propertyGroup('def').set('soundspeed', 'cs(T[1/K])[m/s]');
model.material('mat1').propertyGroup('def').func.create('eta', 'Piecewise');
model.material('mat1').propertyGroup('def').func('eta').set('funcname', 'eta');
model.material('mat1').propertyGroup('def').func('eta').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('eta').set('extrap', 'constant');
model.material('mat1').propertyGroup('def').func('eta').set('pieces', {'200.0' '1600.0' '-8.38278E-7+8.35717342E-8*T^1-7.69429583E-11*T^2+4.6437266E-14*T^3-1.06585607E-17*T^4'});
model.material('mat1').propertyGroup('def').func.create('Cp', 'Piecewise');
model.material('mat1').propertyGroup('def').func('Cp').set('funcname', 'Cp');
model.material('mat1').propertyGroup('def').func('Cp').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('Cp').set('extrap', 'constant');
model.material('mat1').propertyGroup('def').func('Cp').set('pieces', {'200.0' '1600.0' '1047.63657-0.372589265*T^1+9.45304214E-4*T^2-6.02409443E-7*T^3+1.2858961E-10*T^4'});
model.material('mat1').propertyGroup('def').func.create('rho', 'Analytic');
model.material('mat1').propertyGroup('def').func('rho').set('funcname', 'rho');
model.material('mat1').propertyGroup('def').func('rho').set('args', {'pA' 'T'});
model.material('mat1').propertyGroup('def').func('rho').set('expr', 'pA*0.02897/8.314/T');
model.material('mat1').propertyGroup('def').func('rho').set('dermethod', 'manual');
model.material('mat1').propertyGroup('def').func('rho').set('argders', {'pA' 'd(pA*0.02897/8.314/T,pA)'; 'T' 'd(pA*0.02897/8.314/T,T)'});
model.material('mat1').propertyGroup('def').func.create('k', 'Piecewise');
model.material('mat1').propertyGroup('def').func('k').set('funcname', 'k');
model.material('mat1').propertyGroup('def').func('k').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('k').set('extrap', 'constant');
model.material('mat1').propertyGroup('def').func('k').set('pieces', {'200.0' '1600.0' '-0.00227583562+1.15480022E-4*T^1-7.90252856E-8*T^2+4.11702505E-11*T^3-7.43864331E-15*T^4'});
model.material('mat1').propertyGroup('def').func.create('cs', 'Analytic');
model.material('mat1').propertyGroup('def').func('cs').set('funcname', 'cs');
model.material('mat1').propertyGroup('def').func('cs').set('args', {'T'});
model.material('mat1').propertyGroup('def').func('cs').set('expr', 'sqrt(1.4*287*T)');
model.material('mat1').propertyGroup('def').func('cs').set('dermethod', 'manual');
model.material('mat1').propertyGroup('def').func('cs').set('argders', {'T' 'd(sqrt(1.4*287*T),T)'});
model.material('mat1').propertyGroup('def').addInput('temperature');
model.material('mat1').propertyGroup('def').addInput('pressure');

model.physics('spf').feature('inl1').set('LaminarInflowOption', 1, 'Uav');

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.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('itrestart', 50);
model.sol('sol1').feature('s1').feature('i1').set('rhob', 20);
model.sol('sol1').feature('s1').feature('i1').set('maxlinit', 200);
model.sol('sol1').feature('s1').feature('i1').set('nlinnormuse', 'on');
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('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
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.remove('fcDef');
model.sol('sol1').attach('std1');

model.result('pg1').run;

model.physics('spf').feature('out1').set('BoundaryCondition', 1, 'LaminarOutflow');
model.physics('spf').feature('out1').set('LaminarOutflowOption', 1, 'p0_exit');

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.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('itrestart', 50);
model.sol('sol1').feature('s1').feature('i1').set('rhob', 20);
model.sol('sol1').feature('s1').feature('i1').set('maxlinit', 200);
model.sol('sol1').feature('s1').feature('i1').set('nlinnormuse', 'on');
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('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
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.remove('fcDef');
model.sol('sol1').attach('std1');

model.result('pg1').run;

model.physics('spf').feature('inl1').selection.set([7]);
model.physics('spf').feature('out1').selection.set([1]);

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.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('itrestart', 50);
model.sol('sol1').feature('s1').feature('i1').set('rhob', 20);
model.sol('sol1').feature('s1').feature('i1').set('maxlinit', 200);
model.sol('sol1').feature('s1').feature('i1').set('nlinnormuse', 'on');
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('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
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.remove('fcDef');
model.sol('sol1').attach('std1');

model.result('pg1').run;

model.physics('spf').feature('inl1').selection.set([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.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('itrestart', 50);
model.sol('sol1').feature('s1').feature('i1').set('rhob', 20);
model.sol('sol1').feature('s1').feature('i1').set('maxlinit', 200);
model.sol('sol1').feature('s1').feature('i1').set('nlinnormuse', 'on');
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('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('pr').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('va1', 'Vanka');
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('iter', 1);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('seconditer', 3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('relax', 0.3);
model.sol('sol1').feature('s1').feature('i1').feature('mg1').feature('po').feature('va1').set('vankavars', {'mod1_spf_Pinlinl1'});
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.remove('fcDef');
model.sol('sol1').attach('std1');

model.result('pg1').run;

model.physics('spf').feature('inl1').selection.set([]);
model.physics('spf').feature('out1').selection.set([]);
model.physics('spf').feature('out1').set('LaminarOutflowOption', 1, 'Uav');

model.geom('geom1').feature.remove('blk1');
model.geom('geom1').feature.remove('spl1');
model.geom('geom1').feature.remove('ext1');
model.geom('geom1').runAll;
model.geom('geom1').feature('wp1').geom.run('r1');
model.geom('geom1').feature('wp1').geom.feature.create('r2', 'Rectangle');
model.geom('geom1').feature('wp1').geom.feature('r2').setIndex('size', '2', 0);
model.geom('geom1').feature('wp1').geom.feature('r2').setIndex('size', '2', 1);
model.geom('geom1').feature('wp1').geom.runAll;
model.geom('geom1').feature('wp1').geom.runAll;
model.geom('geom1').feature('wp1').geom.feature.remove('r2');
model.geom('geom1').run('wp1');
model.geom('geom1').feature.create('ext1', 'Extrude');
model.geom('geom1').feature('ext1').setIndex('distance', '10', 0);
model.geom('geom1').runAll;
model.geom('geom1').run;

model.material('mat1').selection.set([1]);

model.physics('spf').feature('inl1').selection.set([3]);
model.physics('spf').feature('inl1').set('BoundaryCondition', 1, 'PressureNoViscousStress');
model.physics('spf').feature('inl1').set('BoundaryCondition', 1, 'Velocity');
model.physics('spf').feature('inl1').set('U0in', 1, '1');
model.physics('spf').feature('out1').set('BoundaryCondition', 1, 'Pressure');
model.physics('spf').feature('out1').selection.set([4]);

model.sol('sol1').runAll;

model.result('pg1').run;
model.result('pg1').run;
model.result('pg1').run;
model.result('pg2').run;
model.result('pg2').run;

model.geom('geom1').run('ext1');
model.geom('geom1').feature.create('wp2', 'WorkPlane');
model.geom('geom1').feature('wp2').geom.feature.create('sq1', 'Square');
model.geom('geom1').feature('wp2').geom.feature('sq1').setIndex('pos', '2', 0);
model.geom('geom1').feature('wp2').geom.feature('sq1').setIndex('pos', '1', 0);
model.geom('geom1').feature('wp2').geom.feature('sq1').setIndex('pos', '1', 1);
model.geom('geom1').feature('wp2').geom.run('sq1');
model.geom('geom1').feature('wp2').geom.feature('sq1').setIndex('pos', '.5', 0);
model.geom('geom1').feature('wp2').geom.feature('sq1').setIndex('pos', '.5', 1);
model.geom('geom1').feature('wp2').geom.run('sq1');
model.geom('geom1').feature('wp2').geom.feature('sq1').setIndex('pos', '0', 0);
model.geom('geom1').feature('wp2').geom.feature('sq1').setIndex('pos', '0', 1);
model.geom('geom1').feature('wp2').geom.runAll;
model.geom('geom1').feature('wp2').geom.runAll;
model.geom('geom1').feature('wp2').geom.run('sq1');
model.geom('geom1').run('wp2');
model.geom('geom1').feature.create('spl1', 'Split');
model.geom('geom1').feature('spl1').selection('input').set({'ext1' 'wp2'});
model.geom('geom1').run('spl1');
model.geom('geom1').run;

model.physics('spf').feature('inl1').selection.set([3]);

model.sol('sol1').runAll;

model.result('pg1').run;

model.name('hydrogen.mph');

model.result('pg1').run;

model.geom('geom1').run('spl1');
model.geom('geom1').feature.create('wp3', 'WorkPlane');
model.geom('geom1').feature('wp3').geom.feature.create('r1', 'Rectangle');
model.geom('geom1').feature('wp3').geom.feature('r1').set('base', 'corner');
model.geom('geom1').feature('wp3').geom.feature('r1').set('pos', {'3' '2'});
model.geom('geom1').feature('wp3').geom.feature('r1').set('lx', '1');
model.geom('geom1').feature('wp3').geom.feature('r1').set('ly', '1');
model.geom('geom1').feature('wp3').geom.run('r1');
model.geom('geom1').run;

model.material.remove('mat1');
model.material.create('mat1');
model.material('mat1').name('Air');
model.material('mat1').propertyGroup('def').set('relpermeability', '1');
model.material('mat1').propertyGroup('def').set('relpermittivity', '1');
model.material('mat1').propertyGroup('def').set('dynamicviscosity', 'eta(T[1/K])[Pa*s]');
model.material('mat1').propertyGroup('def').set('ratioofspecificheat', '1.4');
model.material('mat1').propertyGroup('def').set('electricconductivity', '0[S/m]');
model.material('mat1').propertyGroup('def').set('heatcapacity', 'Cp(T[1/K])[J/(kg*K)]');
model.material('mat1').propertyGroup('def').set('density', 'rho(pA[1/Pa],T[1/K])[kg/m^3]');
model.material('mat1').propertyGroup('def').set('thermalconductivity', 'k(T[1/K])[W/(m*K)]');
model.material('mat1').propertyGroup('def').set('soundspeed', 'cs(T[1/K])[m/s]');
model.material('mat1').propertyGroup('def').func.create('eta', 'Piecewise');
model.material('mat1').propertyGroup('def').func('eta').set('funcname', 'eta');
model.material('mat1').propertyGroup('def').func('eta').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('eta').set('extrap', 'constant');
model.material('mat1').propertyGroup('def').func('eta').set('pieces', {'200.0' '1600.0' '-8.38278E-7+8.35717342E-8*T^1-7.69429583E-11*T^2+4.6437266E-14*T^3-1.06585607E-17*T^4'});
model.material('mat1').propertyGroup('def').func.create('Cp', 'Piecewise');
model.material('mat1').propertyGroup('def').func('Cp').set('funcname', 'Cp');
model.material('mat1').propertyGroup('def').func('Cp').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('Cp').set('extrap', 'constant');
model.material('mat1').propertyGroup('def').func('Cp').set('pieces', {'200.0' '1600.0' '1047.63657-0.372589265*T^1+9.45304214E-4*T^2-6.02409443E-7*T^3+1.2858961E-10*T^4'});
model.material('mat1').propertyGroup('def').func.create('rho', 'Analytic');
model.material('mat1').propertyGroup('def').func('rho').set('funcname', 'rho');
model.material('mat1').propertyGroup('def').func('rho').set('args', {'pA' 'T'});
model.material('mat1').propertyGroup('def').func('rho').set('expr', 'pA*0.02897/8.314/T');
model.material('mat1').propertyGroup('def').func('rho').set('dermethod', 'manual');
model.material('mat1').propertyGroup('def').func('rho').set('argders', {'pA' 'd(pA*0.02897/8.314/T,pA)'; 'T' 'd(pA*0.02897/8.314/T,T)'});
model.material('mat1').propertyGroup('def').func.create('k', 'Piecewise');
model.material('mat1').propertyGroup('def').func('k').set('funcname', 'k');
model.material('mat1').propertyGroup('def').func('k').set('arg', 'T');
model.material('mat1').propertyGroup('def').func('k').set('extrap', 'constant');
model.material('mat1').propertyGroup('def').func('k').set('pieces', {'200.0' '1600.0' '-0.00227583562+1.15480022E-4*T^1-7.90252856E-8*T^2+4.11702505E-11*T^3-7.43864331E-15*T^4'});
model.material('mat1').propertyGroup('def').func.create('cs', 'Analytic');
model.material('mat1').propertyGroup('def').func('cs').set('funcname', 'cs');
model.material('mat1').propertyGroup('def').func('cs').set('args', {'T'});
model.material('mat1').propertyGroup('def').func('cs').set('expr', 'sqrt(1.4*287*T)');
model.material('mat1').propertyGroup('def').func('cs').set('dermethod', 'manual');
model.material('mat1').propertyGroup('def').func('cs').set('argders', {'T' 'd(sqrt(1.4*287*T),T)'});
model.material('mat1').propertyGroup('def').addInput('temperature');
model.material('mat1').propertyGroup('def').addInput('pressure');
model.material('mat1').selection.geom('geom1', 2);
model.material('mat1').selection.set([3]);
model.material.create('mat2');
model.material('mat2').name('Hydrogen');
model.material('mat2').propertyGroup('def').set('dynamicviscosity', 'eta(T[1/K])[Pa*s]');
model.material('mat2').propertyGroup('def').set('ratioofspecificheat', '1.41');
model.material('mat2').propertyGroup('def').set('heatcapacity', 'Cp(T[1/K])[J/(kg*K)]');
model.material('mat2').propertyGroup('def').set('density', 'rho(pA[1/Pa],T[1/K])[kg/m^3]');
model.material('mat2').propertyGroup('def').set('thermalconductivity', 'k(T[1/K])[W/(m*K)]');
model.material('mat2').propertyGroup('def').func.create('eta', 'Piecewise');
model.material('mat2').propertyGroup('def').func('eta').set('funcname', 'eta');
model.material('mat2').propertyGroup('def').func('eta').set('arg', 'T');
model.material('mat2').propertyGroup('def').func('eta').set('extrap', 'constant');
model.material('mat2').propertyGroup('def').func('eta').set('pieces', {'200.0' '1300.0' '2.14524642E-6+2.54245E-8*T^1-1.0235587E-11*T^2+2.80895021E-15*T^3'});
model.material('mat2').propertyGroup('def').func.create('Cp', 'Piecewise');
model.material('mat2').propertyGroup('def').func('Cp').set('funcname', 'Cp');
model.material('mat2').propertyGroup('def').func('Cp').set('arg', 'T');
model.material('mat2').propertyGroup('def').func('Cp').set('extrap', 'constant');
model.material('mat2').propertyGroup('def').func('Cp').set('pieces', {'200.0' '1300.0' '10808.501+21.5799904*T^1-0.0444720318*T^2+3.85401176E-5*T^3-1.14979447E-8*T^4'});
model.material('mat2').propertyGroup('def').func.create('rho', 'Analytic');
model.material('mat2').propertyGroup('def').func('rho').set('funcname', 'rho');
model.material('mat2').propertyGroup('def').func('rho').set('args', {'pA' 'T'});
model.material('mat2').propertyGroup('def').func('rho').set('expr', 'pA*0.002016/8.314/T');
model.material('mat2').propertyGroup('def').func('rho').set('dermethod', 'manual');
model.material('mat2').propertyGroup('def').func('rho').set('argders', {'pA' 'd(pA*0.002016/8.314/T,pA)'; 'T' 'd(pA*0.002016/8.314/T,T)'});
model.material('mat2').propertyGroup('def').func.create('k', 'Piecewise');
model.material('mat2').propertyGroup('def').func('k').set('funcname', 'k');
model.material('mat2').propertyGroup('def').func('k').set('arg', 'T');
model.material('mat2').propertyGroup('def').func('k').set('extrap', 'constant');
model.material('mat2').propertyGroup('def').func('k').set('pieces', {'200.0' '1300.0' '0.00517975922+6.72778E-4*T^1-3.0388973E-7*T^2+6.58874687E-11*T^3'});
model.material('mat2').propertyGroup('def').addInput('temperature');
model.material('mat2').propertyGroup('def').addInput('pressure');
model.material('mat2').selection.geom('geom1', 2);
model.material('mat2').selection.set([7]);

model.physics('spf').feature.create('inl2', 'Inlet', 2);
model.physics('spf').feature('inl2').selection.set([7]);
model.physics('spf').feature('inl2').set('U0in', 1, '2');

model.mesh('mesh1').feature.remove('ftet3');
model.mesh('mesh1').feature.create('ftet1', 'FreeTet');
model.mesh('mesh1').run;

model.result('pg1').run;

model.material('mat1').selection.geom('geom1', 3);
model.material('mat1').selection.geom('geom1', 2);
model.material('mat1').selection.set([3]);

model.physics('spf').feature.remove('inl2');

model.result('pg1').run;
model.result('pg1').set('renderdatacached', true);
model.result.table.create('evl3', 'Table');
model.result.table('evl3').comments('Interactive 3D values');
model.result.table('evl3').name('Evaluation 3D');
model.result.table('evl3').addRow([3.56117056911689 2.1407101218288354 4.440892098500626E-16 0]);

model.physics('spf').feature('out1').selection.set([4]);

model.result('pg1').run;

model.material('mat1').selection.geom('geom1', 3);
model.material('mat1').selection.all;
model.material('mat1').selection.set([1]);
model.material('mat1').selection.all;

model.sol('sol1').runAll;

model.result('pg1').run;

model.physics('spf').feature.create('inl2', 'Inlet', 2);
model.physics('spf').feature('inl2').selection.set([7]);
model.physics('spf').feature('inl2').set('U0in', 1, '4');

model.sol('sol1').runAll;

model.result('pg1').run;

model.material('mat2').selection.geom('geom1', 3);
model.material('mat2').selection.all;
model.material('mat2').selection.set([1]);
model.material('mat2').selection.all;
model.material('mat2').selection.geom('geom1', 0);
model.material('mat2').selection.set([8]);
model.material('mat2').selection.geom('geom1', 3);
model.material('mat2').selection.all;

model.physics.create('spf2', 'LaminarFlow', 'geom1');

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

model.physics('spf2').feature.create('inl1', 'Inlet', 2);
model.physics('spf2').feature('inl1').selection.set([7]);
model.physics('spf2').feature('inl1').set('U0in', 1, '5');
model.physics('spf').feature.remove('inl2');

model.material('mat2').selection.set([1]);

model.result('pg1').run;

model.physics('spf2').feature('inl1').set('U0in', 1, '3');
model.physics('spf2').feature.create('out1', 'Outlet', 2);
model.physics('spf2').feature('out1').selection.set([4]);

model.sol('sol1').runAll;

model.result('pg1').run;
model.result('pg1').set('renderdatacached', true);
model.result.table('evl3').addRow([3.2907501462926287 1.4979659418667985 -4.440892098500626E-16 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', 'stat');
model.sol('sol2').feature.create('v1', 'Variables');
model.sol('sol2').feature.create('s1', 'Stationary');
model.sol('sol2').feature('s1').feature.create('se1', 'Segregated');
model.sol('sol2').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol2').feature('s1').feature('se1').feature.create('ss1', 'SegregatedStep');
model.sol('sol2').feature('s1').feature('se1').feature('ss1').set('segvar', {'mod1_u' 'mod1_p'});
model.sol('sol2').feature('s1').feature.create('i1', 'Iterative');
model.sol('sol2').feature('s1').feature('i1').set('linsolver', 'gmres');
model.sol('sol2').feature('s1').feature('i1').set('prefuntype', 'left');
model.sol('sol2').feature('s1').feature('i1').set('itrestart', 50);
model.sol('sol2').feature('s1').feature('i1').set('rhob', 20);
model.sol('sol2').feature('s1').feature('i1').set('maxlinit', 200);
model.sol('sol2').feature('s1').feature('i1').set('nlinnormuse', 'on');
model.sol('sol2').feature('s1').feature('se1').feature('ss1').set('linsolver', 'i1');
model.sol('sol2').feature('s1').feature('i1').feature.create('mg1', 'Multigrid');
model.sol('sol2').feature('s1').feature('i1').feature('mg1').set('prefun', 'gmg');
model.sol('sol2').feature('s1').feature('i1').feature('mg1').set('mcasegen', 'any');
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('pr').feature.create('sl1', 'SORLine');
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('iter', 2);
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('linerelax', 0.2);
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('linealgorithm', 'matrix');
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('linemethodmatrix', 'coupled');
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('linevar', 'mod1_u');
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('seconditer', 2);
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('relax', 0.3);
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('sl1', 'SORLine');
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('iter', 2);
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('linerelax', 0.2);
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('linealgorithm', 'matrix');
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('linemethodmatrix', 'coupled');
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('linevar', 'mod1_u');
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('seconditer', 2);
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('relax', 0.3);
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('cs').feature.create('d1', 'Direct');
model.sol('sol2').feature('s1').feature('i1').feature('mg1').feature('cs').feature('d1').set('linsolver', 'pardiso');
model.sol('sol2').feature('s1').feature('se1').feature.create('ss2', 'SegregatedStep');
model.sol('sol2').feature('s1').feature('se1').feature('ss2').set('segvar', {'mod1_u2' 'mod1_p2'});
model.sol('sol2').feature('s1').feature.create('i2', 'Iterative');
model.sol('sol2').feature('s1').feature('i2').set('linsolver', 'gmres');
model.sol('sol2').feature('s1').feature('i2').set('prefuntype', 'left');
model.sol('sol2').feature('s1').feature('i2').set('itrestart', 50);
model.sol('sol2').feature('s1').feature('i2').set('rhob', 20);
model.sol('sol2').feature('s1').feature('i2').set('maxlinit', 200);
model.sol('sol2').feature('s1').feature('i2').set('nlinnormuse', 'on');
model.sol('sol2').feature('s1').feature('se1').feature('ss2').set('linsolver', 'i2');
model.sol('sol2').feature('s1').feature('i2').feature.create('mg1', 'Multigrid');
model.sol('sol2').feature('s1').feature('i2').feature('mg1').set('prefun', 'gmg');
model.sol('sol2').feature('s1').feature('i2').feature('mg1').set('mcasegen', 'any');
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('pr').feature.create('sl1', 'SORLine');
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('iter', 2);
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('linerelax', 0.2);
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('linealgorithm', 'matrix');
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('linemethodmatrix', 'coupled');
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('linevar', 'mod1_u2');
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('seconditer', 2);
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('relax', 0.3);
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('po').feature.create('sl1', 'SORLine');
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('iter', 2);
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('linerelax', 0.2);
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('linealgorithm', 'matrix');
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('linemethodmatrix', 'coupled');
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('linevar', 'mod1_u2');
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('seconditer', 2);
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('relax', 0.3);
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('cs').feature.create('d1', 'Direct');
model.sol('sol2').feature('s1').feature('i2').feature('mg1').feature('cs').feature('d1').set('linsolver', 'pardiso');
model.sol('sol2').feature('s1').feature.remove('fcDef');
model.sol('sol2').attach('std2');

model.result.create('pg3', 3);
model.result('pg3').set('data', 'dset2');
model.result('pg3').feature.create('surf1', 'Surface');
model.result('pg3').feature('surf1').set('expr', 'spf.U');
model.result('pg3').feature('surf1').set('descr', 'Velocity magnitude');
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', 'spf.U');
model.result('pg4').feature('slc1').set('descr', 'Velocity magnitude');

model.sol('sol2').runAll;

model.result('pg3').run;

model.name('Hydrogen2.mph');

model.result('pg3').run;

model.physics.create('chcs', 'ConcentratedSpecies', 'geom1', {'w31' 'w32'});

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

model.physics('chcs').prop('DiffusionModel').set('DiffusionModel', 1, 'FicksLaw');
model.physics('chcs').feature('cdm1').set('u', {'u' 'v' 'w'});
model.physics('chcs').feature.create('in1', 'Inflow', 2);
model.physics('chcs').feature('in1').selection.set([3]);
model.physics('chcs').feature.create('in2', 'Inflow', 2);
model.physics('chcs').feature('in2').selection.set([7]);
model.physics('chcs').feature('in1').set('w0', 2, '0.5');
model.physics('chcs').feature('in2').set('w0', 2, '0.5');
model.physics('chcs').feature.create('out1', 'Outflow', 2);
model.physics('chcs').feature('out1').selection.set([4]);

model.sol.create('sol3');
model.sol('sol3').study('std3');
model.sol('sol3').feature.create('st1', 'StudyStep');
model.sol('sol3').feature('st1').set('study', 'std3');
model.sol('sol3').feature('st1').set('studystep', 'stat');
model.sol('sol3').feature.create('v1', 'Variables');
model.sol('sol3').feature.create('s1', 'Stationary');
model.sol('sol3').feature('s1').feature.create('se1', 'Segregated');
model.sol('sol3').feature('s1').feature('se1').feature.remove('ssDef');
model.sol('sol3').feature('s1').feature('se1').feature.create('ss1', 'SegregatedStep');
model.sol('sol3').feature('s1').feature('se1').feature('ss1').set('segvar', {'mod1_w32'});
model.sol('sol3').feature('s1').feature('se1').feature('ss1').set('subdamp', 0.7);
model.sol('sol3').feature('s1').feature.create('i1', 'Iterative');
model.sol('sol3').feature('s1').feature('i1').set('linsolver', 'gmres');
model.sol('sol3').feature('s1').feature('i1').set('prefuntype', 'left');
model.sol('sol3').feature('s1').feature('i1').set('rhob', 20);
model.sol('sol3').feature('s1').feature('i1').set('itrestart', 50);
model.sol('sol3').feature('s1').feature('se1').feature('ss1').set('linsolver', 'i1');
model.sol('sol3').feature('s1').feature('i1').feature.create('mg1', 'Multigrid');
model.sol('sol3').feature('s1').feature('i1').feature('mg1').set('prefun', 'gmg');
model.sol('sol3').feature('s1').feature('i1').feature('mg1').set('mcasegen', 'any');
model.sol('sol3').feature('s1').feature('i1').feature('mg1').feature('pr').feature.create('sl1', 'SORLine');
model.sol('sol3').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('iter', 2);
model.sol('sol3').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('linerelax', 0.4);
model.sol('sol3').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('seconditer', 2);
model.sol('sol3').feature('s1').feature('i1').feature('mg1').feature('pr').feature('sl1').set('relax', 0.7);
model.sol('sol3').feature('s1').feature('i1').feature('mg1').feature('po').feature.create('sl1', 'SORLine');
model.sol('sol3').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('iter', 2);
model.sol('sol3').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('linerelax', 0.4);
model.sol('sol3').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('seconditer', 2);
model.sol('sol3').feature('s1').feature('i1').feature('mg1').feature('po').feature('sl1').set('relax', 0.7);
model.sol('sol3').feature('s1').feature('i1').feature('mg1').feature('cs').feature.create('d1', 'Direct');
model.sol('sol3').feature('s1').feature('i1').feature('mg1').feature('cs').feature('d1').set('linsolver', 'pardiso');
model.sol('sol3').feature('s1').feature('se1').feature.create('ss2', 'SegregatedStep');
model.sol('sol3').feature('s1').feature('se1').feature('ss2').set('segvar', {'mod1_u' 'mod1_p'});
model.sol('sol3').feature('s1').feature.create('i2', 'Iterative');
model.sol('sol3').feature('s1').feature('i2').set('linsolver', 'gmres');
model.sol('sol3').feature('s1').feature('i2').set('prefuntype', 'left');
model.sol('sol3').feature('s1').feature('i2').set('itrestart', 50);
model.sol('sol3').feature('s1').feature('i2').set('rhob', 20);
model.sol('sol3').feature('s1').feature('i2').set('maxlinit', 200);
model.sol('sol3').feature('s1').feature('i2').set('nlinnormuse', 'on');
model.sol('sol3').feature('s1').feature('se1').feature('ss2').set('linsolver', 'i2');
model.sol('sol3').feature('s1').feature('i2').feature.create('mg1', 'Multigrid');
model.sol('sol3').feature('s1').feature('i2').feature('mg1').set('prefun', 'gmg');
model.sol('sol3').feature('s1').feature('i2').feature('mg1').set('mcasegen', 'any');
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('pr').feature.create('sl1', 'SORLine');
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('iter', 2);
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('linerelax', 0.3);
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('linealgorithm', 'matrix');
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('linemethodmatrix', 'coupled');
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('linevar', 'mod1_u');
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('maxline', 15);
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('seconditer', 2);
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('pr').feature('sl1').set('relax', 0.7);
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('po').feature.create('sl1', 'SORLine');
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('iter', 2);
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('linerelax', 0.3);
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('linealgorithm', 'matrix');
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('linemethodmatrix', 'coupled');
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('linevar', 'mod1_u');
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('maxline', 15);
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('seconditer', 2);
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('po').feature('sl1').set('relax', 0.7);
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('cs').feature.create('d1', 'Direct');
model.sol('sol3').feature('s1').feature('i2').feature('mg1').feature('cs').feature('d1').set('linsolver', 'pardiso');
model.sol('sol3').feature('s1').feature('se1').feature.create('ss3', 'SegregatedStep');
model.sol('sol3').feature('s1').feature('se1').feature('ss3').set('segvar', {'mod1_u2' 'mod1_p2'});
model.sol('sol3').feature('s1').feature.create('i3', 'Iterative');
model.sol('sol3').feature('s1').feature('i3').set('linsolver', 'gmres');
model.sol('sol3').feature('s1').feature('i3').set('prefuntype', 'left');
model.sol('sol3').feature('s1').feature('i3').set('itrestart', 50);
model.sol('sol3').feature('s1').feature('i3').set('rhob', 20);
model.sol('sol3').feature('s1').feature('i3').set('maxlinit', 200);
model.sol('sol3').feature('s1').feature('i3').set('nlinnormuse', 'on');
model.sol('sol3').feature('s1').feature('se1').feature('ss3').set('linsolver', 'i3');
model.sol('sol3').feature('s1').feature('i3').feature.create('mg1', 'Multigrid');
model.sol('sol3').feature('s1').feature('i3').feature('mg1').set('prefun', 'gmg');
model.sol('sol3').feature('s1').feature('i3').feature('mg1').set('mcasegen', 'any');
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('pr').feature.create('sl1', 'SORLine');
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('pr').feature('sl1').set('iter', 2);
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('pr').feature('sl1').set('linerelax', 0.3);
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('pr').feature('sl1').set('linealgorithm', 'matrix');
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('pr').feature('sl1').set('linemethodmatrix', 'coupled');
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('pr').feature('sl1').set('linevar', 'mod1_u2');
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('pr').feature('sl1').set('maxline', 15);
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('pr').feature('sl1').set('seconditer', 2);
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('pr').feature('sl1').set('relax', 0.7);
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('po').feature.create('sl1', 'SORLine');
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('po').feature('sl1').set('iter', 2);
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('po').feature('sl1').set('linerelax', 0.3);
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('po').feature('sl1').set('linealgorithm', 'matrix');
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('po').feature('sl1').set('linemethodmatrix', 'coupled');
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('po').feature('sl1').set('linevar', 'mod1_u2');
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('po').feature('sl1').set('maxline', 15);
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('po').feature('sl1').set('seconditer', 2);
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('po').feature('sl1').set('relax', 0.7);
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('cs').feature.create('d1', 'Direct');
model.sol('sol3').feature('s1').feature('i3').feature('mg1').feature('cs').feature('d1').set('linsolver', 'pardiso');
model.sol('sol3').feature('s1').feature.remove('fcDef');
model.sol('sol3').attach('std3');

model.result.create('pg5', 3);
model.result('pg5').set('data', 'dset3');
model.result('pg5').feature.create('slc1', 'Slice');
model.result('pg5').feature('slc1').set('expr', {'spf.U'});
model.result('pg5').set('frametype', 'spatial');
model.result('pg5').name('Velocity (spf)');
model.result.dataset.create('surf1', 'Surface');
model.result.dataset('surf1').set('data', 'dset3');
model.result.create('pg6', 3);
model.result('pg6').set('data', 'surf1');
model.result('pg6').set('frametype', 'spatial');
model.result('pg6').feature.create('surf1', 'Surface');
model.result('pg6').feature('surf1').set('expr', {'1'});
model.result('pg6').feature('surf1').set('coloring', 'uniform');
model.result('pg6').feature('surf1').set('color', 'gray');
model.result('pg6').feature.create('con1', 'Contour');
model.result('pg6').feature('con1').set('expr', {'p'});
model.result('pg6').feature('con1').set('number', 40);
model.result('pg6').name('Pressure (spf)');
model.result.create('pg7', 3);
model.result('pg7').set('data', 'dset3');
model.result('pg7').feature.create('slc1', 'Slice');
model.result('pg7').feature('slc1').set('expr', {'spf2.U'});
model.result('pg7').set('frametype', 'spatial');
model.result('pg7').name('Velocity (spf2)');
model.result.dataset.create('surf2', 'Surface');
model.result.dataset('surf2').set('data', 'dset3');
model.result.create('pg8', 3);
model.result('pg8').set('data', 'surf2');
model.result('pg8').set('frametype', 'spatial');
model.result('pg8').feature.create('surf1', 'Surface');
model.result('pg8').feature('surf1').set('expr', {'1'});
model.result('pg8').feature('surf1').set('coloring', 'uniform');
model.result('pg8').feature('surf1').set('color', 'gray');
model.result('pg8').feature.create('con1', 'Contour');
model.result('pg8').feature('con1').set('expr', {'p2'});
model.result('pg8').feature('con1').set('number', 40);
model.result('pg8').name('Pressure (spf2)');
model.result.create('pg9', 3);
model.result('pg9').set('data', 'dset3');
model.result('pg9').feature.create('slc1', 'Slice');
model.result('pg9').feature('slc1').set('expr', 'w31');
model.result('pg9').feature('slc1').set('descr', 'Mass fraction');
model.result('pg9').name('Mass Fraction (chcs)');
model.result.create('pg10', 3);
model.result('pg10').set('data', 'dset3');
model.result('pg10').feature.create('surf1', 'Surface');
model.result('pg10').feature('surf1').set('expr', 'w31');
model.result('pg10').feature('surf1').set('descr', 'Mass fraction');
model.result('pg10').name('Mass Fraction (chcs) 1');
model.result('pg5').run;
model.result('pg9').run;
model.result('pg10').run;

model.physics('spf').feature('inl1').set('U0in', 1, '0.1');
model.physics('spf2').feature('inl1').set('U0in', 1, '0.3');

model.sol('sol3').runAll;

model.result('pg5').run;
model.result('pg9').run;
model.result('pg10').run;

model.name('Hydrogen2.mph');

model.result('pg10').run;

out = model;