Index: /issm/trunk-jpl/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 16043)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 16044)
@@ -148,4 +148,6 @@
 	case 509, name='PigSteaHO';
 	case 510, name='PigSteaFS';
+	case 530, name='PigBalVel1';
+	case 531, name='PigBalVel2';
 	case 511, name='PigCMBFS';
 	case 512, name='PigCMDragHO';
Index: /issm/trunk-jpl/test/NightlyRun/IdToName.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 16043)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 16044)
@@ -158,4 +158,6 @@
 	515  : 'PigTherTranSUPG',
 	516  : 'PigTherSteaSUPG',
+	530  : 'PigBalVel1',
+	531  : 'PigBalVel2',
 	601  : '79NorthMasstransp2d',
 	602  : '79NorthMasstransp2dDG',
Index: /issm/trunk-jpl/test/NightlyRun/test530.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test530.m	(revision 16044)
+++ /issm/trunk-jpl/test/NightlyRun/test530.m	(revision 16044)
@@ -0,0 +1,15 @@
+md=triangle(model(),'../Exp/Pig.exp',20000.);
+md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
+md=parameterize(md,'../Par/Pig.par');
+md=setflowequation(md,'SSA','all');
+md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,BalancevelocitySolutionEnum());
+
+%Fields and tolerances to track changes
+field_names     ={'SurfaceSlopeX','SurfaceSlopeX','Vel'};
+field_tolerances={1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.BalancevelocitySolution.SurfaceSlopeX),...
+	(md.results.BalancevelocitySolution.SurfaceSlopeY),...
+	(md.results.BalancevelocitySolution.Vel),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test530.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test530.py	(revision 16044)
+++ /issm/trunk-jpl/test/NightlyRun/test530.py	(revision 16044)
@@ -0,0 +1,24 @@
+from model import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from EnumDefinitions import *
+from solve import *
+from MatlabFuncs import *
+
+md=triangle(model(),'../Exp/Pig.exp',20000.)
+md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp')
+md=parameterize(md,'../Par/Pig.py')
+md=setflowequation(md,'SSA','all')
+md.cluster=generic('name',oshostname(),'np',3)
+md=solve(md,BalancevelocitySolutionEnum())
+
+# Fields and tolerances to track changes
+field_names     =['SurfaceSlopeX','SurfaceSlopeY','Vel']
+field_tolerances=[1e-13,1e-13,1e-13]
+field_values=[\
+		md.results.BalancevelocitySolution.SurfaceSlopeX,\
+		md.results.BalancevelocitySolution.SurfaceSlopeY,\
+		md.results.BalancevelocitySolution.Vel,\
+		]
Index: /issm/trunk-jpl/test/NightlyRun/test531.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test531.m	(revision 16044)
+++ /issm/trunk-jpl/test/NightlyRun/test531.m	(revision 16044)
@@ -0,0 +1,17 @@
+md=triangle(model(),'../Exp/Pig.exp',20000.);
+md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
+md=parameterize(md,'../Par/Pig.par');
+md.initialization.vx(:)=0;
+md.initialization.vy(:)=0;
+md=setflowequation(md,'SSA','all');
+md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,BalancevelocitySolutionEnum());
+
+%Fields and tolerances to track changes
+field_names     ={'SurfaceSlopeX','SurfaceSlopeX','Vel'};
+field_tolerances={1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.BalancevelocitySolution.SurfaceSlopeX),...
+	(md.results.BalancevelocitySolution.SurfaceSlopeY),...
+	(md.results.BalancevelocitySolution.Vel),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test531.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test531.py	(revision 16044)
+++ /issm/trunk-jpl/test/NightlyRun/test531.py	(revision 16044)
@@ -0,0 +1,26 @@
+from model import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from EnumDefinitions import *
+from solve import *
+from MatlabFuncs import *
+
+md=triangle(model(),'../Exp/Pig.exp',20000.)
+md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp')
+md=parameterize(md,'../Par/Pig.py')
+md.initialization.vx[:]=0.
+md.initialization.vy[:]=0.
+md=setflowequation(md,'SSA','all')
+md.cluster=generic('name',oshostname(),'np',3)
+md=solve(md,BalancevelocitySolutionEnum())
+
+# Fields and tolerances to track changes
+field_names     =['SurfaceSlopeX','SurfaceSlopeY','Vel']
+field_tolerances=[1e-13,1e-13,1e-13]
+field_values=[\
+		md.results.BalancevelocitySolution.SurfaceSlopeX,\
+		md.results.BalancevelocitySolution.SurfaceSlopeY,\
+		md.results.BalancevelocitySolution.Vel,\
+		]