Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 21599)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 21600)
@@ -245,4 +245,5 @@
 	if(frictionlaw==4 || frictionlaw==6) parameters->AddObject(iomodel->CopyConstantObject("md.friction.gamma",FrictionGammaEnum));
 	if(frictionlaw==3) parameters->AddObject(iomodel->CopyConstantObject("md.friction.coupling",FrictionCouplingEnum));
+	if(frictionlaw==9) parameters->AddObject(iomodel->CopyConstantObject("md.friction.gamma",FrictionGammaEnum));
 }/*}}}*/
 
Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 21599)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 21600)
@@ -895,4 +895,5 @@
 	if(frictionlaw==3) parameters->AddObject(iomodel->CopyConstantObject("md.friction.coupling",FrictionCouplingEnum));
 	if(frictionlaw==5) parameters->AddObject(iomodel->CopyConstantObject("md.friction.f",FrictionFEnum));
+	if(frictionlaw==9) parameters->AddObject(iomodel->CopyConstantObject("md.friction.gamma",FrictionGammaEnum));
 
 }/*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 21599)
+++ /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 21600)
@@ -493,11 +493,12 @@
 	/*Compute gamma*/
 	element->parameters->FindParam(&time,TimeEnum);
-	if(time<25e3*yts){
-		gamma = 10.;
-	}
-	else{
-		gamma = 5.;
-	}
-	gamma = 5.;
+	element->parameters->FindParam(&gamma,FrictionGammaEnum);
+	//if(time<25e3*yts){
+	//	gamma = 10.;
+	//}
+	//else{
+	//	gamma = 5.;
+	//}
+	//gamma = 5.;
 
 	/*Compute scaling parameter*/
Index: /issm/trunk-jpl/src/m/classes/frictionjosh.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/frictionjosh.m	(revision 21599)
+++ /issm/trunk-jpl/src/m/classes/frictionjosh.m	(revision 21600)
@@ -8,4 +8,5 @@
 		coefficient = NaN;
 		pressure_adjusted_temperature = NaN;
+		gamma      = 0.;
 	end
 	methods
@@ -26,4 +27,7 @@
 		function self = setdefaultparameters(self) % {{{
 
+			%Default gamma: 1
+			self.gamma = 1.;
+
 		end % }}}
 		function md = checkconsistency(self,md,solution,analyses) % {{{
@@ -34,4 +38,5 @@
 			md = checkfield(md,'fieldname','friction.coefficient','NaN',1,'Inf',1);
 			md = checkfield(md,'fieldname','friction.pressure_adjusted_temperature','NaN',1,'Inf',1);
+			md = checkfield(md,'fieldname','friction.gamma','numel',1,'NaN',1,'Inf',1,'>',0.);
 
 			%Check that temperature is provided
@@ -42,4 +47,5 @@
 			fielddisplay(self,'coefficient','friction coefficient [SI]');
 			fielddisplay(self,'pressure_adjusted_temperature','friction pressure_adjusted_temperature (T - Tpmp) [K]');
+			fielddisplay(self,'gamma','(T - Tpmp)/gamma [K]');
 		end % }}}
 		function marshall(self,prefix,md,fid) % {{{
@@ -48,4 +54,5 @@
 			WriteData(fid,prefix,'class','friction','object',self,'fieldname','coefficient','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
 			WriteData(fid,prefix,'class','friction','object',self,'fieldname','pressure_adjusted_temperature','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
+			WriteData(fid,prefix,'class','friction','object',self,'fieldname','gamma','format','Double');
 		end % }}}
 	end