tseparate files for calculator and calculator usage - 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 bfba4f99f00ee5b5c3e717c9d2f8dd6a7868f537
(DIR) parent 92ea73312f1004063e84f628cce811977e1b228e
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Wed, 13 Aug 2014 11:20:50 +0200
separate files for calculator and calculator usage
Diffstat:
M python/permeability-calculator.py | 7 ++-----
A python/permeability-results.py | 15 +++++++++++++++
2 files changed, 17 insertions(+), 5 deletions(-)
---
(DIR) diff --git a/python/permeability-calculator.py b/python/permeability-calculator.py
t@@ -13,6 +13,8 @@ class PermeabilityCalc:
self.findPermeability()
self.findConductivity()
self.findMeanPorosity()
+ self.printResults()
+ self.plotEvolution()
def readfile(self):
self.sim = sphere.sim(self.sid, fluid=True)
t@@ -171,8 +173,3 @@ class PermeabilityCalc:
plt.savefig(filename)
print('Figure saved as "' + filename + '"')
plt.show()
-
-# Simulation ID
-pc = PermeabilityCalc('permeability-dp=1000.0')
-pc.printResults()
-pc.plotEvolution()
(DIR) diff --git a/python/permeability-results.py b/python/permeability-results.py
t@@ -0,0 +1,15 @@
+#!/usr/bin/env python
+from permeability-calculator import *
+
+pc = PermeabilityCalc('permeability-dp=1000.0')
+pc = PermeabilityCalc('permeability-dp=1000.0-c_phi=1.0-c_grad_p=0.01')
+pc = PermeabilityCalc('permeability-dp=1000.0-c_phi=1.0-c_grad_p=0.5')
+
+pc = PermeabilityCalc('permeability-dp=20000.0-c_phi=1.0-c_grad_p=0.01')
+pc = PermeabilityCalc('permeability-dp=20000.0-c_phi=1.0-c_grad_p=0.1')
+pc = PermeabilityCalc('permeability-dp=20000.0-c_phi=1.0-c_grad_p=0.5')
+
+pc = PermeabilityCalc('permeability-dp=4000.0')
+pc = PermeabilityCalc('permeability-dp=4000.0-c_phi=1.0-c_grad_p=0.01')
+pc = PermeabilityCalc('permeability-dp=4000.0-c_phi=1.0-c_grad_p=0.1')
+pc = PermeabilityCalc('permeability-dp=4000.0-c_phi=1.0-c_grad_p=0.5')