Index: ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 20623)
+++ ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 20624)
@@ -1492,10 +1492,10 @@
 				for(i=0;i<numnodes;i++) B[i]=Cuffey(temperature[i]);
 				element->AddInput(MaterialsRheologyBEnum,&B[0],element->GetElementType());
 				break;
-			case CuffeyTemperateEnum:
-				for(i=0;i<numnodes;i++) B[i]=CuffeyTemperate(temperature[i], waterfraction[i], element->GetMaterialParameter(MaterialsRheologyNEnum));
-				element->AddInput(MaterialsRheologyBEnum,&B[0],element->GetElementType());
-				break;
+// 			case CuffeyTemperateEnum:
+// 				for(i=0;i<numnodes;i++) B[i]=CuffeyTemperate(temperature[i], waterfraction[i], element->GetMaterialParameter(MaterialsRheologyNEnum));
+// 				element->AddInput(MaterialsRheologyBEnum,&B[0],element->GetElementType());
+// 				break;
 			case PatersonEnum:
 				for(i=0;i<numnodes;i++) B[i]=Paterson(temperature[i]);
 				element->AddInput(MaterialsRheologyBEnum,&B[0],element->GetElementType());
Index: ../trunk-jpl/src/c/shared/Elements/elements.h
===================================================================
--- ../trunk-jpl/src/c/shared/Elements/elements.h	(revision 20623)
+++ ../trunk-jpl/src/c/shared/Elements/elements.h	(revision 20624)
@@ -8,7 +8,7 @@
 #include "../Numerics/types.h"
 
 IssmDouble Cuffey(IssmDouble temperature);
-IssmDouble CuffeyTemperate(IssmDouble temperature, IssmDouble waterfraction, IssmDouble stressexp);
+// IssmDouble CuffeyTemperate(IssmDouble temperature, IssmDouble waterfraction, IssmDouble stressexp);
 IssmDouble Paterson(IssmDouble temperature);
 IssmDouble Arrhenius(IssmDouble temperature,IssmDouble depth,IssmDouble n);
 IssmDouble LliboutryDuval(IssmDouble enthalpy, IssmDouble pressure, IssmDouble n, IssmDouble betaCC, IssmDouble referencetemperature, IssmDouble heatcapacity, IssmDouble latentheat);