Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 15275)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 15276)
@@ -6526,5 +6526,5 @@
 			GetInputListOnVertices(&sed_head[0],SedimentHeadEnum); 
 			GetInputListOnVertices(&epl_head[0],EplHeadEnum);
-
+			
 			this->parameters->FindParam(&leakage,HydrologydcLeakageFactorEnum);
 
@@ -6534,8 +6534,8 @@
 			for(int i=0;i<numdof;i++){
 				if(activeEpl[i]==0.0)continue;
-
+				
 				eplstoring[i]=matpar->GetEplStoring();		
 				sedstoring[i]=matpar->GetSedimentStoring();
-
+				
 				if(sed_head[i]>epl_head[i]){
 					wh_trans=sed_trans*sedstoring[i]*(epl_head[i]-sed_head[i])/(leakage*sed_thick);				
@@ -6547,7 +6547,4 @@
 						wh_trans=0.0;
 					}
-					//					if((sed_head[i]+(wh_trans/sedstoring[i]))>h_max){
-					//(h_max-sed_head[i])*eplstoring[i];
-					//}
 				}
 				/*Assign output pointer*/
Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp	(revision 15275)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp	(revision 15276)
@@ -44,5 +44,5 @@
 	hydrocount=1;
 
-	/*Iteration on the two layers + transfer*/
+	/*Iteration on the two layers*/
 	femmodel->SetCurrentConfiguration(HydrologyDCInefficientAnalysisEnum);
 	GetSolutionFromInputsx(&ug_sed, femmodel->elements, femmodel->nodes, femmodel->vertices,femmodel->loads, femmodel->materials, femmodel->parameters);
@@ -115,5 +115,5 @@
 			/*Iteration on the EPL layer*/
 			for(;;){
-				//				femmodel->HydrologyTransferx();
+				femmodel->HydrologyTransferx();
 				femmodel->SystemMatricesx(&Kff, &Kfs, &pf,&df,NULL);
 				CreateNodalConstraintsx(&ys,femmodel->nodes,HydrologyDCEfficientAnalysisEnum);
@@ -127,5 +127,5 @@
 				Mergesolutionfromftogx(&ug_epl,uf,ys,femmodel->nodes,femmodel->parameters); delete ys;
 				InputUpdateFromSolutionx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,ug_epl);
-				femmodel->HydrologyEPLupdateDomainx();
+				//femmodel->HydrologyEPLupdateDomainx();
 				ConstraintsStatex(&constraints_converged,&num_unstable_constraints, femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);