tfix dirichlet bc at dynamic 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 b0ec0f769e28d02c5dc5402ad463e99092e89e1c
 (DIR) parent 8ffec2ec616f479cbae25398373996c4b28f788a
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Mon, 19 Jan 2015 09:55:00 +0100
       
       fix dirichlet bc at dynamic top wall
       
       Diffstat:
         M python/halfshear-darcy-stress-star… |       1 -
         M src/device.cu                       |       3 ++-
       
       2 files changed, 2 insertions(+), 2 deletions(-)
       ---
 (DIR) diff --git a/python/halfshear-darcy-stress-starter.py b/python/halfshear-darcy-stress-starter.py
       t@@ -52,7 +52,6 @@ if fluid:
            sim.initFluid(mu = mu, p = 0.0, cfd_solver = 1)
            sim.setFluidBottomNoFlow()
            sim.setFluidTopFixedPressure()
       -    sim.wmode[0] = 3
            #sim.setDEMstepsPerCFDstep(10)
            sim.setMaxIterations(2e5)
            sim.setPermeabilityPrefactor(k_c)
 (DIR) diff --git a/src/device.cu b/src/device.cu
       t@@ -1259,7 +1259,8 @@ __host__ void DEM::startTime()
        #endif
        
                            // find cell containing top wall
       -                    if (walls.nw > 0 && walls.wmode[0] == 1) {
       +                    if (walls.nw > 0 &&
       +                            (walls.wmode[0] == 1 || walls.wmode[0] == 3)) {
                                wall0_iz = walls.nx->w/(grid.L[2]/grid.num[2]);
                                setNSepsilonAtTopWall<<<dimGridFluid, dimBlockFluid>>>(
                                        dev_ns_epsilon,