Index: /issm/trunk-jpl/src/m/classes/frictionweertman.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/frictionweertman.m	(revision 27492)
+++ /issm/trunk-jpl/src/m/classes/frictionweertman.m	(revision 27493)
@@ -8,4 +8,5 @@
 		C = NaN;
 		m = NaN;
+		linearize  = 0;
 	end
    methods (Static)
@@ -29,4 +30,6 @@
 		function self = setdefaultparameters(self) % {{{
 
+			self.linearize = 0;
+
 		end % }}}
 		function md = checkconsistency(self,md,solution,analyses) % {{{
@@ -36,4 +39,5 @@
 			md = checkfield(md,'fieldname','friction.C','timeseries',1,'NaN',1,'Inf',1);
 			md = checkfield(md,'fieldname','friction.m','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
+			md = checkfield(md,'fieldname','friction.linearize','numel',[1],'values',[0:2]);
 		end % }}}
 		function disp(self) % {{{
@@ -47,4 +51,5 @@
 			fielddisplay(self,'C','friction coefficient [SI]');
 			fielddisplay(self,'m','m exponent');
+			fielddisplay(self,'linearize','0: not linearized, 1: interpolated linearly, 2: constant per element (default is 0)');
 		end % }}}
 		function marshall(self,prefix,md,fid) % {{{
@@ -54,6 +59,5 @@
 			WriteData(fid,prefix,'class','friction','object',self,'fieldname','C','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
 			WriteData(fid,prefix,'class','friction','object',self,'fieldname','m','format','DoubleMat','mattype',2);
-			
-
+			WriteData(fid,prefix,'class','friction','object',self,'fieldname','linearize','format','Integer');
 		end % }}}
 	end
Index: /issm/trunk-jpl/src/m/classes/frictionweertman.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/frictionweertman.py	(revision 27492)
+++ /issm/trunk-jpl/src/m/classes/frictionweertman.py	(revision 27493)
@@ -15,4 +15,5 @@
         self.C = float('NaN')
         self.m = float('NaN')
+        self.linearize = 0
 
     #set defaults
@@ -26,8 +27,10 @@
         string = "%s\n%s" % (string, fielddisplay(self, "C", "friction coefficient [SI]"))
         string = "%s\n%s" % (string, fielddisplay(self, "m", "m exponent"))
+        string = "%s\n%s" % (string, fielddisplay(self, "linearize", "0: not linearized, 1: interpolated linearly, 2: constant per element (default is 0)"))
         return string
     # }}}
 
     def setdefaultparameters(self):  # {{{
+        self.linearize = 0
         return self
     # }}}
@@ -41,4 +44,5 @@
         md = checkfield(md, 'fieldname', 'friction.C', 'timeseries', 1, 'NaN', 1, 'Inf', 1)
         md = checkfield(md, 'fieldname', 'friction.m', 'NaN', 1, 'Inf', 1, 'size', [md.mesh.numberofelements])
+        md = checkfield(md, 'fieldname', 'friction.linearize', 'numel', [1], 'values', [0, 1, 2])
 
         return md
@@ -49,3 +53,4 @@
         WriteData(fid, prefix, 'class', 'friction', 'object', self, 'fieldname', 'C', 'format', 'DoubleMat', 'mattype', 1)
         WriteData(fid, prefix, 'class', 'friction', 'object', self, 'fieldname', 'm', 'format', 'DoubleMat', 'mattype', 2)
+        WriteData(fid, prefix, 'class', 'friction', 'object', self, 'fieldname', 'linearize', 'format', 'Integer')
     # }}}