tmove dp=0 one iz upwards, may change previous results - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
  HTML git clone git://src.adamsgaard.dk/sphere
   DIR Log
   DIR Files
   DIR Refs
   DIR LICENSE
       ---
   DIR commit 9768e6e712769d16a8919dae571519040fcbd4fc
   DIR parent b8293f1e95dc5c366c98e0e282e0c7fa1f33d0e4
  HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu, 22 Jan 2015 10:00:58 +0100
       
       move dp=0 one iz upwards, may change previous results
       
       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 - 1 && 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 - 1 && bc_top == 0))
       +                //|| (z >= wall0_iz - 1 && bc_top == 0))
       +                || (z >= wall0_iz && bc_top == 0))
                    dp_impl = 0.0;
                    //p_new = p;