Changeset 27417 for issm/trunk-jpl/test/NightlyRun/test355.py

Timestamp:

11/24/22 00:49:18 (2 years ago)

Author:

jdquinn

Message:

CHG: Python translation of debris class; standardized tolerances for tests; varied cleanup

File:

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

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

@@ -10 +10 @@
 from generic import generic
 
-#create model:
+# Create model
 md = triangle(model(), '../Exp/Square.exp', 50000.)
 md.mesh.x = md.mesh.x / 100
 md.mesh.y = md.mesh.y / 100
-#miscellaneous
-md = setmask(md, '', '')   #everywhere grounded
+md.miscellaneous.name = 'testChannels'
+
+# Miscellaneous
+md = setmask(md, '', '') # Everywhere grounded
 md = setflowequation(md, 'SSA', 'all')
-md.miscellaneous.name = 'testChannels'
-md.stressbalance.maxiter = 2  #Make sure it runs quickly...
+md.stressbalance.maxiter = 2 # Make sure it runs quickly...
 
-#Some constants
+# Some constants
 md.constants.g = 9.8
 md.materials.rho_ice = 910
 
-#Geometry
+# Geometry
 md.geometry.surface = -0.02 * md.mesh.x + 320
 md.geometry.bed = np.zeros((md.mesh.numberofvertices))
@@ -30 +31 @@
 md.geometry.thickness = md.geometry.surface - md.geometry.bed
 
-#Define initial conditions
+# Define initial conditions
 md.initialization.vx = 1.0e-6 * md.constants.yts * np.ones((md.mesh.numberofvertices))
 md.initialization.vy = np.zeros((md.mesh.numberofvertices))
@@ -37 +38 @@
 md.initialization.hydraulic_potential = md.materials.rho_ice * md.constants.g * md.geometry.thickness
 
-#Materials
+#cMaterials
 md.materials.rheology_B = (5e-25)**(-1./3.) * np.ones((md.mesh.numberofvertices))
 md.materials.rheology_n = 3. * np.ones((md.mesh.numberofelements))
 
-#Friction
+#cFriction
 md.friction.coefficient = np.zeros((md.mesh.numberofvertices))
 md.friction.p = np.ones((md.mesh.numberofelements))
@@ -47 +48 @@
 #md.friction.coupling = 0
 
-#Boundary conditions:
+#Bcoundary conditions:
 md = SetIceSheetBC(md)
 
@@ -74 +75 @@
 
 md.cluster = generic('np', 2)
-md = solve(md, 'Transient')   #or 'tr'
+md = solve(md, 'Transient') # Or 'tr'
 
-#Fields and tolerances to track changes
+# Fields and tolerances to track changes
 field_names = ['HydrologySheetThickness1', 'HydraulicPotential1', 'ChannelArea1',
                'HydrologySheetThickness2', 'HydraulicPotential2', 'ChannelArea2',