Index: /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp	(revision 14803)
+++ /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp	(revision 14804)
@@ -6348,7 +6348,5 @@
 		/*Get the flag to the limitation method*/
 		this->parameters->FindParam(&hmax_flag,HydrologydcSedimentlimitFlagEnum);
-
-		if(this->id==1) printf("flag = %i\n",hmax_flag);
-
+		
 		/*Switch between the different cases*/
 		switch(hmax_flag){
Index: /issm/trunk-jpl/src/c/classes/objects/Loads/Pengrid.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Loads/Pengrid.cpp	(revision 14803)
+++ /issm/trunk-jpl/src/c/classes/objects/Loads/Pengrid.cpp	(revision 14804)
@@ -708,5 +708,4 @@
 	element->GetInputValue(&h,node,SedimentHeadEnum);
 	element->GetHydrologyDCInefficientHmax(&h_max,node);
-	/*printf("hmax %g \n",h_max);*/
 
 	if (h>h_max)
@@ -754,7 +753,6 @@
 	ElementVector* pe=new ElementVector(&node,1,this->parameters);
 
-	/*Get h_max*/
-	h_max = 10000.;
-
+	/*Get h_max and compute penalty*/
+	element->GetHydrologyDCInefficientHmax(&h_max,node);
 	pe->values[0]=kmax*pow(10.,penalty_factor)*h_max;