Changeset 25836 for issm/trunk/test/NightlyRun/test445.py

Timestamp:

12/08/20 08:45:53 (4 years ago)

Author:

Mathieu Morlighem

Message:

merged trunk-jpl and trunk for revision 25834

Location:

issm/trunk

Files:

• . (modified) (1 prop)
• test (modified) (1 prop)
• test/NightlyRun (modified) (1 prop)
• test/NightlyRun/test445.py (modified) (5 diffs)

issm/trunk/test/NightlyRun

@@ -18 +18 @@
 run.old
 run_matlab
+test218.qmu.in
+test218.qmu.out
+test218.qmu.err

@@ -18 +18 @@
 run.old
 run_matlab
+test218.qmu.in
+test218.qmu.out
+test218.qmu.err

issm/trunk/test/NightlyRun/test445.py

@@ -1 +1 @@
 #Test Name: SquareSheetShelfSteaEnthalpyHO3dDakotaSampNeff
+
+# TODO:
+# - Figure out why test fails intermittently on Mac with "IndexError: list 
+# index out of range"
+#
+
+from os import getcwd
+from socket import gethostname
+
 import numpy as np
-from os import getcwd
+
+from ContourToMesh import *
+from dmeth_params_set import *
 from model import *
-from socket import gethostname
+from parameterize import *
+from partitioner import *
+from setflowequation import *
+from setmask import *
+from solve import *
 from triangle import *
-from setmask import *
-from parameterize import *
-from setflowequation import *
-from solve import *
-from partitioner import *
-from dmeth_params_set import *
-from ContourToMesh import *
+
 
 #model not consistent:  equality thickness = surface-base violated
@@ -24 +33 @@
 md.thermal.isenthalpy = 1
 md.thermal.isdynamicbasalspc = 1
-md.initialization.waterfraction = np.zeros((md.mesh.numberofvertices))
-md.initialization.watercolumn = np.zeros((md.mesh.numberofvertices))
+md.initialization.waterfraction = np.zeros((md.mesh.numberofvertices, 1))
+md.initialization.watercolumn = np.zeros((md.mesh.numberofvertices, 1))
 
 md.friction.coupling = 3
@@ -34 +43 @@
 version = float(version[0])
 
+#partitioning
+npart = 10
+partition, md = partitioner(md, 'package', 'chaco', 'npart', npart, 'weighting', 'on')
+partition -= 1
+md.qmu.isdakota = 1
+
 #variables
-md.qmu.variables.neff = normal_uncertain.normal_uncertain('scaled_FrictionEffectivePressure', 1, .05)
-md.qmu.variables.geoflux = normal_uncertain.normal_uncertain('scaled_BasalforcingsGeothermalflux', 1, .05)
+md.qmu.variables.neff = normal_uncertain.normal_uncertain(
+    'descriptor', 'scaled_FrictionEffectivePressure',
+    'mean', np.ones((npart, 1)),
+    'stddev', .05 * np.ones((npart, 1)),
+    'partition', partition
+    )
+md.qmu.variables.geoflux = normal_uncertain.normal_uncertain(
+    'descriptor', 'scaled_BasalforcingsGeothermalflux',
+    'mean', np.ones((npart, 1)),
+    'stddev', .05 * np.ones((npart, 1)),
+    'partition', partition
+    )
 
 #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.MassFlux1 = response_function.response_function('indexed_MassFlux_1', [], [0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999])
-md.qmu.responses.MassFlux2 = response_function.response_function('indexed_MassFlux_2', [], [0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999])
-md.qmu.responses.MassFlux3 = response_function.response_function('indexed_MassFlux_3', [], [0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999])
-md.qmu.responses.MassFlux4 = response_function.response_function('indexed_MassFlux_4', [], [0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999])
-md.qmu.responses.MassFlux5 = response_function.response_function('indexed_MassFlux_5', [], [0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999])
-md.qmu.responses.MassFlux6 = response_function.response_function('indexed_MassFlux_6', [], [0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999])
-md.qmu.responses.MassFlux7 = response_function.response_function('indexed_MassFlux_7', [], [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('descriptor','MaxVel')
+md.qmu.responses.MassFlux1 = response_function.response_function('descriptor','indexed_MassFlux_1')
+md.qmu.responses.MassFlux2 = response_function.response_function('descriptor','indexed_MassFlux_2')
+md.qmu.responses.MassFlux3 = response_function.response_function('descriptor','indexed_MassFlux_3')
+md.qmu.responses.MassFlux4 = response_function.response_function('descriptor','indexed_MassFlux_4')
+md.qmu.responses.MassFlux5 = response_function.response_function('descriptor','indexed_MassFlux_5')
+md.qmu.responses.massFlux6 = response_function.response_function('descriptor','indexed_MassFlux_6')
+md.qmu.responses.massFlux7 = response_function.response_function('descriptor','indexed_MassFlux_7')
 
 #mass flux profiles
@@ -54 +79 @@
 #method
 md.qmu.method = dakota_method.dakota_method('nond_samp')
-md.qmu.method = dmeth_params_set(md.qmu.method, 'seed', 1234, 'samples', 20, 'sample_type', 'random')
+md.qmu.method = dmeth_params_set(
+    md.qmu.method,
+    'seed', 1234,
+    'samples', 20,
+    'sample_type', 'random'
+    )
 
 #parameters
@@ -69 +99 @@
     md.qmu.params.evaluation_concurrency = 1
 
-#partitioning
-md.qmu.numberofpartitions = 10
-md = partitioner(md, 'package', 'chaco', 'npart', md.qmu.numberofpartitions, 'weighting', 'on')
-md.qmu.vpartition = md.qmu.vpartition - 1
-md.qmu.isdakota = 1
-
 md.stressbalance.reltol = 10**-5  #tighten for qmu analyses
-
 
 #solve

@@ -18 +18 @@
 run.old
 run_matlab
+test218.qmu.in
+test218.qmu.out
+test218.qmu.err