Index: /issm/trunk-jpl/src/c/analyses/hydrology_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/hydrology_core.cpp	(revision 15005)
+++ /issm/trunk-jpl/src/c/analyses/hydrology_core.cpp	(revision 15006)
@@ -63,4 +63,5 @@
 
 			/*transfer water column thickness to old water column thickness: */
+			
 			InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,WatercolumnEnum,WaterColumnOldEnum);
 
@@ -78,5 +79,7 @@
 
 		else if (hydrology_model==HydrologydcEnum){
+				InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,SedimentHeadEnum,SedimentHeadOldEnum);
 			femmodel->parameters->FindParam(&isefficientlayer,HydrologydcIsefficientlayerEnum);
+			
 			if(VerboseSolution()) _pprintLine_("   computing water head");
 			solutionsequence_hydro_nonlinear(femmodel);
Index: /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp	(revision 15005)
+++ /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp	(revision 15006)
@@ -2075,5 +2075,7 @@
 				name==GradientEnum ||
 				name==OldGradientEnum ||
-	         name==ConvergedEnum ||
+				name==ConvergedEnum ||
+				name==SedimentHeadOldEnum ||
+				name==SedimentHeadEnum ||
 				name==BasisIntegralEnum ||
 				name==QmuVxEnum ||
@@ -6023,5 +6025,4 @@
 		}
 	}
-
 	/*Clean up and return*/
 	delete gauss;
@@ -6165,9 +6166,10 @@
 
 	if(reCast<bool,IssmDouble>(dt)){
-		old_wh_input=inputs->GetInput(SedimentHeadEnum); _assert_(old_wh_input);
+		old_wh_input=inputs->GetInput(SedimentHeadOldEnum); _assert_(old_wh_input);
 	}
 
 	/* Start  looping on the number of gaussian points: */
 	gauss=new GaussTria(2);
+	
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 
@@ -6190,5 +6192,4 @@
 		}
 	}
-
 	/*Clean up and return*/
 	delete gauss;
@@ -6340,5 +6341,5 @@
 	this->inputs->GetInputValue(&converged,ConvergedEnum);
 	GetInputListOnVertices(&intbasis[0],BasisIntegralEnum);
-	
+
 	if(converged){
 		this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
@@ -6357,6 +6358,8 @@
 
 	/*Add input to the element: */
+
 	this->inputs->AddInput(new TriaP1Input(SedimentHeadEnum,values));
 	this->inputs->AddInput(new TriaP1Input(SedimentHeadResidualEnum,residual));
+	if(converged)this->inputs->AddInput(new TriaP1Input(SedimentHeadOldEnum,values));
 
 	/*Free ressources:*/
Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCInefficient/UpdateElementsHydrologyDCInefficient.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCInefficient/UpdateElementsHydrologyDCInefficient.cpp	(revision 15005)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCInefficient/UpdateElementsHydrologyDCInefficient.cpp	(revision 15006)
@@ -43,4 +43,5 @@
 	iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum);
 	iomodel->FetchDataToInput(elements,SedimentHeadEnum);
+	/*	iomodel->FetchDataToInput(elements,SedimentHeadOldEnum);*/
 
 	/*Free data: */
Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp	(revision 15005)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp	(revision 15006)
@@ -15,5 +15,5 @@
 	Vector<IssmDouble>* uf_old=NULL; 
 	Vector<IssmDouble>* ys=NULL; 
-	IssmDouble sediment_kmax;
+	IssmDouble sediment_kmax,time;
 
 	/*intermediary: */
@@ -37,13 +37,6 @@
 	converged=false;
 
+	femmodel->parameters->FindParam(&time,TimeEnum);
 	for(;;){
-		
-		/*Computing the transfer term
-			
-
-		 */
-
-		
-
 		/*First layer*/
 		femmodel->SetCurrentConfiguration(HydrologyDCInefficientAnalysisEnum);
@@ -53,4 +46,5 @@
 		femmodel->parameters->SetParam(HydrologySedimentEnum,HydrologyLayerEnum);
 		for(;;){
+
 			femmodel->SystemMatricesx(&Kff, &Kfs, &pf,&df, &sediment_kmax);
 			CreateNodalConstraintsx(&ys,femmodel->nodes,HydrologyDCInefficientAnalysisEnum);
@@ -75,5 +69,4 @@
 			if(converged){
 				femmodel->parameters->SetParam(sediment_kmax,HydrologySedimentKmaxEnum);
-
 				InputUpdateFromConstantx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,converged,ConvergedEnum);
 				InputUpdateFromSolutionx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,ug);