Changeset 23793 for issm/trunk-jpl/test/NightlyRun/test218.py

Timestamp:

03/13/19 03:17:46 (6 years ago)

Author:

bdef

Message:

pep8 compliance of NTs

File:

• issm/trunk-jpl/test/NightlyRun/test218.py (modified) (5 diffs)

issm/trunk-jpl/test/NightlyRun/test218.py

@@ -22 +22 @@
 from verbose import *
 
-md = squaremesh(model(),1000000,1000000,5,5)
-md = setmask(md,'all','')
-md = parameterize(md,'../Par/SquareShelf2.py')
-md = setflowequation(md,'SSA','all')
-md.cluster = generic('name',gethostname(),'np',3)
+md = squaremesh(model(), 1000000, 1000000, 5, 5)
+md = setmask(md, 'all', '')
+md = parameterize(md, '../Par/SquareShelf2.py')
+md = setflowequation(md, 'SSA', 'all')
+md.cluster = generic('name', gethostname(), 'np', 3)
 
-#redo the parameter file for this special shelf. 
-#constant thickness, constrained (vy=0) flow into an icefront, 
+#redo the parameter file for this special shelf.
+#constant thickness, constrained (vy = 0) flow into an icefront,
 #from 0 m/yr at the grounding line.
 
@@ -67 +67 @@
 
 #constrain grounding line to 0 velocity
-pos = np.where(md.mesh.y ==ymin)
+pos = np.where(md.mesh.y == ymin)
 md.stressbalance.spcvx[pos] = 0.
 md.stressbalance.spcvy[pos] = 0.
@@ -73 +73 @@
 #partitioning
 md.qmu.numberofpartitions = md.mesh.numberofvertices
-md = partitioner(md,'package','linear','npart',md.qmu.numberofpartitions)
-md.qmu.partition = (md.qmu.partition - 1.).reshape(-1,1).T
+md = partitioner(md, 'package', 'linear', 'npart', md.qmu.numberofpartitions)
+md.qmu.partition = (md.qmu.partition - 1.).reshape(-1, 1).T
 
 #Dakota options
@@ -81 +81 @@
 version = IssmConfig('_DAKOTA_VERSION_')
 # returns tuple "(u'6.2',)" -> unicode string '6.2', convert to float
-version=float(version[0])
+version = float(version[0])
 
 #variables
-md.qmu.variables.rheology_B = normal_uncertain.normal_uncertain('scaled_MaterialsRheologyB',1,.05)
+md.qmu.variables.rheology_B = normal_uncertain.normal_uncertain('scaled_MaterialsRheologyB', 1, .05)
 
 #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
@@ -97 +97 @@
 
 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**-10 #tighten for qmu analysis
+#imperative!
+md.stressbalance.reltol = 10**-10  #tighten for qmu analysis
 md.qmu.isdakota = 1
 
 #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_MaterialsRheologyB','MaxVel').reshape(-1,1)
+md.results.dakota.importancefactors = importancefactors(md, 'scaled_MaterialsRheologyB', 'MaxVel').reshape(-1, 1)
 field_names = ['importancefactors']
 field_tolerances = [1e-10]
-field_values=[md.results.dakota.importancefactors]
+field_values = [md.results.dakota.importancefactors]