Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 21278)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 21279)
@@ -498,7 +498,9 @@
 		//Computing continuum function to apply to transfer term, transfer is null only if
 		// epl_head>sediment_head AND sediment_head>h_max
-		continuum=((1.0/(1.0+exp(-20.0*(sediment_head-epl_head)))))+(1.0/(1.0+exp(-20.0*(hmax-sediment_head))));
-		factor=min(continuum,1.0);
-		transfer=leakage*factor;
+		//let's try without that for a while
+		/* continuum=((1.0/(1.0+exp(-20.0*(sediment_head-epl_head)))))+(1.0/(1.0+exp(-20.0*(hmax-sediment_head)))); */
+		/* factor=min(continuum,1.0); */
+		/* transfer=leakage*factor; */
+		transfer=leakage;
 		break;
 	default:
@@ -540,7 +542,9 @@
 		//Computing continuum function to apply to transfer term, transfer is null only if
 		// epl_head>sediment_head AND sediment_head>h_max
-		continuum=((1.0/(1.0+exp(-20.0*(sediment_head-epl_head)))))+(1.0/(1.0+exp(-20.0*(hmax-sediment_head))));
-		factor=min(continuum,1.0);
-		transfer=sediment_head*leakage*factor;
+		//let's try without that for a while
+		/* continuum=((1.0/(1.0+exp(-20.0*(sediment_head-epl_head)))))+(1.0/(1.0+exp(-20.0*(hmax-sediment_head)))); */
+		/* factor=min(continuum,1.0); */
+		/* transfer=sediment_head*leakage*factor; */
+		transfer=sediment_head*leakage;
 
 		break;
@@ -625,5 +629,7 @@
 				
 				/*Compute first the effective pressure in the EPL*/
-				EPL_N=gravity*((rho_ice*ice_thickness[i])-(rho_water*max(0.0,(eplhead[i]-bed[i]))));
+				//EPL_N=gravity*((rho_ice*ice_thickness[i])-(rho_water*max(0.0,(eplhead[i]-bed[i]))));
+				//allowing EPLN larger than Pi
+				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 the gradient of EPL heads*/
@@ -631,7 +637,10 @@
 				
 				/*And proceed to the real thing*/
-				thickness[i] = old_thickness[i]*(1.0+((rho_water*gravity*dt)/(rho_ice*latentheat))*epl_conductivity*EPLgrad2-
-																				 (2.0*(A*dt/(pow(n,n)))*(pow(EPL_N,n))));
-				
+				/* thickness[i] = old_thickness[i]*(1.0+((rho_water*gravity*dt)/(rho_ice*latentheat))*epl_conductivity*EPLgrad2- */
+				/* 																 (2.0*(A*dt/(pow(n,n)))*(pow(EPL_N,n)))); */
+				thickness[i] = old_thickness[i]*
+					(2.0+((rho_water*gravity*dt)/(rho_ice*latentheat))*epl_conductivity*EPLgrad2-(2.0*(A*dt/(pow(n,n)))*(pow(EPL_N,n))))/
+					(2.0-((rho_water*gravity*dt)/(rho_ice*latentheat))*epl_conductivity*EPLgrad2+(2.0*(A*dt/(pow(n,n)))*(pow(EPL_N,n))));
+				//thickness[i] = 0.5*(thickness[i]+old_thickness[i]);
 				/*Take care of otherthikening*/
 				if(thickness[i]>max_thick){
Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 21278)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 21279)
@@ -629,7 +629,9 @@
 		//Computing continuum function to apply to transfer term, transfer is null only if
 		//epl_head>sediment_head AND sediment_head>h_max
-		continuum=((1.0/(1.0+exp(-20.0*(sediment_head-epl_head)))))+(1.0/(1.0+exp(-20.0*(hmax-sediment_head))));
-		factor=min(continuum,1.0);
-		transfer=leakage*factor;
+		//let's try without that for a while
+		/* continuum=((1.0/(1.0+exp(-20.0*(sediment_head-epl_head)))))+(1.0/(1.0+exp(-20.0*(hmax-sediment_head)))); */
+		/* factor=min(continuum,1.0); */
+		/* transfer=leakage*factor; */
+		transfer=leakage; 
 
 		break;
@@ -668,7 +670,11 @@
 		//Computing continuum function to apply to transfer term, transfer is null only if
 		//epl_head>sediment_head AND sediment_head>h_max
-		continuum=((1.0/(1.0+exp(-20.0*(sediment_head-epl_head)))))+(1.0/(1.0+exp(-20.0*(hmax-sediment_head))));
-		factor=min(continuum,1.0);
-		transfer=epl_head*leakage*factor;
+		//let's try without that for a while
+		/* continuum=((1.0/(1.0+exp(-20.0*(sediment_head-epl_head)))))+(1.0/(1.0+exp(-20.0*(hmax-sediment_head)))); */
+		/* factor=min(continuum,1.0); */
+		/* transfer=epl_head*leakage*factor; */
+
+		transfer=epl_head*leakage;
+
 		break;
 	default:
Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp	(revision 21278)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp	(revision 21279)
@@ -216,5 +216,5 @@
 				//If there is some colapse go through sediment again
 				/* if(ThickCount<L2Count)eplconverged=true; */
-				/* eplcount++; */
+				eplcount++;
 				
 				delete ug_epl_sub_iter;