[cig-commits] r5737 - in long/3D/Gale/trunk: . src/Gale src/PICellerator src/PICellerator/Utils/src src/StGermain src/StgFEM src/Underworld

[walter at geodynamics.org]

Author: walter
Date: 2007-01-09 10:37:16 -0800 (Tue, 09 Jan 2007)
New Revision: 5737

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/src/Gale/
   long/3D/Gale/trunk/src/PICellerator/
   long/3D/Gale/trunk/src/PICellerator/Utils/src/StressBC.c
   long/3D/Gale/trunk/src/StGermain/
   long/3D/Gale/trunk/src/StgFEM/
   long/3D/Gale/trunk/src/Underworld/
Log:
 r1346 at earth:  boo | 2007-01-05 13:49:39 -0800
 Remove all externals



Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
   - 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1345
   + 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1346


Property changes on: long/3D/Gale/trunk/src/Gale
___________________________________________________________________
Name: svn:externals
   - VMake			https://csd.vpac.org/svn/VMake/branches/decomp3d



Property changes on: long/3D/Gale/trunk/src/PICellerator
___________________________________________________________________
Name: svn:externals
   - 

Modified: long/3D/Gale/trunk/src/PICellerator/Utils/src/StressBC.c
===================================================================
--- long/3D/Gale/trunk/src/PICellerator/Utils/src/StressBC.c	2007-01-09 18:37:14 UTC (rev 5736)
+++ long/3D/Gale/trunk/src/PICellerator/Utils/src/StressBC.c	2007-01-09 18:37:16 UTC (rev 5737)
@@ -378,21 +378,53 @@
             coord1=Mesh_GetVertex(mesh,elementNodes[lower]);
             coord2=Mesh_GetVertex(mesh,elementNodes[upper]);
             area=coord2[direction]-coord1[direction];
-
-/*             printf("area %d %d %d %d %lf %lf %lf %lf %lf %lf %lf %lf\n",lElement_I, lower, upper, direction, */
-/*                  Mesh_CoordAt(mesh,mesh->elementNodeTbl[lElement_I][0])[0], */
-/*                  Mesh_CoordAt(mesh,mesh->elementNodeTbl[lElement_I][0])[1], */
-/*                  Mesh_CoordAt(mesh,mesh->elementNodeTbl[lElement_I][1])[0], */
-/*                  Mesh_CoordAt(mesh,mesh->elementNodeTbl[lElement_I][1])[1], */
-/*                  Mesh_CoordAt(mesh,mesh->elementNodeTbl[lElement_I][2])[0], */
-/*                  Mesh_CoordAt(mesh,mesh->elementNodeTbl[lElement_I][2])[1], */
-/*                  Mesh_CoordAt(mesh,mesh->elementNodeTbl[lElement_I][3])[0], */
-/*                  Mesh_CoordAt(mesh,mesh->elementNodeTbl[lElement_I][3])[1], */
-/*                  area); */
-
           }
         else
           {
+            double *coord1, *coord2, *coord3, *coord4;
+            int lower,upper,direction;
+            switch(self->_wall)
+              {
+              case StressBC_Wall_Left:
+                lower=0;
+                upper=3;
+                direction=1;
+                break;
+              case StressBC_Wall_Right:
+                lower=1;
+                upper=2;
+                direction=1;
+                break;
+              case StressBC_Wall_Bottom:
+                lower=0;
+                upper=1;
+                direction=0;
+                break;
+              case StressBC_Wall_Top:
+                lower=3;
+                upper=2;
+                direction=0;
+                break;
+              case StressBC_Wall_Front:
+                lower=0;
+                upper=3;
+                direction=1;
+                break;
+              case StressBC_Wall_Back:
+                lower=0;
+                upper=3;
+                direction=1;
+                break;
+              }
+
+            Mesh_GetIncidence(mesh, Mesh_GetDimSize(mesh), lElement_I,
+                              MT_VERTEX,&elementNodeCount, &elementNodes);
+                              
+            coord1=Mesh_GetVertex(mesh,elementNodes[lower]);
+            coord2=Mesh_GetVertex(mesh,elementNodes[upper]);
+            area=coord2[direction]-coord1[direction];
+
+
             printf("Need to implement area in StressBC.c for dim==3\n");
             assert(0);
           }


Property changes on: long/3D/Gale/trunk/src/StGermain
___________________________________________________________________
Name: svn:externals
   - 


Property changes on: long/3D/Gale/trunk/src/StgFEM
___________________________________________________________________
Name: svn:externals
   - 


Property changes on: long/3D/Gale/trunk/src/Underworld
___________________________________________________________________
Name: svn:externals
   -