Changeset 15806

Timestamp:

08/12/13 16:47:12 (12 years ago)

Author:

seroussi

Message:

BUG: Set up mask.icelevelset values before creating dirichlets

File:

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

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

@@ -38 +38 @@
         md.stressbalance.loadingforce=0*numpy.ones((md.mesh.numberofvertices,3))
 
+        #Icefront position
+        pos=numpy.nonzero(nodeonicefront)[0]
+        md.mask.icelevelset[pos]=0
+
         #First find segments that are not completely on the front
         if md.mesh.dimension==2:
@@ -64 +68 @@
         else:
                 print "      boundary conditions for stressbalance model: spc set as zero"
-
-        #Icefront position
-        pos=numpy.nonzero(nodeonicefront)[0]
-        md.mask.icelevelset[pos]=0
 
         #Create zeros basalforcings and surfaceforcings