Changeset 17590 for issm/trunk-jpl/test/Par/SquareSheetShelf.py

Timestamp:

03/28/14 09:15:19 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: WARNING: md.geometry.bed is now md.geometry.base and bathymetry is bed

File:

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

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

@@ -18 +18 @@
 xmax=max(md.mesh.x)
 md.geometry.thickness=hmax+(hmin-hmax)*(md.mesh.y.reshape(-1,1)-ymin)/(ymax-ymin)+0.1*(hmin-hmax)*(md.mesh.x.reshape(-1,1)-xmin)/(xmax-xmin)
-md.geometry.bed=-md.materials.rho_ice/md.materials.rho_water*md.geometry.thickness
+md.geometry.base=-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.)
-md.geometry.bed[pos]=bed_sheet
-md.geometry.surface=md.geometry.bed+md.geometry.thickness
+md.geometry.base[pos]=bed_sheet
+md.geometry.surface=md.geometry.base+md.geometry.thickness
 
 #Initial velocity