URI: 

       tporosity function seems to work - 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 6be8b36e5ce5d77511a239d32c52e486d037209d
   DIR parent 2e13e8b98b7b64cc9915b87bf4785d0a45b26295
  HTML 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;