tporosity function seems to work - 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 6be8b36e5ce5d77511a239d32c52e486d037209d
 (DIR) parent 2e13e8b98b7b64cc9915b87bf4785d0a45b26295
 (HTM) Author: Anders Damsgaard <adc@geo.au.dk>
       Date:   Mon, 17 Dec 2012 15:47:02 +0100
       
       porosity function seems to work
       
       Diffstat:
         M src/sphere.cpp                      |       3 +++
       
       1 file changed, 3 insertions(+), 0 deletions(-)
       ---
 (DIR) diff --git a/src/sphere.cpp b/src/sphere.cpp
       t@@ -309,6 +309,9 @@ Float sphericalCap(const Float h, const Float r)
        // Calculate the porosity with depth, and write to file in output directory
        void DEM::porosity(const int z_slices)
        {
       +    // The porosity value is higher at the boundaries due 
       +    // to the no-flux BCs.
       +    
            Float top;
            if (walls.nw > 0)
                top = walls.nx->w;