Index: ../trunk-jpl/test/Par/SquareSheetShelf.py
===================================================================
--- ../trunk-jpl/test/Par/SquareSheetShelf.py	(revision 13677)
+++ ../trunk-jpl/test/Par/SquareSheetShelf.py	(revision 13678)
@@ -14,7 +14,7 @@
 hmax=1000.
 ymin=min(md.mesh.y)
 ymax=max(md.mesh.y)
-md.geometry.thickness=hmax+(hmin-hmax)*(md.mesh.y-ymin)/(ymax-ymin)
+md.geometry.thickness=hmax+(hmin-hmax)*(md.mesh.y.reshape(-1,1)-ymin)/(ymax-ymin)
 md.geometry.bed=-md.materials.rho_ice/md.materials.rho_water*md.geometry.thickness
 bed_sheet=-md.materials.rho_ice/md.materials.rho_water*(hmax+(hmin-hmax)*(ymax/2-ymin)/(ymax-ymin))
 pos=numpy.nonzero(md.mesh.y<=ymax/2.)