Changeset 17704

Timestamp:

04/09/14 19:18:18 (11 years ago)

Author:

hongjuy

Message:

CHG: Assign 1 or -1 for vertices not on bed

File:

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

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

@@ -64 +64 @@
         IssmDouble*  serial_vertexfloating = NULL;
         IssmDouble*  phi                   = NULL;
+        int index;
+        int* dof=NULL;
+        IssmDouble* xyz_list=NULL;
 
         /*Initialize vector with number of vertices*/
-        int numberofvertices=vertices->NumberOfVertices();
-        vertexgrounded=new Vector<IssmDouble>(numberofvertices);
-        vertexfloating=new Vector<IssmDouble>(numberofvertices);
-        phi           =xNew<IssmDouble>(numberofvertices);
+        int numberofvertices            = vertices->NumberOfVertices();
+        vertexgrounded                  = new Vector<IssmDouble>(numberofvertices);
+        vertexfloating                  = new Vector<IssmDouble>(numberofvertices);
+        phi                             = xNew<IssmDouble>(numberofvertices);
+        IssmDouble x[numberofvertices];
+        
+        /*Get x coordinate for each vertex*/
+        for(int i = 0;i<elements->Size();i++){
+                Element* element = dynamic_cast<Element*>(elements->GetObjectByOffset(i));
+                element->GetVerticesCoordinates(&xyz_list);
+                element->GetDofList(&dof,FSvelocityEnum,GsetEnum);      
+                for (int j=0;j<3;j++)   x[dof[j]]=xyz_list[3*j];
+        }
 
         /*Fill vector vertices_potentially_floating: */
-        for(int i=0;i<elements->Size();i++){
+        for(int i = 0;i<elements->Size();i++){
                 Element* element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
                 element->FSContactMigration(vertexgrounded,vertexfloating);
@@ -80 +92 @@
         vertexgrounded->Assemble();
         vertexfloating->Assemble();
-        serial_vertexgrounded=vertexgrounded->ToMPISerial();
-        serial_vertexfloating=vertexfloating->ToMPISerial();
-
-        for(int i=0;i<numberofvertices;i++){
+        serial_vertexgrounded = vertexgrounded->ToMPISerial();
+        serial_vertexfloating = vertexfloating->ToMPISerial();
+
+        /*find vertex on GL */
+        for(int i = 0;i<numberofvertices;i++){
                 if(serial_vertexgrounded[i]==1. && serial_vertexfloating[i]==1.){
-                        phi[i]=0.;
-                }
-                else if(serial_vertexgrounded[i]==1.){
-                        phi[i]=1.;
-                }
-                else if(serial_vertexfloating[i]==1.){
-                        phi[i]=-1.;
-                }
-                else{
-                        _error_("not supported (vertexfloating="<<serial_vertexfloating[i]<<" vertexgrounded="<<serial_vertexgrounded[i]<<")");
+                        phi[i] = 0.;
+                        index  = i;
+                }
+        }
+
+        /*assign value for phi, grounded if >0*/
+        for(int i = 0;i<numberofvertices;i++){
+                if (x[i] < x[index]) {
+                        phi[i] = 1.;
+                }
+                else if (x[i] > x[index]) {
+                        phi[i] = -1.;
+                }
+                else if (x[i]==x[index]) {
+                        phi[i] = 0;
                 }
         }