Changeset 24256 for issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py

Timestamp:

10/18/19 06:52:19 (5 years ago)

Author:

bdef

Message:

BUG: fixing some negative number syntax

File:

• issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py (modified) (2 diffs)

issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py ¶

@@ -28 +28 @@
             raise IOError("SetIceShelfBC error message: ice front file '%s' not found." % icefrontfile)
         nodeinsideicefront = ContourToMesh(md.mesh.elements, md.mesh.x, md.mesh.y, icefrontfile, 'node', 2)
-        nodeonicefront = np.logical_and(md.mesh.vertexonboundary, nodeinsideicefront.reshape(- 1))
+        nodeonicefront = np.logical_and(md.mesh.vertexonboundary, nodeinsideicefront.reshape(-1))
     else:
         nodeonicefront = np.zeros((md.mesh.numberofvertices), bool)
@@ -68 +68 @@
         #reshape to rank - 2 if necessary to match spc arrays
         if np.ndim(md.inversion.vx_obs) == 1:
-            md.inversion.vx_obs = md.inversion.vx_obs.reshape(- 1, )
+            md.inversion.vx_obs = md.inversion.vx_obs.reshape(-1, )
         if np.ndim(md.inversion.vy_obs) == 1:
-            md.inversion.vy_obs = md.inversion.vy_obs.reshape(- 1, )
+            md.inversion.vy_obs = md.inversion.vy_obs.reshape(-1, )
         print("      boundary conditions for stressbalance model: spc set as observed velocities")
         md.stressbalance.spcvx[pos] = md.inversion.vx_obs[pos]