Index: ../trunk-jpl/src/m/classes/masstransport.py
===================================================================
--- ../trunk-jpl/src/m/classes/masstransport.py	(revision 18345)
+++ ../trunk-jpl/src/m/classes/masstransport.py	(revision 18346)
@@ -34,7 +34,7 @@
 		string="%s\n%s"%(string,fielddisplay(self,'isfreesurface','do we use free surfaces (FS only) are mass conservation'))
 		string="%s\n%s"%(string,fielddisplay(self,'min_thickness','minimum ice thickness allowed [m]'))
 		string="%s\n%s"%(string,fielddisplay(self,'hydrostatic_adjustment','adjustment of ice shelves surface and bed elevations: ''Incremental'' or ''Absolute'' '))
-		string="%s\n%s"%(string,fielddisplay(self,'stabilization','0: no, 1: artificial_diffusivity, 2: streamline upwinding, 3: discontinuous Galerkin'))
+		string="%s\n%s"%(string,fielddisplay(self,'stabilization','0: no, 1: artificial_diffusivity, 2: streamline upwinding, 3: discontinuous Galerkin, 4: Flux Correction Transport'))
 		string="%s\n%s"%(string,fielddisplay(self,'requested_outputs','additional outputs requested'))
 
 		return string
@@ -73,7 +73,7 @@
 			md = checkfield(md,'fieldname','masstransport.calvingrate','NaN',1,'size',[md.mesh.numberofvertices],'>=',0)
 		md = checkfield(md,'fieldname','masstransport.isfreesurface','values',[0,1])
 		md = checkfield(md,'fieldname','masstransport.hydrostatic_adjustment','values',['Absolute','Incremental'])
-		md = checkfield(md,'fieldname','masstransport.stabilization','values',[0,1,2,3])
+		md = checkfield(md,'fieldname','masstransport.stabilization','values',[0,1,2,3,4])
 		md = checkfield(md,'fieldname','masstransport.min_thickness','>',0)
 		md = checkfield(md,'fieldname','masstransport.requested_outputs','stringrow',1)
 
Index: ../trunk-jpl/src/m/classes/masstransport.m
===================================================================
--- ../trunk-jpl/src/m/classes/masstransport.m	(revision 18345)
+++ ../trunk-jpl/src/m/classes/masstransport.m	(revision 18346)
@@ -126,7 +126,7 @@
 
 			md = checkfield(md,'fieldname','masstransport.isfreesurface','values',[0 1]);
 			md = checkfield(md,'fieldname','masstransport.hydrostatic_adjustment','values',{'Absolute' 'Incremental'});
-			md = checkfield(md,'fieldname','masstransport.stabilization','values',[0 1 2 3]);
+			md = checkfield(md,'fieldname','masstransport.stabilization','values',[0 1 2 3 4]);
 			md = checkfield(md,'fieldname','masstransport.min_thickness','>',0);
 			md = checkfield(md,'fieldname','masstransport.requested_outputs','stringrow',1);
 
@@ -138,7 +138,7 @@
 			fielddisplay(obj,'isfreesurface','do we use free surfaces (FS only) are mass conservation');
 			fielddisplay(obj,'min_thickness','minimum ice thickness allowed [m]');
 			fielddisplay(obj,'hydrostatic_adjustment','adjustment of ice shelves surface and bed elevations: ''Incremental'' or ''Absolute'' ');
-			fielddisplay(obj,'stabilization','0: no, 1: artificial_diffusivity, 2: streamline upwinding, 3: discontinuous Galerkin');
+			fielddisplay(obj,'stabilization','0: no, 1: artificial_diffusivity, 2: streamline upwinding, 3: discontinuous Galerkin, 4: Flux Correction Transport');
 
 			disp(sprintf('\n      %s','Penalty options:'));
 			fielddisplay(obj,'penalty_factor','offset used by penalties: penalty = Kmax*10^offset');
