tcreate vtk files in separate loop - 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 e72f1e870c239ec2e23ecfbc0d8e1f46def5d8f0
 (DIR) parent 69ca3d1d16dafe98bc5c455ce2f8a72840b7c0fa
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Wed, 13 Aug 2014 13:32:43 +0200
       
       create vtk files in separate loop
       
       Diffstat:
         M python/permeability-results.py      |       9 ++++++++-
         M python/permeabilitycalculator.py    |       2 +-
       
       2 files changed, 9 insertions(+), 2 deletions(-)
       ---
 (DIR) diff --git a/python/permeability-results.py b/python/permeability-results.py
       t@@ -1,7 +1,8 @@
        #!/usr/bin/env python
        import numpy
       -import matplotlib.pyplot as plt
       +import sphere
        from permeabilitycalculator import *
       +import matplotlib.pyplot as plt
        
        sids = [
            'permeability-dp=1000.0-c_phi=1.0-c_grad_p=0.01',
       t@@ -25,6 +26,12 @@ for sid in sids:
            c_grad_p[i] = pc.c_grad_p()
            i += 1
                
       +
       +# produce VTK files
       +for sid in sids:
       +    sim = sphere.sim(sid, fluid=True)
       +    sim.writeVTKall()
       +
        fig = plt.figure()
        plt.xlabel('Pressure gradient coefficient $c$ [-]')
        plt.ylabel('Hydraulic conductivity $K$ [m/s]')
 (DIR) diff --git a/python/permeabilitycalculator.py b/python/permeabilitycalculator.py
       t@@ -108,7 +108,7 @@ class PermeabilityCalc:
        
            def printResults(self):
                print('\n### Permeability resuts for "' + self.sid + '" ###')
       -        print('Pressure gradient: dP = ' + str(self.dPdL) + ' Pa/m')
       +        print('Pressure gradient: dPdL = ' + str(self.dPdL) + ' Pa/m')
                print('Flux: Q = ' + str(self.Q) + ' m^3/s')
                print('Intrinsic permeability: k = ' + str(self.k) + ' m^2')
                print('Saturated hydraulic conductivity: K = ' + str(self.K) + ' m/s')