[cig-commits] r16962 - in short/3D/PyLith/trunk/tests/2d: . slipdir
brad at geodynamics.org
brad at geodynamics.org
Wed Jun 9 16:01:15 PDT 2010
Author: brad
Date: 2010-06-09 16:01:14 -0700 (Wed, 09 Jun 2010)
New Revision: 16962
Added:
short/3D/PyLith/trunk/tests/2d/slipdir/
short/3D/PyLith/trunk/tests/2d/slipdir/fault_x.cfg
short/3D/PyLith/trunk/tests/2d/slipdir/fault_x.exo
short/3D/PyLith/trunk/tests/2d/slipdir/fault_x.jou
short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy1.cfg
short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy1.exo
short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy1.jou
short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy2.cfg
short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy2.exo
short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy2.jou
short/3D/PyLith/trunk/tests/2d/slipdir/fault_y.cfg
short/3D/PyLith/trunk/tests/2d/slipdir/fault_y.exo
short/3D/PyLith/trunk/tests/2d/slipdir/fault_y.jou
short/3D/PyLith/trunk/tests/2d/slipdir/matprops.spatialdb
short/3D/PyLith/trunk/tests/2d/slipdir/pylithapp.cfg
short/3D/PyLith/trunk/tests/2d/slipdir/slip.spatialdb
short/3D/PyLith/trunk/tests/2d/slipdir/sliptime.spatialdb
Log:
Added manual test for slip direction.
Added: short/3D/PyLith/trunk/tests/2d/slipdir/fault_x.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/slipdir/fault_x.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/slipdir/fault_x.cfg 2010-06-09 23:01:14 UTC (rev 16962)
@@ -0,0 +1,17 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = fault_x.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = fault_x.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = fault_x-fault.vtk
Added: short/3D/PyLith/trunk/tests/2d/slipdir/fault_x.exo
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/trunk/tests/2d/slipdir/fault_x.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/trunk/tests/2d/slipdir/fault_x.jou
===================================================================
--- short/3D/PyLith/trunk/tests/2d/slipdir/fault_x.jou (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/slipdir/fault_x.jou 2010-06-09 23:01:14 UTC (rev 16962)
@@ -0,0 +1,44 @@
+# Create mesh with +x fault normal.
+
+reset
+
+brick x 200 y 200 z 200
+move volume 1 x 0 y 0 z 100
+
+create surface extended from surface 6
+move surface 7 x -100 y 0 z 0
+webcut volume 1 with plane surface 7
+delete surface 7
+
+create midsurface volume 3 surface 15 16
+create midsurface volume 1 surface 9 11
+
+delete volume 1
+delete volume 3
+
+imprint all
+merge all
+
+surface all size 20
+surface all scheme pave
+mesh surface all
+
+# Create blocks and nodesets
+
+block 1 surface all
+
+group "fault" add node in curve 33
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+group "boundary_nofault" add node in curve 34
+group "boundary_nofault" add node in curve 35
+group "boundary_nofault" add node in curve 36
+group "boundary_nofault" add node in curve 40
+group "boundary_nofault" add node in curve 39
+group "boundary_nofault" add node in curve 38
+group "boundary_nofault" remove node in group fault
+nodeset 1 group boundary_nofault
+nodeset 1 name "boundary_nofault"
+
+export mesh "fault_x.exo" dimension 2 overwrite
Added: short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy1.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy1.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy1.cfg 2010-06-09 23:01:14 UTC (rev 16962)
@@ -0,0 +1,24 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = fault_xy1.exo
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+
+up_dir = [0.0,1.0,0.0]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = fault_xy1.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = fault_xy1-fault.vtk
Added: short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy1.exo
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy1.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy1.jou
===================================================================
--- short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy1.jou (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy1.jou 2010-06-09 23:01:14 UTC (rev 16962)
@@ -0,0 +1,43 @@
+# Create mesh with +xy fault normal.
+
+reset
+
+brick x 200 y 200 z 200
+move volume 1 x 0 y 0 z 100
+
+create surface extended from surface 6
+move surface 7 x -100 y 0 z 0
+surface 7 rotate +45 about z
+webcut volume 1 with plane surface 7
+delete surface 7
+
+create midsurface volume 3 surface 12 13
+create midsurface volume 1 surface 9 10
+
+delete volume 1
+delete volume 3
+
+imprint all
+merge all
+
+surface all size 20
+surface all scheme pave
+mesh surface all
+
+# Create blocks and nodesets
+
+block 1 surface all
+
+group "fault" add node in curve 23
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+group "boundary_nofault" add node in curve 24
+group "boundary_nofault" add node in curve 25
+group "boundary_nofault" add node in curve 27
+group "boundary_nofault" add node in curve 28
+group "boundary_nofault" remove node in group fault
+nodeset 1 group boundary_nofault
+nodeset 1 name "boundary_nofault"
+
+export mesh "fault_xy1.exo" dimension 2 overwrite
Added: short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy2.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy2.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy2.cfg 2010-06-09 23:01:14 UTC (rev 16962)
@@ -0,0 +1,24 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = fault_xy2.exo
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+
+up_dir = [0.0,1.0,0.0]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = fault_xy2.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = fault_xy2-fault.vtk
Added: short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy2.exo
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy2.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy2.jou
===================================================================
--- short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy2.jou (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/slipdir/fault_xy2.jou 2010-06-09 23:01:14 UTC (rev 16962)
@@ -0,0 +1,43 @@
+# Create mesh with +xy fault normal.
+
+reset
+
+brick x 200 y 200 z 200
+move volume 1 x 0 y 0 z 100
+
+create surface extended from surface 6
+move surface 7 x -100 y 0 z 0
+surface 7 rotate -45 about z
+webcut volume 1 with plane surface 7
+delete surface 7
+
+create midsurface volume 3 surface 12 13
+create midsurface volume 1 surface 9 10
+
+delete volume 1
+delete volume 3
+
+imprint all
+merge all
+
+surface all size 20
+surface all scheme pave
+mesh surface all
+
+# Create blocks and nodesets
+
+block 1 surface all
+
+group "fault" add node in curve 23
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+group "boundary_nofault" add node in curve 24
+group "boundary_nofault" add node in curve 25
+group "boundary_nofault" add node in curve 27
+group "boundary_nofault" add node in curve 28
+group "boundary_nofault" remove node in group fault
+nodeset 1 group boundary_nofault
+nodeset 1 name "boundary_nofault"
+
+export mesh "fault_xy2.exo" dimension 2 overwrite
Added: short/3D/PyLith/trunk/tests/2d/slipdir/fault_y.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/slipdir/fault_y.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/slipdir/fault_y.cfg 2010-06-09 23:01:14 UTC (rev 16962)
@@ -0,0 +1,17 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = fault_y.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = fault_y.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = fault_y-fault.vtk
Added: short/3D/PyLith/trunk/tests/2d/slipdir/fault_y.exo
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/trunk/tests/2d/slipdir/fault_y.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/trunk/tests/2d/slipdir/fault_y.jou
===================================================================
--- short/3D/PyLith/trunk/tests/2d/slipdir/fault_y.jou (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/slipdir/fault_y.jou 2010-06-09 23:01:14 UTC (rev 16962)
@@ -0,0 +1,45 @@
+# Create mesh with +y fault normal.
+
+reset
+
+brick x 200 y 200 z 200
+move volume 1 x 0 y 0 z 100
+
+create surface extended from surface 5
+move surface 7 x 0 y -100 z 0
+webcut volume 1 with plane surface 7
+delete surface 7
+
+create midsurface volume 3 surface 15 16
+create midsurface volume 1 surface 9 11
+
+delete volume 1
+delete volume 3
+
+imprint all
+merge all
+
+surface all size 20
+surface all scheme pave
+mesh surface all
+
+# Create blocks and nodesets
+
+block 1 surface all
+
+group "fault" add node in curve 33
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+group "boundary_nofault" add node in curve 35
+group "boundary_nofault" add node in curve 36
+group "boundary_nofault" add node in curve 40
+group "boundary_nofault" add node in curve 39
+group "boundary_nofault" add node in curve 38
+group "boundary_nofault" add node in curve 34
+group "boundary_nofault" remove node in group fault
+nodeset 1 group boundary_nofault
+nodeset 1 name "boundary_nofault"
+
+export mesh "fault_y.exo" dimension 2 overwrite
+
Added: short/3D/PyLith/trunk/tests/2d/slipdir/matprops.spatialdb
===================================================================
--- short/3D/PyLith/trunk/tests/2d/slipdir/matprops.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/slipdir/matprops.spatialdb 2010-06-09 23:01:14 UTC (rev 16962)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = density vs vp
+ value-units = kg/m**3 m/s m/s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 2
+ }
+}
+0.0 0.0 2500.0 1000.0 1732.0
Added: short/3D/PyLith/trunk/tests/2d/slipdir/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/slipdir/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/slipdir/pylithapp.cfg 2010-06-09 23:01:14 UTC (rev 16962)
@@ -0,0 +1,119 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+meshiocubit = 1
+faultcohesivekin = 1
+
+[pylithapp.journal.debug]
+pylithapp = 1
+problem = 1
+explicit = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+coordsys.space_dim = 2
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 2
+
+bc = [boundary]
+interfaces = [fault]
+
+materials = [elastic]
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+dt = 1.0*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+elastic = pylith.materials.ElasticPlaneStrain
+
+[pylithapp.timedependent.materials.elastic]
+
+label = Elastic material
+id = 1
+
+db_properties.iohandler.filename = matprops.spatialdb
+
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.boundary]
+
+bc_dof = [0,1]
+label = boundary_nofault
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+
+id = 100
+label = fault
+
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+eq_srcs.rupture.slip_function = pylith.faults.StepSlipFn
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+
+slip.iohandler.filename = slip.spatialdb
+slip_time.iohandler.filename = sliptime.spatialdb
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_type = gmres
+pc_type = asm
+
+sub_pc_factor_shift_type = nonzero
+
+ksp_rtol = 1.0e-8
+ksp_atol = 1.0e-10
+ksp_max_it = 50
+ksp_gmres_restart = 25
+
+ksp_monitor = true
+#ksp_view = true
+log_summary = true
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+vertex_data_fields = [displacement]
+writer.time_format = %05.2f
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.time_format = %05.2f
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = null.vtk
+cell_info_fields = []
+cell_data_fields = []
Added: short/3D/PyLith/trunk/tests/2d/slipdir/slip.spatialdb
===================================================================
--- short/3D/PyLith/trunk/tests/2d/slipdir/slip.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/slipdir/slip.spatialdb 2010-06-09 23:01:14 UTC (rev 16962)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 2
+ value-names = left-lateral-slip fault-opening
+ value-units = m m
+ num-locs = 1
+ data-dim = 0
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 2
+ }
+}
+0.0 0.0 1.0 0.0
Added: short/3D/PyLith/trunk/tests/2d/slipdir/sliptime.spatialdb
===================================================================
--- short/3D/PyLith/trunk/tests/2d/slipdir/sliptime.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/slipdir/sliptime.spatialdb 2010-06-09 23:01:14 UTC (rev 16962)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 1
+ value-names = slip-time
+ value-units = s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 2
+ }
+}
+0.0 0.0 -1.0
More information about the CIG-COMMITS
mailing list