Index: /issm/trunk-jpl/test/NightlyRun/test101.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test101.js	(revision 20591)
+++ /issm/trunk-jpl/test/NightlyRun/test101.js	(revision 20591)
@@ -0,0 +1,15 @@
+console.log('creating model'); 
+var md = new model();
+
+console.log('meshing');
+triangle(md,square[0],40000); 
+
+console.log('parameterization');
+setmask(md,'all','');
+parameterize(md);
+setflowequation(md,'SSA','all');
+md.verbose.solution=2;  md.verbose.convergence=0;
+
+md=solve(md,StressbalanceSolutionEnum(),'checkconsistency','no');
+
+console.log(md.results['StressbalanceSolution'][0]['Vel']);
Index: /issm/trunk-jpl/test/NightlyRun/test102.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test102.js	(revision 20591)
+++ /issm/trunk-jpl/test/NightlyRun/test102.js	(revision 20591)
@@ -0,0 +1,20 @@
+var md = new model();
+triangle(md,square[0],180000.);
+setmask(md, 'all', '');
+parameterize(md);
+md.extrude(md, 3, 2.);
+setflowequation(md, 'SSA', 'all');
+//md.cluster = new generic('name', oshostname(), 'np', 3); //oshostname
+md=solve(md, StressbalanceSolutionEnum(), 'checkconsistency', 'no');
+
+var field_names      = ['Vx','Vy','Vz','Vel','Pressure'];
+var field_tolerances = [1e-13,1e-13,1e-13,1e-13,1e-13];
+var field_values = [
+    md.results.StressbalanceSolution.Vx,
+    md.results.StressbalanceSolution.Vy,
+    md.results.StressbalanceSolution.Vz,
+    md.results.StressbalanceSolution.Vel,
+    md.results.StressbalanceSolution.Pressure
+];
+
+plotmodel(md, 'data', md.results.StressbalanceSolution[0].Vel);
Index: /issm/trunk-jpl/test/NightlyRun/test104.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test104.js	(revision 20591)
+++ /issm/trunk-jpl/test/NightlyRun/test104.js	(revision 20591)
@@ -0,0 +1,20 @@
+//Test Name: SquareShelfConstrainedStressFS
+var md = new model();
+triangle(md, square[0], 180000.);
+setmask(md, 'all', '');
+parameterize(md);
+md.extrude(md, 3, 2.);
+setflowequation(md, 'FS', 'all');
+md.cluster = new generic('name', oshostname(), 'np', 3);
+md = solve(md, StressbalanceSolutionEnum());
+
+//Fields and tolerances to track changes
+var field_names      = ['Vx','Vy','Vz','Vel','Pressure'];
+var field_tolerances = [1e-08,1e-08,1e-06,1e-08,1e-08];
+var field_values = [
+    md.results.StressbalanceSolution.Vx,
+    md.results.StressbalanceSolution.Vy,
+    md.results.StressbalanceSolution.Vz,
+    md.results.StressbalanceSolution.Vel,
+    md.results.StressbalanceSolution.Pressure
+];
Index: /issm/trunk-jpl/test/NightlyRun/test106.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test106.js	(revision 20591)
+++ /issm/trunk-jpl/test/NightlyRun/test106.js	(revision 20591)
@@ -0,0 +1,15 @@
+//Test Name: SquareShelfConstrainedMasstransp2dDG
+var md = new model();
+triangle(md,square[0],150000.);
+meshconvert(md); // meshconvert not implemented yet @TODO
+setmask(md, 'all', '');
+parameterize(md);
+setflowequation(md, 'SSA', 'all');
+md.cluster=generic('name', oshostname(), 'np', 3);
+md.masstransport.stabilization=3;
+md.masstransport.spcthickness=md.geometry.thickness;
+md=solve(md, MasstransportSolutionEnum());
+
+var field_names = ['Thickness'];
+var field_tolerances = [1e-13];
+var field_values = [md.results.MasstransportSolution.Thickness];