Index: /issm/trunk-jpl/test/NightlyRun/test333.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test333.m	(revision 22286)
+++ /issm/trunk-jpl/test/NightlyRun/test333.m	(revision 22287)
@@ -34,4 +34,10 @@
 
 md=solve(md,'Transient');
+%re-run with no substeps
+mdfine=md;
+mdfine.result=struct();
+mdfine.hydrology.steps_per_step=1;
+mdfine.timestepping.time_step=0.2;
+mdfine=solve(md,'Transient');
 
 %store=md.constants.g*md.hydrology.sediment_porosity*md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility);
@@ -43,22 +49,25 @@
 %Fields and tolerances to track changes
 field_names     ={'SedimentWaterHead1','EplWaterHead1','SedimentHeadResidual1',...
-						'SedimentWaterHead4','EplWaterHead4','SedimentHeadResidual4',...
-						'SedimentWaterHead5','EplWaterHead5','SedimentHeadResidual5',...
-						'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9'};
+                  'SedimentWaterHead4','EplWaterHead4','SedimentHeadResidual4',...
+                  'SedimentWaterHead5','EplWaterHead5','SedimentHeadResidual5',...
+                  'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9',...
+                  'EplWaterHead9'};
 field_tolerances={...
 	1e-13, 1e-13, 1e-13,...
 	1e-13, 1e-13, 1e-13,...
 	1e-13, 5e-12, 1e-11,...
-	1e-13, 5e-12, 1e-11};
-field_values={md.results.TransientSolution(1).SedimentHead, ...
-							md.results.TransientSolution(1).EplHead,...
-							md.results.TransientSolution(1).SedimentHeadResidual,...
-							md.results.TransientSolution(4).SedimentHead,...
-							md.results.TransientSolution(4).EplHead,...
-							md.results.TransientSolution(4).SedimentHeadResidual, ...
-							md.results.TransientSolution(5).SedimentHead,...
-							md.results.TransientSolution(5).EplHead,...
-							md.results.TransientSolution(5).SedimentHeadResidual, ...
-							md.results.TransientSolution(9).SedimentHead,...
-							md.results.TransientSolution(9).EplHead,...
-							md.results.TransientSolution(9).SedimentHeadResidual};
+	1e-13, 5e-12, 1e-11,...
+  1e-13};
+field_values={mdfine.results.TransientSolution(1).SedimentHead, ...
+							mdfine.results.TransientSolution(1).EplHead,...
+							mdfine.results.TransientSolution(1).SedimentHeadResidual,...
+							mdfine.results.TransientSolution(4).SedimentHead,...
+							mdfine.results.TransientSolution(4).EplHead,...
+							mdfine.results.TransientSolution(4).SedimentHeadResidual, ...
+							mdfine.results.TransientSolution(5).SedimentHead,...
+							mdfine.results.TransientSolution(5).EplHead,...
+							mdfine.results.TransientSolution(5).SedimentHeadResidual, ...
+							mdfine.results.TransientSolution(9).SedimentHead,...
+							mdfine.results.TransientSolution(9).EplHead,...
+							mdfine.results.TransientSolution(9).SedimentHeadResidual,...
+              md.results.TransientSolution(0).EplHead};
Index: /issm/trunk-jpl/test/NightlyRun/test333.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test333.py	(revision 22286)
+++ /issm/trunk-jpl/test/NightlyRun/test333.py	(revision 22287)
@@ -42,8 +42,17 @@
 md.hydrology.epl_thick_comp=1
 md.hydrology.epl_max_thickness=1
-md.timestepping.time_step=0.2
+md.hydrology.steps_per_step=10
+md.timestepping.time_step=2.0
 md.timestepping.final_time=2.0
 
 md=solve(md,'Transient')
+
+#re-run with no substeps
+mdfine=md
+mdfine.result=[]
+mdfine.hydrology.steps_per_step=1
+mdfine.timestepping.time_step=0.2
+mdfine=solve(md,'Transient')
+
 
 # sedvol=np.zeros(10)
@@ -59,19 +68,23 @@
 						 'SedimentWaterHead4','EplWaterHead4','SedimentHeadResidual4',
 						 'SedimentWaterHead5','EplWaterHead5','SedimentHeadResidual5',
-						 'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9']
+						 'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9', 
+						 'EplWaterHead9']
 field_tolerances=[1e-13, 1e-13, 1e-13,
 									1e-13, 1e-13, 1e-13,
 									1e-13, 5e-12, 1e-11,
-									1e-13, 5e-12, 1e-11]
-field_values=[md.results.TransientSolution[0].SedimentHead,
-							md.results.TransientSolution[0].EplHead,
-							md.results.TransientSolution[0].SedimentHeadResidual,
-							md.results.TransientSolution[3].SedimentHead,
-							md.results.TransientSolution[3].EplHead,
-							md.results.TransientSolution[3].SedimentHeadResidual,
-							md.results.TransientSolution[4].SedimentHead,
-							md.results.TransientSolution[4].EplHead,
-							md.results.TransientSolution[4].SedimentHeadResidual,
-							md.results.TransientSolution[8].SedimentHead,
-							md.results.TransientSolution[8].EplHead,
-							md.results.TransientSolution[8].SedimentHeadResidual]
+									1e-13, 5e-12, 1e-11,
+									1e-13]
+field_values=[mdfine.results.TransientSolution[0].SedimentHead,
+							mdfine.results.TransientSolution[0].EplHead,
+							mdfine.results.TransientSolution[0].SedimentHeadResidual,
+							mdfine.results.TransientSolution[3].SedimentHead,
+							mdfine.results.TransientSolution[3].EplHead,
+							mdfine.results.TransientSolution[3].SedimentHeadResidual,
+							mdfine.results.TransientSolution[4].SedimentHead,
+							mdfine.results.TransientSolution[4].EplHead,
+							mdfine.results.TransientSolution[4].SedimentHeadResidual,
+							mdfine.results.TransientSolution[8].SedimentHead,
+							mdfine.results.TransientSolution[8].EplHead,
+							mdfine.results.TransientSolution[8].SedimentHeadResidual,
+							md.results.TransientSolution[0].EplHead]
+