Index: ../trunk-jpl/src/m/classes/hydrologydc.py
===================================================================
--- ../trunk-jpl/src/m/classes/hydrologydc.py	(revision 18748)
+++ ../trunk-jpl/src/m/classes/hydrologydc.py	(revision 18749)
@@ -40,6 +40,7 @@
 		self.epl_colapse_thickness    = 0
 		self.epl_max_thickness        = 0
 		self.epl_conductivity         = 0
+		self.eplflip_lock             = 0
 				 
 		#set defaults
 		self.setdefaultparameters()
@@ -89,6 +90,7 @@
 			string="%s\n%s"%(string,fielddisplay(self,'epl_initial_thickness','epl initial thickness [m]'))			
 			string="%s\n%s"%(string,fielddisplay(self,'epl_colapse_thickness','epl colapsing thickness [m]'))
 			string="%s\n%s"%(string,fielddisplay(self,'epl_conductivity','epl conductivity [m^2/s]'))
+			string="%s\n%s"%(string,fielddisplay(self,'eplflip_lock','lock epl activity to avoid flip-floping (default is 0, no stabilization)'))
 		return string
 #}}}
 	def setdefaultparameters(self): #{{{ 
@@ -97,6 +99,7 @@
 		self.water_compressibility    = 5.04e-10
 		self.isefficientlayer         = 1
 		self.penalty_factor           = 3
+		self.penalty_lock             = 0
 		self.rel_tol                  = 1.0e-06
 		self.max_iter                 = 100
 		self.sedimentlimit_flag       = 0
@@ -115,7 +118,8 @@
 		self.epl_colapse_thickness    = 1.0e-3
 		self.epl_max_thickness        = 5.0
 		self.epl_conductivity         = 8.0e-02
-
+		self.eplflip_lock             = 0
+		
 		return self
 	# }}}
 
@@ -135,6 +139,7 @@
 		md = checkfield(md,'fieldname','hydrology.water_compressibility','numel',[1],'>',0.)
 		md = checkfield(md,'fieldname','hydrology.isefficientlayer','numel',[1],'values',[0,1])
 		md = checkfield(md,'fieldname','hydrology.penalty_factor','>',0.,'numel',[1])
+		md = checkfield(md,'fieldname','hydrology.penalty_lock','>=',0.,'numel',[1])
 		md = checkfield(md,'fieldname','hydrology.rel_tol','>',0.,'numel',[1])
 		md = checkfield(md,'fieldname','hydrology.max_iter','>',0.,'numel',[1])
 		md = checkfield(md,'fieldname','hydrology.sedimentlimit_flag','numel',[1],'values',[0,1,2,3])
@@ -160,6 +165,7 @@
 			md = checkfield(md,'fieldname','hydrology.epl_max_thickness','numel',[1],'>',0.)
 			md = checkfield(md,'fieldname','hydrology.epl_initial_thickness','numel',[1],'>',0.)
 			md = checkfield(md,'fieldname','hydrology.epl_colapse_thickness','numel',[1],'>',0.)
+			md = checkfield(md,'fieldname','hydrology.eplflip_lock','>=',0.,'numel',[1])
 			if self.epl_colapse_thickness > self.epl_initial_thickness:
 				md.checkmessage('Colapsing thickness for EPL larger than initial thickness')
 			md = checkfield(md,'fieldname','hydrology.epl_conductivity','numel',[1],'>',0.)
@@ -196,4 +202,5 @@
 			WriteData(fid,'object',self,'fieldname','epl_initial_thickness','format','Double')			
 			WriteData(fid,'object',self,'fieldname','epl_colapse_thickness','format','Double')			
 			WriteData(fid,'object',self,'fieldname','epl_conductivity','format','Double')
+			WriteData(fid,'object',self,'fieldname','eplflip_lock','format','Integer')
 	# }}}
