Index: /issm/trunk/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk/test/NightlyRun/IdToName.m	(revision 7022)
+++ /issm/trunk/test/NightlyRun/IdToName.m	(revision 7023)
@@ -86,4 +86,12 @@
 	case 241, name='SquareShelfDiagMS3dTilingSerial';
 	case 242, name='SquareShelfDiagMS3dTilingParallel';
+	case 243, name='SquareShelfDiagM2dTransientIncrHydroSerial';
+	case 244, name='SquareShelfDiagM2dTransientIncrHydroParallel';
+	case 245, name='SquareShelfDiagM2dTransientIncrNonHydroSerial';
+	case 246, name='SquareShelfDiagM2dTransientIncrNonHydroParallel';
+	case 247, name='SquareShelfDiagP3dTransientIncrHydroSerial';
+	case 248, name='SquareShelfDiagP3dTransientIncrHydroParallel';
+	case 249, name='SquareShelfDiagP3dTransientIncrNonHydroSerial';
+	case 250, name='SquareShelfDiagP3dTransientIncrNonHydroParallel';
 	case 301, name='SquareSheetConstrainedDiagM2dSerial';
 	case 302, name='SquareSheetConstrainedDiagM2dParallel';
Index: /issm/trunk/test/NightlyRun/test243.m
===================================================================
--- /issm/trunk/test/NightlyRun/test243.m	(revision 7023)
+++ /issm/trunk/test/NightlyRun/test243.m	(revision 7023)
@@ -0,0 +1,34 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=none;
+md.hydrostatic_adjustment=IncrementalEnum;
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3'};
+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,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.Transient2DSolution(1).Vx),...
+	PatchToVec(md.results.Transient2DSolution(1).Vy),...
+	PatchToVec(md.results.Transient2DSolution(1).Vel),...
+	PatchToVec(md.results.Transient2DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(1).Bed),...
+	PatchToVec(md.results.Transient2DSolution(1).Surface),...
+	PatchToVec(md.results.Transient2DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient2DSolution(2).Vx),...
+	PatchToVec(md.results.Transient2DSolution(2).Vy),...
+	PatchToVec(md.results.Transient2DSolution(2).Vel),...
+	PatchToVec(md.results.Transient2DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(2).Bed),...
+	PatchToVec(md.results.Transient2DSolution(2).Surface),...
+	PatchToVec(md.results.Transient2DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient2DSolution(3).Vx),...
+	PatchToVec(md.results.Transient2DSolution(3).Vy),...
+	PatchToVec(md.results.Transient2DSolution(3).Vel),...
+	PatchToVec(md.results.Transient2DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(3).Bed),...
+	PatchToVec(md.results.Transient2DSolution(3).Surface),...
+	PatchToVec(md.results.Transient2DSolution(3).Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test244.m
===================================================================
--- /issm/trunk/test/NightlyRun/test244.m	(revision 7023)
+++ /issm/trunk/test/NightlyRun/test244.m	(revision 7023)
@@ -0,0 +1,34 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=setelementstype(md,'macayeal','all');
+md=SetParallel(md,3);
+md.hydrostatic_adjustment=IncrementalEnum;
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3'};
+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,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.Transient2DSolution(1).Vx),...
+	PatchToVec(md.results.Transient2DSolution(1).Vy),...
+	PatchToVec(md.results.Transient2DSolution(1).Vel),...
+	PatchToVec(md.results.Transient2DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(1).Bed),...
+	PatchToVec(md.results.Transient2DSolution(1).Surface),...
+	PatchToVec(md.results.Transient2DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient2DSolution(2).Vx),...
+	PatchToVec(md.results.Transient2DSolution(2).Vy),...
+	PatchToVec(md.results.Transient2DSolution(2).Vel),...
+	PatchToVec(md.results.Transient2DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(2).Bed),...
+	PatchToVec(md.results.Transient2DSolution(2).Surface),...
+	PatchToVec(md.results.Transient2DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient2DSolution(3).Vx),...
+	PatchToVec(md.results.Transient2DSolution(3).Vy),...
+	PatchToVec(md.results.Transient2DSolution(3).Vel),...
+	PatchToVec(md.results.Transient2DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(3).Bed),...
+	PatchToVec(md.results.Transient2DSolution(3).Surface),...
+	PatchToVec(md.results.Transient2DSolution(3).Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test245.m
===================================================================
--- /issm/trunk/test/NightlyRun/test245.m	(revision 7023)
+++ /issm/trunk/test/NightlyRun/test245.m	(revision 7023)
@@ -0,0 +1,35 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=setelementstype(md,'macayeal','all');
+md.bed=md.bed+50; md.surface=md.surface+50;
+md.cluster=none;
+md.hydrostatic_adjustment=IncrementalEnum;
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3'};
+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,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.Transient2DSolution(1).Vx),...
+	PatchToVec(md.results.Transient2DSolution(1).Vy),...
+	PatchToVec(md.results.Transient2DSolution(1).Vel),...
+	PatchToVec(md.results.Transient2DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(1).Bed),...
+	PatchToVec(md.results.Transient2DSolution(1).Surface),...
+	PatchToVec(md.results.Transient2DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient2DSolution(2).Vx),...
+	PatchToVec(md.results.Transient2DSolution(2).Vy),...
+	PatchToVec(md.results.Transient2DSolution(2).Vel),...
+	PatchToVec(md.results.Transient2DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(2).Bed),...
+	PatchToVec(md.results.Transient2DSolution(2).Surface),...
+	PatchToVec(md.results.Transient2DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient2DSolution(3).Vx),...
+	PatchToVec(md.results.Transient2DSolution(3).Vy),...
+	PatchToVec(md.results.Transient2DSolution(3).Vel),...
+	PatchToVec(md.results.Transient2DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(3).Bed),...
+	PatchToVec(md.results.Transient2DSolution(3).Surface),...
+	PatchToVec(md.results.Transient2DSolution(3).Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test246.m
===================================================================
--- /issm/trunk/test/NightlyRun/test246.m	(revision 7023)
+++ /issm/trunk/test/NightlyRun/test246.m	(revision 7023)
@@ -0,0 +1,35 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=setelementstype(md,'macayeal','all');
+md.bed=md.bed+50; md.surface=md.surface+50;
+md=SetParallel(md,3);
+md.hydrostatic_adjustment=IncrementalEnum;
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3'};
+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,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.Transient2DSolution(1).Vx),...
+	PatchToVec(md.results.Transient2DSolution(1).Vy),...
+	PatchToVec(md.results.Transient2DSolution(1).Vel),...
+	PatchToVec(md.results.Transient2DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(1).Bed),...
+	PatchToVec(md.results.Transient2DSolution(1).Surface),...
+	PatchToVec(md.results.Transient2DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient2DSolution(2).Vx),...
+	PatchToVec(md.results.Transient2DSolution(2).Vy),...
+	PatchToVec(md.results.Transient2DSolution(2).Vel),...
+	PatchToVec(md.results.Transient2DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(2).Bed),...
+	PatchToVec(md.results.Transient2DSolution(2).Surface),...
+	PatchToVec(md.results.Transient2DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient2DSolution(3).Vx),...
+	PatchToVec(md.results.Transient2DSolution(3).Vy),...
+	PatchToVec(md.results.Transient2DSolution(3).Vel),...
+	PatchToVec(md.results.Transient2DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(3).Bed),...
+	PatchToVec(md.results.Transient2DSolution(3).Surface),...
+	PatchToVec(md.results.Transient2DSolution(3).Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test247.m
===================================================================
--- /issm/trunk/test/NightlyRun/test247.m	(revision 7023)
+++ /issm/trunk/test/NightlyRun/test247.m	(revision 7023)
@@ -0,0 +1,46 @@
+md=mesh(model,'../Exp/Square.exp',200000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster=none;
+md.hydrostatic_adjustment=IncrementalEnum;
+md=solve(md,'analysis_type',Transient3DSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','MeltingRate3'};
+field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06};
+field_values={...
+	PatchToVec(md.results.Transient3DSolution(1).Vx),...
+	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vz),...
+	PatchToVec(md.results.Transient3DSolution(1).Vel),...
+	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(1).Bed),...
+	PatchToVec(md.results.Transient3DSolution(1).Surface),...
+	PatchToVec(md.results.Transient3DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(1).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(1).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(2).Vx),...
+	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
+	PatchToVec(md.results.Transient3DSolution(2).Vel),...
+	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(2).Bed),...
+	PatchToVec(md.results.Transient3DSolution(2).Surface),...
+	PatchToVec(md.results.Transient3DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(2).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(2).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(3).Vx),...
+	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
+	PatchToVec(md.results.Transient3DSolution(3).Vel),...
+	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(3).Bed),...
+	PatchToVec(md.results.Transient3DSolution(3).Surface),...
+	PatchToVec(md.results.Transient3DSolution(3).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(3).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(3).MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test248.m
===================================================================
--- /issm/trunk/test/NightlyRun/test248.m	(revision 7023)
+++ /issm/trunk/test/NightlyRun/test248.m	(revision 7023)
@@ -0,0 +1,49 @@
+md=mesh(model,'../Exp/Square.exp',200000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all');
+md=SetParallel(md,3);
+md.hydrostatic_adjustment=IncrementalEnum;
+md=solve(md,'analysis_type',Transient3DSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','MeltingRate3'};
+field_tolerances={...
+	1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,...
+	1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,...
+	1e-09,1e-09,1e-08,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05};
+field_values={...
+	PatchToVec(md.results.Transient3DSolution(1).Vx),...
+	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vz),...
+	PatchToVec(md.results.Transient3DSolution(1).Vel),...
+	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(1).Bed),...
+	PatchToVec(md.results.Transient3DSolution(1).Surface),...
+	PatchToVec(md.results.Transient3DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(1).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(1).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(2).Vx),...
+	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
+	PatchToVec(md.results.Transient3DSolution(2).Vel),...
+	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(2).Bed),...
+	PatchToVec(md.results.Transient3DSolution(2).Surface),...
+	PatchToVec(md.results.Transient3DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(2).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(2).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(3).Vx),...
+	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
+	PatchToVec(md.results.Transient3DSolution(3).Vel),...
+	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(3).Bed),...
+	PatchToVec(md.results.Transient3DSolution(3).Surface),...
+	PatchToVec(md.results.Transient3DSolution(3).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(3).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(3).MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test249.m
===================================================================
--- /issm/trunk/test/NightlyRun/test249.m	(revision 7023)
+++ /issm/trunk/test/NightlyRun/test249.m	(revision 7023)
@@ -0,0 +1,47 @@
+md=mesh(model,'../Exp/Square.exp',200000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all');
+md.bed=md.bed+50; md.surface=md.surface+50;
+md.cluster=none;
+md.hydrostatic_adjustment=IncrementalEnum;
+md=solve(md,'analysis_type',Transient3DSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','MeltingRate3'};
+field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06};
+field_values={...
+	PatchToVec(md.results.Transient3DSolution(1).Vx),...
+	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vz),...
+	PatchToVec(md.results.Transient3DSolution(1).Vel),...
+	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(1).Bed),...
+	PatchToVec(md.results.Transient3DSolution(1).Surface),...
+	PatchToVec(md.results.Transient3DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(1).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(1).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(2).Vx),...
+	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
+	PatchToVec(md.results.Transient3DSolution(2).Vel),...
+	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(2).Bed),...
+	PatchToVec(md.results.Transient3DSolution(2).Surface),...
+	PatchToVec(md.results.Transient3DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(2).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(2).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(3).Vx),...
+	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
+	PatchToVec(md.results.Transient3DSolution(3).Vel),...
+	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(3).Bed),...
+	PatchToVec(md.results.Transient3DSolution(3).Surface),...
+	PatchToVec(md.results.Transient3DSolution(3).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(3).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(3).MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test250.m
===================================================================
--- /issm/trunk/test/NightlyRun/test250.m	(revision 7023)
+++ /issm/trunk/test/NightlyRun/test250.m	(revision 7023)
@@ -0,0 +1,50 @@
+md=mesh(model,'../Exp/Square.exp',200000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all');
+md.bed=md.bed+50; md.surface=md.surface+50;
+md=SetParallel(md,3);
+md.hydrostatic_adjustment=IncrementalEnum;
+md=solve(md,'analysis_type',Transient3DSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','MeltingRate3'};
+field_tolerances={...
+	1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,...
+	1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,...
+	1e-09,1e-09,1e-08,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05};
+field_values={...
+	PatchToVec(md.results.Transient3DSolution(1).Vx),...
+	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vz),...
+	PatchToVec(md.results.Transient3DSolution(1).Vel),...
+	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(1).Bed),...
+	PatchToVec(md.results.Transient3DSolution(1).Surface),...
+	PatchToVec(md.results.Transient3DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(1).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(1).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(2).Vx),...
+	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
+	PatchToVec(md.results.Transient3DSolution(2).Vel),...
+	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(2).Bed),...
+	PatchToVec(md.results.Transient3DSolution(2).Surface),...
+	PatchToVec(md.results.Transient3DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(2).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(2).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(3).Vx),...
+	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
+	PatchToVec(md.results.Transient3DSolution(3).Vel),...
+	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(3).Bed),...
+	PatchToVec(md.results.Transient3DSolution(3).Surface),...
+	PatchToVec(md.results.Transient3DSolution(3).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(3).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(3).MeltingRate),...
+	};
