Index: /issm/trunk-jpl/test/Par/RoundSheetShelf.par
===================================================================
--- /issm/trunk-jpl/test/Par/RoundSheetShelf.par	(revision 15609)
+++ /issm/trunk-jpl/test/Par/RoundSheetShelf.par	(revision 15610)
@@ -104,4 +104,5 @@
 pressureload=[pressureload 1*md.mask.elementonfloatingice(pressureload(:,end)) + 0*md.mask.elementongroundedice(pressureload(:,end))];
 md.diagnostic.icefront=pressureload;
+md.mask.icelevelset(find(md.mesh.vertexonboundary))=0;
 md.balancethickness.spcthickness=NaN*ones(md.mesh.numberofvertices,1);
 md.prognostic.spcthickness=NaN*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/Par/RoundSheetShelf.py
===================================================================
--- /issm/trunk-jpl/test/Par/RoundSheetShelf.py	(revision 15609)
+++ /issm/trunk-jpl/test/Par/RoundSheetShelf.py	(revision 15610)
@@ -113,4 +113,6 @@
 pressureload=numpy.hstack((pressureload,(1*md.mask.elementonfloatingice[pressureload[:,-1].astype(int)-1] + 0*md.mask.elementongroundedice[pressureload[:,-1].astype(int)-1]).reshape(-1,1)))
 md.diagnostic.icefront=pressureload
+pos=numpy.nonzero(md.mesh.vertexonboundary)
+md.mask.icelevelset[pos]=0
 md.balancethickness.spcthickness=float('nan')*numpy.ones((md.mesh.numberofvertices,1))
 md.prognostic.spcthickness=float('nan')*numpy.ones((md.mesh.numberofvertices,1))