from model import * from triangle import * from setmask import * from parameterize import * from setflowequation import * from EnumDefinitions import * from solve import * from MatlabFuncs import * md=triangle(model(),'../Exp/Square.exp',180000) md=setmask(md,'all','') md=parameterize(md,'../Par/SquareShelfConstrained.py') md.extrude(3,1) md=setflowequation(md,'macayeal','all') md.cluster=generic('name',oshostname(),'np',3); md.initialization.waterfraction=numpy.zeros(md.mesh.numberofvertices) md.transient.isdiagnostic=0 md.transient.isprognostic=0 md.transient.isthermal=1 md.transient.isgroundingline=0 md.thermal.isenthalpy=1 md=solve(md,TransientSolutionEnum()) #Fields and tolerances to track changes field_names =['Enthalpy1','Waterfraction1','Temperature1',\ 'Enthalpy2','Waterfraction2','Temperature2',\ 'Enthalpy3','Waterfraction3','Temperature3'] field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-10,1e-13] field_values=[\ md.results['TransientSolution'][1]['Enthalpy'],\ md.results['TransientSolution'][1]['Waterfraction'],\ md.results['TransientSolution'][1]['Temperature'],\ md.results['TransientSolution'][2]['Enthalpy'],\ md.results['TransientSolution'][2]['Waterfraction'],\ md.results['TransientSolution'][2]['Temperature'],\ md.results['TransientSolution'][3]['Enthalpy'],\ md.results['TransientSolution'][3]['Waterfraction'],\ md.results['TransientSolution'][3]['Temperature'],\ ]