Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 17044)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 17045)
@@ -774,5 +774,5 @@
 			element->EnthalpyToThermal(&temperature[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");
+			//if(waterfraction[i]>1.) _error_("Water fraction >1 found in solution vector");
 		}
 		element->AddInput(EnthalpyEnum,values,P1Enum);
@@ -948,5 +948,5 @@
 
 	/******** DRAINAGE *****************************************/
-	IssmDouble* drainrate_column  = xNew<IssmDouble>(numsegments); //TODO: xDelete?
+	IssmDouble* drainrate_column  = xNew<IssmDouble>(numsegments);
 	IssmDouble* drainrate_element = xNew<IssmDouble>(numsegments);
 	for(is=0;is<numsegments;is++)	drainrate_column[is]=0.;
@@ -1116,8 +1116,8 @@
 			setspc=((enthalpyup<PureIceEnthalpy(element,pressureup)) && (watercolumn>=0.))?true:false;
 		}
-		else if(watercolumn>0.) // case H<Hpmp && watercolumn>0.
-			setspc=true;
 		else
 			setspc = false;
+
+
 
 		node=element->GetNode(indices[i]);
Index: /issm/trunk-jpl/src/c/cores/thermal_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/thermal_core.cpp	(revision 17044)
+++ /issm/trunk-jpl/src/c/cores/thermal_core.cpp	(revision 17045)
@@ -48,8 +48,10 @@
 		InputDuplicatex(femmodel,EnthalpyEnum,EnthalpyPicardEnum);
 
-		/*Post process*/
-		enthalpy_analysis = new EnthalpyAnalysis();
-		enthalpy_analysis->PostProcessing(femmodel);
-		delete enthalpy_analysis;
+		if(solution_type!=SteadystateSolutionEnum){
+			/*Post process*/
+			enthalpy_analysis = new EnthalpyAnalysis();
+			enthalpy_analysis->PostProcessing(femmodel);
+			delete enthalpy_analysis;
+		}
 	}
 	else{