md=triangle(model,'../Exp/Square.exp',150000); md=setmask(md,'../Exp/SquareShelf.exp',''); md=parameterize(md,'../Par/SquareSheetShelf.par'); md=extrude(md,3,2); md=setflowequation(md,'pattyn','all'); md.cluster=generic('name',oshostname(),'np',3); md.timestepping.time_step=0; md.thermal.isenthalpy=1; md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1); md=solve(md,SteadystateSolutionEnum); %Fields and tolerances to track changes field_names ={'Vx','Vy','Vz','Vel','Pressure','Temperature','Water fraction','Enthalpy'}; field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-13,1e-10,1e-10,1e-10}; field_values={... (md.results.SteadystateSolution.Vx),... (md.results.SteadystateSolution.Vy),... (md.results.SteadystateSolution.Vz),... (md.results.SteadystateSolution.Vel),... (md.results.SteadystateSolution.Pressure),... (md.results.SteadystateSolution.Temperature),... (md.results.SteadystateSolution.Waterfraction),... (md.results.SteadystateSolution.Enthalpy),... };