Changeset 23831 for issm/trunk-jpl/src

Timestamp:

04/08/19 16:08:02 (6 years ago)

Author:

kruegern

Message:

CHG: added a few minor changes, fixed minor runme bug

File:

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

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

@@ -57 +57 @@
                         if solution=='TransientSolution' and md.transient.isgroundingline:
                                 md = checkfield(md,'fieldname','geometry.bed','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
+                                if np.any(self.bed > self.base):
+                                    md.checkmessage('base<bed on one or more vertex')
+                                pos = np.where(md.mask.groundedice_levelset > 0)
+                                if np.any(np.abs(self.bed[pos]-self.base[pos]>10**-9):
+                                    md.checkmessage('equality base=bed on grounded ice violated')
+                                md = checkfield(md,'fieldname','geometry.bed','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
 
                 return md