Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 18590)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 18591)
@@ -4,4 +4,5 @@
 #include "../shared/shared.h"
 #include "../modules/modules.h"
+#include "../solutionsequences/solutionsequences.h"
 
 /*Model processing*/
@@ -199,5 +200,17 @@
 /*Finite Element Analysis*/
 void EnthalpyAnalysis::Core(FemModel* femmodel){/*{{{*/
-	_error_("not implemented");
+	if(VerboseSolution()) _printf0_("   computing enthalpy\n");
+	femmodel->SetCurrentConfiguration(EnthalpyAnalysisEnum);
+	solutionsequence_thermal_nonlinear(femmodel);
+
+	/*transfer enthalpy to enthalpy picard for the next step: */
+	InputDuplicatex(femmodel,EnthalpyEnum,EnthalpyPicardEnum);
+
+	int solution_type;
+	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
+	if(solution_type!=SteadystateSolutionEnum){
+		PostProcessing(femmodel);
+	}
+
 }/*}}}*/
 ElementVector* EnthalpyAnalysis::CreateDVector(Element* element){/*{{{*/
Index: /issm/trunk-jpl/src/c/cores/thermal_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/thermal_core.cpp	(revision 18590)
+++ /issm/trunk-jpl/src/c/cores/thermal_core.cpp	(revision 18591)
@@ -28,18 +28,7 @@
 
 	if(isenthalpy){
-
-		if(VerboseSolution()) _printf0_("   computing enthalpy\n");
-		femmodel->SetCurrentConfiguration(EnthalpyAnalysisEnum);
-		solutionsequence_thermal_nonlinear(femmodel);
-
-		/*transfer enthalpy to enthalpy picard for the next step: */
-		InputDuplicatex(femmodel,EnthalpyEnum,EnthalpyPicardEnum);
-
-		if(solution_type!=SteadystateSolutionEnum){
-			/*Post process*/
-			enthalpy_analysis = new EnthalpyAnalysis();
-			enthalpy_analysis->PostProcessing(femmodel);
-			delete enthalpy_analysis;
-		}
+		enthalpy_analysis = new EnthalpyAnalysis();
+		enthalpy_analysis->Core(femmodel);
+		delete enthalpy_analysis;
 	}
 	else{