[cig-commits] r16653 - in short/3D/PyLith/trunk/examples/3d/hex8: . common
willic3 at geodynamics.org
willic3 at geodynamics.org
Wed May 5 19:31:47 PDT 2010
Author: willic3
Date: 2010-05-05 19:31:46 -0700 (Wed, 05 May 2010)
New Revision: 16653
Added:
short/3D/PyLith/trunk/examples/3d/hex8/common/box_hex8_1000m.exo
short/3D/PyLith/trunk/examples/3d/hex8/common/finalslip.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/common/finalslip_rupture.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/common/fixeddisp_axial.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/common/fixeddisp_shear.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/common/geometry.jou
short/3D/PyLith/trunk/examples/3d/hex8/common/initial_stress.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/common/mat_elastic.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/common/mat_genmaxwell.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/common/mat_maxwell.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/common/mat_powerlaw.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/common/mesh/
short/3D/PyLith/trunk/examples/3d/hex8/common/mesh_hex8_1000m.jou
short/3D/PyLith/trunk/examples/3d/hex8/common/powerlaw_properties.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/common/pylithapp.cfg
short/3D/PyLith/trunk/examples/3d/hex8/common/sliprate_creep.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/common/sliptime.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/common/spatialdb/
short/3D/PyLith/trunk/examples/3d/hex8/step01/
short/3D/PyLith/trunk/examples/3d/hex8/step02/
short/3D/PyLith/trunk/examples/3d/hex8/step03/
short/3D/PyLith/trunk/examples/3d/hex8/step04/
short/3D/PyLith/trunk/examples/3d/hex8/step05/
short/3D/PyLith/trunk/examples/3d/hex8/step06/
short/3D/PyLith/trunk/examples/3d/hex8/step07/
short/3D/PyLith/trunk/examples/3d/hex8/step08/
short/3D/PyLith/trunk/examples/3d/hex8/step09/
Removed:
short/3D/PyLith/trunk/examples/3d/hex8/box_hex8_1000m.exo
short/3D/PyLith/trunk/examples/3d/hex8/finalslip.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/finalslip_rupture.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/fixeddisp_axial.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/fixeddisp_shear.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/geometry.jou
short/3D/PyLith/trunk/examples/3d/hex8/initial_stress.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/mat_elastic.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/mat_genmaxwell.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/mat_maxwell.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/mat_powerlaw.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/mesh_hex8_1000m.jou
short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
short/3D/PyLith/trunk/examples/3d/hex8/sliprate_creep.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/sliptime.spatialdb
Log:
Made some more changes to directory structure.
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/box_hex8_1000m.exo
===================================================================
(Binary files differ)
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/box_hex8_1000m.exo (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/box_hex8_1000m.exo)
===================================================================
(Binary files differ)
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/finalslip.spatialdb (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/finalslip.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/finalslip.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/finalslip.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,31 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of slip on the
+// fault surface. In this case we prescribe a piecewise linear, depth
+// dependent distribution of slip. The slip is 2.0 m right-lateral
+// with 0.25 m of reverse slip at the surface with a linear taper from
+// 2.0 m to 0.0 m from -2 km to -4 km.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = left-lateral-slip reverse-slip fault-opening
+ value-units = m m m
+ num-locs = 3
+ data-dim = 1 // Locations of data points form a line.
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // Specify coordinates in km for convenience.
+ space-dim = 3
+ } // cs-data
+} // SimpleDB
+// Columns are
+// (1) x coordinate (km)
+// (2) y coordinate (km)
+// (3) z coordinate (km)
+// (4) left-lateral-slip (m) (right-lateral is negative)
+// (5) reverse-slip (m)
+// (6) fault-opening (m)
+0.0 0.0 0.0 -2.00 0.25 0.00
+0.0 0.0 -2.0 -2.00 0.00 0.00
+0.0 0.0 -4.0 0.00 0.00 0.00
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/finalslip_rupture.spatialdb (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/finalslip_rupture.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/finalslip_rupture.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/finalslip_rupture.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,31 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of slip on the
+// fault surface. In this case we impose 4.0 m of right-lateral slip
+// for z >= -2.0 km with zero slip for z < -2.0 km. There is no reverse
+// slip or fault opening.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = left-lateral-slip reverse-slip fault-opening
+ value-units = m m m
+ num-locs = 4
+ data-dim = 1 // Data is specified along a line.
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // Specify coordinates in km for convenience.
+ space-dim = 3
+ } // cs-data
+} // SimpleDB
+// Columns are
+// (1) x coordinate (km)
+// (2) y coordinate (km)
+// (3) z coordinate (km)
+// (4) left-lateral-slip (m) (right-lateral is negative)
+// (5) reverse-slip (m)
+// (6) fault-opening (m)
+0.0 0.0 0.00 -4.0 0.0 0.0
+0.0 0.0 -2.00 -4.0 0.0 0.0
+0.0 0.0 -2.01 0.0 0.0 0.0
+0.0 0.0 -4.00 0.0 0.0 0.0
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/fixeddisp_axial.spatialdb (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/fixeddisp_axial.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/fixeddisp_axial.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/fixeddisp_axial.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,32 @@
+// -*- 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
+//
+// dof-0: Ux(x) = 0
+// dof-1: Uy(x) = 1.0 m * x / 3 km
+// dof-2: Uz(x) = 0
+//
+#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
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/fixeddisp_shear.spatialdb (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/fixeddisp_shear.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/fixeddisp_shear.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/fixeddisp_shear.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,32 @@
+// -*- 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.
+//
+// dof-0: Ux(x) = 0
+// dof-1: Uy(x) = 1.0 m * x / 3 km
+// dof-2: Uz(x) = 0
+//
+#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
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/geometry.jou (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/geometry.jou)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/geometry.jou (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/geometry.jou 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,47 @@
+## /tools/common/cubit-10.2/bin/clarox
+## Cubit Version 10.2
+## Cubit Build 24
+## Revised 12/15/2006 16:09:40 MST
+## Running 06/18/2007 10:26:50 AM
+## Command Options:
+## -warning = On
+## -information = On
+
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+
+# Block is 6 km x 6 km x 4 km
+# -3 km <= x <= 3 km
+# -3 km <= y <= 3 km
+# -4 km <= z <= 0 km
+reset
+brick x 6000 y 6000 z 4000
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -2000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+surface 7 name "fault_surface"
+create planar surface with plane zplane offset -2000
+surface 8 name "material_interface"
+
+# ----------------------------------------------------------------------
+# Divide volumes using interface surfaces
+# ----------------------------------------------------------------------
+webcut volume 1 with plane surface fault_surface
+webcut volume 1 with plane surface material_interface
+webcut volume 4 with plane surface material_interface
+volume 1 name "elastic_xpos"
+volume 4 name "elastic_xneg"
+volume 5 name "visco_xpos"
+volume 6 name "visco_xneg"
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/initial_stress.spatialdb (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/initial_stress.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/initial_stress.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/initial_stress.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,36 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the initial stresses for the model
+// The stresses increase linearly with depth.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6 // number of stress components
+
+ // names of the stress components
+ value-names = stress-xx stress-yy stress-zz stress-xy stress-yz stress-xz
+
+ value-units = Pa Pa Pa Pa Pa Pa // units
+
+ num-locs = 2 // number of locations: 1 value each at the upper and lower
+ // surfaces. Linear interpolation should be used.
+
+ data-dim = 1
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) stress-xx (Pa)
+// (5) stress-yy (Pa)
+// (6) stress-zz (Pa)
+// (7) stress-xy (Pa)
+// (8) stress-yz (Pa)
+// (9) stress-xz (Pa)
+0.0 0.0 -4000.0 -98066500 -98066500 -98066500 0.0 0.0 0.0
+0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/mat_elastic.spatialdb (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/mat_elastic.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/mat_elastic.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/mat_elastic.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,26 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties. In this case, the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3 // number of material property values
+ value-names = density vs vp // names of the material property values
+ value-units = kg/m**3 m/s m/s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ 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)
+0.0 0.0 0.0 2500.0 3000.0 5291.502622129181
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/mat_genmaxwell.spatialdb (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/mat_genmaxwell.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/mat_genmaxwell.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/mat_genmaxwell.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,33 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a Generalized Maxwell viscoelastic material. In
+// this case, the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 9 // number of material property values
+ value-names = density vs vp shear-ratio-1 shear-ratio-2 shear-ratio-3 viscosity-1 viscosity-2 viscosity-3 // names of the material property values
+ value-units = kg/m**3 m/s m/s None None None Pa*s Pa*s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ 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) shear ratio 1 (dimensionless)
+// (8) shear ratio 2 (dimensionless)
+// (9) shear ratio 3 (dimensionless)
+// (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
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/mat_maxwell.spatialdb (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/mat_maxwell.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/mat_maxwell.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/mat_maxwell.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a Maxwell viscoelastic material. In this case,
+// the material properties are uniform.
+//
+// The viscosity corresponds to a Maxwell time of 100 years.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ 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) viscosity (Pa-s)
+0.0 0.0 0.0 2500.0 3000.0 5291.502622129181 7.10046e+19
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/mat_powerlaw.spatialdb (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/mat_powerlaw.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/mat_powerlaw.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/mat_powerlaw.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,30 @@
+// -*- 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
+ cs-data = cartesian {
+ to-meters = 1.0
+ 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
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/mesh_hex8_1000m.jou (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/mesh_hex8_1000m.jou)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/mesh_hex8_1000m.jou (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/mesh_hex8_1000m.jou 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,104 @@
+## /tools/common/cubit-10.2/bin/clarox
+## Cubit Version 10.2
+## Cubit Build 24
+## Revised 12/15/2006 16:09:40 MST
+## Running 06/18/2007 10:26:50 AM
+## Command Options:
+## -warning = On
+## -information = On
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size 1000
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 4
+block 1 name "elastic"
+block 2 volume 5 6
+block 2 name "viscoelastic"
+
+# ----------------------------------------------------------------------
+# Create nodeset for fault
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at A
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 20
+group "face_xpos" add node in surface 28
+nodeset 11 group face_xpos
+nodeset 11 name "face xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 30
+group "face_xneg" add node in surface 38
+nodeset 12 group face_xneg
+nodeset 12 name "face xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 21
+group "face_ypos" add node in surface 27
+group "face_ypos" add node in surface 33
+group "face_ypos" add node in surface 35
+nodeset 13 group face_ypos
+nodeset 13 name "face ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 23
+group "face_yneg" add node in surface 25
+group "face_yneg" add node in surface 31
+group "face_yneg" add node in surface 37
+nodeset 14 group face_yneg
+nodeset 14 name "face yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 12
+group "face_zneg" add node in surface 16
+nodeset 15 group face_zneg
+nodeset 15 name "face zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face w/o fault
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 16 group face_zneg_nofault
+nodeset 16 name "face zneg nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 10
+group "face_zpos" add node in surface 17
+nodeset 17 group face_zpos
+nodeset 17 name "face zpos"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "box_hex8_1000m.exo" dimension 3 overwrite
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/powerlaw_properties.spatialdb (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/powerlaw_properties.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/powerlaw_properties.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,23 @@
+#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
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/pylithapp.cfg (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/pylithapp.cfg 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,94 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+# Turn on some journals to show progress.
+[pylithapp.journal.info]
+timedependent = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshiocubit = 1
+implicitelasticity = 1
+faultcohesivekin = 1
+fiatlagrange = 1
+pylithapp = 1
+materials = 1
+powerlaw3d = 1
+##implicit = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1 ; uncomment to get very verbose mesh information
+
+# Change the default mesh reader to the CUBIT reader.
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+# Set filename of mesh to import.
+filename = box_hex8_1000m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+# Define the total time for the simulation and the default time step size.
+total_time = 0.0*s ; total time of simulation
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+
+# Set materials to an array of 2 materials:
+# 'elastic' and 'viscoelastic'
+materials = [elastic,viscoelastic]
+
+# Change from default material (elastic, isotropic, 3-D material) to
+# Maxwell viscoelastic, isotropic 3-D material.
+materials.viscoelastic = pylith.materials.MaxwellIsotropic3D
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+db_properties.iohandler.filename = mat_elastic.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+# We do not define the database info for the viscoelastic material, since
+# a different database type may be used for some examples.
+[pylithapp.timedependent.materials.viscoelastic]
+label = Viscoelastic material
+id = 2
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+# Set the solver options.
+
+[pylithapp.petsc]
+pc_type = asm
+
+# Change the preconditioner settings.
+sub_pc_factor_shift_type = none
+
+ksp_rtol = 1.0e-8
+ksp_max_it = 500
+ksp_gmres_restart = 50
+
+ksp_monitor = true
+ksp_view = true
+ksp_converged_reason = true
+
+snes_monitor = true
+snes_view = true
+snes_converged_reason = true
+
+log_summary = true
+# start_in_debugger = true
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/sliprate_creep.spatialdb (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/sliprate_creep.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/sliprate_creep.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/sliprate_creep.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,32 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of slip rate on
+// the fault surface for a constant slip rate time function. In this
+// case we specify a uniform slip rate of -2.0 cm/yr
+// (6.3376175628057904e-10 m/s) for z < -2.0 km and zero slip for z >=
+// -2.0 km. There is no reverse-slip or fault opening.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = left-lateral-slip reverse-slip fault-opening
+ value-units = cm/year cm/year cm/year
+ num-locs = 4
+ data-dim = 1 // Values are specified along a line.
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // Specify coordinates in km for convenience.
+ space-dim = 3
+ } // cs-data
+} // SimpleDB
+// Columns are
+// (1) x coordinate (km)
+// (2) y coordinate (km)
+// (3) z coordinate (km)
+// (4) left-lateral-slip rate (m/s) (right-lateral is negative)
+// (5) reverse-slip rate (m/s)
+// (6) fault-opening rate (m/s)
+0.0 0.0 0.00 0.00 0.0 0.0
+0.0 0.0 -2.00 0.00 0.0 0.0
+0.0 0.0 -2.01 -2.00 0.0 0.0
+0.0 0.0 -4.00 -2.00 0.0 0.0
Copied: short/3D/PyLith/trunk/examples/3d/hex8/common/sliptime.spatialdb (from rev 16651, short/3D/PyLith/trunk/examples/3d/hex8/sliptime.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/common/sliptime.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/common/sliptime.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -0,0 +1,26 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of the slip
+// initiation time over the fault surface for the kinematic fault
+// rupture. We specify a slip time of 0.0 s so that the slip occurs at
+// the origin time of the earthquake rupture.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 1
+ value-names = slip-time
+ value-units = s
+ num-locs = 1
+ data-dim = 0 // data is uniform (1 location)
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) Slip initiation time (s)
+0.0 0.0 0.0 0.0
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/finalslip.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/finalslip.spatialdb 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/finalslip.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,31 +0,0 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the distribution of slip on the
-// fault surface. In this case we prescribe a piecewise linear, depth
-// dependent distribution of slip. The slip is 2.0 m right-lateral
-// with 0.25 m of reverse slip at the surface with a linear taper from
-// 2.0 m to 0.0 m from -2 km to -4 km.
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 3
- value-names = left-lateral-slip reverse-slip fault-opening
- value-units = m m m
- num-locs = 3
- data-dim = 1 // Locations of data points form a line.
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0e+3 // Specify coordinates in km for convenience.
- space-dim = 3
- } // cs-data
-} // SimpleDB
-// Columns are
-// (1) x coordinate (km)
-// (2) y coordinate (km)
-// (3) z coordinate (km)
-// (4) left-lateral-slip (m) (right-lateral is negative)
-// (5) reverse-slip (m)
-// (6) fault-opening (m)
-0.0 0.0 0.0 -2.00 0.25 0.00
-0.0 0.0 -2.0 -2.00 0.00 0.00
-0.0 0.0 -4.0 0.00 0.00 0.00
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/finalslip_rupture.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/finalslip_rupture.spatialdb 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/finalslip_rupture.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,31 +0,0 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the distribution of slip on the
-// fault surface. In this case we impose 4.0 m of right-lateral slip
-// for z >= -2.0 km with zero slip for z < -2.0 km. There is no reverse
-// slip or fault opening.
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 3
- value-names = left-lateral-slip reverse-slip fault-opening
- value-units = m m m
- num-locs = 4
- data-dim = 1 // Data is specified along a line.
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0e+3 // Specify coordinates in km for convenience.
- space-dim = 3
- } // cs-data
-} // SimpleDB
-// Columns are
-// (1) x coordinate (km)
-// (2) y coordinate (km)
-// (3) z coordinate (km)
-// (4) left-lateral-slip (m) (right-lateral is negative)
-// (5) reverse-slip (m)
-// (6) fault-opening (m)
-0.0 0.0 0.00 -4.0 0.0 0.0
-0.0 0.0 -2.00 -4.0 0.0 0.0
-0.0 0.0 -2.01 0.0 0.0 0.0
-0.0 0.0 -4.00 0.0 0.0 0.0
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/fixeddisp_axial.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/fixeddisp_axial.spatialdb 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/fixeddisp_axial.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,32 +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
-//
-// dof-0: Ux(x) = 0
-// dof-1: Uy(x) = 1.0 m * x / 3 km
-// dof-2: Uz(x) = 0
-//
-#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/fixeddisp_shear.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/fixeddisp_shear.spatialdb 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/fixeddisp_shear.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,32 +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.
-//
-// dof-0: Ux(x) = 0
-// dof-1: Uy(x) = 1.0 m * x / 3 km
-// dof-2: Uz(x) = 0
-//
-#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
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/geometry.jou
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/geometry.jou 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/geometry.jou 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,47 +0,0 @@
-## /tools/common/cubit-10.2/bin/clarox
-## Cubit Version 10.2
-## Cubit Build 24
-## Revised 12/15/2006 16:09:40 MST
-## Running 06/18/2007 10:26:50 AM
-## Command Options:
-## -warning = On
-## -information = On
-
-# ----------------------------------------------------------------------
-# Create block
-# ----------------------------------------------------------------------
-
-# Block is 6 km x 6 km x 4 km
-# -3 km <= x <= 3 km
-# -3 km <= y <= 3 km
-# -4 km <= z <= 0 km
-reset
-brick x 6000 y 6000 z 4000
-
-# Translate block so the top is at z=0
-volume 1 move x 0 y 0 z -2000
-
-# ----------------------------------------------------------------------
-# Create interface surfaces
-# ----------------------------------------------------------------------
-create planar surface with plane xplane offset 0
-surface 7 name "fault_surface"
-create planar surface with plane zplane offset -2000
-surface 8 name "material_interface"
-
-# ----------------------------------------------------------------------
-# Divide volumes using interface surfaces
-# ----------------------------------------------------------------------
-webcut volume 1 with plane surface fault_surface
-webcut volume 1 with plane surface material_interface
-webcut volume 4 with plane surface material_interface
-volume 1 name "elastic_xpos"
-volume 4 name "elastic_xneg"
-volume 5 name "visco_xpos"
-volume 6 name "visco_xneg"
-
-# ----------------------------------------------------------------------
-# Imprint all volumes, merging surfaces
-# ----------------------------------------------------------------------
-imprint all with volume all
-merge all
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/initial_stress.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/initial_stress.spatialdb 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/initial_stress.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,36 +0,0 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the initial stresses for the model
-// The stresses increase linearly with depth.
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 6 // number of stress components
-
- // names of the stress components
- value-names = stress-xx stress-yy stress-zz stress-xy stress-yz stress-xz
-
- value-units = Pa Pa Pa Pa Pa Pa // units
-
- num-locs = 2 // number of locations: 1 value each at the upper and lower
- // surfaces. Linear interpolation should be used.
-
- data-dim = 1
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0
- space-dim = 3
- }
-}
-// Columns are
-// (1) x coordinate (m)
-// (2) y coordinate (m)
-// (3) z coordinate (m)
-// (4) stress-xx (Pa)
-// (5) stress-yy (Pa)
-// (6) stress-zz (Pa)
-// (7) stress-xy (Pa)
-// (8) stress-yz (Pa)
-// (9) stress-xz (Pa)
-0.0 0.0 -4000.0 -98066500 -98066500 -98066500 0.0 0.0 0.0
-0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/mat_elastic.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/mat_elastic.spatialdb 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/mat_elastic.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,26 +0,0 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the distribution of material
-// properties. In this case, the material properties are uniform.
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 3 // number of material property values
- value-names = density vs vp // names of the material property values
- value-units = kg/m**3 m/s m/s // units
- num-locs = 1 // number of locations
- data-dim = 0
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0
- 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)
-0.0 0.0 0.0 2500.0 3000.0 5291.502622129181
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/mat_genmaxwell.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/mat_genmaxwell.spatialdb 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/mat_genmaxwell.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,33 +0,0 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the distribution of material
-// properties for a Generalized Maxwell viscoelastic material. In
-// this case, the material properties are uniform.
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 9 // number of material property values
- value-names = density vs vp shear-ratio-1 shear-ratio-2 shear-ratio-3 viscosity-1 viscosity-2 viscosity-3 // names of the material property values
- value-units = kg/m**3 m/s m/s None None None Pa*s Pa*s Pa*s // units
- num-locs = 1 // number of locations
- data-dim = 0
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0
- 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) shear ratio 1 (dimensionless)
-// (8) shear ratio 2 (dimensionless)
-// (9) shear ratio 3 (dimensionless)
-// (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
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/mat_maxwell.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/mat_maxwell.spatialdb 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/mat_maxwell.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,30 +0,0 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the distribution of material
-// properties for a Maxwell viscoelastic material. In this case,
-// the material properties are uniform.
-//
-// The viscosity corresponds to a Maxwell time of 100 years.
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 4 // number of material property values
- value-names = density vs vp viscosity // names of the material property values
- value-units = kg/m**3 m/s m/s Pa*s // units
- num-locs = 1 // number of locations
- data-dim = 0
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0
- 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) viscosity (Pa-s)
-0.0 0.0 0.0 2500.0 3000.0 5291.502622129181 7.10046e+19
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/mat_powerlaw.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/mat_powerlaw.spatialdb 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/mat_powerlaw.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,30 +0,0 @@
-// -*- 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
- cs-data = cartesian {
- to-meters = 1.0
- 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
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/mesh_hex8_1000m.jou
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/mesh_hex8_1000m.jou 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/mesh_hex8_1000m.jou 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,104 +0,0 @@
-## /tools/common/cubit-10.2/bin/clarox
-## Cubit Version 10.2
-## Cubit Build 24
-## Revised 12/15/2006 16:09:40 MST
-## Running 06/18/2007 10:26:50 AM
-## Command Options:
-## -warning = On
-## -information = On
-
-# ----------------------------------------------------------------------
-# Generate geometry
-# ----------------------------------------------------------------------
-playback 'geometry.jou'
-
-# ----------------------------------------------------------------------
-# Set discretization size
-# ----------------------------------------------------------------------
-volume all size 1000
-
-# ----------------------------------------------------------------------
-# Generate the mesh
-# ----------------------------------------------------------------------
-mesh volume all
-
-# ----------------------------------------------------------------------
-# Create blocks for materials
-# ----------------------------------------------------------------------
-block 1 volume 1 4
-block 1 name "elastic"
-block 2 volume 5 6
-block 2 name "viscoelastic"
-
-# ----------------------------------------------------------------------
-# Create nodeset for fault
-# ----------------------------------------------------------------------
-group "fault" add node in fault_surface
-group "fault" add node in fault_surface at A
-nodeset 10 group fault
-nodeset 10 name "fault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +x face
-# ----------------------------------------------------------------------
-group "face_xpos" add node in surface 20
-group "face_xpos" add node in surface 28
-nodeset 11 group face_xpos
-nodeset 11 name "face xpos"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -x face
-# ----------------------------------------------------------------------
-group "face_xneg" add node in surface 30
-group "face_xneg" add node in surface 38
-nodeset 12 group face_xneg
-nodeset 12 name "face xneg"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face
-# ----------------------------------------------------------------------
-group "face_ypos" add node in surface 21
-group "face_ypos" add node in surface 27
-group "face_ypos" add node in surface 33
-group "face_ypos" add node in surface 35
-nodeset 13 group face_ypos
-nodeset 13 name "face ypos"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face
-# ----------------------------------------------------------------------
-group "face_yneg" add node in surface 23
-group "face_yneg" add node in surface 25
-group "face_yneg" add node in surface 31
-group "face_yneg" add node in surface 37
-nodeset 14 group face_yneg
-nodeset 14 name "face yneg"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -z face
-# ----------------------------------------------------------------------
-group "face_zneg" add node in surface 12
-group "face_zneg" add node in surface 16
-nodeset 15 group face_zneg
-nodeset 15 name "face zneg"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -z face w/o fault
-# ----------------------------------------------------------------------
-group "face_zneg_nofault" add node in face_zneg
-group "face_zneg_nofault" remove node in fault
-nodeset 16 group face_zneg_nofault
-nodeset 16 name "face zneg nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +z face
-# ----------------------------------------------------------------------
-group "face_zpos" add node in surface 10
-group "face_zpos" add node in surface 17
-nodeset 17 group face_zpos
-nodeset 17 name "face zpos"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "box_hex8_1000m.exo" dimension 3 overwrite
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/powerlaw_properties.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -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
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,94 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-# Turn on some journals to show progress.
-[pylithapp.journal.info]
-timedependent = 1
-implicit = 1
-petsc = 1
-solverlinear = 1
-meshiocubit = 1
-implicitelasticity = 1
-faultcohesivekin = 1
-fiatlagrange = 1
-pylithapp = 1
-materials = 1
-powerlaw3d = 1
-##implicit = 1
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-#debug = 1 ; uncomment to get very verbose mesh information
-
-# Change the default mesh reader to the CUBIT reader.
-reader = pylith.meshio.MeshIOCubit
-
-[pylithapp.mesh_generator.reader]
-# Set filename of mesh to import.
-filename = box_hex8_1000m.exo
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.formulation.time_step]
-# Define the total time for the simulation and the default time step size.
-total_time = 0.0*s ; total time of simulation
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-# Set materials to an array of 2 materials:
-# 'elastic' and 'viscoelastic'
-materials = [elastic,viscoelastic]
-
-# Change from default material (elastic, isotropic, 3-D material) to
-# Maxwell viscoelastic, isotropic 3-D material.
-materials.viscoelastic = pylith.materials.MaxwellIsotropic3D
-
-[pylithapp.timedependent.materials.elastic]
-label = Elastic material
-id = 1
-db_properties.iohandler.filename = mat_elastic.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-# We do not define the database info for the viscoelastic material, since
-# a different database type may be used for some examples.
-[pylithapp.timedependent.materials.viscoelastic]
-label = Viscoelastic material
-id = 2
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-# ----------------------------------------------------------------------
-# PETSc
-# ----------------------------------------------------------------------
-# Set the solver options.
-
-[pylithapp.petsc]
-pc_type = asm
-
-# Change the preconditioner settings.
-sub_pc_factor_shift_type = none
-
-ksp_rtol = 1.0e-8
-ksp_max_it = 500
-ksp_gmres_restart = 50
-
-ksp_monitor = true
-ksp_view = true
-ksp_converged_reason = true
-
-snes_monitor = true
-snes_view = true
-snes_converged_reason = true
-
-log_summary = true
-# start_in_debugger = true
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/sliprate_creep.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/sliprate_creep.spatialdb 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/sliprate_creep.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,32 +0,0 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the distribution of slip rate on
-// the fault surface for a constant slip rate time function. In this
-// case we specify a uniform slip rate of -2.0 cm/yr
-// (6.3376175628057904e-10 m/s) for z < -2.0 km and zero slip for z >=
-// -2.0 km. There is no reverse-slip or fault opening.
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 3
- value-names = left-lateral-slip reverse-slip fault-opening
- value-units = cm/year cm/year cm/year
- num-locs = 4
- data-dim = 1 // Values are specified along a line.
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0e+3 // Specify coordinates in km for convenience.
- space-dim = 3
- } // cs-data
-} // SimpleDB
-// Columns are
-// (1) x coordinate (km)
-// (2) y coordinate (km)
-// (3) z coordinate (km)
-// (4) left-lateral-slip rate (m/s) (right-lateral is negative)
-// (5) reverse-slip rate (m/s)
-// (6) fault-opening rate (m/s)
-0.0 0.0 0.00 0.00 0.0 0.0
-0.0 0.0 -2.00 0.00 0.0 0.0
-0.0 0.0 -2.01 -2.00 0.0 0.0
-0.0 0.0 -4.00 -2.00 0.0 0.0
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/sliptime.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/sliptime.spatialdb 2010-05-06 02:23:35 UTC (rev 16652)
+++ short/3D/PyLith/trunk/examples/3d/hex8/sliptime.spatialdb 2010-05-06 02:31:46 UTC (rev 16653)
@@ -1,26 +0,0 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the distribution of the slip
-// initiation time over the fault surface for the kinematic fault
-// rupture. We specify a slip time of 0.0 s so that the slip occurs at
-// the origin time of the earthquake rupture.
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 1
- value-names = slip-time
- value-units = s
- num-locs = 1
- data-dim = 0 // data is uniform (1 location)
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0
- space-dim = 3
- }
-}
-// Columns are
-// (1) x coordinate (m)
-// (2) y coordinate (m)
-// (3) z coordinate (m)
-// (4) Slip initiation time (s)
-0.0 0.0 0.0 0.0
More information about the CIG-COMMITS
mailing list