Index: ../trunk-jpl/src/m/classes/hydrologydc.py
===================================================================
--- ../trunk-jpl/src/m/classes/hydrologydc.py	(revision 18584)
+++ ../trunk-jpl/src/m/classes/hydrologydc.py	(revision 18585)
@@ -22,6 +22,7 @@
 		self.sedimentlimit            = 0
 		self.transfer_flag            = 0
 		self.leakage_factor           = 0
+                self.basal_moulin_input       = float('NaN')
 
 		self.spcsediment_head         = float('NaN')
 		self.sediment_transmitivity   = float('NaN')
@@ -35,10 +36,12 @@
 		self.epl_compressibility      = 0
 		self.epl_porosity             = 0
 		self.epl_initial_thickness    = 0
+		self.epl_max_thickness        = 0
 		self.epl_conductivity         = 0
                 
 		#set defaults
 		self.setdefaultparameters()
+
                 #}}}
                 
 	def __repr__(self): # {{{
@@ -58,7 +61,8 @@
 
             if self.sedimentlimit_flag==1:
                 string="%s\n%s"%(string,fielddisplay(self,'sedimentlimit','user defined upper limit for the inefficient layer [m]'))
-                string="%s\n%s"%(string,fielddisplay(self,'transfer_flag',['what kind of transfer method is applied between the layers']))
+                string="%s\n%s"%(string,fielddisplay(self,'basal_moulin_input','water flux at a given point [m3 s-1]'))
+                string="%s\n%s"%(string,fielddisplay(self,'transfer_flag','what kind of transfer method is applied between the layers'))
                 string='%55s  0: no transfer',' '
                 string='%55s  1: constant leakage factor: %s',' ','leakage_factor'
                 
@@ -102,11 +106,20 @@
         self.epl_compressibility      = 1.0e-08
         self.epl_porosity             = 0.4
         self.epl_initial_thickness    = 1.0
+        self.epl_max_thickness        = 5.0
         self.epl_conductivity         = 8.0e-02
         
         return self
     # }}}
     
+    def initialize(self): # {{{
+        if numpy.all(numpy.isnan(self.basal_moulin_input):
+            self.basal_moulin_input=numpy.zeros(md.mesh.numberofvertices,1)
+            print"      no hydrology.basal_moulin_input specified: values set as zero"
+
+        return self
+
+    # }}}
     def checkconsistency(self,md,solution,analyses): #{{{ 
         
         #Early return
@@ -126,7 +139,8 @@
     
         if self.transfer_flag==1:
             md = checkfield(md,'fieldname','hydrology.leakage_factor','>',0,'numel',1)
-            
+
+        md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,'forcing',1)
         md = checkfield(md,'fieldname','hydrology.spcsediment_head','forcing',1)
         md = checkfield(md,'fieldname','hydrology.sediment_compressibility','>',0,'numel',1)
         md = checkfield(md,'fieldname','hydrology.sediment_porosity','>',0,'numel',1)
@@ -158,6 +172,7 @@
         if self.transfer_flag==1:
             WriteData(fid,'object',self,'fieldname','leakage_factor','format','Double')
 
+	WriteData(fid,'object',self,'fieldname','basal_moulin_input','format','DoubleMat','mattype',1,'forcinglength',md.mesh.numberofvertices+1)
         WriteData(fid,'object',self,'fieldname','spcsediment_head','format','DoubleMat','mattype',1,'forcinglength',md.mesh.numberofvertices+1)
         WriteData(fid,'object',self,'fieldname','sediment_compressibility','format','Double')
         WriteData(fid,'object',self,'fieldname','sediment_porosity','format','Double')			
