Index: /issm/trunk-jpl/src/c/Makefile.am
===================================================================
--- /issm/trunk-jpl/src/c/Makefile.am	(revision 20624)
+++ /issm/trunk-jpl/src/c/Makefile.am	(revision 20625)
@@ -116,4 +116,5 @@
 					./shared/Sorting/binary_search.cpp\
 					./shared/Elements/Cuffey.cpp\
+					./shared/Elements/CuffeyTemperate.cpp\
 					./shared/Elements/StressIntensityIntegralWeight.cpp\
 					./shared/Elements/Paterson.cpp\
Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 20624)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 20625)
@@ -1493,8 +1493,8 @@
 				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]);
Index: /issm/trunk-jpl/src/c/shared/Elements/elements.h
===================================================================
--- /issm/trunk-jpl/src/c/shared/Elements/elements.h	(revision 20624)
+++ /issm/trunk-jpl/src/c/shared/Elements/elements.h	(revision 20625)
@@ -9,5 +9,5 @@
 
 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);