source: issm/trunk-jpl/test/NightlyRun/test213.py@ 21049

Last change on this file since 21049 was 21049, checked in by agscott1, 9 years ago

CHG: Replaced Enums with Strings in matlab and python. Updated corresponding cpp code.

File size: 1.7 KB
Line 
1#Test Name: SquareShelfCMBSSA3d
2from MatlabFuncs import *
3from model import *
4from EnumDefinitions import *
5from numpy import *
6from triangle import *
7from setmask import *
8from parameterize import *
9from setflowequation import *
10from solve import *
11from generic import generic
12
13md=triangle(model(),'../Exp/Square.exp',200000)
14md=setmask(md,'all','')
15md=parameterize(md,'../Par/SquareShelf.py')
16md.extrude(3,1.)
17md=setflowequation(md,'SSA','all')
18
19
20# control parameters
21
22md.inversion.iscontrol=1
23md.inversion.control_parameters=['MaterialsRheologyBbar']
24md.inversion.min_parameters=10**6*ones((md.mesh.numberofvertices,1))
25md.inversion.max_parameters=2.*10**9*ones((md.mesh.numberofvertices,1))
26md.inversion.nsteps=2
27md.inversion.cost_functions=[101]
28md.inversion.cost_functions_coefficients=ones((md.mesh.numberofvertices,1))
29md.inversion.gradient_scaling=10**7*ones((md.inversion.nsteps,1))
30md.inversion.maxiter_per_step=2.*ones((md.inversion.nsteps,1))
31md.inversion.step_threshold=0.3*ones((md.inversion.nsteps,1))
32md.inversion.vx_obs=md.initialization.vx
33md.inversion.vy_obs=md.initialization.vy
34
35
36md.cluster=generic('name',oshostname(),'np',3)
37md=solve(md,'StressbalanceSolution')
38
39
40# Fields and tolerances to track changes
41
42field_names =['Gradient','Misfits','MaterialsRheologyBbar','Pressure','Vel','Vx','Vy']
43field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
44#field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
45field_values=[\
46 md.results.StressbalanceSolution.Gradient1,\
47 md.results.StressbalanceSolution.J,\
48 md.results.StressbalanceSolution.MaterialsRheologyBbar,\
49 md.results.StressbalanceSolution.Pressure,\
50 md.results.StressbalanceSolution.Vel,\
51 md.results.StressbalanceSolution.Vx,\
52 md.results.StressbalanceSolution.Vy,\
53]
Note: See TracBrowser for help on using the repository browser.