[cig-commits] r16913 - in short/3D/PyLith/trunk/examples/3d/hex8: . spatialdb spatialdb/powerlaw
willic3 at geodynamics.org
willic3 at geodynamics.org
Sun Jun 6 20:36:13 PDT 2010
Author: willic3
Date: 2010-06-06 20:36:12 -0700 (Sun, 06 Jun 2010)
New Revision: 16913
Removed:
short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/fixeddisp_axial.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/fixeddisp_shear.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw_properties.spatialdb
Modified:
short/3D/PyLith/trunk/examples/3d/hex8/README
short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_genmaxwell.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_powerlaw.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw/powerlaw_gendb.cfg
short/3D/PyLith/trunk/examples/3d/hex8/step15.cfg
short/3D/PyLith/trunk/examples/3d/hex8/step16.cfg
short/3D/PyLith/trunk/examples/3d/hex8/step17.cfg
Log:
Removed unused databases, changed power-law to use the powerlaw_gendb.py
utility (rather than the simplified spatialdb), and changed the gravity
examples to use the generalized Maxwell model.
Modified: short/3D/PyLith/trunk/examples/3d/hex8/README
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/README 2010-06-07 00:22:20 UTC (rev 16912)
+++ short/3D/PyLith/trunk/examples/3d/hex8/README 2010-06-07 03:36:12 UTC (rev 16913)
@@ -1,7 +1,7 @@
-The examples in this directory include (1) simple shearing of a box
-and (2) slip on a through-going, vertical, strike-slip fault in a box
-that is 6 km x 6 km x 4 km with linear hexahedral cells that have
-edges 1.0 km long.
+The examples in this directory form a step-by-step sequence of 17 problems,
+each building on the one before (for the most part). All of the examples
+use the same mesh, which was created by Cubit. The mesh is 6 km x 6 km x 4
+km with linear hexahedral cells that have edges 1.0 km long.
The box spans the volume:
@@ -9,9 +9,9 @@
-3 km <= y <= +3 km
-6 km <= z <= 0 km.
-The mesh is generated using CUBIT. Journal files are included and are
-annotated to guide you through the GUI to replicate the commands in
-the journal files should you prefer to use the GUI.
+The mesh is generated using CUBIT. Journal files are included in the mesh
+directory and are annotated to guide you through the GUI to replicate the
+commands in the journal files should you prefer to use the GUI.
NOTE: Importing Exodus files into PyLith requires the netcdf
library. This is included in the PyLith binary distribution. If you
@@ -23,59 +23,99 @@
You can examine the Exodus file exported from CUBIT using the ncdump
command.
-See the configuration files shearxy.cfg and dislocation.cfg on
-instructions on how to run these simulations. The configuration files
-also contain information about the simulation parameters.
+The examples are all described in the PyLith manual. Results of each
+example are put in the stepxx directories after running the simulation.
+Each simulation may be run by typing:
+pylith stepxx.cfg
+where xx is the number of the desired example.
+
+
DESCRIPTION OF FILES
-box_hex8_1000m.exo - Exodus file containing mesh exported from CUBIT
+README - This file.
-dislocation.cfg - Simulation parameters for fault slip problem
+pylithapp.cfg - PyLith configuration file containing parameters common to
+all simulations.
-finalslip.spatialdb - Spatial database for spatial variation of final slip
+figures/hex8-mesh.jpg - JPEG image showing mesh (generated with ParaView).
-finalslip_rupture.spatialdb - Spatial database for spatial variation
-of final slip for earthquake rupture in Savage-Prescot type problem.
+mesh/box_hex8_1000m.exo - Exodus file containing mesh exported from Cubit.
-fixeddisp_shear.spatialdb - Spatial database for spatial variation of
-the displacement field in the Dirichlet (prescribed displacement)
-boundary conditions for the shear problem
+mesh/geometry.jou - Cubit journal file (script) to generate solid model
+geometry
-geometry.jou - CUBIT journal file (script) to generate solid model geometry
+mesh/mesh_hex8_1000m.jou - Cubit journal file (script) to run geometry
+script, generate the mesh, create blocks and nodesets, and export the
+mesh to an Exodus file
-gravity.cfg - Simulation parameters for gravitational body forces problem
+spatialdb/finalslip.spatialdb - Spatial database defining fault slip for
+example step03.
-mat_elastic.spatialdb - Spatial database specifying the parameters for
-the physical properties of a 3-D elastic, isotropic material.
+spatialdb/finalslip_rupture.spatialdb - Spatial database defining fault
+rupture for examples step06, step07, step08, and step09.
-mat_maxwell.spatialdb - Spatial database specifying the parameters for
-the physical properties of a 3-D Maxwell viscoelastic, isotropic
-material.
+spatialdb/initial_stress.spatialdb - Spatial database defining initial
+stresses for example step16.
-mat_genmaxwell.spatialdb - Spatial database specifying the parameters for
-the physical properties of a 3-D generalized Maxwell viscoelastic, isotropic
-material.
+spatialdb/mat_druckerprager.spatialdb - Spatial database defining
+Drucker-Prager elastoplastic material properties for example step09.
-mat_powerlaw.spatialdb - Spatial database specifying the parameters for
-the physical properties of a 3-D power-law viscoelastic, isotropic
-material.
+spatialdb/mat_elastic.spatialdb - Spatial database defining elastic
+material properties for all examples.
-mesh_hex8_1000m.cub - Native CUBIT file containing finite-element mesh
+spatialdb/mat_genmaxwell.spatialdb - Spatial database defining
+generalized Maxwell viscoelastic material properties. This spatial database
+is not used by any examples, but could easily be exchanged for the material
+models being used in one of the other examples (an 'exercise for the
+reader').
-mesh_hex8_1000m.jou - CUBIT journal file (script) to run geometry
-script, generate the mesh, create blocks and nodesets, and export the
-mesh to an Exodus file
+spatialdb/mat_maxwell.spatialdb - Spatial database defining Maxwell
+viscoelastic material properties for examples step04, step05, step06,
+step07, step15, step16, and step17.
-pylithapp.cfg - General simulation parameters for the mesh
+spatialdb/mat_powerlaw.spatialdb - Spatial database defining power-law
+viscoelastic material properties for example step08.
-savageprescott.cfg - Simulation parameters for Savage-Prescott type problem
+spatialdb/sliprate_creep.spatialdb - Spatial database defining slip on the
+creeping portion of the fault for examples step06, step07, step08, and
+step09.
-shearxy.cfg - Simulation parameters for the shear problem
+spatialdb/sliptime.spatialdb - Spatial database defining the distribution
+of slip initial time over the fault for examples step03, step06, step07,
+step08, and step09.
-sliprate_creep.spatialdb - Spatial database for spatial variation
-of slip rate for creepin portion of fault in Savage-Prescot type problem.
+spatialdb/tractions_axial_shear.spatialdb - Spatial database defining
+axial and shear tractions for example step02.
-sliptime.spatialdb - Spatial database for spatial variation of slip
-initiation time
+spatialdb/powerlaw - This directory contains files needed to generate
+power-law properties for PyLith, given the power-law parameters, a
+temperature distribution, and a set of points for which properties are
+desired. The powerlaw_gendb.py utility code is used from this directory
+to create spatialdb/mat_powerlaw.spatialdb.
+
+stepxx.cfg - PyLith configuration file containing parameters for a
+particular example. The example problems are briefly described below:
+
+step01: Dirichlet BC (static)
+step02: Neumann BC (static)
+step03: Earthquake rupture (static)
+step04: Dirichlet BC (quasi-static)
+step05: Dirichlet + Neumann (quasi-static)
+step06: Multiple earthquake rupture + creep (quasi-static)
+step07: Earthquake rupture + creep + Dirichlet BC (quasi-static)
+step08: Same as step07 with power-law rheology (quasi-static)
+step09: Same as step08 with Drucker-Prager elastoplastic rheology
+ (quasi-static)
+step10: Static friction (stick) + Dirichlet BC (static)
+step11: Static friction (slip) + Dirichlet BC (static)
+step12: Static friction + Dirichlet BC (quasi-static)
+step13: Slip-weakening friction + Dirichlet BC (quasi-static)
+step14: Rate- and state-friction + Dirichlet BC (quasi-static)
+step15: Gravitational body forces with generalized Maxwell rheology
+ (quasi-static)
+step16: Gravitational body forces + initial stresses with generalized
+ Maxwell rheology (quasi-static)
+step17: Gravitational body forces + finite strain with generalized Maxwell
+ rheology (quasi-static)
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/fixeddisp_axial.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/fixeddisp_axial.spatialdb 2010-06-07 00:22:20 UTC (rev 16912)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/fixeddisp_axial.spatialdb 2010-06-07 03:36:12 UTC (rev 16913)
@@ -1,29 +0,0 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the distribution of the
-// displacement field for Dirichlet boundary conditions associated
-// with axial displacements on -x and +x
-//
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 3
- value-names = dof-0 dof-1 dof-2
- value-units = m m m
- num-locs = 2
- data-dim = 1 // locations form a line
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0e+3 // specify coordinates in km
- space-dim = 3
- }
-}
-// Columns are
-// (1) x coordinate (km)
-// (2) y coordinate (km)
-// (3) z coordinate (km)
-// (4) Ux (m)
-// (5) Uy (m)
-// (6) Uz (m)
--3.0 0.0 0.0 +1.0 0.0 0.0
-+3.0 0.0 0.0 -1.0 0.0 0.0
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/fixeddisp_shear.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/fixeddisp_shear.spatialdb 2010-06-07 00:22:20 UTC (rev 16912)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/fixeddisp_shear.spatialdb 2010-06-07 03:36:12 UTC (rev 16913)
@@ -1,28 +0,0 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the distribution of the
-// displacement field for Dirichlet boundary conditions associated
-// with shear in the xy plane.
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 3
- value-names = displacement-x displacement-y displacement-z
- value-units = m m m
- num-locs = 2
- data-dim = 1 // locations form a line
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0e+3 // specify coordinates in km
- space-dim = 3
- }
-}
-// Columns are
-// (1) x coordinate (km)
-// (2) y coordinate (km)
-// (3) z coordinate (km)
-// (4) Ux (m)
-// (5) Uy (m)
-// (6) Uz (m)
--3.0 0.0 0.0 0.0 -1.0 0.0
-+3.0 0.0 0.0 0.0 +1.0 0.0
Modified: short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_genmaxwell.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_genmaxwell.spatialdb 2010-06-07 00:22:20 UTC (rev 16912)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_genmaxwell.spatialdb 2010-06-07 03:36:12 UTC (rev 16913)
@@ -4,6 +4,15 @@
// properties for a Generalized Maxwell viscoelastic material. In
// this case, the material properties are uniform.
//
+// The first viscosity corresponds to a Maxwell time of 10 years, the
+// second corresponds to a Maxwell time of 100 years, and the third
+// corresponds to a Maxwell time of 1000 years. The shear ratio is set so
+// that the first two models are weighted by 0.4 and the third is weighted
+// by 0.2. Since these add up to 1.0, the material is still a fluid. If
+// the sum of the shear ratios was less than one, the remainder goes into a
+// spring in parallel with the Maxwell models (and the material is thus a
+// solid).
+//
#SPATIAL.ascii 1
SimpleDB {
num-values = 9 // number of material property values
@@ -30,4 +39,4 @@
// (10) viscosity 1 (Pa-s)
// (11) viscosity 2 (Pa-S)
// (12) viscosity 3 (Pa-s)
-0.0 0.0 0.0 2700.0 1054.09255338946 1825.741858350554 0.5 0.5 0.0 1.0e19 1.0e18 1.0e20
+0.0 0.0 0.0 2500.0 3000.0 5291.502622129181 0.4 0.4 0.2 7.10046e18 7.10046e19 7.10046e20
Modified: short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_powerlaw.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_powerlaw.spatialdb 2010-06-07 00:22:20 UTC (rev 16912)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_powerlaw.spatialdb 2010-06-07 03:36:12 UTC (rev 16913)
@@ -1,30 +1,23 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the distribution of material
-// properties for a power-law viscoelastic material. In this case,
-// the material properties are uniform.
-//
#SPATIAL.ascii 1
SimpleDB {
- num-values = 6 // number of material property values
- value-names = density vs vp reference-strain-rate reference-stress power-law-exponent // names of the material property values
- value-units = kg/m**3 m/s m/s 1/s Pa None // units
- num-locs = 1 // number of locations
- data-dim = 0
- space-dim = 3
+ num-values = 3
+ value-names = reference-stress reference-strain-rate power-law-exponent
+ value-units = Pa 1/s none
+ num-locs = 10
+ data-dim = 1
+ space-dim = 3
cs-data = cartesian {
- to-meters = 1.0
- space-dim = 3
- }
+ to-meters = 1
+ space-dim = 3
}
-// Columns are
-// (1) x coordinate (m)
-// (2) y coordinate (m)
-// (3) z coordinate (m)
-// (4) density (kg/m^3)
-// (5) vs (m/s)
-// (6) vp (m/s)
-// (7) reference-strain-rate (1/s)
-// (8) reference-stress (Pa)
-// (8) power-law-exponent (no units)
-0.0 0.0 0.0 2500.0 3000.0 5291.502622129181 1.0e-6 1.25992105e+10 3.0
+}
+ 0.000000e+00 0.000000e+00 0.000000e+00 1.818610e+16 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -5.000000e+02 1.977176e+13 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -1.000000e+03 2.801905e+11 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -1.500000e+03 1.528518e+10 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -1.900000e+03 2.687566e+09 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -2.100000e+03 1.201397e+12 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -2.500000e+03 1.483187e+11 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -3.000000e+03 1.798919e+10 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -3.500000e+03 3.284422e+09 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -4.000000e+03 8.099233e+08 1.000000e-06 3.500000e+00
Modified: short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw/powerlaw_gendb.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw/powerlaw_gendb.cfg 2010-06-07 00:22:20 UTC (rev 16912)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw/powerlaw_gendb.cfg 2010-06-07 03:36:12 UTC (rev 16913)
@@ -26,4 +26,4 @@
reader.filename = powerlaw_points.txt
[powerlaw_gendb.iohandler]
-filename = powerlaw_properties.spatialdb
+filename = ../mat_powerlaw.spatialdb
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw_properties.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw_properties.spatialdb 2010-06-07 00:22:20 UTC (rev 16912)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw_properties.spatialdb 2010-06-07 03:36:12 UTC (rev 16913)
@@ -1,23 +0,0 @@
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 3
- value-names = reference-stress reference-strain-rate power-law-exponent
- value-units = Pa 1/s none
- num-locs = 10
- data-dim = 1
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1
- space-dim = 3
-}
-}
- 0.000000e+00 0.000000e+00 0.000000e+00 1.818610e+16 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -5.000000e+02 1.977176e+13 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -1.000000e+03 2.801905e+11 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -1.500000e+03 1.528518e+10 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -1.900000e+03 2.687566e+09 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -2.100000e+03 1.201397e+12 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -2.500000e+03 1.483187e+11 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -3.000000e+03 1.798919e+10 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -3.500000e+03 3.284422e+09 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -4.000000e+03 8.099233e+08 1.000000e-06 3.500000e+00
Modified: short/3D/PyLith/trunk/examples/3d/hex8/step15.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step15.cfg 2010-06-07 00:22:20 UTC (rev 16912)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step15.cfg 2010-06-07 03:36:12 UTC (rev 16913)
@@ -17,9 +17,9 @@
# The resulting stress field should be rho*g*h for Stress_zz, and should
# be rho*g*h*pr/(1-pr), for Stress_xx and Stress_yy, where pr is Poisson's
# ratio.
-# We make the lower_crust material Maxwell viscoelastic, and allow the
-# stresses to relax for 200 years, so they begin approaching and isotropic
-# stress state in the lower crust.
+# We make the lower_crust material generalized Maxwell viscoelastic, and
+# allow the stresses to relax for 200 years, so they begin approaching and
+# isotropic stress state in the lower crust.
# ----------------------------------------------------------------------
# RUNNING THE SIMULATON
@@ -68,18 +68,19 @@
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
-# Change material type of lower crust to Maxwell viscoelastic.
+# Change material type of lower crust to generalized Maxwell viscoelastic.
[pylithapp.timedependent]
-materials.lower_crust = pylith.materials.MaxwellIsotropic3D
+materials.lower_crust = pylith.materials.GenMaxwellIsotropic3D
# Provide a spatial database from which to obtain property values.
-# Since there are additional properties and state variables for the Maxwell
-# model, we explicitly request that they be output. Properties are named in
-# cell_info_fields and state variables are named in cell_data_fields.
+# Since there are additional properties and state variables for the
+# generalized Maxwell model, we explicitly request that they be output.
+# Properties are named in cell_info_fields and state variables are named in
+# cell_data_fields.
[pylithapp.timedependent.materials.lower_crust]
-db_properties.iohandler.filename = spatialdb/mat_maxwell.spatialdb
-output.cell_info_fields = [density,mu,lambda,maxwell_time]
-output.cell_data_fields = [total_strain,stress,viscous_strain]
+db_properties.iohandler.filename = spatialdb/mat_genmaxwell.spatialdb
+output.cell_info_fields = [density,mu,lambda,shear_ratio,maxwell_time]
+output.cell_data_fields = [total_strain,stress,viscous_strain_1,viscous_strain_2,viscous_strain_3]
# ----------------------------------------------------------------------
# boundary conditions
Modified: short/3D/PyLith/trunk/examples/3d/hex8/step16.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step16.cfg 2010-06-07 00:22:20 UTC (rev 16912)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step16.cfg 2010-06-07 03:36:12 UTC (rev 16913)
@@ -68,18 +68,19 @@
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
-# Change material type of lower crust to Maxwell viscoelastic.
+# Change material type of lower crust to generalized Maxwell viscoelastic.
[pylithapp.timedependent]
-materials.lower_crust = pylith.materials.MaxwellIsotropic3D
+materials.lower_crust = pylith.materials.GenMaxwellIsotropic3D
# Provide a spatial database from which to obtain property values.
-# Since there are additional properties and state variables for the Maxwell
-# model, we explicitly request that they be output. Properties are named in
-# cell_info_fields and state variables are named in cell_data_fields.
+# Since there are additional properties and state variables for the
+# generalized Maxwell model, we explicitly request that they be output.
+# Properties are named in cell_info_fields and state variables are named in
+# cell_data_fields.
[pylithapp.timedependent.materials.lower_crust]
-db_properties.iohandler.filename = spatialdb/mat_maxwell.spatialdb
-output.cell_info_fields = [density,mu,lambda,maxwell_time]
-output.cell_data_fields = [total_strain,stress,viscous_strain]
+db_properties.iohandler.filename = spatialdb/mat_genmaxwell.spatialdb
+output.cell_info_fields = [density,mu,lambda,shear_ratio,maxwell_time]
+output.cell_data_fields = [total_strain,stress,viscous_strain_1,viscous_strain_2,viscous_strain_3]
# ----------------------------------------------------------------------
# initial stresses
Modified: short/3D/PyLith/trunk/examples/3d/hex8/step17.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step17.cfg 2010-06-07 00:22:20 UTC (rev 16912)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step17.cfg 2010-06-07 03:36:12 UTC (rev 16913)
@@ -65,18 +65,19 @@
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
-# Change material type of lower crust to Maxwell viscoelastic.
+# Change material type of lower crust to generalized Maxwell viscoelastic.
[pylithapp.timedependent]
-materials.lower_crust = pylith.materials.MaxwellIsotropic3D
+materials.lower_crust = pylith.materials.GenMaxwellIsotropic3D
# Provide a spatial database from which to obtain property values.
-# Since there are additional properties and state variables for the Maxwell
-# model, we explicitly request that they be output. Properties are named in
-# cell_info_fields and state variables are named in cell_data_fields.
+# Since there are additional properties and state variables for the
+# generalized Maxwell model, we explicitly request that they be output.
+# Properties are named in cell_info_fields and state variables are named in
+# cell_data_fields.
[pylithapp.timedependent.materials.lower_crust]
-db_properties.iohandler.filename = spatialdb/mat_maxwell.spatialdb
-output.cell_info_fields = [density,mu,lambda,maxwell_time]
-output.cell_data_fields = [total_strain,stress,viscous_strain]
+db_properties.iohandler.filename = spatialdb/mat_genmaxwell.spatialdb
+output.cell_info_fields = [density,mu,lambda,shear_ratio,maxwell_time]
+output.cell_data_fields = [total_strain,stress,viscous_strain_1,viscous_strain_2,viscous_strain_3]
# ----------------------------------------------------------------------
# boundary conditions
More information about the CIG-COMMITS
mailing list