[cig-commits] r6644 - mc/3D/ellipsis3d/trunk

[walter at geodynamics.org]

Author: walter
Date: 2007-04-23 23:03:52 -0700 (Mon, 23 Apr 2007)
New Revision: 6644

Modified:
   mc/3D/ellipsis3d/trunk/Allocation_of_memory.c
   mc/3D/ellipsis3d/trunk/Boundary_conditions_velocity.c
Log:
Initialize some uninitialized variables

Modified: mc/3D/ellipsis3d/trunk/Allocation_of_memory.c
===================================================================
--- mc/3D/ellipsis3d/trunk/Allocation_of_memory.c	2007-04-24 06:02:26 UTC (rev 6643)
+++ mc/3D/ellipsis3d/trunk/Allocation_of_memory.c	2007-04-24 06:03:52 UTC (rev 6644)
@@ -20,6 +20,7 @@
 */
 
 
+#include <string.h>
 #include <math.h>
 #include "element_definitions.h"
 #include "global_defs.h"
@@ -66,6 +67,15 @@
   E->depl = (standard_precision *) Malloc0((E->mesh.nno+1)*sizeof(standard_precision)); 
   /*RAA: 24/09/02, depl is melting stuff by C. O'Neill */
 
+  memset(E->T,0,(E->mesh.nno+1)*sizeof(standard_precision));
+  memset(E->nQ,0,(E->mesh.nno+1)*sizeof(standard_precision));
+  memset(E->edot,0,(E->mesh.nno+1)*sizeof(standard_precision));
+  memset(E->strs,0,(E->mesh.nno+1)*sizeof(standard_precision));
+  memset(E->strd,0,(E->mesh.nno+1)*sizeof(standard_precision));
+  memset(E->strd1,0,(E->mesh.nno+1)*sizeof(standard_precision));
+  memset(E->depl,0,(E->mesh.nno+1)*sizeof(standard_precision));
+
+
   if(2==E->mesh.nsd)
     E->Psi=(standard_precision *) Malloc0((E->mesh.nno+1)*sizeof(standard_precision));
   

Modified: mc/3D/ellipsis3d/trunk/Boundary_conditions_velocity.c
===================================================================
--- mc/3D/ellipsis3d/trunk/Boundary_conditions_velocity.c	2007-04-24 06:02:26 UTC (rev 6643)
+++ mc/3D/ellipsis3d/trunk/Boundary_conditions_velocity.c	2007-04-24 06:03:52 UTC (rev 6644)
@@ -352,7 +352,10 @@
         fprintf(E->fp1," nno Vbs1,2  %d   %g  %g\n",i,E->Vb[1][0][i],E->Vb[2][0][i]);
      } */
       
-    
+
+    /* Set the boundary to non-skew by default. */
+    E->control.HAVE_SKEWBCS = 0;
+
     /* Skew boundary conditions */
 
 #if 0