tadd lines noting D0, D50 and D100 - 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 3b2f81307fd4d1eb1dda6c012c9795cc89e4486a
(DIR) parent fc0101e370e16bcd388383a32f38e23a6894f3f5
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Thu, 14 Aug 2014 13:54:50 +0200
add lines noting D0, D50 and D100
Diffstat:
M python/sphere.py | 8 ++++++++
1 file changed, 8 insertions(+), 0 deletions(-)
---
(DIR) diff --git a/python/sphere.py b/python/sphere.py
t@@ -4164,10 +4164,18 @@ class sim:
t[i-1] = sim.time_current[0]
dh[i-1] = h - sim.w_x[0]
+ # save consolidation variables
+ self.D0 = h
+ self.D100 = h - dh[-1]
+ self.D50 = (self.D0 + self.D100)/2.0
+
fig = plt.figure()
plt.xlabel('Time [s]')
plt.ylabel('Consolidation [m]')
plt.semilogx(t, dh)
+ plt.axhline(y = self.D0)
+ plt.axhline(y = self.D50)
+ plt.axhline(y = self.D100)
plt.grid()
plt.savefig(self.sid + '-loadcurve.' + graphics_format)
plt.clf()