[cig-commits] r6633 - mc/3D/ellipsis3d/trunk

[walter at geodynamics.org]

Author: walter
Date: 2007-04-23 09:18:40 -0700 (Mon, 23 Apr 2007)
New Revision: 6633

Modified:
   mc/3D/ellipsis3d/trunk/Tracer_read.c
Log:
handle the case where a particle is just on the edge of a box

Modified: mc/3D/ellipsis3d/trunk/Tracer_read.c
===================================================================
--- mc/3D/ellipsis3d/trunk/Tracer_read.c	2007-04-23 03:13:24 UTC (rev 6632)
+++ mc/3D/ellipsis3d/trunk/Tracer_read.c	2007-04-23 16:18:40 UTC (rev 6633)
@@ -707,15 +707,15 @@
     
       for(m=1;m<=E->tracer.NUM_TRACERS;m++) { /* RAA --> added this distinction  */
 	if(2==dims) {  
-	  if((E->tracer.tx[m] > rect_tx1[i]) && (E->tracer.tx[m] <= rect_tx2[i]) &&
-	     (E->tracer.tz[m] > rect_tz1[i]) && (E->tracer.tz[m] <= rect_tz2[i]) ) {
+	  if((E->tracer.tx[m] >= rect_tx1[i]) && (E->tracer.tx[m] <= rect_tx2[i]) &&
+	     (E->tracer.tz[m] >= rect_tz1[i]) && (E->tracer.tz[m] <= rect_tz2[i]) ) {
 	      E->tracer.property_group[m] = rect_tracer_colour[i];
 	  }
 	}
 	else if(3==dims) {	/* RAA --> added this bit  */
-	  if((E->tracer.tx[m] > rect_tx1[i]) && (E->tracer.tx[m] <= rect_tx2[i]) &&
-	     (E->tracer.tz[m] > rect_tz1[i]) && (E->tracer.tz[m] <= rect_tz2[i]) &&
-	     (E->tracer.ty[m] > rect_ty1[i]) && (E->tracer.ty[m] <= rect_ty2[i]) ) {
+	  if((E->tracer.tx[m] >= rect_tx1[i]) && (E->tracer.tx[m] <= rect_tx2[i]) &&
+	     (E->tracer.tz[m] >= rect_tz1[i]) && (E->tracer.tz[m] <= rect_tz2[i]) &&
+	     (E->tracer.ty[m] >= rect_ty1[i]) && (E->tracer.ty[m] <= rect_ty2[i]) ) {
 	      E->tracer.property_group[m] = rect_tracer_colour[i];
 	  }
         }