[cig-commits] r7246 - mc/3D/CitcomS/branches/v2.2/lib

[tan2 at geodynamics.org]

Author: tan2
Date: 2007-06-14 14:39:17 -0700 (Thu, 14 Jun 2007)
New Revision: 7246

Modified:
   mc/3D/CitcomS/branches/v2.2/lib/Topo_gravity.c
Log:
Porting r7245 from trunk -- Zeroing geoid arrays before every geoid calculation. This fixes issue114

Modified: mc/3D/CitcomS/branches/v2.2/lib/Topo_gravity.c
===================================================================
--- mc/3D/CitcomS/branches/v2.2/lib/Topo_gravity.c	2007-06-14 21:38:19 UTC (rev 7245)
+++ mc/3D/CitcomS/branches/v2.2/lib/Topo_gravity.c	2007-06-14 21:39:17 UTC (rev 7246)
@@ -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++)  {