Index: /issm/trunk-jpl/test/Par/RoundSheetShelf.py
===================================================================
--- /issm/trunk-jpl/test/Par/RoundSheetShelf.py	(revision 20991)
+++ /issm/trunk-jpl/test/Par/RoundSheetShelf.py	(revision 20992)
@@ -14,5 +14,5 @@
 hmin=300.
 hmax=1000.
-radius=numpy.sqrt(md.mesh.x*md.mesh.x+md.mesh.y*md.mesh.y)
+radius=numpy.sqrt(md.mesh.x.reshape(-1,1)*md.mesh.x.reshape(-1,1)+md.mesh.y.reshape(-1,1)*md.mesh.y.reshape(-1,1))
 ymin=numpy.min(radius)
 ymax=numpy.max(radius)