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

[tan2 at geodynamics.org]

Author: tan2
Date: 2008-11-04 16:19:05 -0800 (Tue, 04 Nov 2008)
New Revision: 13260

Modified:
   mc/3D/CitcomS/trunk/lib/Full_version_dependent.c
Log:
Minor fix on node counting


Modified: mc/3D/CitcomS/trunk/lib/Full_version_dependent.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Full_version_dependent.c	2008-11-04 22:50:02 UTC (rev 13259)
+++ mc/3D/CitcomS/trunk/lib/Full_version_dependent.c	2008-11-05 00:19:05 UTC (rev 13260)
@@ -77,6 +77,8 @@
   E->mesh.elz = E->mesh.noz-1;
 
   /* number of nodes, excluding overlaping nodes between processors */
+  /* each cap has one row of nox and one row of noy overlapped, exclude these nodes.
+   * nodes at north/south poles are exclued by all caps, include them by 2*noz*/
   E->mesh.nno = E->sphere.caps*(E->mesh.nox-1)*(E->mesh.noy-1)*E->mesh.noz + 2*E->mesh.noz;
 
   E->mesh.nel = E->sphere.caps*E->mesh.elx*E->mesh.elz*E->mesh.ely;
@@ -106,7 +108,7 @@
     E->mesh.ELX[i] = nox-1;
     E->mesh.ELY[i] = noy-1;
     E->mesh.ELZ[i] = noz-1;
-    E->mesh.NNO[i] = E->sphere.caps * (nox-1) * (noy-1) * noz;
+    E->mesh.NNO[i] = E->sphere.caps * (nox-1) * (noy-1) * noz + 2 * noz;
     E->mesh.NEL[i] = E->sphere.caps * (nox-1) * (noz-1) * (noy-1);
     E->mesh.NPNO[i] = E->mesh.NEL[i] ;
     E->mesh.NOX[i] = nox;