thermal_core.cpp File Reference

: core of the thermal solution More...

#include "./cores.h"
#include "../toolkits/toolkits.h"
#include "../classes/classes.h"
#include "../analyses/analyses.h"
#include "../shared/shared.h"
#include "../modules/modules.h"
#include "../solutionsequences/solutionsequences.h"

Go to the source code of this file.

Functions
void	thermal_core (FemModel *femmodel)

Detailed Description

: core of the thermal solution

Definition in file thermal_core.cpp.

Function Documentation

thermal_core()

void thermal_core

(

FemModel *

femmodel

)

Definition at line 13 of file thermal_core.cpp.

                                     {
 
   /*Start profiler*/
   femmodel->profiler->Start(THERMALCORE);
 
   /*intermediary*/
   bool   save_results,isenthalpy;
   bool   dakota_analysis;
   int    solution_type,numoutputs;
   char** requested_outputs = NULL;
   EnthalpyAnalysis * enthalpy_analysis = NULL;
 
   /*first recover parameters common to all solutions*/
   femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
   femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
   femmodel->parameters->FindParam(&isenthalpy,ThermalIsenthalpyEnum);
   femmodel->parameters->FindParam(&numoutputs,ThermalNumRequestedOutputsEnum);
   if(numoutputs) femmodel->parameters->FindParam(&requested_outputs,&numoutputs,ThermalRequestedOutputsEnum);
 
   /*Calculate geothermalflux*/
   GeothermalFluxx(femmodel);
 
   if(isenthalpy){
      femmodel->InputMakeDiscontinuous(BasalforcingsGroundediceMeltingRateEnum);
      enthalpy_analysis = new EnthalpyAnalysis();
      enthalpy_analysis->Core(femmodel);
      delete enthalpy_analysis;
   }
   else{
      if(VerboseSolution()) _printf0_("   computing temperatures\n");
      femmodel->SetCurrentConfiguration(ThermalAnalysisEnum);
      solutionsequence_thermal_nonlinear(femmodel);
 
      if(VerboseSolution()) _printf0_("   computing melting\n");
      femmodel->SetCurrentConfiguration(MeltingAnalysisEnum);
      solutionsequence_linear(femmodel);
   }
 
   if(save_results){
      if(VerboseSolution()) _printf0_("   saving thermal results\n");
      femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs);
   }
 
   /*Free ressources:*/
   if(numoutputs){for(int i=0;i<numoutputs;i++){xDelete<char>(requested_outputs[i]);} xDelete<char*>(requested_outputs);}
 
   /*End profiler*/
        femmodel->profiler->Stop(THERMALCORE);
}