Index: /issm/trunk-jpl/src/c/modules/PostprocessingEnthalpyx/PostprocessingEnthalpyx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/PostprocessingEnthalpyx/PostprocessingEnthalpyx.cpp	(revision 16027)
+++ /issm/trunk-jpl/src/c/modules/PostprocessingEnthalpyx/PostprocessingEnthalpyx.cpp	(revision 16027)
@@ -0,0 +1,31 @@
+/*!\file PostprocessingEnthalpyx.cpp
+ * \brief: execute processing of variables after computation: water drainage, subglacial melting rates, update of thermal BCs
+ */
+
+#include "./PostprocessingEnthalpyx.h"
+#include "../../shared/shared.h"
+#include "../../toolkits/toolkits.h"
+
+void PostprocessingEnthalpyx(FemModel* femmodel){
+    /*Intermediary*/
+	int i;
+	Element* element=NULL;
+
+    /*drain excess water fraction: */
+for (i=0;i<femmodel->elements->Size();i++){
+        element=dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));
+        element->DrainWaterfraction();
+	}
+
+    /*Compute basal melting rates: */
+for (i=0;i<femmodel->elements->Size();i++){
+        element=dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));
+        element->ComputeBasalMeltingrate();
+	}
+
+    /*Update basal dirichlet BCs for enthalpy: */
+for (i=0;i<femmodel->elements->Size();i++){
+        element=dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));
+        element->UpdateThermalBasalConstraints();
+	}
+}
Index: /issm/trunk-jpl/src/c/modules/PostprocessingEnthalpyx/PostprocessingEnthalpyx.h
===================================================================
--- /issm/trunk-jpl/src/c/modules/PostprocessingEnthalpyx/PostprocessingEnthalpyx.h	(revision 16027)
+++ /issm/trunk-jpl/src/c/modules/PostprocessingEnthalpyx/PostprocessingEnthalpyx.h	(revision 16027)
@@ -0,0 +1,13 @@
+/*!\file PostprocessingEnthalpy.h
+ * \brief: header file for PostprocessingEnthalpy
+ */
+
+#ifndef _POSTPROCESSINGENTHALPYX_H
+#define _POSTPROCESSINGENTHALPYX_H
+
+#include "../../classes/classes.h"
+#include "../../toolkits/toolkits.h"
+
+void PostprocessingEnthalpyx(FemModel* femmodel);
+
+#endif  /* _POSTPROCESSINGENTHALPYX_H */