Index: /issm/trunk-jpl/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/objects/Elements/Penta.cpp	(revision 11313)
+++ /issm/trunk-jpl/src/c/objects/Elements/Penta.cpp	(revision 11314)
@@ -4271,5 +4271,9 @@
 	if(converged){
 		/*Convert enthalpy into temperature and water fraction*/
-		for(i=0;i<numdof;i++) matpar->EnthalpyToThermal(&temperatures[i],&waterfraction[i],values[i],pressure[i]);
+		for(i=0;i<numdof;i++){
+			matpar->EnthalpyToThermal(&temperatures[i],&waterfraction[i],values[i],pressure[i]);
+			if(waterfraction[i]<0) _error_("Negative water fraction found in solution vector");
+			if(waterfraction[i]>1) _error_("Water fraction >1 found in solution vector");
+		}
 			
 		this->inputs->AddInput(new PentaP1Input(EnthalpyEnum,values));