Index: /issm/trunk-jpl/test/NightlyRun/yellowstone.lsf
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/yellowstone.lsf	(revision 22884)
+++ /issm/trunk-jpl/test/NightlyRun/yellowstone.lsf	(revision 22884)
@@ -0,0 +1,36 @@
+!/bin/tcsh
+#BSUB -P P93300301
+#BSUB -W 12:0
+#BSUB -n 1
+#BSUB -J test101
+#BSUB -o /glade/u/home/elarour/trunk-jpl/exeuction/test101-03-06-2017-16-12-56-77385/test101.out 
+#BSUB -o /glade/u/home/elarour/trunk-jpl/exeuction/test101-03-06-2017-16-12-56-77385/test101.err 
+#BSUB -q premium
+module load ncarenv/1.0
+module load ncarbinlibs/1.1
+module load perlmods/5.0
+module load gmake/4.1
+module load python/2.7.7
+module load all-python-libs
+module load git/2.3.0
+module load intel/15.0.3
+module load mkl/11.1.2
+module load esmf
+module load esmf-6.3.0rp1-defio-mpi-O
+module load netcdf-mpi/4.3.3.1
+module load pnetcdf/1.6.1
+module load ncarcompilers/1.0
+module load cmake/3.0.2
+module load matlab/R2015b
+module load issm
+export OMP_STACKSIZE=256M
+export MP_LABELIO=yes
+export MP_INFOLEVEL=2
+export MP_SHARED_MEMORY=yes
+export MP_EUILIB=us
+export MP_STDOUTMODE=unordered
+export MP_RC_USE_LMC=yes
+export MP_MPILIB=mpich2
+cd /glade/u/home/elarour/trunk-jpl/exeuction/test101-03-06-2017-16-12-56-77385/
+
+mpirun.lsf -np 1 /glade/u/home/elarour/trunk-jpl/build-fw/bin/issm.exe StressbalanceSolution /glade/u/home/elarour/trunk-jpl/exeuction/test101-03-06-2017-16-12-56-77385 test101
Index: /issm/trunk-jpl/test/NightlyRun/yellowstone.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/yellowstone.m	(revision 22884)
+++ /issm/trunk-jpl/test/NightlyRun/yellowstone.m	(revision 22884)
@@ -0,0 +1,46 @@
+%Test Name: SquareShelfConstrainedStressSSA2d
+md=triangle(model(),'../Exp/Square.exp',50000.);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setflowequation(md,'SSA','all');
+
+%md.cluster=generic('name',oshostname(),'np',2);
+md.cluster=yellowstone;
+md.cluster.login='elarour';
+md.cluster.codepath='/glade/u/home/elarour/trunk-jpl/build-fw/bin';
+md.cluster.executionpath='/glade/u/home/elarour/trunk-jpl/execution';
+
+%output
+md.stressbalance.requested_outputs={'default','DeviatoricStressxx','DeviatoricStressyy','DeviatoricStressxy','MassFlux1','MassFlux2','MassFlux3','MassFlux4','MassFlux5','MassFlux6'};
+md.outputdefinition.definitions={...
+	massfluxatgate('name','MassFlux1','profilename',['../Exp/MassFlux1.exp'],'definitionstring','Outputdefinition1'),...
+	massfluxatgate('name','MassFlux2','profilename',['../Exp/MassFlux2.exp'],'definitionstring','Outputdefinition2'),...
+	massfluxatgate('name','MassFlux3','profilename',['../Exp/MassFlux3.exp'],'definitionstring','Outputdefinition3'),...
+	massfluxatgate('name','MassFlux4','profilename',['../Exp/MassFlux4.exp'],'definitionstring','Outputdefinition4'),...
+	massfluxatgate('name','MassFlux5','profilename',['../Exp/MassFlux5.exp'],'definitionstring','Outputdefinition5'),...
+	massfluxatgate('name','MassFlux6','profilename',['../Exp/MassFlux6.exp'],'definitionstring','Outputdefinition6')...
+	};
+
+md=solve(md,'Stressbalance');
+
+%Fields and tolerances to track changes
+field_names     ={'Vx','Vy','Vel','Pressure',...
+	'DeviatoricStressxx','DeviatoricStressyy','DeviatoricStressxy','MassFlux1','MassFlux2','MassFlux3','MassFlux4','MassFlux5','MassFlux6'};
+field_tolerances={2e-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),...
+	(md.results.StressbalanceSolution.Vy),...
+	(md.results.StressbalanceSolution.Vel),...
+	(md.results.StressbalanceSolution.Pressure),...
+	(md.results.StressbalanceSolution.DeviatoricStressxx),...
+	(md.results.StressbalanceSolution.DeviatoricStressyy),...
+	(md.results.StressbalanceSolution.DeviatoricStressxy),...
+	(md.results.StressbalanceSolution.MassFlux1),...
+	(md.results.StressbalanceSolution.MassFlux2),...
+	(md.results.StressbalanceSolution.MassFlux3),...
+	(md.results.StressbalanceSolution.MassFlux4),...
+	(md.results.StressbalanceSolution.MassFlux5),...
+	(md.results.StressbalanceSolution.MassFlux6)...
+	};