Index: /issm/trunk-jpl/test/NightlyRun/test801.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test801.py	(revision 17447)
+++ /issm/trunk-jpl/test/NightlyRun/test801.py	(revision 17448)
@@ -9,5 +9,5 @@
 from MatlabFuncs import *
 
-md=triangle(model(),'../Exp/Square.exp',30000)
+md=triangle(model(),'../Exp/Square.exp',50000)
 md=setmask(md,'','')
 md=parameterize(md,'../Par/ValleyGlacierShelf.py')
Index: /issm/trunk-jpl/test/Par/ValleyGlacierShelf.py
===================================================================
--- /issm/trunk-jpl/test/Par/ValleyGlacierShelf.py	(revision 17448)
+++ /issm/trunk-jpl/test/Par/ValleyGlacierShelf.py	(revision 17448)
@@ -0,0 +1,98 @@
+import os.path
+import netCDF4
+import numpy
+import inspect
+import math
+from verbose import verbose
+from InterpFromMeshToMesh2d import InterpFromMeshToMesh2d
+from paterson import paterson
+from SetIceShelfBC import SetIceShelfBC
+
+#Start defining model parameters here
+x=md.mesh.x.reshape(-1,1)
+y=md.mesh.y.reshape(-1,1)
+xmin, xmax = min(x), max(x)
+ymin, ymax = min(y), max(y)
+Lx=(xmax-xmin)
+Ly=(ymax-ymin)
+xm,ym = (xmin+xmax)/2., (ymin+ymax)/2.
+
+#Geometry: U-shaped valley in y direction
+thk_center = 1000.
+thk_margin = 0.5*thk_center
+bmax=0.
+bmin=-thk_center*md.materials.rho_ice/md.materials.rho_water
+
+alpha=2./3.
+slope = (bmin-bmax)*(x-xmin)/(Lx*alpha) + bmax
+md.geometry.surface= (thk_center+bmax) + slope 
+md.geometry.bed=bmax + slope + 4./Ly**2*(thk_center-thk_margin)*(numpy.power(y-ym,2))
+md.geometry.thickness=md.geometry.surface - md.geometry.bed
+md.geometry.bathymetry = md.geometry.bed
+
+#Mask
+md.mask.ice_levelset=x - alpha*Lx
+md.mask.groundedice_levelset= numpy.ones((md.mesh.numberofvertices,1))
+
+#Initial velocity 
+md.initialization.vx=numpy.zeros((md.mesh.numberofvertices,1))
+md.initialization.vy=numpy.zeros((md.mesh.numberofvertices,1))
+md.initialization.vz=numpy.zeros((md.mesh.numberofvertices,1))
+md.initialization.pressure=numpy.zeros((md.mesh.numberofvertices,1))
+
+#Materials
+md.initialization.temperature=(273.15-5.)*numpy.ones((md.mesh.numberofvertices,1))
+md.materials.rheology_B=paterson(md.initialization.temperature)
+md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
+
+#Groundingline
+md.groundingline.migration='SubelementMigration'
+
+#Damage
+md.damage.D=numpy.zeros((md.mesh.numberofvertices,1))
+
+#Surface mass balance and basal melting
+md.surfaceforcings.mass_balance=0.3*numpy.ones((md.mesh.numberofvertices,1))
+md.basalforcings.melting_rate=md.surfaceforcings.mass_balance
+
+#Friction
+md.friction.coefficient=20.*numpy.ones((md.mesh.numberofvertices,1))
+md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0.
+md.friction.p=numpy.ones((md.mesh.numberofelements,1))
+md.friction.q=numpy.ones((md.mesh.numberofelements,1))
+
+#Transient
+md.transient.isstressbalance=True
+md.transient.islevelset=True
+md.transient.ismasstransport=False
+md.transient.isthermal=False
+md.transient.isgroundingline=True
+md.transient.isgia=False
+
+#Stressbalance
+md.stressbalance.maxiter=100
+md.stressbalance.viscosity_overshoot=0.0
+md.stressbalance.restol=0.05
+md.stressbalance.reltol=0.05
+md.stressbalance.abstol=float('nan')
+
+#Masstransport
+md.masstransport.calvingrate=0.*numpy.ones((md.mesh.numberofvertices,1))
+md.masstransport.stabilization=1.
+
+#Numerical parameters
+md.thermal.stabilization=1.
+md.settings.waitonlock=30
+md.steadystate.reltol=0.05
+md.timestepping.time_step=1.
+md.timestepping.final_time=3.
+
+#Verbose
+md.verbose = verbose(0)
+
+#Deal with boundary conditions:
+md = SetIceShelfBC(md)
+
+#Change name so that no tests have the same name
+if len(inspect.stack()) > 2:
+	md.miscellaneous.name = os.path.basename(inspect.stack()[2][1]).split('.')[0]