Changeset 24313 for issm/trunk/test/NightlyRun/test412.py

Timestamp:

11/01/19 12:01:57 (5 years ago)

Author:

Mathieu Morlighem

Message:

merged trunk-jpl and trunk for revision 24310

Location:

issm/trunk

Files:

• . (modified) (1 prop)
• test (modified) (1 prop)
• test/NightlyRun/test412.py (modified) (4 diffs)

issm/trunk/test/NightlyRun/test412.py

@@ -1 +1 @@
 #Test Name: SquareSheetShelfDiadSSA3dDakota
-import numpy as np
 from model import *
 from socket import gethostname
@@ -10 +9 @@
 from partitioner import *
 from importancefactors import *
-
 from normal_uncertain import *
 from response_function import *
 
-md = triangle(model(),'../Exp/Square.exp',300000.)
-md = setmask(md,'../Exp/SquareShelf.exp','')
-md = parameterize(md,'../Par/SquareSheetShelf.py')
-md = setflowequation(md,'SSA','all')
-md.cluster = generic('name',gethostname(),'np',3)
+md = triangle(model(), '../Exp/Square.exp', 300000.)
+md = setmask(md, '../Exp/SquareShelf.exp', '')
+md = parameterize(md, '../Par/SquareSheetShelf.py')
+md = setflowequation(md, 'SSA', 'all')
+md.cluster = generic('name', gethostname(), 'np', 3)
 
 #partitioning
 md.qmu.numberofpartitions = md.mesh.numberofvertices
-md = partitioner(md,'package','linear','npart',md.qmu.numberofpartitions)
-md.qmu.partition = md.qmu.partition-1
+md = partitioner(md, 'package', 'linear', 'npart', md.qmu.numberofpartitions)
+md.qmu.vpartition = md.qmu.vpartition - 1
 md.qmu.isdakota = 1
 
@@ -33 +31 @@
 
 #variables
-md.qmu.variables.rho_ice = normal_uncertain.normal_uncertain('MaterialsRhoIce',md.materials.rho_ice,0.01)
-md.qmu.variables.drag_coefficient = normal_uncertain.normal_uncertain('scaled_FrictionCoefficient',1,0.01)
+md.qmu.variables.rho_ice = normal_uncertain.normal_uncertain('MaterialsRhoIce', md.materials.rho_ice, 0.01)
+md.qmu.variables.drag_coefficient = normal_uncertain.normal_uncertain('scaled_FrictionCoefficient', 1, 0.01)
 
 #responses
-md.qmu.responses.MaxVel = response_function.response_function('MaxVel',[],[0.0001,0.001,0.01,0.25,0.5,0.75,0.99,0.999,0.9999])
+md.qmu.responses.MaxVel = response_function.response_function('MaxVel', [], [0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999])
 
 #method
@@ -47 +45 @@
 
 if version >= 6:
-        md.qmu.params.analysis_driver = 'matlab'
-        md.qmu.params.evaluation_scheduling = 'master'
-        md.qmu.params.processors_per_evaluation = 2
+    md.qmu.params.analysis_driver = 'matlab'
+    md.qmu.params.evaluation_scheduling = 'master'
+    md.qmu.params.processors_per_evaluation = 2
 else:
-        md.qmu.params.analysis_driver = 'stressbalance'
-        md.qmu.params.evaluation_concurrency = 1
+    md.qmu.params.analysis_driver = 'stressbalance'
+    md.qmu.params.evaluation_concurrency = 1
 
-#imperative! 
-md.stressbalance.reltol = 10**-5 #tighten for qmu analyses
+#imperative!
+md.stressbalance.reltol = 10**-5  #tighten for qmu analyses
 
 #solve
-md.verbose = verbose('000000000')       # this line is recommended
-md = solve(md,'Stressbalance','overwrite','y')
+md.verbose = verbose('000000000')  # this line is recommended
+md = solve(md, 'Stressbalance', 'overwrite', 'y')
 
 #Fields and tolerances to track changes
 md.qmu.results = md.results.dakota
-md.results.dakota.importancefactors = importancefactors(md,'scaled_FrictionCoefficient','MaxVel').T
+md.results.dakota.importancefactors = importancefactors(md, 'scaled_FrictionCoefficient', 'MaxVel').T
 field_names = ['importancefactors']
 field_tolerances = [1e-10]