Index: /issm/trunk/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk/test/NightlyRun/IdToName.m	(revision 6222)
+++ /issm/trunk/test/NightlyRun/IdToName.m	(revision 6223)
@@ -193,4 +193,6 @@
 elseif (id==526), name='PigCMDragSteaM3dparallel';
 elseif (id==527), name='PigBamgMesh';
+elseif (id==529), name='PigTherTranStabilizedSerial';
+elseif (id==530), name='PigTherTranStabilizedParallel';
 elseif (id==601), name='79NorthProg2dSerial';
 elseif (id==602), name='79NorthProg2dParallel';
Index: /issm/trunk/test/NightlyRun/test529.m
===================================================================
--- /issm/trunk/test/NightlyRun/test529.m	(revision 6223)
+++ /issm/trunk/test/NightlyRun/test529.m	(revision 6223)
@@ -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=solve(md,'analysis_type',ThermalSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Temperature1','MeltingRate1', ...
+				      'Temperature2','MeltingRate2'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.ThermalSolution(1).Temperature),...
+	PatchToVec(md.results.ThermalSolution(1).MeltingRate),...
+	PatchToVec(md.results.ThermalSolution(2).Temperature),...
+	PatchToVec(md.results.ThermalSolution(2).MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test530.m
===================================================================
--- /issm/trunk/test/NightlyRun/test530.m	(revision 6223)
+++ /issm/trunk/test/NightlyRun/test530.m	(revision 6223)
@@ -0,0 +1,19 @@
+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=solve(md,'analysis_type',ThermalSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Temperature1','MeltingRate1', ...
+				      'Temperature2','MeltingRate2'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.ThermalSolution(1).Temperature),...
+	PatchToVec(md.results.ThermalSolution(1).MeltingRate),...
+	PatchToVec(md.results.ThermalSolution(2).Temperature),...
+	PatchToVec(md.results.ThermalSolution(2).MeltingRate),...
+	};
