[cig-commits] r16338 - in short/3D/PyLith/trunk/examples/3d/hex8: . powerlaw
willic3 at geodynamics.org
willic3 at geodynamics.org
Wed Feb 24 21:02:28 PST 2010
Author: willic3
Date: 2010-02-24 21:02:27 -0800 (Wed, 24 Feb 2010)
New Revision: 16338
Modified:
short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg
short/3D/PyLith/trunk/examples/3d/hex8/gravity.cfg
short/3D/PyLith/trunk/examples/3d/hex8/gravity_istress.cfg
short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg
short/3D/PyLith/trunk/examples/3d/hex8/shearxy.cfg
Log:
Increased the temperature to get slightly more interesting power-law
behavior. Also added output.cell_data_fields for all examples.
Modified: short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg 2010-02-25 04:03:53 UTC (rev 16337)
+++ short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg 2010-02-25 05:02:27 UTC (rev 16338)
@@ -110,6 +110,7 @@
[pylithapp.timedependent.materials.viscoelastic.output]
cell_info_fields = [density,mu,lambda,maxwell_time]
+cell_data_fields = [total_strain,stress,viscous_strain]
cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = dislocation-viscoelastic.vtk
Modified: short/3D/PyLith/trunk/examples/3d/hex8/gravity.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/gravity.cfg 2010-02-25 04:03:53 UTC (rev 16337)
+++ short/3D/PyLith/trunk/examples/3d/hex8/gravity.cfg 2010-02-25 05:02:27 UTC (rev 16338)
@@ -98,5 +98,6 @@
[pylithapp.timedependent.materials.viscoelastic.output]
cell_info_fields = [density,mu,lambda,maxwell_time]
+cell_data_fields = [total_strain,stress,viscous_strain]
cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = gravity-viscoelastic.vtk
Modified: short/3D/PyLith/trunk/examples/3d/hex8/gravity_istress.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/gravity_istress.cfg 2010-02-25 04:03:53 UTC (rev 16337)
+++ short/3D/PyLith/trunk/examples/3d/hex8/gravity_istress.cfg 2010-02-25 05:02:27 UTC (rev 16338)
@@ -111,5 +111,6 @@
[pylithapp.timedependent.materials.viscoelastic.output]
cell_info_fields = [density,mu,lambda,maxwell_time]
+cell_data_fields = [total_strain,stress,viscous_strain]
cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = gravity_istress-viscoelastic.vtk
Modified: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb 2010-02-25 04:03:53 UTC (rev 16337)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw/temperature.spatialdb 2010-02-25 05:02:27 UTC (rev 16338)
@@ -1,8 +1,8 @@
// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
//
// This spatial database defines a simple linear 1D temperature profile,
-// assuming a temperature of 0 degrees C (273 degrees K) at the surface,
-// increasing to 400 degrees C (773 degrees K) at 4 km depth.
+// assuming a temperature of 100 degrees C (373 degrees K) at the surface,
+// increasing to 650 degrees C (923 degrees K) at 4 km depth.
//
#SPATIAL.ascii 1
SimpleDB {
@@ -22,5 +22,5 @@
// (2) y coordinate (m)
// (3) z coordinate (m)
// (4) temperature (K)
-0.0 0.0 0.0 273.0
-0.0 0.0 -4.0 673.0
+0.0 0.0 0.0 373.0
+0.0 0.0 -4.0 923.0
Modified: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb 2010-02-25 04:03:53 UTC (rev 16337)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb 2010-02-25 05:02:27 UTC (rev 16338)
@@ -11,12 +11,12 @@
space-dim = 3
}
}
- 0.000000e+00 0.000000e+00 0.000000e+00 8.804741e+20 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -5.000000e+02 8.804741e+20 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -1.000000e+03 8.804741e+20 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -1.500000e+03 8.804741e+20 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -2.000000e+03 8.804741e+20 1.000000e-06 1.500000e+00
- 0.000000e+00 0.000000e+00 -2.500000e+03 3.163259e+14 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -3.000000e+03 3.163259e+14 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -3.500000e+03 3.163259e+14 1.000000e-06 3.500000e+00
- 0.000000e+00 0.000000e+00 -4.000000e+03 3.163259e+14 1.000000e-06 3.500000e+00
+ 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.818610e+16 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -1.000000e+03 1.818610e+16 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -1.500000e+03 1.818610e+16 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -2.000000e+03 1.818610e+16 1.000000e-06 1.500000e+00
+ 0.000000e+00 0.000000e+00 -2.500000e+03 1.935443e+11 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -3.000000e+03 1.935443e+11 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -3.500000e+03 1.935443e+11 1.000000e-06 3.500000e+00
+ 0.000000e+00 0.000000e+00 -4.000000e+03 1.935443e+11 1.000000e-06 3.500000e+00
Modified: short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg 2010-02-25 04:03:53 UTC (rev 16337)
+++ short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg 2010-02-25 05:02:27 UTC (rev 16338)
@@ -55,6 +55,7 @@
# [pylithapp.timedependent.materials.viscoelastic]
# db_properties.iohandler.filename = mat_maxwell.spatialdb
# output.cell_info_fields = [density,mu,lambda,maxwell_time]
+# output.cell_data_fields = [total_strain,stress,viscous_strain]
# *** END MAXWELL SECTION ***
# ----------------------------------------------------------------------
@@ -72,6 +73,7 @@
db_properties.db_A = spatialdata.spatialdb.SimpleDB
db_properties.db_B = spatialdata.spatialdb.SimpleDB
output.cell_info_fields = [density,mu,lambda,reference_strain_rate,reference_stress,power_law_exponent]
+output.cell_data_fields = [total_strain,stress,viscous_strain]
[pylithapp.timedependent.materials.viscoelastic.db_properties]
values_A = [density,vs,vp]
@@ -166,7 +168,8 @@
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.domain]
-skip = 9
+output_freq = time_step
+time_step = 5.0*year
writer.filename = savageprescott.vtk
writer.time_format = %04.0f
writer.time_constant = 1.0*year
@@ -188,14 +191,16 @@
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
-skip = 9
+output_freq = time_step
+time_step = 5.0*year
cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = savageprescott-elastic.vtk
writer.time_format = %04.0f
writer.time_constant = 1.0*year
[pylithapp.timedependent.materials.viscoelastic.output]
-skip = 9
+output_freq = time_step
+time_step = 5.0*year
cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = savageprescott-viscoelastic.vtk
writer.time_format = %04.0f
Modified: short/3D/PyLith/trunk/examples/3d/hex8/shearxy.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/shearxy.cfg 2010-02-25 04:03:53 UTC (rev 16337)
+++ short/3D/PyLith/trunk/examples/3d/hex8/shearxy.cfg 2010-02-25 05:02:27 UTC (rev 16338)
@@ -90,5 +90,6 @@
[pylithapp.timedependent.materials.viscoelastic.output]
cell_info_fields = [density,mu,lambda,maxwell_time]
+cell_data_fields = [total_strain,stress,viscous_strain]
cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = shearxy-viscoelastic.vtk
More information about the CIG-COMMITS
mailing list