Changeset 10440 for issm/trunk/src/c/objects/Elements/Penta.cpp

Timestamp:

11/02/11 16:51:01 (13 years ago)

Author:

Mathieu Morlighem

Message:

Temporary commit for spc in Stokes under grounded ice, still not perfect

File:

• issm/trunk/src/c/objects/Elements/Penta.cpp (modified) (1 diff)

issm/trunk/src/c/objects/Elements/Penta.cpp ¶

@@ -2350 +2350 @@
         }
 
+}
+/*}}}*/
+/*FUNCTION Penta::ResetCoordinateSystem{{{1*/
+void  Penta::ResetCoordinateSystem(void){
+
+        int    approximation;
+        bool   stokes_basal_penalty;
+        double slopex[NUMVERTICES];
+        double slopey[NUMVERTICES];
+        double xz_plane[6];
+
+        /*Do not change CS if stokes_basal_penalty is activated*/
+        parameters->FindParam(&stokes_basal_penalty,DiagnosticStokesBasalPenaltyEnum);
+        if(stokes_basal_penalty) return;
+
+        /*For Stokes only: we want the CS to be tangential to the bedrock*/
+        inputs->GetInputValue(&approximation,ApproximationEnum);
+        if(IsFloating() || !IsOnBed() || (approximation!=StokesApproximationEnum && approximation!=MacAyealStokesApproximationEnum &&  approximation!=PattynStokesApproximationEnum)) return;
+
+        /*Get slope on each node*/
+        GetInputListOnVertices(&slopex[0],BedSlopeXEnum);
+        GetInputListOnVertices(&slopey[0],BedSlopeYEnum);
+
+        /*Loop over basal nodes (first 3) and update their CS*/
+        for(int i=0;i<NUMVERTICES2D;i++){
+
+                /*New X axis             New Z axis*/
+                xz_plane[0]=1.;          xz_plane[3]=-slopex[i];  
+                xz_plane[1]=0.;          xz_plane[4]=-slopey[i];  
+                xz_plane[2]=slopex[i];   xz_plane[5]=1.;          
+
+                XZvectorsToCoordinateSystem(&this->nodes[i]->coord_system[0][0],&xz_plane[0]);
+        }
 }
 /*}}}*/