tallow propagation of top wall fluid pressures downwards - 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 aec1f5c5adb57d50565a0bf015fd7408566c4cbf
 (DIR) parent 9215e58a3f6d52a2db23665514e82283ed256265
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu, 22 Jan 2015 15:12:05 +0100
       
       allow propagation of top wall fluid pressures downwards
       
       Diffstat:
         M src/darcy.cuh                       |       4 ++--
       
       1 file changed, 2 insertions(+), 2 deletions(-)
       ---
 (DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -976,7 +976,7 @@ __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 && bc_top == 0))
                    dp_expl = 0.0;
        
        #ifdef REPORT_FORCING_TERMS
       t@@ -1128,7 +1128,7 @@ __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 && bc_top == 0))
                    dp_impl = 0.0;
                    //p_new = p;