[cig-commits] r16656 - in short/3D/PyLith/trunk/examples/3d/hex8: . mesh spatialdb

brad at geodynamics.org brad at geodynamics.org
Thu May 6 10:05:56 PDT 2010


Author: brad
Date: 2010-05-06 10:05:56 -0700 (Thu, 06 May 2010)
New Revision: 16656

Added:
   short/3D/PyLith/trunk/examples/3d/hex8/mesh/
   short/3D/PyLith/trunk/examples/3d/hex8/mesh/box_hex8_1000m.exo
   short/3D/PyLith/trunk/examples/3d/hex8/mesh/geometry.jou
   short/3D/PyLith/trunk/examples/3d/hex8/mesh/mesh_hex8_1000m.jou
   short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/
   short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/finalslip.spatialdb
   short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/finalslip_rupture.spatialdb
   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/initial_stress.spatialdb
   short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_elastic.spatialdb
   short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_genmaxwell.spatialdb
   short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_maxwell.spatialdb
   short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_powerlaw.spatialdb
   short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw_properties.spatialdb
   short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/sliprate_creep.spatialdb
   short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/sliptime.spatialdb
Removed:
   short/3D/PyLith/trunk/examples/3d/hex8/common/
Log:
Refactored common into mesh and spatialdb directories.

Copied: short/3D/PyLith/trunk/examples/3d/hex8/mesh/box_hex8_1000m.exo (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/box_hex8_1000m.exo)
===================================================================
(Binary files differ)

Copied: short/3D/PyLith/trunk/examples/3d/hex8/mesh/geometry.jou (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/geometry.jou)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/mesh/geometry.jou	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/mesh/geometry.jou	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/mesh/mesh_hex8_1000m.jou (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/mesh_hex8_1000m.jou)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/mesh/mesh_hex8_1000m.jou	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/mesh/mesh_hex8_1000m.jou	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/spatialdb/finalslip.spatialdb (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/finalslip.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/finalslip.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/finalslip.spatialdb	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/spatialdb/finalslip_rupture.spatialdb (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/finalslip_rupture.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/finalslip_rupture.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/finalslip_rupture.spatialdb	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/spatialdb/fixeddisp_axial.spatialdb (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/fixeddisp_axial.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/fixeddisp_axial.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/fixeddisp_axial.spatialdb	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/spatialdb/fixeddisp_shear.spatialdb (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/fixeddisp_shear.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/fixeddisp_shear.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/fixeddisp_shear.spatialdb	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/spatialdb/initial_stress.spatialdb (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/initial_stress.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/initial_stress.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/initial_stress.spatialdb	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/spatialdb/mat_elastic.spatialdb (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/mat_elastic.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_elastic.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_elastic.spatialdb	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/spatialdb/mat_genmaxwell.spatialdb (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/mat_genmaxwell.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_genmaxwell.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_genmaxwell.spatialdb	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/spatialdb/mat_maxwell.spatialdb (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/mat_maxwell.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_maxwell.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_maxwell.spatialdb	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/spatialdb/mat_powerlaw.spatialdb (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/mat_powerlaw.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_powerlaw.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/mat_powerlaw.spatialdb	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/spatialdb/powerlaw_properties.spatialdb (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/powerlaw_properties.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw_properties.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw_properties.spatialdb	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/spatialdb/sliprate_creep.spatialdb (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/sliprate_creep.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/sliprate_creep.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/sliprate_creep.spatialdb	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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/spatialdb/sliptime.spatialdb (from rev 16654, short/3D/PyLith/trunk/examples/3d/hex8/common/sliptime.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/sliptime.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/sliptime.spatialdb	2010-05-06 17:05:56 UTC (rev 16656)
@@ -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



More information about the CIG-COMMITS mailing list