Changeset 17431

Timestamp:

03/14/14 11:22:55 (11 years ago)

Author:

jbondzio

Message:

ADD: coupling LSM - SIA flow model. Use this as a template for coupling LSM to an analysis.

File:

• issm/trunk-jpl/src/c/analyses/StressbalanceSIAAnalysis.cpp (modified) (9 diffs)

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

@@ -19 +19 @@
         /*Fetch data needed: */
         bool   isSIA;
+        bool   islevelset;
         iomodel->Constant(&isSIA,FlowequationIsSIAEnum);
+        iomodel->Constant(&islevelset,TransientIslevelsetEnum);
 
         /*Now, is the flag SIA on? otherwise, do nothing: */
@@ -38 +40 @@
         iomodel->FetchDataToInput(elements,ThicknessEnum);
         iomodel->FetchDataToInput(elements,FrictionCoefficientEnum);
+        if(islevelset){
+                iomodel->FetchDataToInput(elements,IceMaskNodeActivationEnum);
+                if(iomodel->meshtype!=Mesh2DhorizontalEnum)
+                        iomodel->FetchDataToInput(elements,MeshVertexonbedEnum); // required for updating active nodes
+        }
 
         /*Free data: */
@@ -150 +157 @@
 }/*}}}*/
 ElementMatrix* StressbalanceSIAAnalysis::CreateKMatrix(Element* element){/*{{{*/
+
+        /* Check if ice in element */
+        if(!element->IsIceInElement()) return NULL;
+
         int meshtype;
         element->FindParam(&meshtype,MeshTypeEnum);
@@ -162 +173 @@
 ElementMatrix* StressbalanceSIAAnalysis::CreateKMatrix2D(Element* element){/*{{{*/
 
+        /* Check if ice in element */
+        if(!element->IsIceInElement()) return NULL;
+
         /*Intermediaries */
         IssmDouble connectivity;
@@ -180 +194 @@
 }/*}}}*/
 ElementMatrix* StressbalanceSIAAnalysis::CreateKMatrix3D(Element* element){/*{{{*/
+
+        /* Check if ice in element */
+        if(!element->IsIceInElement()) return NULL;
 
         /*Intermediaries */
@@ -244 +261 @@
 ElementVector* StressbalanceSIAAnalysis::CreatePVector(Element* element){/*{{{*/
 
+        /* Check if ice in element */
+        if(!element->IsIceInElement()) return NULL;
+
         int meshtype;
         element->FindParam(&meshtype,MeshTypeEnum);
@@ -255 +275 @@
 }/*}}}*/
 ElementVector* StressbalanceSIAAnalysis::CreatePVector2D(Element* element){/*{{{*/
+
+        /* Check if ice in element */
+        if(!element->IsIceInElement()) return NULL;
 
         /*Intermediaries */
@@ -326 +349 @@
 }/*}}}*/
 ElementVector* StressbalanceSIAAnalysis::CreatePVector3D(Element* element){/*{{{*/
+
+        /* Check if ice in element */
+        if(!element->IsIceInElement()) return NULL;
 
         /*Intermediaries */
@@ -550 +576 @@
 }/*}}}*/
 void StressbalanceSIAAnalysis::UpdateConstraints(FemModel* femmodel){/*{{{*/
-        /*Default, do nothing*/
+
+        bool islevelset;
+        femmodel->parameters->FindParam(&islevelset,TransientIslevelsetEnum);
+        if(islevelset){
+                _printf0_("   Updating active and non-active nodes for StressbalanceSIAAnalysis \n");
+                SetActiveNodesLSMx(femmodel);
+        }
         return;
 }/*}}}*/