Index: /issm/trunk-jpl/src/m/classes/hydrologydc.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/hydrologydc.m	(revision 21965)
+++ /issm/trunk-jpl/src/m/classes/hydrologydc.m	(revision 21966)
@@ -15,4 +15,5 @@
 		sedimentlimit            = 0;
 		transfer_flag            = 0;
+		unconfined_flag          = 0;
 		leakage_factor           = 0;
 		basal_moulin_input       = NaN;
@@ -79,4 +80,5 @@
 			self.sedimentlimit            = 0;
 			self.transfer_flag            = 0;
+			self.unconfined_flag          = 0;
 			self.leakage_factor           = 10.0;
 
@@ -113,4 +115,5 @@
 			md = checkfield(md,'fieldname','hydrology.sedimentlimit_flag','numel',[1],'values',[0 1 2 3]);
 			md = checkfield(md,'fieldname','hydrology.transfer_flag','numel',[1],'values',[0 1]);
+			md = checkfield(md,'fieldname','hydrology.unconfined_flag','numel',[1],'values',[0 1]);
 			if self.sedimentlimit_flag==1,
 				md = checkfield(md,'fieldname','hydrology.sedimentlimit','>',0,'numel',1);
@@ -170,4 +173,7 @@
 				fielddisplay(self,'leakage_factor','user defined leakage factor [m]');
 			end
+			fielddisplay(self,'unconfined_flag',['Do you want unconfined scheme to be used (transitory)']);
+			disp(sprintf('%55s  0: confined only',' '));
+			disp(sprintf('%55s  1: confined unconfined'));
 			fielddisplay(self,'basal_moulin_input','water flux at a given point [m3 s-1]');
 			disp(sprintf('   - for the sediment layer'));
@@ -207,4 +213,5 @@
 			WriteData(fid,prefix,'object',self,'fieldname','sedimentlimit_flag','format','Integer');
 			WriteData(fid,prefix,'object',self,'fieldname','transfer_flag','format','Integer');
+			WriteData(fid,prefix,'object',self,'fieldname','unconfined_flag','format','Integer');
 			if self.sedimentlimit_flag==1,
 				WriteData(fid,prefix,'object',self,'fieldname','sedimentlimit','format','Double');
Index: /issm/trunk-jpl/src/m/classes/hydrologydc.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/hydrologydc.py	(revision 21965)
+++ /issm/trunk-jpl/src/m/classes/hydrologydc.py	(revision 21966)
@@ -23,4 +23,5 @@
 		self.sedimentlimit            = 0
 		self.transfer_flag            = 0
+		self.unconfined_flag          = 0
 		self.leakage_factor           = 0
 		self.basal_moulin_input       = float('NaN')
@@ -73,4 +74,8 @@
 			string="%s\n%s"%(string,fielddisplay(self,'leakage_factor','user defined leakage factor [m]'))
 
+		string="%s\n%s"%(string,fielddisplay(self,'unconfined_flag','using an unconfined scheme or not (transitory)'))
+		string="%s\n\t\t%s"%(string,'0: Confined only')
+		string="%s\n\t\t%s"%(string,'1: Confined-Unconfined')
+			 
 		string="%s\n%s"%(string,'   - for the sediment layer')
 		string="%s\n%s"%(string,fielddisplay(self,'spcsediment_head','sediment water head constraints (NaN means no constraint) [m above MSL]'))
@@ -115,4 +120,5 @@
 		self.sedimentlimit            = 0
 		self.transfer_flag            = 0
+		self.unconfined_flag          = 0
 		self.leakage_factor           = 10.0
 
@@ -154,4 +160,5 @@
 		md = checkfield(md,'fieldname','hydrology.sedimentlimit_flag','numel',[1],'values',[0,1,2,3])
 		md = checkfield(md,'fieldname','hydrology.transfer_flag','numel',[1],'values',[0,1])
+		md = checkfield(md,'fieldname','hydrology.unconfined_flag','numel',[1],'values',[0,1])
 
 		if self.sedimentlimit_flag==1:
@@ -191,4 +198,5 @@
 		WriteData(fid,prefix,'object',self,'fieldname','sedimentlimit_flag','format','Integer')
 		WriteData(fid,prefix,'object',self,'fieldname','transfer_flag','format','Integer')
+		WriteData(fid,prefix,'object',self,'fieldname','unconfined_flag','format','Integer')
 		if self.sedimentlimit_flag==1:
 			WriteData(fid,prefix,'object',self,'fieldname','sedimentlimit','format','Double')