#Test Name: SquareSheetConstrainedTherTran 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,'','') md=parameterize(md,'../Par/SquareSheetConstrained.py') md.extrude(3,1.) md=setflowequation(md,'SSA','all') md.cluster=generic('name',oshostname(),'np',3) md.verbose=verbose('convergence',True,'solution',True) md.transient.isstressbalance=False md.transient.ismasstransport=False md.transient.issmb=True md.transient.isthermal=True md.transient.isgroundingline=False md=solve(md,'Transient') #Fields and tolerances to track changes field_names =['Temperature','BasalforcingsGroundediceMeltingRate'] field_tolerances=[1e-13,1e-13] field_values=[\ md.results.TransientSolution[0].Temperature,\ md.results.TransientSolution[0].BasalforcingsGroundediceMeltingRate,\ ]