Changeset 15967 for issm/trunk-jpl/src/c/modules/GroundinglineMigrationx/GroundinglineMigrationx.cpp

Timestamp:

08/27/13 08:35:37 (12 years ago)

Author:

seroussi

Message:

CHG: trying to finish Grounding line changes

File:

• issm/trunk-jpl/src/c/modules/GroundinglineMigrationx/GroundinglineMigrationx.cpp (modified) (6 diffs)

issm/trunk-jpl/src/c/modules/GroundinglineMigrationx/GroundinglineMigrationx.cpp

@@ -51 +51 @@
 }
 
-/*FUNCTION CreateNodesOnFloatingIce {{{*/
-Vector<IssmDouble>* CreateNodesOnFloatingIce(Nodes* nodes,int configuration_type){ 
-
-        int     i,numnods;
-        Vector<IssmDouble>*   vec_nodes_on_floatingice = NULL;
-        Node *node                     = NULL;
-
-        /*First, initialize nodes_on_floatingice, which will track which nodes have changed status: */
-        numnods=nodes->NumberOfNodes(configuration_type);
-        vec_nodes_on_floatingice=new Vector<IssmDouble>(numnods);
-
-        /*Loop through nodes, and fill vec_nodes_on_floatingice: */
-        for(i=0;i<nodes->Size();i++){
-                node=(Node*)nodes->GetObjectByOffset(i);
-                if(node->InAnalysis(configuration_type)){
-                        if(node->IsFloating()){
-                                vec_nodes_on_floatingice->SetValue(node->Sid(),1.0,INS_VAL);
-                        }
-                }
-        }
-
-        /*Assemble vector: */
-        vec_nodes_on_floatingice->Assemble();
-
-        return vec_nodes_on_floatingice;
-}
-/*%}}}*/
 /*FUNCTION PotentialUngrounding {{{*/
 IssmDouble*    PotentialUngrounding(Elements* elements,Vertices* vertices,Parameters* parameters){ 
@@ -111 +84 @@
         int      numberofvertices;
         int      nflipped,local_nflipped;
-        IssmDouble*  nodes_on_floatingice                  = NULL;
+        IssmDouble*  phi                  = NULL;
         IssmDouble*  elements_neighboring_floatingce      = NULL;
         Vector<IssmDouble>*      vec_elements_neighboring_floatingice = NULL;
-        Vector<IssmDouble>*      vec_nodes_on_floatingice              = NULL;
+        Vector<IssmDouble>*      vec_phi              = NULL;
         Node*    node                                  = NULL;
         Element* element                               = NULL;
@@ -122 +95 @@
         numberofvertices=vertices->NumberOfVertices();
 
-        /*recover vec_nodes_on_floatingice*/
-        vec_nodes_on_floatingice=CreateNodesOnFloatingIce(nodes,analysis_type);
-        nodes_on_floatingice=vec_nodes_on_floatingice->ToMPISerial();
+        /*recover vec_phi*/
+        vec_phi=new Vector<IssmDouble>(vertices->NumberOfVertices());
+        for(i=0;i<elements->Size();i++){
+                Element* element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
+                element->GetVectorFromInputs(vec_phi,MaskGroundediceLevelsetEnum);
+        }
+        phi=vec_phi->ToMPISerial();
 
         nflipped=1; //bootstrap
@@ -135 +112 @@
                 for(i=0;i<elements->Size();i++){
                         element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
-                        vec_elements_neighboring_floatingice->SetValue(element->Sid(),element->IsNodeOnShelfFromFlags(nodes_on_floatingice)?1.0:0.0,INS_VAL);
+                        vec_elements_neighboring_floatingice->SetValue(element->Sid(),element->IsNodeOnShelfFromFlags(phi)?1.0:0.0,INS_VAL);
                 }
 
@@ -147 +124 @@
                         element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
                         if(reCast<int,IssmDouble>(elements_neighboring_floatingce[element->Sid()])){
-                                local_nflipped+=element->UpdatePotentialUngrounding(vertices_potentially_ungrounding,vec_nodes_on_floatingice,nodes_on_floatingice);
+                                local_nflipped+=element->UpdatePotentialUngrounding(vertices_potentially_ungrounding,vec_phi,phi);
                         }
                 }
-                vec_nodes_on_floatingice->Assemble();
+                vec_phi->Assemble();
 
                 ISSM_MPI_Allreduce(&local_nflipped,&nflipped,1,ISSM_MPI_INT,ISSM_MPI_SUM,IssmComm::GetComm());
@@ -157 +134 @@
                 /*Avoid leaks: */
                 xDelete<IssmDouble>(elements_neighboring_floatingce);
-                xDelete<IssmDouble>(nodes_on_floatingice);
+                xDelete<IssmDouble>(phi);
 
                 /*Assemble and serialize:*/
                 delete vec_elements_neighboring_floatingice;
-                nodes_on_floatingice=vec_nodes_on_floatingice->ToMPISerial();
+                phi=vec_phi->ToMPISerial();
         }
 
         /*Free ressources:*/
-        delete vec_nodes_on_floatingice;
+        delete vec_phi;
         xDelete<IssmDouble>(elements_neighboring_floatingce);
 
-        return nodes_on_floatingice;
+        return phi;
 }
 /*}}}*/