Changeset 27167 for issm/trunk-jpl/test/NightlyRun/test2004.py

Timestamp:

08/01/22 06:36:56 (3 years ago)

Author:

bdef

Message:

CHG: minor formating

File:

• issm/trunk-jpl/test/NightlyRun/test2004.py (modified) (12 diffs)

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

@@ -2 +2 @@
 #
 # TODO:
-# - Save boundaries by name to some data structure with tag so that they can be 
-#   copied to basins that use identical shapefiles and projections (will need 
-#   to check if the cost of the additional structure, checks, and copying are 
+# - Save boundaries by name to some data structure with tag so that they can be
+#   copied to basins that use identical shapefiles and projections (will need
+#   to check if the cost of the additional structure, checks, and copying are
 #   greater than the cost of projecting).
 #
@@ -41 +41 @@
 #HemisphereWest #{{{
 sl.addbasin(
-    basin('continent', 'hemispherewest', 'name', 'hemispherewest', 'proj', laea(0, -90), 'boundaries', [ # Peru projection 3587
+    basin('continent', 'hemispherewest', 'name', 'hemispherewest', 'proj', laea(0, -90), 'boundaries', [  # Peru projection 3587
         boundary('shppath', shppath, 'shpfilename', 'HemisphereSplit', 'proj', proj4326, 'orientation', 'reverse'),
         boundary('shppath', shppath, 'shpfilename', 'NorthAntarctica', 'proj', proj3031),
@@ -49 +49 @@
         boundary('shppath', shppath, 'shpfilename', 'RonneWestSummit', 'proj', proj3031),
         boundary('shppath', shppath, 'shpfilename', 'WestAntarctica2', 'proj', proj3031, 'orientation', 'reverse'),
-        boundary('shppath', shppath, 'shpfilename', 'SouthAntarctica', 'proj', proj3031)
-        ]
+        boundary('shppath', shppath, 'shpfilename', 'SouthAntarctica', 'proj', proj3031)]
     )
 )
@@ -59 +58 @@
         boundary('shppath', shppath, 'shpfilename', 'SouthAntarctica', 'proj', proj3031),
         boundary('shppath', shppath, 'shpfilename', 'RossIceShelf', 'proj', proj3031),
-        boundary('shppath', shppath, 'shpfilename', 'RossWestFront','proj', proj3031),
-        boundary('shppath', shppath, 'shpfilename', 'RossFront', 'proj', proj3031, 'orientation', 'reverse')
-        ]
+        boundary('shppath', shppath, 'shpfilename', 'RossWestFront', 'proj', proj3031),
+        boundary('shppath', shppath, 'shpfilename', 'RossFront', 'proj', proj3031, 'orientation', 'reverse')]
     )
 )
@@ -67 +65 @@
 #HemisphereEast: {{{
 sl.addbasin(
-    basin('continent', 'hemisphereeast', 'name', 'hemisphereeast', 'proj', laea(0, +90), 'boundaries', [ #Australian projection lat, long
+    basin('continent', 'hemisphereeast', 'name', 'hemisphereeast', 'proj', laea(0, +90), 'boundaries', [  #Australian projection lat, long
         boundary('shppath', shppath, 'shpfilename', 'HemisphereSplit', 'proj', proj4326),
         boundary('shppath', shppath, 'shpfilename', 'SouthAntarctica', 'proj', proj3031),
-        boundary('shppath', shppath, 'shpfilename','RossFront', 'proj', proj3031),
+        boundary('shppath', shppath, 'shpfilename', 'RossFront', 'proj', proj3031),
         boundary('shppath', shppath, 'shpfilename', 'RossWestFront', 'proj', proj3031),
-        boundary('shppath', shppath, 'shpfilename', 'EastAntarctica2', 'proj', proj3031,'orientation', 'reverse'),
-        boundary('shppath', shppath, 'shpfilename', 'NorthAntarctica', 'proj', proj3031)
-        ]
+        boundary('shppath', shppath, 'shpfilename', 'EastAntarctica2', 'proj', proj3031, 'orientation', 'reverse'),
+        boundary('shppath', shppath, 'shpfilename', 'NorthAntarctica', 'proj', proj3031)]
     )
 )
@@ -90 +87 @@
         boundary('shppath', shppath, 'shpfilename', 'RonneIceShelf', 'proj', proj3031),
         boundary('shppath', shppath, 'shpfilename', 'RonneEastSummit', 'proj', proj3031),
-        boundary('shppath', shppath, 'shpfilename', 'RonneBrunt', 'proj', proj3031)
-        ]
+        boundary('shppath', shppath, 'shpfilename', 'RonneBrunt', 'proj', proj3031)]
     )
 )
@@ -101 +97 @@
         boundary('shppath', shppath, 'shpfilename', 'RonneIceShelf', 'proj', proj3031),
         boundary('shppath', shppath, 'shpfilename', 'RonneEastSummit', 'proj', proj3031),
-        boundary('shppath', shppath, 'shpfilename', 'RonneFront', 'proj', proj3031, 'orientation', 'reverse')
-        ]
+        boundary('shppath', shppath, 'shpfilename', 'RonneFront', 'proj', proj3031, 'orientation', 'reverse')]
     )
 )
@@ -114 +109 @@
 hmin = hmin * 1000
 hmax = hmax * 1000
-tolerance = 100 # tolerance of 100m on Earth position when merging 3d meshes
+tolerance = 100  # tolerance of 100m on Earth position when merging 3d meshes
 threshold = 5
-defaultoptions = [
-    'KeepVertices', 0,
-    'MaxCornerAngle', 0.0000000001,
-    'NoBoundaryRefinement', 1
-]
+defaultoptions = ['KeepVertices', 0,
+                  'MaxCornerAngle', 0.0000000001,
+                  'NoBoundaryRefinement', 1]
 alreadyloaded = 0
 #}}}
@@ -134 +127 @@
 
     # Mesh
-    md = bamg(model(), 'domain', domain, 'subdomains', coastline, 'hmin', hmin, 'hmax', hmax, *defaultoptions) # NOTE: Unpacking defaultoptions with '*'
+    md = bamg(model(), 'domain', domain, 'subdomains', coastline, 'hmin', hmin, 'hmax', hmax, *defaultoptions)  # NOTE: Unpacking defaultoptions with '*'
 
     # Miscellaneous
@@ -254 +247 @@
     landelsconocean = landels[np.nonzero(sumconnectedisonocean)[0]]
 
-    ind1 = np.hstack((
-        md.mesh.elements[landelsconocean, 0],
-        md.mesh.elements[landelsconocean, 1],
-        md.mesh.elements[landelsconocean, 2]
-    ))
-    ind2 = np.hstack((
-        md.mesh.elements[landelsconocean, 1],
-        md.mesh.elements[landelsconocean, 2],
-        md.mesh.elements[landelsconocean, 0]
-    ))
+    ind1 = np.hstack((md.mesh.elements[landelsconocean, 0],
+                      md.mesh.elements[landelsconocean, 1],
+                      md.mesh.elements[landelsconocean, 2]))
+    ind2 = np.hstack((md.mesh.elements[landelsconocean, 1],
+                      md.mesh.elements[landelsconocean, 2],
+                      md.mesh.elements[landelsconocean, 0]))
 
     # Edge ind1 and ind2
@@ -325 +314 @@
         delHGLA = InterpFromMeshToMesh2d(index, longGLA, latGLA, delHGLA, longe, late)
 
-        # NOTE: For some reason, cannot apply pos to multiple arrays in a 
-        #       singlelike we might do in MATLAB. Instead, we iterate over 
+        # NOTE: For some reason, cannot apply pos to multiple arrays in a
+        #       singlelike we might do in MATLAB. Instead, we iterate over
         #       elements of pos.
         #
@@ -461 +450 @@
 md.solidearth.settings.rotation = 0
 md.solidearth.settings.viscous = 0
-md.solidearth.requested_outputs = [
-    'default',
-    'DeltaIceThickness',
-    'Sealevel',
-    'Bed',
-    'SealevelBarystaticIceMask',
-    'SealevelBarystaticOceanMask',
-]
+md.solidearth.requested_outputs = ['default',
+                                   'DeltaIceThickness',
+                                   'Sealevel',
+                                   'Bed',
+                                   'SealevelBarystaticIceMask',
+                                   'SealevelBarystaticOceanMask']
 md = solve(md, 'Transient')
 Seustatic = md.results.TransientSolution.Sealevel