tUse custom ranges for simulations, change to linux path - 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 2c37b16a9453da73ba601a75c2385c80f4e5bab1
(DIR) parent bfeffc1ed018cf598acd283f96978f4e73ae8b7e
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Wed, 7 Jan 2015 11:39:00 +0100
Use custom ranges for simulations, change to linux path
Diffstat:
M python/halfshear-darcy-internals.py | 16 +++++++++++-----
1 file changed, 11 insertions(+), 5 deletions(-)
---
(DIR) diff --git a/python/halfshear-darcy-internals.py b/python/halfshear-darcy-internals.py
t@@ -13,15 +13,21 @@ from permeabilitycalculator import *
import matplotlib.pyplot as plt
from matplotlib.ticker import MaxNLocator
-steps = [10, 50, 100, 1000, 1999]
-nsteps_avg = 1 # no. of steps to average over
-#nsteps_avg = 100 # no. of steps to average over
-
#sigma0 = float(sys.argv[1])
sigma0 = 20000.0
#k_c = 3.5e-13
k_c = float(sys.argv[1])
+if k_c == 3.5e-15:
+ steps = [1232, 1332, 1433, 1534, 1635]
+elif k_c == 3.5e-13:
+ steps = [100, 200, 300, 410, 515]
+else:
+ steps = [10, 50, 100, 1000, 1999]
+nsteps_avg = 1 # no. of steps to average over
+#nsteps_avg = 100 # no. of steps to average over
+
+
sid = 'halfshear-darcy-sigma0=' + str(sigma0) + '-k_c=' + str(k_c) + \
'-mu=1.797e-06-velfac=1.0-shear'
sim = sphere.sim(sid, fluid=True)
t@@ -213,5 +219,5 @@ ax[3].locator_params(nbins=3)
filename = 'halfshear-darcy-internals-k_c=%.0e.pdf' % (k_c)
plt.savefig(filename)
-shutil.copyfile(filename, '/Users/adc/articles/own/2/graphics/' + filename)
+shutil.copyfile(filename, '/home/adc/articles/own/2/graphics/' + filename)
print(filename)