Index: /issm/trunk-jpl/src/c/modules/PostprocessingEnthalpyx/PostprocessingEnthalpyx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/PostprocessingEnthalpyx/PostprocessingEnthalpyx.cpp	(revision 16359)
+++ /issm/trunk-jpl/src/c/modules/PostprocessingEnthalpyx/PostprocessingEnthalpyx.cpp	(revision 16360)
@@ -19,8 +19,9 @@
 
 	/*drain excess water fraction: */
-	for (i=0;i<femmodel->elements->Size();i++){
+	/*for (i=0;i<femmodel->elements->Size();i++){
 		element=dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));
 		element->DrainWaterfraction();
 	}
+	*/
 
 	/*Update basal dirichlet BCs for enthalpy: */
Index: /issm/trunk-jpl/src/c/shared/Elements/DrainageFunctionWaterfraction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Elements/DrainageFunctionWaterfraction.cpp	(revision 16359)
+++ /issm/trunk-jpl/src/c/shared/Elements/DrainageFunctionWaterfraction.cpp	(revision 16360)
@@ -9,37 +9,37 @@
 /*FUNCTION IssmDouble DrainageFunctionWaterfraction()*/
 IssmDouble DrainageFunctionWaterfraction(IssmDouble waterfraction, IssmDouble dt=0.){
-    /* DrainageFunctionWaterfraction returns how much of the waterfraction is drained per year */
-    _assert_(waterfraction>=0.);
-    _assert_(dt>=0.);
+	/* DrainageFunctionWaterfraction returns how much of the waterfraction is drained per year */
+	_assert_(waterfraction>=0.);
+	_assert_(dt>=0.);
 
-    IssmDouble w0=0.01, w1=0.02, w2=0.03;
-    IssmDouble Dret, D0=0, D1=0.005, D2=0.05;
-		IssmDouble yts=365*24*60*60;
-		dt/=yts;
+	IssmDouble w0=0.01, w1=0.02, w2=0.03;
+	IssmDouble Dret, D0=0, D1=0.005, D2=0.05;
+	IssmDouble yts=365*24*60*60;
+	dt/=yts;
 
-    /*get drainage function value*/
-    if((w0==w1)||(w1==w2)||(w0==w2))
-        _error_("Error: equal ordinates in DrainageFunctionWaterfraction -> division by zero. Abort");
-    if(waterfraction<=w0)
-        Dret=D0;
-    if((waterfraction>w0) && (waterfraction<=w1))
-        Dret=(D1-D0)/(w1-w0)*(waterfraction-w0)+D0;
-    if((waterfraction>w1) && (waterfraction<=w2))
-        Dret=(D2-D1)/(w2-w1)*(waterfraction-w1)+D1;
-    else 
-        Dret=D2;
+	/*get drainage function value*/
+	if((w0==w1)||(w1==w2)||(w0==w2))
+		_error_("Error: equal ordinates in DrainageFunctionWaterfraction -> division by zero. Abort");
+	if(waterfraction<=w0)
+		Dret=D0;
+	if((waterfraction>w0) && (waterfraction<=w1))
+		Dret=(D1-D0)/(w1-w0)*(waterfraction-w0)+D0;
+	if((waterfraction>w1) && (waterfraction<=w2))
+		Dret=(D2-D1)/(w2-w1)*(waterfraction-w1)+D1;
+	else 
+		Dret=D2;
 
-    /*check if dt*Dret>waterfraction. If so, drain whole waterfraction*/
-    if(dt==0.){
-      if(Dret>waterfraction)
-        return waterfraction;
-      else
-        return Dret;
-    }
-    else{
-      if(dt*Dret>waterfraction)
-        return waterfraction;
-      else
-        return dt*Dret;
-    }
+	/*check if dt*Dret>waterfraction. If so, drain whole waterfraction*/
+	if(dt==0.){
+		if(Dret>waterfraction)
+			return waterfraction;
+		else
+			return Dret;
+	}
+	else{
+		if(dt*Dret>waterfraction)
+			return waterfraction;
+		else
+			return dt*Dret;
+	}
 }