#Test Name: SquareSheetShelfSteaSSA3d import numpy from model import * from EnumDefinitions import * from MatlabFuncs import * from triangle import * from setmask import * from parameterize import * from setflowequation import * from solve import * md=triangle(model(),'../Exp/Square.exp',150000.) md=setmask(md,'../Exp/SquareShelf.exp','') md=parameterize(md,'../Par/SquareSheetShelf.py') md.extrude(3,2.) md=setflowequation(md,'SSA','all') md.cluster=generic('name',oshostname(),'np',3) md.timestepping.time_step=0 md=solve(md,'SteadystateSolution') #Fields and tolerances to track changes field_names =['Vx','Vy','Vz','Vel','Pressure','Temperature','BasalforcingsGroundediceMeltingRate'] field_tolerances=[1e-09,1e-09,1e-08,1e-09,1e-13,5e-10,1e-06] 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.BasalforcingsGroundediceMeltingRate,\ ]