Changeset 23033

Timestamp:

07/31/18 15:03:26 (7 years ago)

Author:

caronlam

Message:

fixing consitency checks in slr.py

File:

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

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

@@ -161 +161 @@
                 if (size(self.love_h,0) != size(self.love_k,0) | size(self.love_h,0) != size(self.love_l,0)):
                         error('slr error message: love numbers should be provided at the same level of accuracy')
+
+                #cross check that whereever we have an ice load, the mask is <0 on each vertex: 
+                pos=np.where(self.deltathickness)
+                maskpos=md.mask.ice_levelset[md.mesh.elements[pos,:]] 
+                els=np.where(maskpos>0)
+                if len(els[1])>0:
+                        warnings.warn('slr checkconsistency fail: there are elements with ice loads where some vertices are not on the ice!')
+                
+                #check that  if geodetic is requested, we are a mesh3dsurface model (planet), or if we are not, 
+                #a coupler to a planet model is provided. 
+                if self.geodetic:
+                        if md.transient.iscoupler: 
+                                #we are good; 
+                        else
+                                if domaintype(md.mesh)=='mesh3dsurface':
+                                        #we are good
+                                else
+                                        error('model is requesting geodetic computations without being a mesh3dsurface, or being coupled to one!');
                 return md
         # }}}