tnormalize the thickness change - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
  HTML git clone git://src.adamsgaard.dk/sphere
   DIR Log
   DIR Files
   DIR Refs
   DIR LICENSE
       ---
   DIR commit 89fc26ac72d9da35930497454992f7a34a1c842a
   DIR parent 96eeffa6fcc7b60bcc548030d18412780338bfe5
  HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Wed, 27 Aug 2014 10:04:22 +0200
       
       normalize the thickness change
       
       Diffstat:
         M python/consolidation-curve.py       |      11 +++++++++--
       
       1 file changed, 9 insertions(+), 2 deletions(-)
       ---
   DIR diff --git a/python/consolidation-curve.py b/python/consolidation-curve.py
       t@@ -67,12 +67,19 @@ for c_grad_p in c_grad_p_list:
        
            c += 1
        
       +# Normalize the thickness change
       +min_H = 0.0
       +for c in range(len(c_grad_p_list)):
       +    min_H_c = numpy.min(H[c])
       +    if min_H_c < min_H:
       +        min_H = min_H_c
        
        plt.xlabel('Time [s]')
       -plt.ylabel('Thickness change [m]')
       +plt.ylabel('Normalized thickness change [-]')
        #plt.ticklabel_format(style='sci', axis='y', scilimits=(0,0))
        for c in range(len(c_grad_p_list)):
       -    plt.semilogx(t[c], H[c], 'o-', label='$c$ = %.2f' % (c_grad_p_list[c]))
       +    H[c] /= -min_H_c
       +    plt.semilogx(t[c], H[c], '.-', label='$c$ = %.2f' % (c_grad_p_list[c]))
        plt.grid()
        
        plt.legend(loc=0)