Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 17737)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 17738)
@@ -54,4 +54,5 @@
 	iomodel->FetchDataToInput(elements,SedimentHeadEnum);
 	iomodel->FetchDataToInput(elements,HydrologydcEplInitialThicknessEnum);
+	iomodel->FetchDataToInput(elements,HydrologydcEplMaxThicknessEnum);
 	iomodel->FetchDataToInput(elements,HydrologydcSedimentTransmitivityEnum);
 	iomodel->FetchDataToInput(elements,HydrologydcEplThicknessEnum);
Index: /issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp	(revision 17737)
+++ /issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp	(revision 17738)
@@ -79,4 +79,5 @@
 				iomodel->Constant(&this->epl_porosity,HydrologydcEplPorosityEnum);
 				iomodel->Constant(&this->epl_init_thickness,HydrologydcEplInitialThicknessEnum);
+				iomodel->Constant(&this->epl_max_thickness,HydrologydcEplMaxThicknessEnum);
 				iomodel->Constant(&this->epl_conductivity,HydrologydcEplConductivityEnum);
 		}
@@ -263,4 +264,5 @@
 		case HydrologydcEplConductivityEnum:         return this->epl_conductivity;
 		case HydrologydcEplInitialThicknessEnum:     return this->epl_init_thickness;
+		case HydrologydcEplMaxThicknessEnum:         return this->epl_max_thickness;
 		case HydrologydcWaterCompressibilityEnum:    return this->water_compressibility;
 		case HydrologyshreveCREnum:                  return this->hydro_CR;
@@ -401,4 +403,9 @@
 IssmDouble Matpar::GetEplInitialThickness(){
 	return epl_init_thickness;		 
+}		 
+/*}}}*/	
+/*FUNCTION Matpar::GetEplMaxThickness {{{*/
+IssmDouble Matpar::GetEplMaxThickness(){
+	return epl_max_thickness;		 
 }		 
 /*}}}*/	
Index: /issm/trunk-jpl/src/c/classes/Materials/Matpar.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Materials/Matpar.h	(revision 17737)
+++ /issm/trunk-jpl/src/c/classes/Materials/Matpar.h	(revision 17738)
@@ -49,4 +49,5 @@
 		IssmDouble  epl_porosity;
 		IssmDouble  epl_init_thickness;
+		IssmDouble  epl_max_thickness;
 		IssmDouble  epl_conductivity;	 
 
@@ -121,4 +122,5 @@
 		IssmDouble GetEplConductivity();
 		IssmDouble GetEplInitialThickness();
+		IssmDouble GetEplMaxThickness();
 		IssmDouble TMeltingPoint(IssmDouble pressure);
 		IssmDouble PureIceEnthalpy(IssmDouble pressure);
Index: /issm/trunk-jpl/src/m/classes/hydrologydc.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/hydrologydc.m	(revision 17737)
+++ /issm/trunk-jpl/src/m/classes/hydrologydc.m	(revision 17738)
@@ -30,4 +30,5 @@
 		epl_porosity             = 0;
 		epl_initial_thickness    = 0;
+		epl_max_thickness    = 0;
 		epl_conductivity         = 0;
   end
@@ -72,4 +73,5 @@
 			obj.epl_porosity             = 0.4;
 			obj.epl_initial_thickness    = 1.0;
+			obj.epl_max_thickness    = 5.0;
 			obj.epl_conductivity         = 8.0e-02;
 
@@ -110,4 +112,5 @@
 				md = checkfield(md,'fieldname','hydrology.epl_porosity','>',0,'numel',1);
 				md = checkfield(md,'fieldname','hydrology.epl_initial_thickness','>',0,'numel',1);
+				md = checkfield(md,'fieldname','hydrology.epl_max_thickness','>',0,'numel',1);
 				md = checkfield(md,'fieldname','hydrology.epl_conductivity','>',0,'numel',1);
 	    end
@@ -153,4 +156,5 @@
 				fielddisplay(obj,'epl_porosity','epl [dimensionless]');
 				fielddisplay(obj,'epl_initial_thickness','epl initial thickness [m]');
+				fielddisplay(obj,'epl_max_thickness','epl maximal thickness [m]');
 				fielddisplay(obj,'epl_conductivity','epl conductivity [m^2/s]');
 	    end
@@ -189,4 +193,5 @@
 				WriteData(fid,'object',obj,'fieldname','epl_porosity','format','Double');			
 				WriteData(fid,'object',obj,'fieldname','epl_initial_thickness','format','Double');
+				WriteData(fid,'object',obj,'fieldname','epl_max_thickness','format','Double');
 				WriteData(fid,'object',obj,'fieldname','epl_conductivity','format','Double');
 			end
Index: /issm/trunk-jpl/test/NightlyRun/test3300.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3300.m	(revision 17737)
+++ /issm/trunk-jpl/test/NightlyRun/test3300.m	(revision 17738)
@@ -24,4 +24,5 @@
 md.hydrology.epl_conductivity=1.5e-2;
 md.hydrology.epl_initial_thickness=1.0;
+md.hydrology.epl_max_thickness=5.0;
 
 md.hydrology.transfer_flag=1;
Index: /issm/trunk-jpl/test/NightlyRun/test333.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test333.m	(revision 17737)
+++ /issm/trunk-jpl/test/NightlyRun/test333.m	(revision 17738)
@@ -23,4 +23,5 @@
 md.hydrology.epl_conductivity=30;
 md.hydrology.epl_initial_thickness=1;
+md.hydrology.epl_max_thickness=1;
 md.timestepping.time_step=0.2;
 md.timestepping.final_time=2.0;
Index: /issm/trunk-jpl/test/NightlyRun/test335.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test335.m	(revision 17737)
+++ /issm/trunk-jpl/test/NightlyRun/test335.m	(revision 17738)
@@ -23,4 +23,5 @@
 md.hydrology.epl_conductivity=30;
 md.hydrology.epl_initial_thickness=1;
+md.hydrology.epl_max_thickness=1;
 md.timestepping.time_step=0.2;
 md.timestepping.final_time=2.0;