Index: /issm/trunk/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk/test/NightlyRun/IdToName.m	(revision 6265)
+++ /issm/trunk/test/NightlyRun/IdToName.m	(revision 6266)
@@ -197,4 +197,6 @@
 	case 529, name='PigTherTranStabilizedSerial';
 	case 530, name='PigTherTranStabilizedParallel';
+	case 529, name='PigTherSteaStabilizedSerial';
+	case 530, name='PigTherSteaStabilizedParallel';
 	case 601, name='79NorthProg2dSerial';
 	case 602, name='79NorthProg2dParallel';
Index: /issm/trunk/test/NightlyRun/test531.m
===================================================================
--- /issm/trunk/test/NightlyRun/test531.m	(revision 6266)
+++ /issm/trunk/test/NightlyRun/test531.m	(revision 6266)
@@ -0,0 +1,17 @@
+md=mesh(model,'../Exp/Pig.exp',100000);
+md=geography(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
+md=parameterize(md,'../Par/Pig.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all');
+md.artificial_diffusivity=2;
+md.dt=0;
+md.min_thermal_constraints=40;
+md=solve(md,'analysis_type',ThermalSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.ThermalSolution.Temperature),...
+	PatchToVec(md.results.ThermalSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test532.m
===================================================================
--- /issm/trunk/test/NightlyRun/test532.m	(revision 6266)
+++ /issm/trunk/test/NightlyRun/test532.m	(revision 6266)
@@ -0,0 +1,18 @@
+md=mesh(model,'../Exp/Pig.exp',100000);
+md=geography(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
+md=parameterize(md,'../Par/Pig.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all');
+md.artificial_diffusivity=2;
+md=SetParallel(md,3);
+md.dt=0;
+md.min_thermal_constraints=40;
+md=solve(md,'analysis_type',ThermalSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.ThermalSolution.Temperature),...
+	PatchToVec(md.results.ThermalSolution.MeltingRate),...
+	};