Index: /issm/trunk-jpl/test/NightlyRun/test336.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test336.py	(revision 17148)
+++ /issm/trunk-jpl/test/NightlyRun/test336.py	(revision 17149)
@@ -10,7 +10,27 @@
 md=triangle(model(),'../Exp/Square.exp',150000.)
 md=setmask(md,'','')
-md=parameterize(md,'../Par/SquareSheetLevelset.py') 
+md=parameterize(md,'../Par/SquareSheetConstrained.py') 
 md=setflowequation(md,'SSA','all')
 md.cluster=generic('name',oshostname(),'np',3)
+
+# parameters
+md.transient.isstressbalance=True
+md.transient.ismasstransport=False
+md.transient.isthermal=False
+md.transient.isgroundingline=False
+md.transient.isgia=False
+md.transient.islevelset=True
+
+# init levelset function
+ymin=numpy.min(md.mesh.y)
+ymax=numpy.max(md.mesh.y)
+xmin=min(md.mesh.x)
+xmax=max(md.mesh.x)
+xmed=(xmax+xmin)/2
+ymed=(ymax+ymin)/2
+distx=numpy.absolute(md.mesh.x.reshape(-1,1)-xmed)
+disty=numpy.absolute(md.mesh.y.reshape(-1,1)-ymed)
+md.mask.ice_levelset=numpy.maximum(distx,disty)-1.e5
+
 md=solve(md,TransientSolutionEnum())