Index: /issm/trunk-jpl/test/NightlyRun/test333.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test333.m	(revision 17197)
+++ /issm/trunk-jpl/test/NightlyRun/test333.m	(revision 17198)
@@ -26,16 +26,9 @@
 md=solve(md,HydrologySolutionEnum());
 
-%store=md.constants.g*md.hydrology.sediment_porosity* ...
-%			md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility)
-
-%sed=ones(1,size(md.results.HydrologySolution,2));
-%epl=ones(1,size(md.results.HydrologySolution,2));
-%res=ones(1,size(md.results.HydrologySolution,2));
-%input=ones(1,size(md.results.HydrologySolution,2));
-%for i= 1:size(md.results.HydrologySolution,2)
-%	sed(i)=mean(md.results.HydrologySolution(i).SedimentHead);
-%	res(i)=mean(md.results.HydrologySolution(i).SedimentHeadResidual);
-%	epl(i)=mean(md.results.HydrologySolution(i).EplHead);
-%	input(i)=2.0*(i*0.2);
+%store=md.constants.g*md.hydrology.sediment_porosity*md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility);
+%sedstore=20.0*store;
+%for i=1:10
+%diff=(mean(md.results.HydrologySolution(i).EplHead)*store+ ...
+%			mean(md.results.HydrologySolution(i).SedimentHead)*sedstore-0.4*i)
 %end