| 1 | from model import * | 
|---|
| 2 | from triangle import * | 
|---|
| 3 | from setmask import * | 
|---|
| 4 | from parameterize import * | 
|---|
| 5 | from setflowequation import * | 
|---|
| 6 | from EnumDefinitions import * | 
|---|
| 7 | from solve import * | 
|---|
| 8 | from MatlabFuncs import * | 
|---|
| 9 |  | 
|---|
| 10 | md=triangle(model(),'../Exp/Square.exp',180000.) | 
|---|
| 11 | md=setmask(md,'all','') | 
|---|
| 12 | md=parameterize(md,'../Par/SquareShelfConstrained.py') | 
|---|
| 13 | md.extrude(3,2.) | 
|---|
| 14 | md=setflowequation(md,'SSA','all') | 
|---|
| 15 | md.cluster=generic('name',oshostname(),'np',1) | 
|---|
| 16 | md.autodiff.isautodiff=True | 
|---|
| 17 | md.toolkits.DefaultAnalysis=issmgslsolver() | 
|---|
| 18 | md=solve(md,StressbalanceSolutionEnum()) | 
|---|
| 19 |  | 
|---|
| 20 | #Fields and tolerances to track changes | 
|---|
| 21 | field_names     =['Vx','Vy','Vz','Vel','Pressure'] | 
|---|
| 22 | field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13] | 
|---|
| 23 | field_values=[\ | 
|---|
| 24 | md.results.StressbalanceSolution.Vx,\ | 
|---|
| 25 | md.results.StressbalanceSolution.Vy,\ | 
|---|
| 26 | md.results.StressbalanceSolution.Vz,\ | 
|---|
| 27 | md.results.StressbalanceSolution.Vel,\ | 
|---|
| 28 | md.results.StressbalanceSolution.Pressure,\ | 
|---|
| 29 | ] | 
|---|