#Test Name: SquareSheetConstrainedTranCflSIA3d import numpy 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',150000.) md=setmask(md,'','') md=parameterize(md,'../Par/SquareSheetConstrained.py') md.extrude(5,1.2) md=setflowequation(md,'SIA','all') md.cluster=generic('name',oshostname(),'np',3) md.timestepping.time_adapt=1 md.timestepping.final_time=16000. md=solve(md,'Transient') #Fields and tolerances to track changes field_names =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1',\ 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2',\ 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsGroundediceMeltingRate3'] field_tolerances=[1e-13, 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\ 5*1e-11,1e-13,5e-13,1e-13,1e-13,1e-13,5e-13,1e-13,1e-13,1e-13,\ 1e-10, 1e-12,1e-12,1e-12,1e-13,1e-13,1e-12,1e-12,1e-13,1e-13] field_values=[\ md.results.TransientSolution[0].Vx,\ md.results.TransientSolution[0].Vy,\ md.results.TransientSolution[0].Vz,\ md.results.TransientSolution[0].Vel,\ md.results.TransientSolution[0].Pressure,\ md.results.TransientSolution[0].Base,\ md.results.TransientSolution[0].Surface,\ md.results.TransientSolution[0].Thickness,\ md.results.TransientSolution[0].Temperature,\ md.results.TransientSolution[0].BasalforcingsGroundediceMeltingRate,\ md.results.TransientSolution[1].Vx,\ md.results.TransientSolution[1].Vy,\ md.results.TransientSolution[1].Vz,\ md.results.TransientSolution[1].Vel,\ md.results.TransientSolution[1].Pressure,\ md.results.TransientSolution[1].Base,\ md.results.TransientSolution[1].Surface,\ md.results.TransientSolution[1].Thickness,\ md.results.TransientSolution[1].Temperature,\ md.results.TransientSolution[1].BasalforcingsGroundediceMeltingRate,\ md.results.TransientSolution[2].Vx,\ md.results.TransientSolution[2].Vy,\ md.results.TransientSolution[2].Vz,\ md.results.TransientSolution[2].Vel,\ md.results.TransientSolution[2].Pressure,\ md.results.TransientSolution[2].Base,\ md.results.TransientSolution[2].Surface,\ md.results.TransientSolution[2].Thickness,\ md.results.TransientSolution[2].Temperature,\ md.results.TransientSolution[2].BasalforcingsGroundediceMeltingRate,\ ]