Changeset 17002

Timestamp:

12/03/13 15:57:09 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: trying to add EnthalpyPostProcessing back to EnthalpyAnalysis.cpp

Location:

issm/trunk-jpl/src/c/analyses

Files:

• EnthalpyAnalysis.cpp (modified) (1 diff)
• EnthalpyAnalysis.h (modified) (1 diff)

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

@@ -785 +785 @@
 }/*}}}*/
 
+
+/*Modules*/
+void EnthalpyAnalysis::PostProcessing(FemModel* femmodel){/*{{{*/
+
+        /*Compute basal melting rates: */
+        for(int i=0;i<femmodel->elements->Size();i++){
+                Element* element=dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));
+                ComputeBasalMeltingrate(element);
+        }
+
+        /*drain excess water fraction: */
+        //for(int i=0;i<femmodel->elements->Size();i++){
+        //      element=dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));
+        //      element->DrainWaterfraction();
+        //}
+
+        /*Update basal dirichlet BCs for enthalpy: */
+        for(int i=0;i<femmodel->elements->Size();i++){
+                Element* element=dynamic_cast<Element*>(femmodel->elements->GetObjectByOffset(i));
+                UpdateBasalConstraints(element);
+        }
+}/*}}}*/
+void EnthalpyAnalysis::ComputeBasalMeltingrate(Element* element){/*{{{*/
+
+}/*}}}*/
+void EnthalpyAnalysis::DrainWaterfraction(Element* element){/*{{{*/
+
+}/*}}}*/
+void EnthalpyAnalysis::UpdateBasalConstraints(Element* element){/*{{{*/
+
+}/*}}}*/
+
 /*Intermediaries*/
 IssmDouble EnthalpyAnalysis::EnthalpyDiffusionParameter(Element* element,IssmDouble enthalpy,IssmDouble pressure){/*{{{*/

issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.h

@@ -36 +36 @@
                 void InputUpdateFromSolution(IssmDouble* solution,Element* element);
 
+                /*Modules*/
+                static void PostProcessing(FemModel* femmodel);
+                static void ComputeBasalMeltingrate(Element* element);
+                static void DrainWaterfraction(Element* element);
+                static void UpdateBasalConstraints(Element* element);
+
                 /*Intermediaries*/
                 IssmDouble EnthalpyDiffusionParameter(Element* element,IssmDouble enthalpy,IssmDouble pressure);