bmb_core.cpp File Reference

: core of the bmb (Basal mass balance) solution More...

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

Go to the source code of this file.

Functions
void	bmb_core (FemModel *femmodel)

Detailed Description

: core of the bmb (Basal mass balance) solution

Definition in file bmb_core.cpp.

Function Documentation

bmb_core()

void bmb_core

(

FemModel *

femmodel

)

Definition at line 12 of file bmb_core.cpp.

                                 {
 
   /*First, get BMB model from parameters*/
   int  basalforcing_model;
   bool isplume = false;
   femmodel->parameters->FindParam(&basalforcing_model,BasalforcingsEnum);
 
   if(VerboseSolution()) _printf0_("   computing basal mass balance\n");
 
   /*In some cases we need to run additional analyses to get the required input data*/
   if(basalforcing_model==BasalforcingsPicoEnum){
      femmodel->parameters->FindParam(&isplume,BasalforcingsPicoIsplumeEnum);
      if(isplume){
         femmodel->SetCurrentConfiguration(L2ProjectionBaseAnalysisEnum);
         femmodel->parameters->SetParam(BaseSlopeXEnum,InputToL2ProjectEnum);
         solutionsequence_linear(femmodel);
         femmodel->parameters->SetParam(BaseSlopeYEnum,InputToL2ProjectEnum);
         solutionsequence_linear(femmodel);
         femmodel->SetCurrentConfiguration(GLheightadvectionAnalysisEnum);
         solutionsequence_linear(femmodel);
      }
   }
 
   /*Call module now*/
   FloatingiceMeltingRatex(femmodel);
 
   /*Extrude basal melt*/
   femmodel->parameters->SetParam(BasalforcingsFloatingiceMeltingRateEnum,InputToExtrudeEnum);
   extrudefrombase_core(femmodel);
}