Changeset 13678

Timestamp:

10/15/12 15:07:53 (12 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixed SquareSheetShelf.py (needed to reshape thickness)

File:

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

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

@@ -15 +15 @@
 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))