Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 16724)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 16725)
@@ -57,5 +57,5 @@
 	iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum);
 	iomodel->FetchDataToInput(elements,EplHeadEnum);
-	iomodel->FetchDataToInput(elements,HydrologydcEplThicknessEnum);
+	iomodel->FetchDataToInput(elements,HydrologydcEplInitialThicknessEnum);
 	
 }/*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16724)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16725)
@@ -2724,4 +2724,5 @@
 				name==EplHeadOldEnum ||
 				name==HydrologydcEplThicknessOldEnum ||
+				name==HydrologydcEplInitialThicknessEnum ||
 				name==HydrologydcEplThicknessEnum ||
 				name==HydrologydcMaskEplactiveEnum ||
@@ -10873,5 +10874,5 @@
 	const int   numdof2d = NDOF1 *NUMVERTICES2D;
 	bool        isefficientlayer;
-	IssmDouble  n,A,dt;
+	IssmDouble  n,A,dt,init_thick;
 	IssmDouble  rho_water,rho_ice;
 	IssmDouble  gravity,latentheat,EPLgrad;
@@ -10896,7 +10897,7 @@
 		latentheat       = matpar->GetLatentHeat();
 		epl_conductivity = matpar->GetEplConductivity();
+		init_thick       = matpar->GetEplInitialThickness();
 		n                = material->GetN();
 		A                = material->GetA();
-		
 		
 		GetInputListOnVertices(&activeEpl[0],HydrologydcMaskEplactiveEnum);
@@ -10911,5 +10912,5 @@
 			/*Keeping thickness to 1 if EPL is not active*/
 			if(activeEpl[i]==0.0){
-				thickness[i]=1.0;
+				thickness[i]=init_thick;
 				thickness[i+numdof2d]=thickness[i];
 			}
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16724)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16725)
@@ -2166,4 +2166,5 @@
 				name==EplHeadEnum ||
 				name==HydrologydcEplThicknessOldEnum ||
+				name==HydrologydcEplInitialThicknessEnum ||
 				name==HydrologydcEplThicknessEnum ||
 				name==HydrologydcMaskEplactiveEnum ||
@@ -7370,5 +7371,5 @@
 	const int   numdof         = NDOF1 *NUMVERTICES;
 	bool        isefficientlayer;
-	IssmDouble  n,A,dt;
+	IssmDouble  n,A,dt,init_thick;
 	IssmDouble  rho_water,rho_ice;
 	IssmDouble  gravity,latentheat,EPLgrad;
@@ -7386,10 +7387,12 @@
 	if(isefficientlayer){
 		/*For now, assuming just one way to compute EPL thickness*/
-		rho_water  = matpar->GetRhoWater();
-		rho_ice    = matpar->GetRhoIce();
-		gravity    = matpar->GetG();
-		latentheat = matpar->GetLatentHeat();
-		n          = material->GetN();
-		A          = material->GetAbar();
+		rho_water        = matpar->GetRhoWater();
+		rho_ice          = matpar->GetRhoIce();
+		gravity          = matpar->GetG();
+		latentheat       = matpar->GetLatentHeat();
+		epl_conductivity = matpar->GetEplConductivity();
+		init_thick       = matpar->GetEplInitialThickness();
+		n                = material->GetN();
+		A                = material->GetAbar();
 		
 		GetInputListOnVertices(&activeEpl[0],HydrologydcMaskEplactiveEnum);
@@ -7401,10 +7404,8 @@
 		GetInputListOnVertices(&bed[0],BedEnum);
 		
-		epl_conductivity = matpar->GetEplConductivity();
-		
 		for(int i=0;i<numdof;i++){
 			/*Keeping thickness to 1 if EPL is not active*/
 			if(activeEpl[i]==0.0){
-				thickness[i]=1.0;
+				thickness[i]=init_thick;
 			}
 			else{
Index: /issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp	(revision 16724)
+++ /issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp	(revision 16725)
@@ -66,4 +66,5 @@
 				iomodel->Constant(&this->epl_compressibility,HydrologydcEplCompressibilityEnum);
 				iomodel->Constant(&this->epl_porosity,HydrologydcEplPorosityEnum);
+				iomodel->Constant(&this->epl_init_thickness,HydrologydcEplInitialThicknessEnum);
 				iomodel->Constant(&this->epl_conductivity,HydrologydcEplConductivityEnum);
 		}
@@ -392,4 +393,9 @@
 }		 
 /*}}}*/	
+/*FUNCTION Matpar::GetEplInitialThickness {{{*/
+IssmDouble Matpar::GetEplInitialThickness(){
+	return epl_init_thickness;		 
+}		 
+/*}}}*/	
 /*FUNCTION Matpar::GetEplConductivity {{{*/
 IssmDouble Matpar::GetEplConductivity(){
Index: /issm/trunk-jpl/src/c/classes/Materials/Matpar.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Materials/Matpar.h	(revision 16724)
+++ /issm/trunk-jpl/src/c/classes/Materials/Matpar.h	(revision 16725)
@@ -48,5 +48,6 @@
 
 		IssmDouble  epl_compressibility;
-		IssmDouble  epl_porosity;	 
+		IssmDouble  epl_porosity;
+		IssmDouble  epl_init_thickness;
 		IssmDouble  epl_conductivity;	 
 
@@ -125,4 +126,5 @@
 		IssmDouble GetSedimentThickness();
 		IssmDouble GetEplConductivity();
+		IssmDouble GetEplInitialThickness();
 		IssmDouble TMeltingPoint(IssmDouble pressure);
 		IssmDouble PureIceEnthalpy(IssmDouble pressure);
Index: /issm/trunk-jpl/src/c/cores/hydrology_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/hydrology_core.cpp	(revision 16724)
+++ /issm/trunk-jpl/src/c/cores/hydrology_core.cpp	(revision 16725)
@@ -84,5 +84,6 @@
 			if (isefficientlayer){
 				InputDuplicatex(femmodel,EplHeadEnum,EplHeadOldEnum);
-				InputDuplicatex(femmodel,HydrologydcEplThicknessEnum,HydrologydcEplThicknessOldEnum);
+				InputDuplicatex(femmodel,HydrologydcEplInitialThicknessEnum,HydrologydcEplThicknessEnum);
+				InputDuplicatex(femmodel,HydrologydcEplInitialThicknessEnum,HydrologydcEplThicknessOldEnum);
 			}