granular-channel-hydro

commit 4516a5f45c72fc1a57aa0d5da689b2225948ebbf
parent 8983440e41f994aedaf6b34648e690dedb03818a
Author: Anders Damsgaard Christensen <adc@geo.au.dk>
Date:   Fri, 13 Jan 2017 22:37:49 -0800

keep grid points within model domain

Diffstat:
M
granular_channel_drainage/model.py
 | 
4
++--

1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/granular_channel_drainage/model.py b/granular_channel_drainage/model.py
@@ -50,8 +50,8 @@ class model:
                 #numpy.random.normal(0., noise_amplitude, N)
             #numpy.random.shuffle(x)
             #numpy.random.shuffle(y)
-            xPoints = numpy.linspace(0., Lx, Nx)
-            yPoints = numpy.linspace(0., Ly, Ny)
+            xPoints = numpy.linspace(Lx/Nx, Lx - Lx/Nx, Nx)
+            yPoints = numpy.linspace(Ly/Ny, Ly - Ly/Ny, Ny)
             gridPoints = numpy.array([[x,y] for y in yPoints for x in xPoints])
             gridPoints[::2, 1] = gridPoints[::2, 1] + Lx/Nx*0.5
             N = len(gridPoints[:, 0])