Index: /issm/trunk-jpl/test/NightlyRun/test336.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test336.py	(revision 17250)
+++ /issm/trunk-jpl/test/NightlyRun/test336.py	(revision 17251)
@@ -8,9 +8,10 @@
 from MatlabFuncs import *
 
-md=triangle(model(),'../Exp/Square.exp',150000.)
+md=triangle(model(),'../Exp/Square.exp',100000.)
 md=setmask(md,'','')
 md=parameterize(md,'../Par/SquareSheetConstrained.py') 
 md=setflowequation(md,'SSA','all')
-md.cluster=generic('name',oshostname(),'np',3)
+md.cluster=generic('name',oshostname(),'np',1)
+# md.cluster=generic('name',oshostname(),'np',3)
 
 # parameters
@@ -23,13 +24,15 @@
 
 # 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
+xmin, xmax=min(md.mesh.x), max(md.mesh.x)
+ymin, ymax=min(md.mesh.y), max(md.mesh.y)
+xmed, ymed=(xmax+xmin)/2.,(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)-0.2*(xmax-xmin)
+
+D=numpy.sqrt(numpy.power(md.mesh.x.reshape(-1,1)-xmed,2.)+numpy.power(md.mesh.y.reshape(-1,1)-ymed,2.))
+R=0.2*(xmax-xmin)
+md.mask.ice_levelset=D-R
 
 md=solve(md,TransientSolutionEnum())