Changeset 16000

Timestamp:

08/28/13 13:25:52 (12 years ago)

Author:

seroussi

Message:

BUG: fixed problem with RoundSheetShelf python

File:

• issm/trunk-jpl/test/Par/RoundSheetShelf.py (modified) (1 diff)

issm/trunk-jpl/test/Par/RoundSheetShelf.py

@@ -76 +76 @@
 
 #Friction
+radius=1.e6
+shelfextent=2.e5
 md.friction.coefficient=20.*numpy.ones((md.mesh.numberofvertices,1))
-md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
+xelem=numpy.mean(md.mesh.x[md.mesh.elements.astype(int)-1],axis=1)
+yelem=numpy.mean(md.mesh.y[md.mesh.elements.astype(int)-1],axis=1)
+rad=numpy.sqrt(xelem**2+yelem**2)
+flags=numpy.zeros(md.mesh.numberofelements)
+pos=numpy.nonzero(rad>=(radius-shelfextent))
+md.friction.coefficient(md.mesh.elements(pos,:))=0.
 md.friction.p=numpy.ones((md.mesh.numberofelements,1))
 md.friction.q=numpy.ones((md.mesh.numberofelements,1))