Changeset 17414

Timestamp:

03/12/14 08:39:59 (11 years ago)

Author:

seroussi

Message:

CHG: local referential coordinates in tria as well

File:

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

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

@@ -1952 +1952 @@
                         IssmDouble theta = atan(slope);
 
+                        /*New X axis                  New Z axis*/
+                        xz_plane[0]=cos(theta);       xz_plane[3]=0.;  
+                        xz_plane[1]=sin(theta);       xz_plane[4]=0.;  
+                        xz_plane[2]=0.;               xz_plane[5]=1.;          
                         if(groundedice>0){
-                                /*New X axis                  New Z axis*/
-                                xz_plane[0]=cos(theta);       xz_plane[3]=0.;  
-                                xz_plane[1]=sin(theta);       xz_plane[4]=0.;  
-                                xz_plane[2]=0.;               xz_plane[5]=1.;          
                                 this->nodes[i]->DofInSSet(1); //vy
                         }
                         else{
-                                /*New X axis                  New Z axis*/
-                                xz_plane[0]=1.;               xz_plane[3]=0;  
-                                xz_plane[1]=0.;               xz_plane[4]=0;  
-                                xz_plane[2]=0;                xz_plane[5]=1.;          
                                 this->nodes[i]->DofInFSet(1); //vy
                         }