Index: /issm/trunk/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk/test/NightlyRun/IdToName.m	(revision 7044)
+++ /issm/trunk/test/NightlyRun/IdToName.m	(revision 7045)
@@ -94,4 +94,8 @@
 	case 249, name='SquareShelfDiagP3dTransientIncrNonHydroSerial';
 	case 250, name='SquareShelfDiagP3dTransientIncrNonHydroParallel';
+	case 251, name='SquareShelfTranCflM2dSerial';
+	case 252, name='SquareShelfTranCflM2dParallel';
+	case 253, name='SquareShelfTranCflP3dSerial';
+	case 254, name='SquareShelfTranCflP3dParallel';
 	case 301, name='SquareSheetConstrainedDiagM2dSerial';
 	case 302, name='SquareSheetConstrainedDiagM2dParallel';
@@ -138,4 +142,8 @@
 	case 343, name='SquareSheetConstrainedCMDragS3dSerial';
 	case 344, name='SquareSheetConstrainedCMDragS3dParallel';
+	case 345, name='SquareSheetConstrainedTranCflM2dSerial';
+	case 346, name='SquareSheetConstrainedTranCflM2dParallel';
+	case 347, name='SquareSheetConstrainedTranCflH3dSerial';
+	case 348, name='SquareSheetConstrainedTranCflH3dParallel';
 	case 401, name='SquareSheetShelfDiagM2dSerial';
 	case 402, name='SquareSheetShelfDiagM2dParallel';
Index: /issm/trunk/test/NightlyRun/test251.m
===================================================================
--- /issm/trunk/test/NightlyRun/test251.m	(revision 7045)
+++ /issm/trunk/test/NightlyRun/test251.m	(revision 7045)
@@ -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.cluster=none;
+md.time_adapt=1;
+md.ndt=10;
+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/test252.m
===================================================================
--- /issm/trunk/test/NightlyRun/test252.m	(revision 7045)
+++ /issm/trunk/test/NightlyRun/test252.m	(revision 7045)
@@ -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=SetParallel(md,3);
+md.time_adapt=1;
+md.ndt=10;
+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/test253.m
===================================================================
--- /issm/trunk/test/NightlyRun/test253.m	(revision 7045)
+++ /issm/trunk/test/NightlyRun/test253.m	(revision 7045)
@@ -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.cluster=none;
+md.time_adapt=1;
+md.ndt=10;
+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/test254.m
===================================================================
--- /issm/trunk/test/NightlyRun/test254.m	(revision 7045)
+++ /issm/trunk/test/NightlyRun/test254.m	(revision 7045)
@@ -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=SetParallel(md,3);
+md.time_adapt=1;
+md.ndt=10;
+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/test345.m
===================================================================
--- /issm/trunk/test/NightlyRun/test345.m	(revision 7045)
+++ /issm/trunk/test/NightlyRun/test345.m	(revision 7045)
@@ -0,0 +1,35 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=none;
+md.time_adapt=1;
+md.ndt=600;
+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/test346.m
===================================================================
--- /issm/trunk/test/NightlyRun/test346.m	(revision 7045)
+++ /issm/trunk/test/NightlyRun/test346.m	(revision 7045)
@@ -0,0 +1,35 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md=SetParallel(md,3);
+md.time_adapt=1;
+md.ndt=600;
+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/test347.m
===================================================================
--- /issm/trunk/test/NightlyRun/test347.m	(revision 7045)
+++ /issm/trunk/test/NightlyRun/test347.m	(revision 7045)
@@ -0,0 +1,47 @@
+md=mesh(model,'../Exp/Square.exp',200000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,4,1.2);
+md=setelementstype(md,'hutter','all');
+md.cluster=none;
+md.time_adapt=1;
+md.ndt=500;
+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-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,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
+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/test348.m
===================================================================
--- /issm/trunk/test/NightlyRun/test348.m	(revision 7045)
+++ /issm/trunk/test/NightlyRun/test348.m	(revision 7045)
@@ -0,0 +1,49 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,5,1.2);
+md=setelementstype(md,'hutter','all');
+md=SetParallel(md,3);
+md.time_adapt=1;
+md.ndt=500;
+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-13,  1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...,
+						5*1e-11,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
+						1e-10,  1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
+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),...
+	};
