tdisable normalization while simulations aren't done - 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 eeed41b75866f4132dc200378894ba4e40ad97e6
 (DIR) parent 1464c92fe0e7fffaf127a24b567a9fb224370ca4
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Wed, 27 Aug 2014 10:28:09 +0200
       
       disable normalization while simulations aren't done
       
       Diffstat:
         M python/consolidation-curve.py       |      19 ++++++++++---------
       
       1 file changed, 10 insertions(+), 9 deletions(-)
       ---
 (DIR) diff --git a/python/consolidation-curve.py b/python/consolidation-curve.py
       t@@ -24,8 +24,8 @@ for c_grad_p in c_grad_p_list:
                             str(c_phi) + '-c_grad_p=' + str(c_grad_p)
            if os.path.isfile('../output/' + sid + '.status.dat'):
                sim = sphere.sim(sid, fluid=True)
       -        t[c] = numpy.empty(sim.status())
       -        H[c] = numpy.empty(sim.status())
       +        t[c] = numpy.ones(sim.status())
       +        H[c] = numpy.ones(sim.status())
        
                #sim.visualize('walls')
                #sim.writeVTKall()
       t@@ -70,17 +70,18 @@ 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
       +#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('Normalized thickness change [-]')
       +#plt.ylabel('Normalized thickness change [-]')
       +plt.ylabel('Thickness change [m]')
        #plt.ticklabel_format(style='sci', axis='y', scilimits=(0,0))
        for c in range(len(c_grad_p_list)):
       -    H[c] /= -min_H_c
       +    #H[c] /= -min_H_c
            plt.semilogx(t[c], H[c], '.-', label='$c$ = %.2f' % (c_grad_p_list[c]))
        plt.grid()