Index: /issm/trunk-jpl/src/m/classes/debris.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/debris.m	(revision 27852)
+++ /issm/trunk-jpl/src/m/classes/debris.m	(revision 27853)
@@ -12,4 +12,5 @@
 		removalmodel             = 0;
 		displacementmodel        = 0;
+		max_displacementvelocity = 0;
 		removal_slope_threshold  = 0;
 		removal_stress_threshold = 0;
@@ -66,5 +67,5 @@
 		function list = defaultoutputs(self,md) % {{{
 
-			list = {'DebrisThickness'};
+			list = {'DebrisThickness','DebrisMaskNodeActivation','VxDebris','VyDebris'};
 
 		end % }}}
@@ -92,4 +93,7 @@
 			self.removal_stress_threshold=0;
 
+			%Max velocity for displacementmodel (1)
+			self.max_displacementvelocity=0;
+
 			%default output
 			self.requested_outputs={'default'};
@@ -101,5 +105,5 @@
 
 			md = checkfield(md,'fieldname','debris.spcthickness');
-			md = checkfield(md,'fieldname','debris.stabilization','values',[0 1 2 3]);
+			md = checkfield(md,'fieldname','debris.stabilization','values',[0 1 2 3 4 5]);
 			md = checkfield(md,'fieldname','debris.min_thickness','>=',0);
 			md = checkfield(md,'fieldname','debris.packingfraction','>=',0);
@@ -108,4 +112,5 @@
 			md = checkfield(md,'fieldname','debris.removal_slope_threshold','>=',0);
 			md = checkfield(md,'fieldname','debris.removal_stress_threshold','>=',0);
+			md = checkfield(md,'fieldname','debris.max_displacementvelocity','>=',0);
 			md = checkfield(md,'fieldname','debris.requested_outputs','stringrow',1);
 			if ~any(isnan(md.stressbalance.vertex_pairing)),
@@ -120,5 +125,6 @@
 			fielddisplay(self,'stabilization','0: no stabilization, 1: artificial diffusion, 2: streamline upwinding, 3: streamline upwind Petrov-Galerkin (SUPG)');
 			fielddisplay(self,'removalmodel','frontal removal of debris. 0: no removal, 1: Slope-triggered debris removal, 2: driving-stress triggered debris removal');
-			fielddisplay(self,'displacementmodel','debris displacement. 0: no displacement, 1: ...');
+			fielddisplay(self,'max_displacementvelocity','maximum velocity of debris transport (v_ice + v_displacement) (m/a)');
+			fielddisplay(self,'displacementmodel','debris displacement. 0: no displacement, 1: additional debris velocity above the critical slope/stress threshold');
 			fielddisplay(self,'removal_slope_threshold','critical slope (degrees) for removalmodel (1)');
 			fielddisplay(self,'removal_stress_threshold','critical stress (Pa) for removalmodel (2)');
@@ -135,4 +141,5 @@
 			WriteData(fid,prefix,'object',self,'fieldname','removalmodel','format','Integer');
 			WriteData(fid,prefix,'object',self,'fieldname','displacementmodel','format','Integer');
+			WriteData(fid,prefix,'object',self,'fieldname','max_displacementvelocity','format','Double');
 			WriteData(fid,prefix,'object',self,'fieldname','removal_slope_threshold','format','Double');
 			WriteData(fid,prefix,'object',self,'fieldname','removal_stress_threshold','format','Double');
@@ -156,4 +163,5 @@
 			writejsdouble(fid,[modelname '.debris.removalmodel'],self.removalmodel);
 			writejsdouble(fid,[modelname '.debris.displacementmodel'],self.displacementmodel);
+			writejsdouble(fid,[modelname '.debris.max_displacementvelocity'],self.displacementmodel);
 			writejsdouble(fid,[modelname '.debris.removal_slope_threshold'],self.removal_slope_threshold);
 			writejsdouble(fid,[modelname '.debris.removal_stress_threshold'],self.removal_stress_threshold);
Index: /issm/trunk-jpl/src/m/classes/model.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.m	(revision 27852)
+++ /issm/trunk-jpl/src/m/classes/model.m	(revision 27853)
@@ -1289,4 +1289,5 @@
 			if isfield(structmd,'cfl_coefficient'), md.timestepping.cfl_coefficient=structmd.cfl_coefficient; end
 			if isfield(structmd,'spcthickness'), md.masstransport.spcthickness=structmd.spcthickness; end
+			if isfield(structmd,'spcthickness'), md.debris.spcthickness=structmd.spcthickness; end
 			if isfield(structmd,'artificial_diffusivity'), md.masstransport.stabilization=structmd.artificial_diffusivity; end
 			if isfield(structmd,'hydrostatic_adjustment'), md.masstransport.hydrostatic_adjustment=structmd.hydrostatic_adjustment; end
Index: /issm/trunk-jpl/src/m/solve/parseresultsfromdisk.m
===================================================================
--- /issm/trunk-jpl/src/m/solve/parseresultsfromdisk.m	(revision 27852)
+++ /issm/trunk-jpl/src/m/solve/parseresultsfromdisk.m	(revision 27853)
@@ -240,4 +240,8 @@
 		field = field*yts;
 	elseif strcmp(fieldname,'VyAverage'),
+		field = field*yts;
+	elseif strcmp(fieldname,'VxDebris'),
+                field = field*yts;
+        elseif strcmp(fieldname,'VyDebris'),
 		field = field*yts;
 	elseif strcmp(fieldname,'BasalforcingsGroundediceMeltingRate'),