Index: /issm/trunk-jpl/src/c/shared/Elements/CuffeyTemperate.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Elements/CuffeyTemperate.cpp	(revision 20625)
+++ /issm/trunk-jpl/src/c/shared/Elements/CuffeyTemperate.cpp	(revision 20626)
@@ -7,12 +7,9 @@
 #include <math.h>
 
-#include "../Numerics/types.h"
+#include "./elements.h"
 
 IssmDouble CuffeyTemperate(IssmDouble temperature, IssmDouble waterfraction, IssmDouble stressexp){
 
-	/*output: */
-	IssmDouble B;
-	B=Cuffey(temperature)*pow(1+181.25*max(0., min(0.01, waterfraction)), -1./stressexp); 
+	return Cuffey(temperature)*pow(1+181.25*max(0., min(0.01, waterfraction)), -1./stressexp); 
 
-	return B;
 }