Index: /issm/trunk-jpl/src/m/classes/hydrologydc.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/hydrologydc.m	(revision 22880)
+++ /issm/trunk-jpl/src/m/classes/hydrologydc.m	(revision 22881)
@@ -75,4 +75,5 @@
 
 		function self = initialize(self,md) % {{{
+			self.epl_colapse_thickness = self.sediment_transmitivity/self.epl_conductivity;
 			if isnan(self.basal_moulin_input),
 				self.basal_moulin_input=zeros(md.mesh.numberofvertices,1);
@@ -107,10 +108,10 @@
 
 			self.epl_compressibility      = 1.0e-08;
+			self.epl_conductivity         = 8.0e-02;
 			self.epl_porosity             = 0.4;
 			self.epl_initial_thickness    = 1.0;
-			self.epl_colapse_thickness    = 1.0e-3;
+			self.epl_colapse_thickness    = self.sediment_transmitivity/self.epl_conductivity;
 			self.epl_thick_comp           = 1;
 			self.epl_max_thickness        = 5.0;
-			self.epl_conductivity         = 8.0e-02;
 			self.eplflip_lock             = 0;
 		end
Index: /issm/trunk-jpl/src/m/classes/hydrologydc.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/hydrologydc.py	(revision 22880)
+++ /issm/trunk-jpl/src/m/classes/hydrologydc.py	(revision 22881)
@@ -133,10 +133,10 @@
 
 		self.epl_compressibility      = 1.0e-08
+		self.epl_conductivity         = 8.0e-02
 		self.epl_porosity             = 0.4
 		self.epl_initial_thickness    = 1.0
-		self.epl_colapse_thickness    = 1.0e-3
+		self.epl_colapse_thickness    = self.sediment_transmitivity/self.epl_conductivity
 		self.epl_thick_comp           = 1
 		self.epl_max_thickness        = 5.0
-		self.epl_conductivity         = 8.0e-02
 		self.eplflip_lock             = 0
 
@@ -156,4 +156,5 @@
 
 	def initialize(self,md): # {{{
+		self.epl_colapse_thickness = self.sediment_transmitivity/self.epl_conductivity
 		if np.all(np.isnan(self.basal_moulin_input)):
 			self.basal_moulin_input=np.zeros((md.mesh.numberofvertices))