[cig-commits] r7413 - short/3D/PyLith/trunk/examples/twocells/twohex8

willic3 at geodynamics.org willic3 at geodynamics.org
Fri Jun 22 22:09:03 PDT 2007


Author: willic3
Date: 2007-06-22 22:09:03 -0700 (Fri, 22 Jun 2007)
New Revision: 7413

Modified:
   short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation.cfg
   short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_disp.spatialdb
   short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_slip.spatialdb
   short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_sliprate.spatialdb
   short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_sliptime.spatialdb
Log:
Fixed fault files.
All problems now appear to run correctly.



Modified: short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation.cfg	2007-06-23 04:49:03 UTC (rev 7412)
+++ short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation.cfg	2007-06-23 05:09:03 UTC (rev 7413)
@@ -17,7 +17,7 @@
 # Journal settings in addition to those given in 'pylithapp.cfg'
 # ----------------------------------------------------------------------
 [pylithapp.journal.info]
-quadrature1d = 1
+quadrature2d = 1
 faultcohesivekin = 1
 
 # ----------------------------------------------------------------------
@@ -30,8 +30,9 @@
 # We want an implicit formulation.
 formulation = pylith.problems.Implicit
 
-# This is a container for a boundary condition specified at a point.
-bc = pylith.bc.BCFourSides
+# This is a container for a boundary condition specified on six sides
+# of a rectangular solid domain.
+bc = pylith.bc.BCSixSides
 
 # This is a container for a single fault.
 interfaces = pylith.faults.SingleFault
@@ -44,39 +45,38 @@
 # Boundary conditions to be applied to the negative x-side of the mesh.
 [pylithapp.timedependent.bc.x_neg]
 
-# We are fixing the 0 (x) and 1 (y) degrees of freedom.
-fixed_dof = [0, 1]
+# We are fixing all 3 degrees of freedom.
+fixed_dof = [0, 1, 2]
 
 # The nodes associated with this boundary condition have the name
 # 'x_neg' in the mesh file.
 label = x_neg
 
-# We are assigning the label 'Dirichlet BC' to the database.
-db.label = Dirichlet BC -x edge
+# We are assigning the label 'Dirichlet BC -x face' to the database.
+db.label = Dirichlet BC -x face
 
 # The name of the file containing the spatial database for the BC
 # specification.
 db.iohandler.filename = dislocation_disp.spatialdb
-db.query_type = linear
 
 
 # Boundary conditions to be applied to the positive x-side of the mesh.
 [pylithapp.timedependent.bc.x_pos]
 
-# We are fixing the 0 (x) and 1 (y) degrees of freedom.
-fixed_dof = [0, 1]
+# We are fixing all 3 degrees of freedom.
+fixed_dof = [0, 1, 2]
 
 # The nodes associated with this boundary condition have the name
 # 'x_pos' in the mesh file.
 label = x_pos
 
-# We are assigning the label 'Dirichlet BC' to the database.
-db.label = Dirichlet BC +x edge
+# We are assigning the label 'Dirichlet BC +x face' to the database.
+db.label = Dirichlet BC +x face
 
 # The name of the file containing the spatial database for the BC
 # specification.
 db.iohandler.filename = dislocation_disp.spatialdb
-db.query_type = linear
+
 # ----------------------------------------------------------------------
 # faults
 # Provide information on the fault (interface).
@@ -96,9 +96,10 @@
 # system of equations.
 mat_db.iohandler.filename = matprops.spatialdb
 
-# The quadrature for a 2D fault is 1D with a linear shape.
-quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
-quadrature.cell.shape = line
+# The quadrature for a 3D fault is 2D.
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
 
 # Give the spatial databases specifying the fault slip.
 [pylithapp.timedependent.interfaces.fault.eq_src.slip_function]

Modified: short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_disp.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_disp.spatialdb	2007-06-23 04:49:03 UTC (rev 7412)
+++ short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_disp.spatialdb	2007-06-23 05:09:03 UTC (rev 7413)
@@ -6,13 +6,13 @@
 // This follows the format for a Simple DB (the only type presently available).
 SimpleDB {
 
-  // There are two values specified in the database, corresponding to the
-  // constraint values for the x (dof-0) and y (dof-1) degrees of freedom.
-  num-values = 2
-  value-names =  dof-0 dof-1
+  // There are 3 values specified in the database, corresponding to the
+  // constraint values for the x, y, and z degrees of freedom.
+  num-values = 3
+  value-names =  dof-0 dof-1 dof-2
 
   // The constraint values (displacements) have units of meters.
-  value-units =  m  m
+  value-units =  m  m  m
 
   // The values are specified at two spatial locations.
   num-locs = 2
@@ -21,8 +21,8 @@
   // is being specified at points.
   data-dim = 1
 
-  // The spatial dimension of the database is 2.
-  space-dim = 2
+  // The spatial dimension of the database is 3.
+  space-dim = 3
 
   // We are specifying the data in a Cartesian coordinate system.
   cs-data = cartesian {
@@ -30,15 +30,15 @@
     // Our units are already in meters, so we can just multiply by one.
     to-meters = 1.0
 
-    // We are using a 2D Cartesian coordinate system.
-    space-dim = 2
+    // We are using a 3D Cartesian coordinate system.
+    space-dim = 3
   }
 }
 
 // This is where the data is specified.
 // As described in dislocation.cfg, this database will be used to specify
-// data along an edge (pylith.bc.BCFourSides).
+// data along faces (pylith.bc.BCSixSides).
 // The entries are:
-// X-coord, Y-coord, dof-0 (x constraint), dof-1 (y constraint).
--2.0  0.0     0.0  0.0
- 2.0  0.0     0.0  0.0
+// X-coord, Y-coord, Z-coord, dof-0, dof-1, dof-2.
+-1.0  0.0  0.0     0.0  0.0  0.0
+ 1.0  0.0  0.0     0.0  0.0  0.0

Modified: short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_slip.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_slip.spatialdb	2007-06-23 04:49:03 UTC (rev 7412)
+++ short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_slip.spatialdb	2007-06-23 05:09:03 UTC (rev 7413)
@@ -4,13 +4,14 @@
 // dislocation (kinematic fault) example.
 SimpleDB {
 
-  // There are two values specified in the database, corresponding to the
-  // slip values for the left-lateral and fault-opening components.
-  num-values = 2
-  value-names =  left-lateral-slip fault-opening
+  // There are 3 values specified in the database, corresponding to the
+  // slip values for the left-lateral, reverse-slip, and fault-opening
+  // components.
+  num-values = 3
+  value-names =  left-lateral-slip  reverse-slip  fault-opening
 
   // The fault slip values have units of meters.
-  value-units =  m  m
+  value-units =  m  m  m
 
   // The values are specified at one spatial location.
   num-locs = 1
@@ -19,8 +20,8 @@
   // is specified for all points in the set.
   data-dim = 0
 
-  // The spatial dimension of the database is 2.
-  space-dim = 2
+  // The spatial dimension of the database is 3.
+  space-dim = 3
 
   // We are specifying the data in a Cartesian coordinate system.
   cs-data = cartesian {
@@ -28,8 +29,8 @@
     // Our units are already in meters, so we can just multiply by one.
     to-meters = 1.0
 
-    // We are using a 2D Cartesian coordinate system.
-    space-dim = 2
+    // We are using a 3D Cartesian coordinate system.
+    space-dim = 3
   }
 }
 
@@ -37,5 +38,5 @@
 // As described in dislocation.cfg, this database will be used to specify
 // uniform data for a fault (pylith.faults.SingleFault).
 // The entries are:
-// X-Coord, Y-Coord, left-lateral slip, fault-opening slip
-0.0  0.0   0.01  0.0
+// X-Coord, Y-Coord, X-Coord, left-lateral slip, reverse slip, fault-opening slip
+0.0  0.0  0.0   0.01  0.0  0.0

Modified: short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_sliprate.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_sliprate.spatialdb	2007-06-23 04:49:03 UTC (rev 7412)
+++ short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_sliprate.spatialdb	2007-06-23 05:09:03 UTC (rev 7413)
@@ -19,8 +19,8 @@
   // is specified for all points in the set.
   data-dim = 0
 
-  // The spatial dimension of the database is 2.
-  space-dim = 2
+  // The spatial dimension of the database is 3.
+  space-dim = 3
 
   // We are specifying the data in a Cartesian coordinate system.
   cs-data = cartesian {
@@ -28,8 +28,8 @@
     // Our units are already in meters, so we can just multiply by one.
     to-meters = 1.0
 
-    // We are using a 2D Cartesian coordinate system.
-    space-dim = 2
+    // We are using a 3D Cartesian coordinate system.
+    space-dim = 3
   }
 }
 
@@ -37,5 +37,5 @@
 // As described in dislocation.cfg, this database will be used to specify
 // uniform data for a fault (pylith.faults.SingleFault).
 // The entries are:
-// X-Coord, Y-Coord, slip rate
-0.0  0.0  1.0e+6
+// X-Coord, Y-Coord, Z-Coord, slip rate
+0.0  0.0  0.0  1.0e+6

Modified: short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_sliptime.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_sliptime.spatialdb	2007-06-23 04:49:03 UTC (rev 7412)
+++ short/3D/PyLith/trunk/examples/twocells/twohex8/dislocation_sliptime.spatialdb	2007-06-23 05:09:03 UTC (rev 7413)
@@ -19,8 +19,8 @@
   // is specified for all points in the set.
   data-dim = 0
 
-  // The spatial dimension of the database is 2.
-  space-dim = 2
+  // The spatial dimension of the database is 3.
+  space-dim = 3
 
   // We are specifying the data in a Cartesian coordinate system.
   cs-data = cartesian {
@@ -28,8 +28,8 @@
     // Our units are already in meters, so we can just multiply by one.
     to-meters = 1.0
 
-    // We are using a 2D Cartesian coordinate system.
-    space-dim = 2
+    // We are using a 3D Cartesian coordinate system.
+    space-dim = 3
   }
 }
 
@@ -37,5 +37,5 @@
 // As described in dislocation.cfg, this database will be used to specify
 // uniform data for a fault (pylith.faults.SingleFault).
 // The entries are:
-// X-Coord, Y-Coord, slip time
-0.0  0.0   -1.0
+// X-Coord, Y-Coord, Z-Coord, slip time
+0.0  0.0  0.0   -1.0



More information about the cig-commits mailing list