Index: /issm/trunk-jpl/test/NightlyRun/test540.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test540.m	(revision 23828)
+++ /issm/trunk-jpl/test/NightlyRun/test540.m	(revision 23829)
@@ -12,4 +12,5 @@
 md.frontalforcings.meltingrate = zeros(md.mesh.numberofvertices,1);
 md.transient.ismovingfront = 1;
+md.transient.isgroundingline = 1;
 md.levelset.spclevelset = NaN(md.mesh.numberofvertices,1);
 pos = find(md.mesh.vertexonboundary);
Index: /issm/trunk-jpl/test/NightlyRun/test540.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test540.py	(revision 23828)
+++ /issm/trunk-jpl/test/NightlyRun/test540.py	(revision 23829)
@@ -21,4 +21,5 @@
 md.frontalforcings.meltingrate = np.zeros((md.mesh.numberofvertices,))
 md.transient.ismovingfront = 1
+md.transient.isgroundingline = 1
 md.levelset.spclevelset = float('NaN') * np.ones((md.mesh.numberofvertices,))
 pos = np.where(md.mesh.vertexonboundary)
Index: /issm/trunk-jpl/test/Par/Pig.par
===================================================================
--- /issm/trunk-jpl/test/Par/Pig.par	(revision 23828)
+++ /issm/trunk-jpl/test/Par/Pig.par	(revision 23829)
@@ -9,4 +9,5 @@
 surface   = archread('../Data/Pig.arch','surface');
 thickness = archread('../Data/Pig.arch','thickness');
+bed       = archread('../Data/Pig.arch','bed');
 
 x			 = x{1};
@@ -17,4 +18,5 @@
 surface	 = surface{1};
 thickness = thickness{1};
+bed       = bed{1};
 
 md.inversion.vx_obs   =InterpFromMeshToMesh2d(index,x,y,vx_obs,md.mesh.x,md.mesh.y);
@@ -23,4 +25,7 @@
 md.geometry.thickness=InterpFromMeshToMesh2d(index,x,y,thickness,md.mesh.x,md.mesh.y);
 md.geometry.base=md.geometry.surface-md.geometry.thickness;
+md.geometry.bed =md.geometry.base;
+pos = find(md.mask.groundedice_levelset<0);
+md.geometry.bed(pos) =InterpFromMeshToMesh2d(index,x,y,bed,md.mesh.x(pos),md.mesh.y(pos));
 clear surface thickness vx_obs vy_obs x y index;
 md.initialization.vx=md.inversion.vx_obs;