Fix node addressing in sc_capillary.

darthez · Dec 24, 2013 · 8d9b353 · 8d9b353
1 parent e53083d
commit 8d9b353
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/examples/binary_fluid/sc_capillary.py b/examples/binary_fluid/sc_capillary.py
@@ -21,8 +21,8 @@ def boundary_conditions(self, hx, hy):
         rem_y = (self.gy - chan_diam) / 2
 
         geometry = np.zeros(hx.shape, dtype=np.bool)
-        geometry[0,:] = True
-        geometry[self.gy-1,:] = True
+        geometry[hy == 0] = True
+        geometry[hy == self.gy - 1] = True
         geometry[np.logical_and(
                     hy < rem_y,
                     hy < rem_y - (np.abs((hx - self.gx/2)) - chan_len/2)