Changeset 17144

Timestamp:

01/17/14 16:31:39 (11 years ago)

Author:

jbondzio

Message:

ADD: setting spcs on ice domain for extrapolation

Location:

issm/trunk-jpl/src/c/analyses

Files:

• ExtrapolationAnalysis.cpp (modified) (2 diffs)
• ExtrapolationAnalysis.h (modified) (1 diff)

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

@@ -56 +56 @@
 
         if(VerboseSolution()) _printf0_("extrapolation: call computational core:\n");
+        UpdateConstraints(femmodel->elements);
         solutionsequence_linear(femmodel);
 
@@ -230 +231 @@
 
 }/*}}}*/
-
+void ExtrapolationAnalysis::UpdateConstraints(Elements* elements){/*{{{*/
+
+        for(int i=0;i<elements->Size();i++){
+                Element* element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
+                this->SetConstraintsOnIce(element);
+        }
+
+}/*}}}*/
+void ExtrapolationAnalysis::SetConstraintsOnIce(Element* element){/*{{{*/
+
+        int numnodes=element->GetNumberOfNodes();       
+
+        /* Intermediaries */
+        int extvar_enum;
+        IssmDouble phi,value;
+        Node* node = NULL;
+
+        /* Get further parameters */
+        element->FindParam(&extvar_enum, ExtrapolationVariableEnum);
+        
+        Input* levelset_input=element->GetInput(MaskIceLevelsetEnum); _assert_(levelset_input);
+        Input* extvar_input=element->GetInput(extvar_enum); _assert_(extvar_input);
+        
+        GaussPenta* gauss=new GaussPenta();
+        for(int in=0;in<numnodes;in++){
+                gauss->GaussNode(element->GetElementType(),in);
+                node=element->GetNode(in);
+                levelset_input->GetInputValue(&phi,gauss);
+                extvar_input->GetInputValue(&value,gauss);
+                if(phi<=0.){
+                        /* if ice, set dirichlet BC */
+                        node->ApplyConstraint(1,value);
+                }
+                else {
+                        /* no ice, set no spc */
+                        node->DofInFSet(0); 
+                }
+        }
+        delete gauss;
+}/*}}}*/

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

@@ -30 +30 @@
         void GetB(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss);
         void GetBprime(IssmDouble* Bprime,Element* element,IssmDouble* xyz_list,Gauss* gauss);
-
+        void UpdateConstraints(Elements* elements);
+        void SetConstraintsOnIce(Element* element);
 };
 #endif