source: issm/oecreview/Archive/13393-13976/ISSM-13517-13518.diff

Last change on this file was 13980, checked in by Mathieu Morlighem, 12 years ago

preparing oecreview for 13393-13976'
File size: 5.4 KB
RevLine 
[13980]1Index: ../trunk-jpl/test/NightlyRun/test212.py
2===================================================================
3--- ../trunk-jpl/test/NightlyRun/test212.py (revision 0)
4+++ ../trunk-jpl/test/NightlyRun/test212.py (revision 13518)
5@@ -0,0 +1,61 @@
6+"""
7+== == == == == == == == == == == == == == == == == == ==
8+Auto generated python script for ISSM: test212.m
9+Created on 2012-09-27 via translateToPy.py Ver 1.0 by mikep
10+== == == == == == == == == == == == == == == == == == ==
11+
12+Matlab script conversion into python
13+translateToPy.py Author: Michael Pellegrin
14+translateToPy.py Date: 09/24/12
15+== == == == == == == == == == == == == == == == == == ==
16+"""
17+
18+from MatlabFuncs import *
19+from model import *
20+from EnumDefinitions import *
21+from numpy import *
22+from triangle import *
23+from setmask import *
24+from parameterize import *
25+from setflowequation import *
26+from solve import *
27+
28+md=triangle(model(),'../Exp/Square.exp',200000)
29+md=setmask(md,'all','')
30+md=parameterize(md,'../Par/SquareShelf.py')
31+md=setflowequation(md,'macayeal','all')
32+
33+
34+# control parameters
35+
36+md.inversion.iscontrol=1
37+md.inversion.control_parameters=['MaterialsRheologyBbar']
38+md.inversion.min_parameters=10**6*ones((md.mesh.numberofvertices,1))
39+md.inversion.max_parameters=2*10**9*ones((md.mesh.numberofvertices,1))
40+md.inversion.nsteps=2
41+md.inversion.cost_functions=101*ones((md.inversion.nsteps,1))
42+md.inversion.cost_functions_coefficients=ones((md.mesh.numberofvertices,1))
43+md.inversion.gradient_scaling=10**7*ones((md.inversion.nsteps,1))
44+md.inversion.maxiter_per_step=2*ones((md.inversion.nsteps,1))
45+md.inversion.step_threshold=0.3*ones((md.inversion.nsteps,1))
46+md.inversion.vx_obs=md.initialization.vx
47+md.inversion.vy_obs=md.initialization.vy
48+
49+
50+md.cluster=generic('name',oshostname(),'np',3)
51+md=solve(md,DiagnosticSolutionEnum())
52+
53+
54+# Fields and tolerances to track changes
55+
56+field_names =['Gradient','Misfits','MaterialsRheologyBbar','Pressure','Vel','Vx','Vy']
57+field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
58+field_values=[\
59+md.results['DiagnosticSolution'][1]['Gradient1'],\
60+md.results['DiagnosticSolution'][1]['J'],\
61+md.results['DiagnosticSolution'][1]['MaterialsRheologyBbar'],\
62+md.results['DiagnosticSolution'][1]['Pressure'],\
63+md.results['DiagnosticSolution'][1]['Vel'],\
64+md.results['DiagnosticSolution'][1]['Vx'],\
65+md.results['DiagnosticSolution'][1]['Vy'],\
66+]
67Index: ../trunk-jpl/test/NightlyRun/test212.m
68===================================================================
69--- ../trunk-jpl/test/NightlyRun/test212.m (revision 13517)
70+++ ../trunk-jpl/test/NightlyRun/test212.m (revision 13518)
71@@ -24,7 +24,7 @@
72 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
73 field_values={...
74 (md.results.DiagnosticSolution.Gradient1),...
75- md.results.DiagnosticSolution.J,...
76+ (md.results.DiagnosticSolution.J),...
77 (md.results.DiagnosticSolution.MaterialsRheologyBbar),...
78 (md.results.DiagnosticSolution.Pressure),...
79 (md.results.DiagnosticSolution.Vel),...
80Index: ../trunk-jpl/test/NightlyRun/InNeedOfDebugging/test212.py
81===================================================================
82--- ../trunk-jpl/test/NightlyRun/InNeedOfDebugging/test212.py (revision 13517)
83+++ ../trunk-jpl/test/NightlyRun/InNeedOfDebugging/test212.py (revision 13518)
84@@ -1,61 +0,0 @@
85-"""
86-== == == == == == == == == == == == == == == == == == ==
87-Auto generated python script for ISSM: test212.m
88-Created on 2012-09-27 via translateToPy.py Ver 1.0 by mikep
89-== == == == == == == == == == == == == == == == == == ==
90-
91-Matlab script conversion into python
92-translateToPy.py Author: Michael Pellegrin
93-translateToPy.py Date: 09/24/12
94-== == == == == == == == == == == == == == == == == == ==
95-"""
96-
97-from MatlabFuncs import *
98-from model import *
99-from EnumDefinitions import *
100-from numpy import *
101-from triangle import *
102-from setmask import *
103-from parameterize import *
104-from setflowequation import *
105-from solve import *
106-
107-md=triangle(model(),'../Exp/Square.exp',200000)
108-md=setmask(md,'all','')
109-md=parameterize(md,'../Par/SquareShelf.py')
110-md=setflowequation(md,'macayeal','all')
111-
112-
113-# control parameters
114-
115-md.inversion.iscontrol=1
116-md.inversion.control_parameters=['MaterialsRheologyBbar']
117-md.inversion.min_parameters=10^6*ones(md.mesh.numberofvertices,1)
118-md.inversion.max_parameters=2*10^9*ones(md.mesh.numberofvertices,1)
119-md.inversion.nsteps=2
120-md.inversion.cost_functions=101*ones(md.inversion.nsteps,1)
121-md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,1)
122-md.inversion.gradient_scaling=10^7*ones(md.inversion.nsteps,1)
123-md.inversion.maxiter_per_step=2*ones(md.inversion.nsteps,1)
124-md.inversion.step_threshold=0.3*ones(md.inversion.nsteps,1)
125-md.inversion.vx_obs=md.initialization.vx
126-md.inversion.vy_obs=md.initialization.vy
127-
128-
129-md.cluster=generic('name',oshostname(),'np',3)
130-md=solve(md,DiagnosticSolutionEnum())
131-
132-
133-# Fields and tolerances to track changes
134-
135-field_names =['Gradient' 'Misfits' 'MaterialsRheologyBbar' 'Pressure' 'Vel' 'Vx' 'Vy']
136-field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
137-field_values=[\
138-md.results['DiagnosticSolution'][1]['Gradient1'],\
139-md.results['DiagnosticSolution.J'],\
140-md.results['DiagnosticSolution'][1]['MaterialsRheologyBbar'],\
141-md.results['DiagnosticSolution'][1]['Pressure'],\
142-md.results['DiagnosticSolution'][1]['Vel'],\
143-md.results['DiagnosticSolution'][1]['Vx'],\
144-md.results['DiagnosticSolution'][1]['Vy'],\
145-]
Note: See TracBrowser for help on using the repository browser.