tx-axis in shear visualization now shear strain - 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 2a487de72d1a16ff8631293498605819772b8f15
 (DIR) parent 386dfba962f5e9bec5930c4b410c40e23ab5faf0
 (HTM) Author: Anders Damsgaard <adc@geo.au.dk>
       Date:   Wed, 19 Dec 2012 11:46:49 +0100
       
       x-axis in shear visualization now shear strain
       
       Diffstat:
         M python/sphere.py                    |      12 ++++++------
       
       1 file changed, 6 insertions(+), 6 deletions(-)
       ---
 (DIR) diff --git a/python/sphere.py b/python/sphere.py
       t@@ -1317,19 +1317,19 @@ def visualize(project, method = 'energy', savefig = True, outformat = 'png'):
                    # Plot stresses
                    if (outformat != 'txt'):
                        ax1 = plt.subplot2grid((2,1),(0,0))
       -                ax1.set_xlabel('Shear distance [m]')
       +                ax1.set_xlabel('Shear strain [-]')
                        ax1.set_ylabel('Stress [Pa]')
       -                ax1.plot(xdisp, sigma_eff, '+-g', label="$\sigma'$")
       -                ax1.plot(xdisp, sigma_def, '+-b', label="$\sigma_0$")
       -                ax1.plot(xdisp, tau, '+-r', label="$\\tau$")
       +                ax1.plot(xdisp/w_x0, sigma_eff, '+-g', label="$\sigma'$")
       +                ax1.plot(xdisp/w_x0, sigma_def, '+-b', label="$\sigma_0$")
       +                ax1.plot(xdisp/w_x0, tau, '+-r', label="$\\tau$")
                        ax1.legend(loc = 4)
                        ax1.grid()
        
                        # Plot dilation
                        ax2 = plt.subplot2grid((2,1),(1,0))
       -                ax2.set_xlabel('Shear distance [m]')
       +                ax2.set_xlabel('Shear strain [-]')
                        ax2.set_ylabel('Dilation [%]')
       -                ax2.plot(xdisp, dilation, '+-')
       +                ax2.plot(xdisp/w_x0, dilation, '+-')
                        ax2.grid()
                else :
                    print("Visualization type '" + method + "' not understood")