[cig-commits] r12773 - short/3D/PyLith/trunk/tests/2d/twofaults_quad4
brad at geodynamics.org
brad at geodynamics.org
Tue Sep 2 10:24:34 PDT 2008
Author: brad
Date: 2008-09-02 10:24:34 -0700 (Tue, 02 Sep 2008)
New Revision: 12773
Added:
short/3D/PyLith/trunk/tests/2d/twofaults_quad4/box.mesh
Removed:
short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_para.mesh
Modified:
short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_para.cfg
short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_perp.cfg
Log:
Updated perp fault test.
Copied: short/3D/PyLith/trunk/tests/2d/twofaults_quad4/box.mesh (from rev 12771, short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_para.mesh)
===================================================================
--- short/3D/PyLith/trunk/tests/2d/twofaults_quad4/box.mesh (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/twofaults_quad4/box.mesh 2008-09-02 17:24:34 UTC (rev 12773)
@@ -0,0 +1,123 @@
+mesh = {
+ dimension = 2
+ use-index-zero = true
+ vertices = {
+ dimension = 2
+ count = 16
+ coordinates = {
+ 0 -3.0 3.0
+ 1 -1.0 3.0
+ 2 1.0 3.0
+ 3 3.0 3.0
+ 4 -3.0 1.0
+ 5 -1.0 1.0
+ 6 1.0 1.0
+ 7 3.0 1.0
+ 8 -3.0 -1.0
+ 9 -1.0 -1.0
+ 10 1.0 -1.0
+ 11 3.0 -1.0
+ 12 -3.0 -3.0
+ 13 -1.0 -3.0
+ 14 1.0 -3.0
+ 15 3.0 -3.0
+ }
+ }
+ cells = {
+ count = 9
+ num-corners = 4
+ simplices = {
+ 0 0 4 5 1
+ 1 1 5 6 2
+ 2 2 6 7 3
+ 3 4 8 9 5
+ 4 5 9 10 6
+ 5 6 10 11 7
+ 6 8 12 13 9
+ 7 9 13 14 10
+ 8 10 14 15 11
+ }
+ material-ids = {
+ 0 10
+ 1 10
+ 2 10
+ 3 10
+ 4 10
+ 5 10
+ 6 10
+ 7 10
+ 8 10
+ }
+ }
+ group = {
+ name = faultL
+ type = vertices
+ count = 4
+ indices = {
+ 1
+ 5
+ 9
+ 13
+ }
+ }
+ group = {
+ name = faultR
+ type = vertices
+ count = 4
+ indices = {
+ 2
+ 6
+ 10
+ 14
+ }
+ }
+ group = {
+ name = faultT
+ type = vertices
+ count = 2
+ indices = {
+ 6
+ 7
+ }
+ }
+ group = {
+ name = faultB
+ type = vertices
+ count = 2
+ indices = {
+ 10
+ 11
+ }
+ }
+ group = {
+ name = edge -x
+ type = vertices
+ count = 4
+ indices = {
+ 0
+ 4
+ 8
+ 12
+ }
+ }
+ group = {
+ name = edge +x
+ type = vertices
+ count = 4
+ indices = {
+ 3
+ 7
+ 11
+ 15
+ }
+ }
+ group = {
+ name = edge +x corners
+ type = vertices
+ count = 2
+ indices = {
+ 3
+ 15
+ }
+ }
+}
Modified: short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_para.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_para.cfg 2008-09-02 16:54:59 UTC (rev 12772)
+++ short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_para.cfg 2008-09-02 17:24:34 UTC (rev 12773)
@@ -23,7 +23,7 @@
#debug = 1
[pylithapp.mesh_generator.importer]
-filename = quad4_para.mesh
+filename = box.mesh
coordsys.space_dim = 2
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_perp.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_perp.cfg 2008-09-02 16:54:59 UTC (rev 12772)
+++ short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_perp.cfg 2008-09-02 17:24:34 UTC (rev 12773)
@@ -6,7 +6,6 @@
# ----------------------------------------------------------------------
[pylithapp.journal.info]
timedependent = 1
-explicit = 1
implicit = 1
petsc = 1
solverlinear = 1
@@ -14,98 +13,95 @@
homogeneous = 1
implicitelasticity = 1
quadrature1d = 1
-fiatsimplex = 1
+fiatlagrange = 1
faultcohesivekin = 1
# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
-debug = 1
+#debug = 1
[pylithapp.mesh_generator.importer]
-filename = quad4h.mesh
+filename = box.mesh
coordsys.space_dim = 2
# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
[pylithapp.timedependent]
-total_time = 0.0*s
-default_dt = 1.0*s
dimension = 2
-formulation = pylith.problems.Implicit
-bc = pylith.bc.BCFourSides
-interfaces = faults
-#interfaces = pylith.faults.SingleFault
-materials = materials
+bc = [neg,pos]
+interfaces = [faultL,faultT,faultB]
+materials = [elastic]
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+dt = 1.0*s
+
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials.one]
-label = elastic material A
+[pylithapp.timedependent.materials]
+elastic = pylith.materials.ElasticPlaneStrain
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
id = 10
-db.iohandler.filename = matprops.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature2D
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 2
-
-[pylithapp.timedependent.materials.two]
-label = elastic material B
-id = 11
db.iohandler.filename = matprops.spatialdb
-quadrature = pylith.feassemble.quadrature.Quadrature2D
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-[pylithapp.timedependent.materials.three]
-label = elastic material C
-id = 12
-db.iohandler.filename = matprops.spatialdb
-quadrature = pylith.feassemble.quadrature.Quadrature2D
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_neg]
+[pylithapp.timedependent.bc.neg]
fixed_dof = [0, 1]
-label = corners
-db.label = Dirichlet BC on corners
-db.iohandler.filename = dislocation_disp.spatialdb
-db.query_type = nearest
+label = edge -x
+db.label = Dirichlet BC on -x edge
+[pylithapp.timedependent.bc.pos]
+fixed_dof = [0, 1]
+label = edge +x corners
+db.label = Dirichlet BC on +x corners
+
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.one]
-id = 1
-label = faultA
+[pylithapp.timedependent.interfaces.faultL]
+id = 100
+label = faultL
quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
quadrature.cell.shape = line
mat_db.iohandler.filename = matprops.spatialdb
-[pylithapp.timedependent.interfaces.one.eq_src.slip_function]
-slip.query_type = nearest
-slip.iohandler.filename = dislocation_slipA.spatialdb
-slip.label = slip for fault A
-slip_rate.iohandler.filename = dislocation_sliprate.spatialdb
-slip_time.iohandler.filename = dislocation_sliptime.spatialdb
+[pylithapp.timedependent.interfaces.faultL.eq_srcs.rupture.slip_function]
+slip.iohandler.filename = slipLL.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
-[pylithapp.timedependent.interfaces.two]
-id = 2
+[pylithapp.timedependent.interfaces.faultT]
+id = 101
+label = faultT
+quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
+quadrature.cell.shape = line
+mat_db.iohandler.filename = matprops.spatialdb
+
+[pylithapp.timedependent.interfaces.faultT.eq_srcs.rupture.slip_function]
+slip.iohandler.filename = slipLL.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
+
+
+[pylithapp.timedependent.interfaces.faultB]
+id = 102
label = faultB
quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
quadrature.cell.shape = line
mat_db.iohandler.filename = matprops.spatialdb
-[pylithapp.timedependent.interfaces.two.eq_src.slip_function]
-slip.label = slip for fault B
-slip.iohandler.filename = dislocation_slipB.spatialdb
-slip_rate.iohandler.filename = dislocation_sliprate.spatialdb
-slip_time.iohandler.filename = dislocation_sliptime.spatialdb
+[pylithapp.timedependent.interfaces.faultB.eq_srcs.rupture.slip_function]
+slip.iohandler.filename = slipRL.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
# ----------------------------------------------------------------------
@@ -120,4 +116,14 @@
# output
# ----------------------------------------------------------------------
[pylithapp.problem.formulation.output.output.writer]
-filename = dislocation.vtk
+filename = dislocation-perp.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = dislocation-perp-statevars-elastic.vtk
+
+[pylithapp.problem.interfaces.faultL.output.writer]
+filename = dislocation-perp-faultL.vtk
+
+[pylithapp.problem.interfaces.faultT.output.writer]
+filename = dislocation-perp-faultT.vtk
Deleted: short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_para.mesh
===================================================================
--- short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_para.mesh 2008-09-02 16:54:59 UTC (rev 12772)
+++ short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_para.mesh 2008-09-02 17:24:34 UTC (rev 12773)
@@ -1,96 +0,0 @@
-mesh = {
- dimension = 2
- use-index-zero = true
- vertices = {
- dimension = 2
- count = 16
- coordinates = {
- 0 -3.0 3.0
- 1 -1.0 3.0
- 2 1.0 3.0
- 3 3.0 3.0
- 4 -3.0 1.0
- 5 -1.0 1.0
- 6 1.0 1.0
- 7 3.0 1.0
- 8 -3.0 -1.0
- 9 -1.0 -1.0
- 10 1.0 -1.0
- 11 3.0 -1.0
- 12 -3.0 -3.0
- 13 -1.0 -3.0
- 14 1.0 -3.0
- 15 3.0 -3.0
- }
- }
- cells = {
- count = 9
- num-corners = 4
- simplices = {
- 0 0 4 5 1
- 1 1 5 6 2
- 2 2 6 7 3
- 3 4 8 9 5
- 4 5 9 10 6
- 5 6 10 11 7
- 6 8 12 13 9
- 7 9 13 14 10
- 8 10 14 15 11
- }
- material-ids = {
- 0 10
- 1 10
- 2 10
- 3 10
- 4 10
- 5 10
- 6 10
- 7 10
- 8 10
- }
- }
- group = {
- name = faultL
- type = vertices
- count = 4
- indices = {
- 1
- 5
- 9
- 13
- }
- }
- group = {
- name = faultR
- type = vertices
- count = 4
- indices = {
- 2
- 6
- 10
- 14
- }
- }
- group = {
- name = edge -x
- type = vertices
- count = 4
- indices = {
- 0
- 4
- 8
- 12
- }
- }
- group = {
- name = edge +x
- type = vertices
- count = 4
- indices = {
- 3
- 7
- 11
- 15
- }
- }
-}
More information about the cig-commits
mailing list