Changeset 27434

Timestamp:

11/28/22 01:31:42 (2 years ago)

Author:

jdquinn

Message:

BUG: Tricky matrix handling in check

File:

• issm/trunk-jpl/src/m/classes/initialization.py (modified) (1 diff)

issm/trunk-jpl/src/m/classes/initialization.py

@@ -82 +82 @@
             md = checkfield(md, 'fieldname', 'initialization.vy', 'NaN', 1, 'Inf', 1, 'size', [md.mesh.numberofvertices])
             # Triangle with zero velocity
-            if np.any(np.logical_and(np.sum(np.abs(md.initialization.vx[md.mesh.elements - 1]), axis=1) == 0, np.sum(np.abs(md.initialization.vy[md.mesh.elements - 1]), axis=1) == 0, np.min(md.mask.ice_levelset[md.mesh.elements - 1], axis=1).reshape(-1, 1) < 0)):
+            if np.any(np.logical_and(np.sum(np.abs(md.initialization.vx[md.mesh.elements - 1]), axis=1).reshape(-1, 1) == 0, np.sum(np.abs(md.initialization.vy[md.mesh.elements - 1]), axis=1).reshape(-1, 1) == 0, np.min(md.mask.ice_levelset[md.mesh.elements - 1], axis=1).reshape(-1, 1) < 0)):
                 md.checkmessage('at least one triangle has all its vertices with a zero velocity')
         if 'ThermalAnalysis' in analyses and not solution == 'TransientSolution' and not md.transient.isthermal: