textend dirichlet BC one cell further down from top wall - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
 (HTM) git clone git://src.adamsgaard.dk/sphere
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) LICENSE
       ---
 (DIR) commit 70bb2b42c839195212ba778c5f498a4ea82b0e8c
 (DIR) parent fb19e2f8cd6c9b146d4a04d34d756b985028ae1a
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Wed,  4 Feb 2015 16:03:49 +0100
       
       extend dirichlet BC one cell further down from top wall
       
       Diffstat:
         M src/darcy.cuh                       |       6 ++++--
       
       1 file changed, 4 insertions(+), 2 deletions(-)
       ---
 (DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -976,7 +976,8 @@ __global__ void firstDarcySolution(
                // Dirichlet BC at dynamic top wall. wall0_iz will be larger than the
                // grid if the wall isn't dynamic
                if ((bc_bot == 0 && z == 0) || (bc_top == 0 && z == nz-1)
       -                || (z >= wall0_iz && bc_top == 0))
       +                || (z >= wall0_iz-1 && bc_top == 0))
       +                //|| (z >= wall0_iz && bc_top == 0))
                    dp_expl = 0.0;
        
        #ifdef REPORT_FORCING_TERMS
       t@@ -1128,7 +1129,8 @@ __global__ void updateDarcySolution(
                // Dirichlet BC at dynamic top wall. wall0_iz will be larger than the
                // grid if the wall isn't dynamic
                if ((bc_bot == 0 && z == 0) || (bc_top == 0 && z == nz-1)
       -                || (z >= wall0_iz && bc_top == 0))
       +                || (z >= wall0_iz-1 && bc_top == 0))
       +                //|| (z >= wall0_iz && bc_top == 0))
                    dp_impl = 0.0;
                    //p_new = p;