Index: /issm/trunk-jpl/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 11280)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 11281)
@@ -44,4 +44,8 @@
 	case 136, name='SquareShelfConstrainedBalThic2dDGParallel';
 	case 137, name='SquareBamgMesh';
+	case 139, name='SquareShelfConstrainedEnthalpySteaSerial';
+	case 140, name='SquareShelfConstrainedEnthalpySteaParallel';
+	case 141, name='SquareShelfConstrainedEnthalpyTranSerial';
+	case 142, name='SquareShelfConstrainedEnthalpyTranParallel';
 	case 201, name='SquareShelfDiagM2dSerial';
 	case 202, name='SquareShelfDiagM2dParallel';
@@ -160,4 +164,8 @@
 	case 347, name='SquareSheetConstrainedTranCflH3dSerial';
 	case 348, name='SquareSheetConstrainedTranCflH3dParallel';
+	case 349, name='SquareSheetConstrainedEnthalpySteaSerial';
+	case 350, name='SquareSheetConstrainedEnthalpySteaParallel';
+	case 351, name='SquareSheetConstrainedEnthalpyTranSerial';
+	case 352, name='SquareSheetConstrainedEnthalpyTranParallel';
 	case 401, name='SquareSheetShelfDiagM2dSerial';
 	case 402, name='SquareSheetShelfDiagM2dParallel';
Index: /issm/trunk-jpl/test/NightlyRun/test139.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test139.m	(revision 11281)
+++ /issm/trunk-jpl/test/NightlyRun/test139.m	(revision 11281)
@@ -0,0 +1,18 @@
+md=triangle(model,'../Exp/Square.exp',180000);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,1);
+md=setflowequation(md,'macayeal','all');
+md.timestepping.time_step=0;
+md.cluster=none;
+md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+md=solve(md,EnthalpySolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Enthalpy','Waterfraction','Temperature'};
+field_tolerances={1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.EnthalpySolution.Enthalpy),...
+	(md.results.EnthalpySolution.Waterfraction),...
+	(md.results.EnthalpySolution.Temperature),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test140.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test140.m	(revision 11281)
+++ /issm/trunk-jpl/test/NightlyRun/test140.m	(revision 11281)
@@ -0,0 +1,19 @@
+md=triangle(model,'../Exp/Square.exp',180000);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,1);
+md=setflowequation(md,'macayeal','all');
+md.timestepping.time_step=0;
+md.cluster=generic('name',oshostname(),'np',3);
+md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+%md=solve(md,ThermalSolutionEnum);
+md=solve(md,EnthalpySolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Enthalpy','Waterfraction','Temperature'};
+field_tolerances={1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.EnthalpySolution.Enthalpy),...
+	(md.results.EnthalpySolution.Waterfraction),...
+	(md.results.EnthalpySolution.Temperature),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test141.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test141.m	(revision 11281)
+++ /issm/trunk-jpl/test/NightlyRun/test141.m	(revision 11281)
@@ -0,0 +1,25 @@
+md=triangle(model,'../Exp/Square.exp',180000);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,1);
+md=setflowequation(md,'macayeal','all');
+md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+md.cluster=none;
+md=solve(md,EnthalpySolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Enthalpy1','Waterfraction1','Temperature1',...
+	'Enthalpy2','Waterfraction2','Temperature2',...
+	'Enthalpy3','Waterfraction3','Temperature3'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.EnthalpySolution(1).Enthalpy),...
+	(md.results.EnthalpySolution(1).Waterfraction),...
+	(md.results.EnthalpySolution(1).Temperature),...
+	(md.results.EnthalpySolution(2).Enthalpy),...
+	(md.results.EnthalpySolution(2).Waterfraction),...
+	(md.results.EnthalpySolution(2).Temperature),...
+	(md.results.EnthalpySolution(3).Enthalpy),...
+	(md.results.EnthalpySolution(3).Waterfraction),...
+	(md.results.EnthalpySolution(3).Temperature),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test142.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test142.m	(revision 11281)
+++ /issm/trunk-jpl/test/NightlyRun/test142.m	(revision 11281)
@@ -0,0 +1,26 @@
+md=triangle(model,'../Exp/Square.exp',180000);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,1);
+md=setflowequation(md,'macayeal','all');
+md.cluster=generic('name',oshostname(),'np',3);
+md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+%md=solve(md,ThermalSolutionEnum);
+md=solve(md,EnthalpySolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Enthalpy1','Waterfraction1','Temperature1',...
+	'Enthalpy2','Waterfraction2','Temperature2',...
+	'Enthalpy3','Waterfraction3','Temperature3'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.EnthalpySolution(1).Enthalpy),...
+	(md.results.EnthalpySolution(1).Waterfraction),...
+	(md.results.EnthalpySolution(1).Temperature),...
+	(md.results.EnthalpySolution(2).Enthalpy),...
+	(md.results.EnthalpySolution(2).Waterfraction),...
+	(md.results.EnthalpySolution(2).Temperature),...
+	(md.results.EnthalpySolution(3).Enthalpy),...
+	(md.results.EnthalpySolution(3).Waterfraction),...
+	(md.results.EnthalpySolution(3).Temperature),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test349.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test349.m	(revision 11281)
+++ /issm/trunk-jpl/test/NightlyRun/test349.m	(revision 11281)
@@ -0,0 +1,18 @@
+md=triangle(model,'../Exp/Square.exp',180000);
+md=setmask(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,3,1);
+md=setflowequation(md,'macayeal','all');
+md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+md.timestepping.time_step=0;
+md.cluster=none;
+md=solve(md,EnthalpySolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Enthalpy','Waterfraction','Temperature'};
+field_tolerances={1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.EnthalpySolution.Enthalpy),...
+	(md.results.EnthalpySolution.Waterfraction),...
+	(md.results.EnthalpySolution.Temperature),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test350.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test350.m	(revision 11281)
+++ /issm/trunk-jpl/test/NightlyRun/test350.m	(revision 11281)
@@ -0,0 +1,20 @@
+md=triangle(model,'../Exp/Square.exp',180000);
+md=setmask(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,3,1);
+md=setflowequation(md,'macayeal','all');
+md.timestepping.time_step=0;
+md.cluster=generic('name',oshostname(),'np',3);
+md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+%md.debug.valgrind=1;
+md=solve(md,EnthalpySolutionEnum);
+%md=solve(md,ThermalSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Enthalpy','Waterfraction','Temperature'};
+field_tolerances={1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.EnthalpySolution.Enthalpy),...
+	(md.results.EnthalpySolution.Waterfraction),...
+	(md.results.EnthalpySolution.Temperature),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test351.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test351.m	(revision 11281)
+++ /issm/trunk-jpl/test/NightlyRun/test351.m	(revision 11281)
@@ -0,0 +1,25 @@
+md=triangle(model,'../Exp/Square.exp',180000);
+md=setmask(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,3,1);
+md=setflowequation(md,'macayeal','all');
+md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+md.cluster=none;
+md=solve(md,EnthalpySolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Enthalpy1','Waterfraction1','Temperature1',...
+	'Enthalpy2','Waterfraction2','Temperature2',...
+	'Enthalpy3','Waterfraction3','Temperature3'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.EnthalpySolution(1).Enthalpy),...
+	(md.results.EnthalpySolution(1).Waterfraction),...
+	(md.results.EnthalpySolution(1).Temperature),...
+	(md.results.EnthalpySolution(2).Enthalpy),...
+	(md.results.EnthalpySolution(2).Waterfraction),...
+	(md.results.EnthalpySolution(2).Temperature),...
+	(md.results.EnthalpySolution(3).Enthalpy),...
+	(md.results.EnthalpySolution(3).Waterfraction),...
+	(md.results.EnthalpySolution(3).Temperature),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test352.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test352.m	(revision 11281)
+++ /issm/trunk-jpl/test/NightlyRun/test352.m	(revision 11281)
@@ -0,0 +1,26 @@
+md=triangle(model,'../Exp/Square.exp',180000);
+md=setmask(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,3,1);
+md=setflowequation(md,'macayeal','all');
+md.cluster=generic('name',oshostname(),'np',3);
+md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+%md=solve(md,ThermalSolutionEnum);
+md=solve(md,EnthalpySolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Enthalpy1','Waterfraction1','Temperature1',...
+	'Enthalpy2','Waterfraction2','Temperature2',...
+	'Enthalpy3','Waterfraction3','Temperature3'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.EnthalpySolution(1).Enthalpy),...
+	(md.results.EnthalpySolution(1).Waterfraction),...
+	(md.results.EnthalpySolution(1).Temperature),...
+	(md.results.EnthalpySolution(2).Enthalpy),...
+	(md.results.EnthalpySolution(2).Waterfraction),...
+	(md.results.EnthalpySolution(2).Temperature),...
+	(md.results.EnthalpySolution(3).Enthalpy),...
+	(md.results.EnthalpySolution(3).Waterfraction),...
+	(md.results.EnthalpySolution(3).Temperature),...
+	};