tdo not zero fluid-pressure gradient unless top is Dirichlet BC - 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 de8f4eb1433514d7729abe95e9e1367161bb1ef6
 (DIR) parent 1a3f65f59a59b92c76e2babaf7aa5c01df5d93f8
 (HTM) Author: Anders Damsgaard Christensen <adc@geo.au.dk>
       Date:   Fri, 26 Aug 2016 14:53:21 -0700
       
       do not zero fluid-pressure gradient unless top is Dirichlet BC
       
       Diffstat:
         M src/darcy.cuh                       |       2 +-
       
       1 file changed, 1 insertion(+), 1 deletion(-)
       ---
 (DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -1232,7 +1232,7 @@ __global__ void findDarcyPressureForceLinear(
                                grad_p_iter = MAKE_FLOAT3(0., 0., 0.);
        
                            // Add Neumann BC at top wall
       -                    if (i_z + d_iz >= wall0_iz - 1)
       +                    if (bc_top == 0 && i_z + d_iz >= wall0_iz - 1)
                                grad_p_iter.z = 0.0;
        
                            n = MAKE_FLOAT3(dx*d_ix, dy*d_iy, dz*d_iz);