Index: /issm/trunk/src/m/classes/model.m
===================================================================
--- /issm/trunk/src/m/classes/model.m	(revision 16561)
+++ /issm/trunk/src/m/classes/model.m	(revision 16562)
@@ -680,5 +680,7 @@
 			if ~isnan(md.initialization.temperature),md.initialization.temperature=project3d(md,'vector',md.initialization.temperature,'type','node');end;
 			if ~isnan(md.initialization.waterfraction),md.initialization.waterfraction=project3d(md,'vector',md.initialization.waterfraction,'type','node');end;
-            if ~isnan(md.initialization.watercolumn),md.initialization.watercolumn=project3d(md,'vector',md.initialization.watercolumn,'type','node','layer',1);end;
+      if ~isnan(md.initialization.watercolumn),md.initialization.watercolumn=project3d(md,'vector',md.initialization.watercolumn,'type','node','layer',1);end;
+      if ~isnan(md.initialization.sediment_head),md.initialization.sediment_head=project3d(md,'vector',md.initialization.sediment_head,'type','node','layer',1);end;
+      if ~isnan(md.initialization.epl_head),md.initialization.epl_head=project3d(md,'vector',md.initialization.epl_head,'type','node','layer',1);end;
 
 			%bedinfo and surface info
@@ -715,4 +717,9 @@
 			md.stressbalance.referential=project3d(md,'vector',md.stressbalance.referential,'type','node');
 			md.stressbalance.loadingforce=project3d(md,'vector',md.stressbalance.loadingforce,'type','node');
+			if isa(md.hydrology,'hydrologydc');
+				md.hydrology.spcsediment_head=project3d(md,'vector',md.hydrology.spcsediment_head,'type','node','layer',1);
+				md.hydrology.spcepl_head=project3d(md,'vector',md.hydrology.spcepl_head,'type','node','layer',1);
+				md.hydrology.mask_eplactive=project3d(md,'vector',md.hydrology.mask_eplactive,'type','node','layer',1);
+	    end
 
 			%connectivity
Index: /issm/trunk/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk/test/NightlyRun/IdToName.m	(revision 16561)
+++ /issm/trunk/test/NightlyRun/IdToName.m	(revision 16562)
@@ -104,4 +104,6 @@
 	case 332, name='SquareSheetConstrainedHydrologyDC';
 	case 333, name='SquareSheetHydrologyDCTwoLayers';
+	case 334, name='SquareSheetConstrainedExtrudedHydrologyDC';
+	case 335, name='SquareSheetExtrudedHydrologyDCTwoLayers';
 	case 401, name='SquareSheetShelfStressSSA2d';
 	case 402, name='SquareSheetShelfStressSSA3d';
Index: /issm/trunk/test/NightlyRun/test334.m
===================================================================
--- /issm/trunk/test/NightlyRun/test334.m	(revision 16562)
+++ /issm/trunk/test/NightlyRun/test334.m	(revision 16562)
@@ -0,0 +1,28 @@
+md=triangle(model(),'../Exp/Square.exp',100000.);
+md=setmask(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setflowequation(md,'SSA','all');
+md.cluster=generic('name',oshostname(),'np',1);
+md.hydrology=(hydrologydc);
+md.hydrology.isefficientlayer=0;
+md.hydrology.sedimentlimit_flag=1;
+md.hydrology.sedimentlimit=8000.0;
+md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices,1);
+md.hydrology.spcsediment_head=NaN*ones(md.mesh.numberofvertices,1);
+pos=find(md.mesh.y==0);
+md.hydrology.spcsediment_head(pos)=0.0;
+md.basalforcings.melting_rate = 2.0*ones(md.mesh.numberofvertices,1);
+md.hydrology.sediment_transmitivity=3;
+md.timestepping.time_step=0;
+md.timestepping.final_time=1.0;
+%md.verbose=verbose('1111111');
+md=extrude(md,3,1.1);
+md=solve(md,HydrologySolutionEnum());
+
+%Fields and tolerances to track changes
+%you can also compare with an analitic solution, but it is exact
+%only if no limits are applied
+%analitic=(md.mesh.y.^2-2*md.mesh.y*1.0e6)*(-2.0/(2*md.constants.yts*md.hydrology.sediment_transmitivity))
+field_names     ={'SedimentWaterHead','SedimentHeadResidual'};
+field_tolerances={1e-13, 2.5e-10};
+field_values={md.results.HydrologySolution.SedimentHead,md.results.HydrologySolution.SedimentHeadResidual}
Index: /issm/trunk/test/NightlyRun/test335.m
===================================================================
--- /issm/trunk/test/NightlyRun/test335.m	(revision 16562)
+++ /issm/trunk/test/NightlyRun/test335.m	(revision 16562)
@@ -0,0 +1,62 @@
+md=triangle(model(),'../Exp/Square.exp',100000.);
+md=setmask(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setflowequation(md,'SSA','all');
+md.cluster=generic('name',oshostname(),'np',1);
+md.hydrology=(hydrologydc);
+md.hydrology.isefficientlayer=1;
+md.hydrology.sedimentlimit_flag=1;
+md.hydrology.sedimentlimit=800.0;
+md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices,1);
+md.hydrology.spcsediment_head=NaN*ones(md.mesh.numberofvertices,1);
+md.basalforcings.melting_rate = 2.0*ones(md.mesh.numberofvertices,1);
+md.hydrology.sediment_transmitivity=3;
+
+md.initialization.epl_head=0.0*ones(md.mesh.numberofvertices,1);
+md.hydrology.spcepl_head=NaN*ones(md.mesh.numberofvertices,1);
+md.hydrology.mask_eplactive=0*ones(md.mesh.numberofvertices,1);
+md.basalforcings.melting_rate = 2.0*ones(md.mesh.numberofvertices,1);
+md.hydrology.epl_transmitivity=30;
+md.hydrology.epl_thickness=1;
+md.timestepping.time_step=0.2;
+md.timestepping.final_time=2.0;
+
+%md.verbose.solution=1;
+md=extrude(md,3,1.1);
+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);
+%end
+
+%Fields and tolerances to track changes
+field_names     ={'SedimentWaterHead1','EplWaterHead1','SedimentHeadResidual1',...
+								 'SedimentWaterHead4','EplWaterHead4','SedimentHeadResidual4',...
+								 'SedimentWaterHead5','EplWaterHead5','SedimentHeadResidual5',...
+								 'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9'};
+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};
+field_values={md.results.HydrologySolution(1).SedimentHead, ...
+							md.results.HydrologySolution(1).EplHead,...
+							md.results.HydrologySolution(1).SedimentHeadResidual,...
+							md.results.HydrologySolution(4).SedimentHead,...
+							md.results.HydrologySolution(4).EplHead,...
+							md.results.HydrologySolution(4).SedimentHeadResidual, ...
+							md.results.HydrologySolution(5).SedimentHead,...
+							md.results.HydrologySolution(5).EplHead,...
+							md.results.HydrologySolution(5).SedimentHeadResidual, ...
+							md.results.HydrologySolution(9).SedimentHead,...
+							md.results.HydrologySolution(9).EplHead,...
+							md.results.HydrologySolution(9).SedimentHeadResidual};