Index: /issm/trunk-jpl/src/m/classes/debris.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/debris.m	(revision 27854)
+++ /issm/trunk-jpl/src/m/classes/debris.m	(revision 27855)
@@ -110,7 +110,7 @@
 			md = checkfield(md,'fieldname','debris.removalmodel','values',[0 1 2]);
 			md = checkfield(md,'fieldname','debris.displacementmodel','values',[0 1 2]);
+			md = checkfield(md,'fieldname','debris.max_displacementvelocity','>=',0);
 			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)),
@@ -125,6 +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: additional debris velocity above the critical slope/stress threshold');
 			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)');
Index: /issm/trunk-jpl/src/m/classes/debris.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/debris.py	(revision 27854)
+++ /issm/trunk-jpl/src/m/classes/debris.py	(revision 27855)
@@ -21,4 +21,5 @@
         self.removalmodel = 0
         self.displacementmodel = 0
+        self.max_displacementvelocity = 0
         self.removal_slope_threshold = 0
         self.removal_stress_threshold = 0
@@ -44,4 +45,5 @@
         s += '{}\n'.format(fielddisplay(self,'removalmodel','frontal removal of debris. 0: no removal, 1: Slope-triggered debris removal, 2: driving-stress triggered debris removal'))
         s += '{}\n'.format(fielddisplay(self,'displacementmodel','debris displacement. 0: no displacement, 1: ...'))
+        s += '{}\n'.format(fielddisplay(self,'max_displacementvelocity','maximum velocity of debris transport (v_ice + v_displacement) (m/a)'))
         s += '{}\n'.format(fielddisplay(self,'removal_slope_threshold','critical slope (degrees) for removalmodel (1)'))
         s += '{}\n'.format(fielddisplay(self,'removal_stress_threshold','critical stress (Pa) for removalmodel (2)'))
@@ -54,5 +56,5 @@
 
     def defaultoutputs(self, md):  # {{{
-        return ['DebrisThickness']
+        return ['DebrisThickness', 'DebrisMaskNodeActivation', 'VxDebris', 'VyDebris']
     # }}}
 
@@ -79,4 +81,7 @@
         self.removal_stress_threshold = 0
 
+        # Max velocity for displacementmodel (1)
+        self.max_displacementvelocity = 0
+
         # Default output
         self.requested_outputs = ['default']
@@ -90,9 +95,10 @@
 
         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)
         md = checkfield(md, 'fieldname', 'debris.removalmodel', 'values', [0, 1, 2])
         md = checkfield(md, 'fieldname', 'debris.displacementmodel', 'values', [0, 1, 2])
+        md = checkfield(md, 'fieldname', 'debris.max_displacementvelocity', '>=', 0)
         md = checkfield(md, 'fieldname', 'debris.removal_slope_threshold', '>=', 0)
         md = checkfield(md, 'fieldname', 'debris.removal_stress_threshold', '>=', 0)
@@ -110,4 +116,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')
Index: /issm/trunk-jpl/src/m/solve/parseresultsfromdisk.m
===================================================================
--- /issm/trunk-jpl/src/m/solve/parseresultsfromdisk.m	(revision 27854)
+++ /issm/trunk-jpl/src/m/solve/parseresultsfromdisk.m	(revision 27855)
@@ -242,6 +242,6 @@
 		field = field*yts;
 	elseif strcmp(fieldname,'VxDebris'),
-                field = field*yts;
-        elseif strcmp(fieldname,'VyDebris'),
+		field = field*yts;
+	elseif strcmp(fieldname,'VyDebris'),
 		field = field*yts;
 	elseif strcmp(fieldname,'BasalforcingsGroundediceMeltingRate'),
Index: /issm/trunk-jpl/src/m/solve/parseresultsfromdisk.py
===================================================================
--- /issm/trunk-jpl/src/m/solve/parseresultsfromdisk.py	(revision 27854)
+++ /issm/trunk-jpl/src/m/solve/parseresultsfromdisk.py	(revision 27855)
@@ -191,4 +191,8 @@
             field = field * yts
         elif fieldname == 'VyAverage':
+            field = field * yts
+        elif fieldname == 'VxDebris':
+            field = field * yts
+        elif fieldname == 'VyDebris':
             field = field * yts
         elif fieldname == 'BasalforcingsGroundediceMeltingRate':