Changeset 26103

Timestamp:

03/16/21 08:46:08 (4 years ago)

Author:

Mathieu Morlighem

Message:

CHG: moveing Element Coordinate outside of sealevel stuff

File:

• issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -1284 +1284 @@
         _printf_("list of vertices on surface: "<<values[0]<<" "<<values[1]<<" "<<values[2]);
         _error_("Could not find 2 vertices on surface");
+}
+/*}}}*/
+void       Tria::ElementCoordinates(Vector<IssmDouble>* vxe,Vector<IssmDouble>* vye,Vector<IssmDouble>* vze, bool spherical){/*{{{*/
+
+        /*Look for x,y,z coordinates:*/
+        IssmDouble xyz_list[NUMVERTICES][3];
+        ::GetVerticesCoordinates(&xyz_list[0][0],this->vertices,NUMVERTICES);
+
+        /*Find centroid:*/
+        IssmDouble xe=(xyz_list[0][0]+xyz_list[1][0]+xyz_list[2][0])/3.0;
+        IssmDouble ye=(xyz_list[0][1]+xyz_list[1][1]+xyz_list[2][1])/3.0;
+        IssmDouble ze=(xyz_list[0][2]+xyz_list[1][2]+xyz_list[2][2])/3.0;
+
+        if(!spherical){
+                vxe->SetValue(this->sid,xe,INS_VAL);
+                vye->SetValue(this->sid,ye,INS_VAL);
+                vze->SetValue(this->sid,ze,INS_VAL);
+        }
+        else _error_("spherical coordinates not supported yet!");
+        return;
 }
 /*}}}*/
@@ -6219 +6239 @@
 }
 /*}}}*/
-void    Tria::ElementCoordinates(Vector<IssmDouble>* vxe,Vector<IssmDouble>* vye,Vector<IssmDouble>* vze, bool spherical){/*{{{*/
-
-        /*Look for x,y,z coordinates:*/
-        IssmDouble xyz_list[NUMVERTICES][3];
-        IssmDouble xe,ye,ze;
-
-        /*Element XYZ list*/
-        ::GetVerticesCoordinates(&xyz_list[0][0],this->vertices,NUMVERTICES);
-
-        /*Find centroid:*/
-        xe=(xyz_list[0][0]+xyz_list[1][0]+xyz_list[2][0])/3.0;
-        ye=(xyz_list[0][1]+xyz_list[1][1]+xyz_list[2][1])/3.0;
-        ze=(xyz_list[0][2]+xyz_list[1][2]+xyz_list[2][2])/3.0;
-
-        if(!spherical){
-                vxe->SetValue(this->sid,xe,INS_VAL);
-                vye->SetValue(this->sid,ye,INS_VAL);
-                vze->SetValue(this->sid,ze,INS_VAL);
-        }
-        else _error_("spherical coordinates not supported yet!");
-        return;
-}
-/*}}}*/
 
 //new logic