source: issm/trunk/src/c/modules/PostprocessingEnthalpyx/PostprocessingEnthalpyx.cpp@ 16137

/*!\file PostprocessingEnthalpyx.cpp
 * \brief: execute processing of variables after computation: water drainage, subglacial melting rates, update of thermal BCs
 */

#include "./PostprocessingEnthalpyx.h"
#include "../../shared/shared.h"
#include "../../toolkits/toolkits.h"

void PostprocessingEnthalpyx(FemModel* femmodel){
    /*Intermediary*/
        int i;
        Element* element=NULL;

        /*drain excess water fraction: */
        for (i=0;i<femmodel->elements->Size();i++){
                element=dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));
                element->DrainWaterfraction();
        }

        /*Compute basal melting rates: */
        for (i=0;i<femmodel->elements->Size();i++){
                element=dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));
                element->ComputeBasalMeltingrate();
        }

        /*Update basal dirichlet BCs for enthalpy: */
        for (i=0;i<femmodel->elements->Size();i++){
                element=dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));
                element->UpdateThermalBasalConstraints();
        }
}