Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 17011)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 17012)
@@ -59,4 +59,6 @@
 	iomodel->FetchDataToInput(elements,HydrologydcEplInitialThicknessEnum);
 	
+	elements->InputDuplicate(HydrologydcEplInitialThicknessEnum,HydrologydcEplThicknessEnum);
+
 }/*}}}*/
 void HydrologyDCEfficientAnalysis::CreateNodes(Nodes* nodes,IoModel* iomodel){/*{{{*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17011)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17012)
@@ -4834,5 +4834,5 @@
 	IssmDouble  n,A,dt,init_thick;
 	IssmDouble  rho_water,rho_ice;
-	IssmDouble  gravity,latentheat,EPLgrad;
+	IssmDouble  gravity,latentheat,EPLgrad2;
 	IssmDouble  EPL_N,epl_conductivity;
 	IssmDouble  activeEpl[numdof],thickness[numdof];
@@ -4876,9 +4876,10 @@
 					EPL_N=gravity*((rho_ice*ice_thickness[i])-(rho_water*(eplhead[i]-bed[i])));
 					if(EPL_N<0.0)EPL_N=0.0;
-					/*Get then the square of th gradient of EPL heads*/
-					EPLgrad = (epl_slopeX[i]*epl_slopeX[i]+epl_slopeY[i]*epl_slopeY[i]);
+					/*Get then the square of the gradient of EPL heads*/
+					EPLgrad2 = (epl_slopeX[i]*epl_slopeX[i]+epl_slopeY[i]*epl_slopeY[i]);
 					
 					/*And proceed to the real thing*/
-					thickness[i] = old_thickness[i]*(1+((rho_water*gravity*dt)/(rho_ice*latentheat))*epl_conductivity*EPLgrad-2.0*(A*dt/(pow(n,n)))*(pow(EPL_N,n)));
+					thickness[i] = old_thickness[i]*(1+((rho_water*gravity*dt)/(rho_ice*latentheat))*epl_conductivity*EPLgrad2-2.0*(A*dt/(pow(n,n)))*(pow(EPL_N,n)));
+
 			}
 		}
Index: /issm/trunk-jpl/src/c/cores/hydrology_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/hydrology_core.cpp	(revision 17011)
+++ /issm/trunk-jpl/src/c/cores/hydrology_core.cpp	(revision 17012)
@@ -84,6 +84,5 @@
 			if (isefficientlayer){
 				InputDuplicatex(femmodel,EplHeadEnum,EplHeadOldEnum);
-				InputDuplicatex(femmodel,HydrologydcEplInitialThicknessEnum,HydrologydcEplThicknessEnum);
-				InputDuplicatex(femmodel,HydrologydcEplInitialThicknessEnum,HydrologydcEplThicknessOldEnum);
+				InputDuplicatex(femmodel,HydrologydcEplThicknessEnum,HydrologydcEplThicknessOldEnum);
 			}
 			
Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp	(revision 17011)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp	(revision 17012)
@@ -178,4 +178,5 @@
 				solutionsequence_linear(femmodel);
 
+				
 				femmodel->HydrologyEPLThicknessx();