Index: /issm/trunk-jpl/test/NightlyRun/test101.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test101.m	(revision 16392)
+++ /issm/trunk-jpl/test/NightlyRun/test101.m	(revision 16393)
@@ -3,13 +3,14 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=setflowequation(md,'SSA','all');
-md.cluster=generic('name',oshostname(),'np',3);
-md.stressbalance.requested_outputs={'StressTensor'};
+md.cluster=generic('name',oshostname(),'np',2);
+md.stressbalance.requested_outputs={'StressTensor','MassFlux1','MassFlux2','MassFlux3','MassFlux4','MassFlux5','MassFlux6'};
 md=solve(md,StressbalanceSolutionEnum());
 
 %Fields and tolerances to track changes
 field_names     ={'Vx','Vy','Vel','Pressure',...
-	'StressTensorxx','StressTensoryy','StressTensorxy'};
+	'StressTensorxx','StressTensoryy','StressTensorxy','MassFlux1','MassFlux2','MassFlux3','MassFlux4','MassFlux5','MassFlux6'};
 field_tolerances={1e-13,1e-13,1e-13,1e-13,...
-	1e-13,1e-13,1e-13};
+	1e-13,1e-13,1e-13,...
+	1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13 };
 field_values={...
 	(md.results.StressbalanceSolution.Vx),...
@@ -20,3 +21,9 @@
 	(md.results.StressbalanceSolution.StressTensoryy),...
 	(md.results.StressbalanceSolution.StressTensorxy),...
+	(md.results.StressbalanceSolution.MassFlux1),...
+	(md.results.StressbalanceSolution.MassFlux2),...
+	(md.results.StressbalanceSolution.MassFlux3),...
+	(md.results.StressbalanceSolution.MassFlux4),...
+	(md.results.StressbalanceSolution.MassFlux5),...
+	(md.results.StressbalanceSolution.MassFlux6)...
 	};
Index: /issm/trunk-jpl/test/NightlyRun/test101.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test101.py	(revision 16392)
+++ /issm/trunk-jpl/test/NightlyRun/test101.py	(revision 16393)
@@ -13,12 +13,13 @@
 md=setflowequation(md,'SSA','all')
 md.cluster=generic('name',oshostname(),'np',3)
-md.stressbalance.requested_outputs=['StressTensor']
+md.stressbalance.requested_outputs=['StressTensor','MassFlux1','MassFlux2','MassFlux3','MassFlux4','MassFlux5','MassFlux6']
 md=solve(md,StressbalanceSolutionEnum())
 
 #Fields and tolerances to track changes
 field_names     =['Vx','Vy','Vel','Pressure',\
-	'StressTensorxx','StressTensoryy','StressTensorxy']
+	'StressTensorxx','StressTensoryy','StressTensorxy','MassFlux1','MassFlux2','MassFlux3','MassFlux4','MassFlux5','MassFlux6']
 field_tolerances=[1e-13,1e-13,1e-13,1e-13,\
-	1e-13,1e-13,1e-13]
+	1e-13,1e-13,1e-13,\
+	1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13 ]
 field_values=[\
 	md.results.StressbalanceSolution.Vx,\
@@ -29,3 +30,8 @@
 	md.results.StressbalanceSolution.StressTensoryy,\
 	md.results.StressbalanceSolution.StressTensorxy,\
-	]
+	md.results.StressbalanceSolution.MassFlux1,\
+	md.results.StressbalanceSolution.MassFlux2,\
+	md.results.StressbalanceSolution.MassFlux3,\
+	md.results.StressbalanceSolution.MassFlux4,\
+	md.results.StressbalanceSolution.MassFlux5,\
+	md.results.StressbalanceSolution.MassFlux6 ]