Index: /issm/trunk-jpl/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/objects/Elements/Penta.cpp	(revision 11362)
+++ /issm/trunk-jpl/src/c/objects/Elements/Penta.cpp	(revision 11363)
@@ -3903,16 +3903,16 @@
 		enthalpy_input->GetInputValue(&enthalpy,gauss);
 		pressure_input->GetInputValue(&pressure,gauss);
-		if(enthalpy>matpar->PureIceEnthalpy(pressure)){
-			enthalpy_input->GetInputValue(&enthalpyup,gaussup);
-			pressure_input->GetInputValue(&pressureup,gaussup);
-			if(enthalpyup>matpar->PureIceEnthalpy(pressureup)){
-				//do nothing, don't add heatflux
-			}
-			else{
-				//need to change spcenthalpy according to Aschwanden 
-				//NEED TO UPDATE
-			}
-		}
-		else{
+//		if(enthalpy>matpar->PureIceEnthalpy(pressure)){
+//			enthalpy_input->GetInputValue(&enthalpyup,gaussup);
+//			pressure_input->GetInputValue(&pressureup,gaussup);
+//			if(enthalpyup>matpar->PureIceEnthalpy(pressureup)){
+//				//do nothing, don't add heatflux
+//			}
+//			else{
+//				//need to change spcenthalpy according to Aschwanden 
+//				//NEED TO UPDATE
+//			}
+//		}
+//		else{
 			geothermalflux_input->GetInputValue(&geothermalflux_value,gauss);
 			friction->GetAlpha2(&alpha2,gauss,VxEnum,VyEnum,VzEnum);
@@ -3925,5 +3925,5 @@
 
 			for(i=0;i<numdof;i++) pe->values[i]+=scalar*basis[i];
-		}
+//		}
 	}