Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 21376)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 21377)
@@ -1497,4 +1497,8 @@
 				/*Do nothing: B is not temperature dependent*/
 				break;
+			case BuddJackaEnum:
+				for(i=0;i<numnodes;i++) B[i]=BuddJacka(temperature[i]);
+				element->AddInput(MaterialsRheologyBEnum,&B[0],element->GetElementType());
+				break;
 			case CuffeyEnum:
 				for(i=0;i<numnodes;i++) B[i]=Cuffey(temperature[i]);
Index: /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp	(revision 21376)
+++ /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp	(revision 21377)
@@ -768,4 +768,8 @@
 				/*Do nothing: B is not temperature dependent*/
 				break;
+			case BuddJackaEnum:
+				for(i=0;i<numnodes;i++) B[i]=BuddJacka(values[i]);
+				element->AddInput(MaterialsRheologyBEnum,&B[0],element->GetElementType());
+				break;
 			case CuffeyEnum:
 				for(i=0;i<numnodes;i++) B[i]=Cuffey(values[i]);
Index: /issm/trunk-jpl/src/m/classes/matestar.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/matestar.m	(revision 21376)
+++ /issm/trunk-jpl/src/m/classes/matestar.m	(revision 21377)
@@ -124,5 +124,5 @@
 			md = checkfield(md,'fieldname','materials.rheology_Ec','>',0,'size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
 			md = checkfield(md,'fieldname','materials.rheology_Es','>',0,'size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
-			md = checkfield(md,'fieldname','materials.rheology_law','values',{'None' 'Cuffey' 'CuffeyTemperate' 'Paterson' 'Arrhenius' 'LliboutryDuval'});
+			md = checkfield(md,'fieldname','materials.rheology_law','values',{'None' 'BuddJacka' 'Cuffey' 'CuffeyTemperate' 'Paterson' 'Arrhenius' 'LliboutryDuval'});
 
 			if ismember('GiaAnalysis',analyses),
@@ -155,5 +155,5 @@
 			fielddisplay(self,'rheology_Ec','compressive enhancement factor');
 			fielddisplay(self,'rheology_Es','shear enhancement factor');
-			fielddisplay(self,'rheology_law',['law for the temperature dependance of the rheology: ''None'', ''Cuffey'', ''CuffeyTemperate'', ''Paterson'', ''Arrhenius'' or ''LliboutryDuval''']);
+			fielddisplay(self,'rheology_law',['law for the temperature dependance of the rheology: ''None'', ''BuddJacka'', ''Cuffey'', ''CuffeyTemperate'', ''Paterson'', ''Arrhenius'' or ''LliboutryDuval''']);
 			fielddisplay(self,'lithosphere_shear_modulus','Lithosphere shear modulus [Pa]');
 			fielddisplay(self,'lithosphere_density','Lithosphere density [g/cm^-3]');
Index: /issm/trunk-jpl/src/m/classes/matice.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/matice.m	(revision 21376)
+++ /issm/trunk-jpl/src/m/classes/matice.m	(revision 21377)
@@ -114,5 +114,5 @@
 			md = checkfield(md,'fieldname','materials.rheology_B','>',0,'timeseries',1,'NaN',1,'Inf',1);
 			md = checkfield(md,'fieldname','materials.rheology_n','>',0,'size',[md.mesh.numberofelements 1]);
-			md = checkfield(md,'fieldname','materials.rheology_law','values',{'None' 'Cuffey' 'CuffeyTemperate' 'Paterson' 'Arrhenius' 'LliboutryDuval'});
+			md = checkfield(md,'fieldname','materials.rheology_law','values',{'None' 'BuddJacka' 'Cuffey' 'CuffeyTemperate' 'Paterson' 'Arrhenius' 'LliboutryDuval'});
 
 			if ismember('GiaAnalysis',analyses),
@@ -144,5 +144,5 @@
 			fielddisplay(self,'rheology_B','flow law parameter [Pa/s^(1/n)]');
 			fielddisplay(self,'rheology_n','Glen''s flow law exponent');
-			fielddisplay(self,'rheology_law',['law for the temperature dependance of the rheology: ''None'', ''Cuffey'', ''CuffeyTemperate'', ''Paterson'', ''Arrhenius'' or ''LliboutryDuval''']);
+			fielddisplay(self,'rheology_law',['law for the temperature dependance of the rheology: ''None'', ''BuddJacka'', Cuffey'', ''CuffeyTemperate'', ''Paterson'', ''Arrhenius'' or ''LliboutryDuval''']);
 			fielddisplay(self,'lithosphere_shear_modulus','Lithosphere shear modulus [Pa]');
 			fielddisplay(self,'lithosphere_density','Lithosphere density [g/cm^-3]');