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

[tan2 at geodynamics.org]

Author: tan2
Date: 2011-11-13 21:06:30 -0800 (Sun, 13 Nov 2011)
New Revision: 19192

Modified:
   mc/3D/CitcomS/trunk/lib/Regional_boundary_conditions.c
Log:
Fixed a bug in reflective TBC where y-index should be used, not z-index.
The bug is harmless though.

I checked the code history. This bug exists since CitcomS 1.0, but not in
CitcomCU nor CitcomT.


Modified: mc/3D/CitcomS/trunk/lib/Regional_boundary_conditions.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Regional_boundary_conditions.c	2011-11-13 05:39:15 UTC (rev 19191)
+++ mc/3D/CitcomS/trunk/lib/Regional_boundary_conditions.c	2011-11-14 05:06:30 UTC (rev 19192)
@@ -383,7 +383,7 @@
           node2 = i +(j-1)*E->lmesh.noz + (E->lmesh.noy-1)*E->lmesh.noz*E->lmesh.nox;
           E->node[m][node2] = E->node[m][node2] & (~TBY);
           E->node[m][node2] = E->node[m][node2] | FBY;
-          E->sphere.cap[m].TB[3][node2] = 0.0;
+          E->sphere.cap[m].TB[2][node2] = 0.0;
           }    /* end loop for i and j */
 
   return;