Index: /issm/trunk-jpl/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 16379)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 16380)
@@ -170,4 +170,5 @@
 	case 612, name='79NorthCMBalThic2dDG';
 	case 613, name='79NorthCMBalThicVxVy';
+	case 701, name='FlowbandFSshelf';
 	case 1101, name='ISMIPAHO';
 	case 1102, name='ISMIPAFS';
Index: /issm/trunk-jpl/test/NightlyRun/test701.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test701.m	(revision 16380)
+++ /issm/trunk-jpl/test/NightlyRun/test701.m	(revision 16380)
@@ -0,0 +1,58 @@
+x =[1:100:3000]';
+h=linspace(1000,300,numel(x))';
+b=-917/1023*h;
+
+md=bamgflowband(model(),x,b+h,b,'hmax',80);
+
+%Geometry
+md.geometry.surface   = interp1(x,b+h,md.mesh.x);
+md.geometry.bed       = interp1(x,b,md.mesh.x);
+md.geometry.thickness = md.geometry.surface-md.geometry.bed;
+
+%mask
+md.mask.ice_levelset  = ones(md.mesh.numberofvertices,1);
+md.mask.ice_levelset(find(vertexflags(md.mesh,2))) = 0;
+md.mask.groundedice_levelset = double(md.mesh.x<0)-.5;
+
+%materials
+md.initialization.temperature=(273-20)*ones(md.mesh.numberofvertices,1);
+md.materials.rheology_B=paterson(md.initialization.temperature);
+md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.damage.D=zeros(md.mesh.numberofvertices,1);
+md.damage.spcdamage=NaN(md.mesh.numberofvertices,1);
+
+%friction
+md.friction.coefficient=zeros(md.mesh.numberofvertices,1);
+md.friction.coefficient(find(vertexflags(md.mesh,1)))=20;
+md.friction.p=ones(md.mesh.numberofelements,1);
+md.friction.q=ones(md.mesh.numberofelements,1);
+
+%Boundary conditions
+md.stressbalance.referential  = NaN*ones(md.mesh.numberofvertices,6);
+md.stressbalance.loadingforce = 0*ones(md.mesh.numberofvertices,3);
+md.stressbalance.spcvx = NaN*ones(md.mesh.numberofvertices,1);
+md.stressbalance.spcvy = NaN*ones(md.mesh.numberofvertices,1);
+md.stressbalance.spcvz = NaN*ones(md.mesh.numberofvertices,1);
+md.stressbalance.spcvx(find(vertexflags(md.mesh,4)))=0;
+md.stressbalance.spcvy(find(vertexflags(md.mesh,4)))=0;
+%md.stressbalance.spcvx(find(vertexflags(md.mesh,4) | vertexflags(md.mesh,2)))=0;
+%md.stressbalance.spcvy(find(vertexflags(md.mesh,4) | vertexflags(md.mesh,2)))=0;
+
+%Misc
+md=setflowequation(md,'FS','all');
+md.flowequation.fe_FS='TaylorHood';
+md.stressbalance.abstol=NaN;
+md.miscellaneous.name = 'flowline';
+
+%Go solve
+md.cluster=generic('np',2);
+md=solve(md,StressbalanceSolutionEnum());
+
+%Fields and tolerances to track changes
+field_names     ={'Vx','Vy','Vel','Pressure'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.StressbalanceSolution.Vx),...
+	(md.results.StressbalanceSolution.Vy),...
+	(md.results.StressbalanceSolution.Vel),...
+	(md.results.StressbalanceSolution.Pressure)};