[cig-commits] r19698 - short/3D/PyLith/trunk/tests/2d/plasticity/cyclic

[brad at geodynamics.org]

Author: brad
Date: 2012-02-28 15:37:54 -0800 (Tue, 28 Feb 2012)
New Revision: 19698

Modified:
   short/3D/PyLith/trunk/tests/2d/plasticity/cyclic/plot_data.py
   short/3D/PyLith/trunk/tests/2d/plasticity/cyclic/pylithapp.cfg
Log:
Correct plotting of 3-D results (adjust indices to components matching 2-D plane strain).

Modified: short/3D/PyLith/trunk/tests/2d/plasticity/cyclic/plot_data.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/plasticity/cyclic/plot_data.py	2012-02-28 23:24:07 UTC (rev 19697)
+++ short/3D/PyLith/trunk/tests/2d/plasticity/cyclic/plot_data.py	2012-02-28 23:37:54 UTC (rev 19698)
@@ -64,8 +64,8 @@
 ax = pylab.subplot(3, 3, 7)
 pylab.plot(quasi2D['time'], quasi2D['total_strain'][:,0,2], 'b-',
            dyn2D['time'], dyn2D['total_strain'][:,0,2], 'r--',
-           quasi3D['time'], quasi3D['total_strain'][:,0,2], 'g-',
-           dyn3D['time'], dyn3D['total_strain'][:,0,2], 'm--',
+           quasi3D['time'], quasi3D['total_strain'][:,0,3], 'g-',
+           dyn3D['time'], dyn3D['total_strain'][:,0,3], 'm--',
     )
 ax.set_xlabel("Time (yr)")
 
@@ -113,8 +113,8 @@
 ax = pylab.subplot(3, 3, 9)
 pylab.plot(quasi2D['time'], quasi2D['stress'][:,0,2], 'b-',
            dyn2D['time'], dyn2D['stress'][:,0,2], 'r--',
-           quasi3D['time'], quasi3D['stress'][:,0,2], 'b-',
-           dyn3D['time'], dyn3D['stress'][:,0,2], 'r--',
+           quasi3D['time'], quasi3D['stress'][:,0,3], 'b-',
+           dyn3D['time'], dyn3D['stress'][:,0,3], 'r--',
     )
 ax.set_xlabel("Time (yr)")
 

Modified: short/3D/PyLith/trunk/tests/2d/plasticity/cyclic/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/plasticity/cyclic/pylithapp.cfg	2012-02-28 23:24:07 UTC (rev 19697)
+++ short/3D/PyLith/trunk/tests/2d/plasticity/cyclic/pylithapp.cfg	2012-02-28 23:37:54 UTC (rev 19698)
@@ -33,7 +33,7 @@
 db_properties = spatialdata.spatialdb.UniformDB
 db_properties.label = DP elastoplastic db
 db_properties.values = [vp,vs,density,friction-angle,cohesion,dilatation-angle]
-db_properties.data = [5291.502622129181*m/s,3000.0*m/s,2700.0*kg/m**3,20.0*degree,1.0*MPa,0.0*degree]
+db_properties.data = [5291.502622129181*m/s,3000.0*m/s,2700.0*kg/m**3,20.0*degree,1.0*MPa,1.0*degree]
 fit_mohr_coulomb = middle
 
 quadrature.cell = pylith.feassemble.FIATLagrange