Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15896)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15897)
@@ -5281,5 +5281,5 @@
 	for(int iv=0;iv<3;iv++){
 		gauss->GaussVertex(iv);
-		gaussup->GaussNode(iv+3);
+		gaussup->GaussVertex(iv+3);
 
         // TODO: make sure that no node is computed twice (insert mask)
@@ -5350,5 +5350,5 @@
     /*Intermediaries*/
     bool isenthalpy;
-    IssmDouble waterfraction_array[NUMVERTICES], temperature[NUMVERTICES];
+    IssmDouble waterfraction[NUMVERTICES], temperature[NUMVERTICES];
     IssmDouble enthalpy[NUMVERTICES], pressure[NUMVERTICES]; 
     IssmDouble latentheat, dt;
@@ -5371,5 +5371,5 @@
         enthalpy_input->GetInputValue(&enthalpy[iv], gauss);
         pressure_input->GetInputValue(&pressure[iv], gauss);
-        matpar->EnthalpyToThermal(&temperature[iv], &waterfraction[iv], enthalpy[iv],pressure[iv]); 
+        matpar->EnthalpyToThermal(&temperature[iv],&waterfraction[iv], enthalpy[iv],pressure[iv]); 
     
         /*drain water fraction*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 15896)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 15897)
@@ -359,5 +359,5 @@
 		void           UpdateThermalBasalConstraints(void);
 		void           ComputeBasalMeltrate(void);
-        void           DrainWaterfraction(void);
+		void           DrainWaterfraction(void);
 		#endif
 
Index: /issm/trunk-jpl/src/c/shared/Elements/elements.h
===================================================================
--- /issm/trunk-jpl/src/c/shared/Elements/elements.h	(revision 15896)
+++ /issm/trunk-jpl/src/c/shared/Elements/elements.h	(revision 15897)
@@ -19,5 +19,5 @@
 				     IssmDouble* PrecipitationsPresentday,
 				     IssmDouble* TemperaturesLgm, IssmDouble* TemperaturesPresentday, 
-					     IssmDouble* monthlytemperaturesout, IssmDouble* monthlyprecout);
+					  IssmDouble* monthlytemperaturesout, IssmDouble* monthlyprecout);
 IssmDouble DrainageFunctionWaterfraction(IssmDouble waterfraction);