[cig-commits] r7245 - mc/3D/CitcomS/trunk/lib

[tan2 at geodynamics.org]

Author: tan2
Date: 2007-06-14 14:38:19 -0700 (Thu, 14 Jun 2007)
New Revision: 7245

Modified:
   mc/3D/CitcomS/trunk/lib/Topo_gravity.c
Log:
Zeroing geoid arrays before every geoid calculation. This fixes issue114

Modified: mc/3D/CitcomS/trunk/lib/Topo_gravity.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Topo_gravity.c	2007-06-14 21:11:24 UTC (rev 7244)
+++ mc/3D/CitcomS/trunk/lib/Topo_gravity.c	2007-06-14 21:38:19 UTC (rev 7245)
@@ -452,6 +452,12 @@
     /* cos coeff */
     geoid[1] = (float*)malloc((E->sphere.hindice+1)*sizeof(float));
 
+    /* reset arrays */
+    for (p = 0; p <= E->sphere.hindice; p++) {
+        harm_geoid[0][p] = 0;
+        harm_geoid[1][p] = 0;
+    }
+
     /* loop over each layer */
     for(k=1;k<E->lmesh.noz;k++)  {