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

[tan2 at geodynamics.org]

Author: tan2
Date: 2011-09-17 20:54:13 -0700 (Sat, 17 Sep 2011)
New Revision: 18932

Modified:
   mc/3D/CitcomS/trunk/lib/Initial_temperature.c
Log:
fixed a bug that when the perturblayer is on the internal boundary, the temperature perturbation wasn't added to the lower processor.

Modified: mc/3D/CitcomS/trunk/lib/Initial_temperature.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Initial_temperature.c	2011-09-18 01:25:44 UTC (rev 18931)
+++ mc/3D/CitcomS/trunk/lib/Initial_temperature.c	2011-09-18 03:54:13 UTC (rev 18932)
@@ -442,10 +442,10 @@
         kk = E->convection.load_depth[p];
         con = E->convection.perturb_mag[p];
 
-        if ( (kk < 1) || (kk >= gnoz) ) continue; /* layer kk is outside domain */
+        if ( (kk < 1) || (kk > gnoz) ) continue; /* layer kk is outside domain */
 
         k = kk - E->lmesh.nzs + 1; /* convert global nz to local nz */
-        if ( (k < 1) || (k >= noz) ) continue; /* layer k is not inside this proc. */
+        if ( (k < 1) || (k > noz) ) continue; /* layer k is not inside this proc. */
         if (E->parallel.me_loc[1] == 0 && E->parallel.me_loc[2] == 0
             && E->sphere.capid[1] == 1 )
             fprintf(stderr,"Initial temperature perturbation:  layer=%d  mag=%g  l=%d  m=%d\n", kk, con, ll, mm);