source: issm/oecreview/Archive/18296-19100/ISSM-18590-18591.diff

../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp

@@ -3 +3 @@
 #include "../classes/classes.h"
 #include "../shared/shared.h"
 #include "../modules/modules.h"
+#include "../solutionsequences/solutionsequences.h"
 
 /*Model processing*/
 int  EnthalpyAnalysis::DofsPerNode(int** doflist,int domaintype,int approximation){/*{{{*/
@@ -198 +199 @@
 
 /*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){/*{{{*/
         /*Default, return NULL*/

../trunk-jpl/src/c/cores/thermal_core.cpp

@@ -27 +27 @@
         if(numoutputs) femmodel->parameters->FindParam(&requested_outputs,&numoutputs,ThermalRequestedOutputsEnum);
 
         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{
                 if(VerboseSolution()) _printf0_("   computing temperatures\n");