Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16899)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16900)
@@ -7167,5 +7167,4 @@
 					/*Assign output pointer*/
 					transfer->SetValue(doflist[i],wh_trans,INS_VAL);
-
 				}
 			}
Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp	(revision 16899)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp	(revision 16900)
@@ -57,4 +57,9 @@
 		femmodel->SetCurrentConfiguration(HydrologyDCEfficientAnalysisEnum);
 		GetSolutionFromInputsx(&ug_epl,femmodel);
+		/*Initialize the transfer input*/
+		HydrologyDCInefficientAnalysis* analysis = new HydrologyDCInefficientAnalysis();
+		analysis->ElementizeEplMask(femmodel);
+		delete analysis;
+		femmodel->HydrologyTransferx();
 	}
 
@@ -83,12 +88,4 @@
 		sedconverged=false;
 		for(;;){
-			if(isefficientlayer){
-				/*Updating Nodal Mask*/
-				HydrologyDCInefficientAnalysis* analysis = new HydrologyDCInefficientAnalysis();
-				analysis->ElementizeEplMask(femmodel);
-				delete analysis;
-			}
-			femmodel->HydrologyTransferx();
-
 			SystemMatricesx(&Kff,&Kfs,&pf,&df,&sediment_kmax,femmodel);
 			CreateNodalConstraintsx(&ys,femmodel->nodes,HydrologyDCInefficientAnalysisEnum);
@@ -104,4 +101,5 @@
 			InputUpdateFromSolutionx(femmodel,ug_sed);
 			ConstraintsStatex(&constraints_converged,&num_unstable_constraints,femmodel);
+
 			if (!sedconverged){
 				if(VerboseConvergence()) _printf0_("   # Sediment unstable constraints = " << num_unstable_constraints << "\n");
@@ -114,4 +112,11 @@
 
 			if(sedconverged){
+				if(isefficientlayer){
+				 	/*Updating Nodal Mask*/
+					HydrologyDCInefficientAnalysis* analysis = new HydrologyDCInefficientAnalysis();
+					analysis->ElementizeEplMask(femmodel);
+					delete analysis;
+					femmodel->HydrologyTransferx();
+				}
 				femmodel->parameters->SetParam(sediment_kmax,HydrologySedimentKmaxEnum);
 				InputUpdateFromConstantx(femmodel,sedconverged,ConvergedEnum);
@@ -147,13 +152,8 @@
 
 				femmodel->SetCurrentConfiguration(HydrologyDCEfficientAnalysisEnum);
-				//				femmodel->HydrologyEPLThicknessx();
+
 				//updating mask after the computation of the epl thickness (Allow to close too thin EPL)
 				femmodel->HydrologyEPLupdateDomainx();
-				/*Updating Nodal Mask*/
-				HydrologyDCInefficientAnalysis* analysis = new HydrologyDCInefficientAnalysis();
-				analysis->ElementizeEplMask(femmodel);
-				delete analysis;
-
-				femmodel->HydrologyTransferx();
+
 				SystemMatricesx(&Kff,&Kfs,&pf,&df,NULL,femmodel);
 				CreateNodalConstraintsx(&ys,femmodel->nodes,HydrologyDCEfficientAnalysisEnum);
@@ -183,4 +183,11 @@
 
 				if(eplconverged){
+
+					Updating Nodal Mask
+					HydrologyDCInefficientAnalysis* analysis = new HydrologyDCInefficientAnalysis();
+					analysis->ElementizeEplMask(femmodel);
+					delete analysis;
+					femmodel->HydrologyTransferx();
+
 					InputUpdateFromConstantx(femmodel,eplconverged,ConvergedEnum);
 					InputUpdateFromConstantx(femmodel,sediment_kmax,MeltingOffsetEnum);
@@ -190,4 +197,6 @@
 			}
 		}
+		
+		
 		/*System convergence check*/
 		if(!hydroconverged){