tdarcy values should only be written if they are computed - 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 35b801d036542c0258011cbfa832189afe081ba3
 (DIR) parent 92486a3dafd6d17fa87078ae20a80b59d1a9a75a
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Mon, 21 Oct 2013 12:54:40 +0200
       
       darcy values should only be written if they are computed
       
       Diffstat:
         M src/device.cu                       |      10 ++++++----
       
       1 file changed, 6 insertions(+), 4 deletions(-)
       ---
 (DIR) diff --git a/src/device.cu b/src/device.cu
       t@@ -1119,10 +1119,12 @@ __host__ void DEM::startTime()
                    writebin(file);
        
                    // Write Darcy arrays
       -            sprintf(file,"output/%s.d_phi.output%05d.bin", sid.c_str(), time.step_count);
       -            writeDarcyArray(d_phi, file);
       -            sprintf(file,"output/%s.d_K.output%05d.bin", sid.c_str(), time.step_count);
       -            writeDarcyArray(d_K, file);
       +            if (params.nu > 0.0 && darcy == 1) {
       +                sprintf(file,"output/%s.d_phi.output%05d.bin", sid.c_str(), time.step_count);
       +                writeDarcyArray(d_phi, file);
       +                sprintf(file,"output/%s.d_K.output%05d.bin", sid.c_str(), time.step_count);
       +                writeDarcyArray(d_K, file);
       +            }
        
                    if (CONTACTINFO == 1) {
                        // Write contact information to stdout