Index: ../trunk-jpl/src/m/classes/hydrologydc.m
===================================================================
--- ../trunk-jpl/src/m/classes/hydrologydc.m	(revision 18747)
+++ ../trunk-jpl/src/m/classes/hydrologydc.m	(revision 18748)
@@ -32,6 +32,7 @@
 		epl_colapse_thickness    = 0;
 		epl_max_thickness        = 0;
 		epl_conductivity         = 0;
+		eplflip_lock             = 0;
   end
 	methods
 		% {{{ function obj = hydrologydc(varargin) 
@@ -58,6 +59,7 @@
 			obj.water_compressibility    = 5.04e-10;
 			obj.isefficientlayer         = 1;
 			obj.penalty_factor           = 3;
+			obj.penalty_lock             = 0;
 			obj.rel_tol                  = 1.0e-06;
 			obj.max_iter                 = 100;
 			obj.sedimentlimit_flag       = 0;
@@ -76,7 +78,7 @@
 			obj.epl_colapse_thickness    = 1.0e-3;
 			obj.epl_max_thickness        = 5.0;
 			obj.epl_conductivity         = 8.0e-02;
-
+			obj.eplflip_lock             = 0;
 		end 
 		% }}}
 		% {{{ function md = checkconsistency(obj,md,solution,analyses) 
@@ -89,6 +91,7 @@
 			md = checkfield(md,'fieldname','hydrology.water_compressibility','>',0,'numel',1);
 			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]);
@@ -116,6 +119,7 @@
 				md = checkfield(md,'fieldname','hydrology.epl_colapse_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);
+				md = checkfield(md,'fieldname','hydrology.eplflip_lock','>=',0,'numel',1);
 				if (obj.epl_colapse_thickness>obj.epl_initial_thickness),
 					md = checkmessage(md,'Colapsing thickness for EPL larger than initial thickness');
 				end 
@@ -164,6 +168,7 @@
 				fielddisplay(obj,'epl_colapse_thickness','epl colapsing thickness [m]');
 				fielddisplay(obj,'epl_max_thickness','epl maximal thickness [m]');
 				fielddisplay(obj,'epl_conductivity','epl conductivity [m^2/s]');
+				fielddisplay(obj,'eplflip_lock','lock the epl activation to avoid fli-floping (default is 0, no stabilization)');
 	    end
 
 		end 
@@ -202,7 +207,8 @@
 				WriteData(fid,'object',obj,'fieldname','epl_colapse_thickness','format','Double');
 				WriteData(fid,'object',obj,'fieldname','epl_max_thickness','format','Double');
 				WriteData(fid,'object',obj,'fieldname','epl_conductivity','format','Double');
-			end
+				WriteData(fid,'object',obj,'fieldname','eplflip_lock','format','Integer');
+	    end
 		end 
 % }}}
   end
Index: ../trunk-jpl/src/m/enum/HydrologydcEplflipLockEnum.m
===================================================================
--- ../trunk-jpl/src/m/enum/HydrologydcEplflipLockEnum.m	(revision 0)
+++ ../trunk-jpl/src/m/enum/HydrologydcEplflipLockEnum.m	(revision 18748)
@@ -0,0 +1,11 @@
+function macro=HydrologydcEplflipLockEnum()
+%HYDROLOGYDCEPLFLIPLOCKENUM - Enum of HydrologydcEplflipLock
+%
+%   WARNING: DO NOT MODIFY THIS FILE
+%            this file has been automatically generated by src/c/shared/Enum/Synchronize.sh
+%            Please read src/c/shared/Enum/README for more information
+%
+%   Usage:
+%      macro=HydrologydcEplflipLockEnum()
+
+macro=StringToEnum('HydrologydcEplflipLock');
Index: ../trunk-jpl/src/m/enum/EnumDefinitions.py
===================================================================
--- ../trunk-jpl/src/m/enum/EnumDefinitions.py	(revision 18747)
+++ ../trunk-jpl/src/m/enum/EnumDefinitions.py	(revision 18748)
@@ -136,6 +136,7 @@
 def HydrologydcLeakageFactorEnum(): return StringToEnum("HydrologydcLeakageFactor")[0]
 def HydrologydcPenaltyFactorEnum(): return StringToEnum("HydrologydcPenaltyFactor")[0]
 def HydrologydcPenaltyLockEnum(): return StringToEnum("HydrologydcPenaltyLock")[0]
+def HydrologydcEplflipLockEnum(): return StringToEnum("HydrologydcEplflipLock")[0]
 def HydrologydcBasalMoulinInputEnum(): return StringToEnum("HydrologydcBasalMoulinInput")[0]
 def HydrologyLayerEnum(): return StringToEnum("HydrologyLayer")[0]
 def HydrologySedimentEnum(): return StringToEnum("HydrologySediment")[0]
