Index: /issm/trunk-jpl/test/NightlyRun/test408.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test408.m	(revision 21603)
+++ /issm/trunk-jpl/test/NightlyRun/test408.m	(revision 21604)
@@ -3,6 +3,15 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
+md.geometry.bed=md.geometry.base;
+pos=find(md.mask.groundedice_levelset<0);
+md.geometry.bed(pos)=md.geometry.base(pos)-10;
+md.friction=frictioncoulomb();
+md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
+md.friction.p=ones(md.mesh.numberofelements,1);
+md.friction.q=ones(md.mesh.numberofelements,1);
+md.friction.coefficientcoulomb=0.02*ones(md.mesh.numberofvertices,1);
 md=extrude(md,3,1);
 md.transient.isthermal=0;
+md.transient.isgroundingline=1;
 md=setflowequation(md,'SSA','all');
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk-jpl/test/NightlyRun/test408.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test408.py	(revision 21603)
+++ /issm/trunk-jpl/test/NightlyRun/test408.py	(revision 21604)
@@ -9,11 +9,21 @@
 from setflowequation import *
 from solve import *
+from frictioncoulomb import frictioncoulomb
 
 md=triangle(model(),'../Exp/Square.exp',150000.)
 md=setmask(md,'../Exp/SquareShelf.exp','')
 md=parameterize(md,'../Par/SquareSheetShelf.py')
+md.geometry.bed=md.geometry.base
+pos=numpy.nonzero(md.mask.groundedice_levelset<0.)
+md.geometry.bed[pos]=md.geometry.base[pos]-10
+md.friction=frictioncoulomb()
+md.friction.coefficient=20*np.ones(md.mesh.numberofvertices)
+md.friction.p=1*np.ones(md.mesh.numberofvertices)
+md.friction.q=1*np.ones(md.mesh.numberofvertices)
+md.friction.coefficientcoulomb=0.02*np.ones(md.mesh.numberofvertices)
+md.transient.isthermal=False
+md.transient.isgroundingline=True
 md.extrude(3,1.)
 md=setflowequation(md,'SSA','all')
-md.transient.isthermal=False
 md.cluster=generic('name',gethostname(),'np',3)
 md.transient.requested_outputs=['default','GroundedArea','FloatingArea','TotalFloatingBmb','TotalGroundedBmb','TotalSmb']