[cig-commits] r13409 - in short/3D/PyLith/benchmarks/trunk/static: . nondimensional
brad at geodynamics.org
brad at geodynamics.org
Tue Nov 25 17:32:15 PST 2008
Author: brad
Date: 2008-11-25 17:32:15 -0800 (Tue, 25 Nov 2008)
New Revision: 13409
Added:
short/3D/PyLith/benchmarks/trunk/static/nondimensional/
short/3D/PyLith/benchmarks/trunk/static/nondimensional/dislocation_slip.spatialdb
short/3D/PyLith/benchmarks/trunk/static/nondimensional/dislocation_slip_norm.spatialdb
short/3D/PyLith/benchmarks/trunk/static/nondimensional/dislocation_sliptime.spatialdb
short/3D/PyLith/benchmarks/trunk/static/nondimensional/matprops.spatialdb
short/3D/PyLith/benchmarks/trunk/static/nondimensional/matprops_norm.spatialdb
short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_dislocation.cfg
short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_dislocation_norm.cfg
short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_geometry.jou
short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_quad4.exo.gz
short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_quad4.jou
short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_quad4_norm.exo.gz
short/3D/PyLith/benchmarks/trunk/static/nondimensional/volume_dislocation.cfg
short/3D/PyLith/benchmarks/trunk/static/nondimensional/volume_geometry.jou
short/3D/PyLith/benchmarks/trunk/static/nondimensional/volume_hex8.jou
Log:
Added nondimensional test problem.
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/dislocation_slip.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/nondimensional/dislocation_slip.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/nondimensional/dislocation_slip.spatialdb 2008-11-26 01:32:15 UTC (rev 13409)
@@ -0,0 +1,13 @@
+#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 {
+ space-dim = 2
+ }
+}
+0.0 0.0 -1.0 0.0
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/dislocation_slip_norm.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/nondimensional/dislocation_slip_norm.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/nondimensional/dislocation_slip_norm.spatialdb 2008-11-26 01:32:15 UTC (rev 13409)
@@ -0,0 +1,13 @@
+#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 {
+ space-dim = 2
+ }
+}
+0.0 0.0 -1.0e-3 0.0
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/dislocation_sliptime.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/nondimensional/dislocation_sliptime.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/nondimensional/dislocation_sliptime.spatialdb 2008-11-26 01:32:15 UTC (rev 13409)
@@ -0,0 +1,13 @@
+#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 {
+ space-dim = 2
+ }
+}
+0.0 0.0 -1.0
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/matprops.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/nondimensional/matprops.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/nondimensional/matprops.spatialdb 2008-11-26 01:32:15 UTC (rev 13409)
@@ -0,0 +1,13 @@
+#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 {
+ space-dim = 2
+ }
+}
+0.0 0.0 2500.0 3000.0 5400.0
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/matprops_norm.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/nondimensional/matprops_norm.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/nondimensional/matprops_norm.spatialdb 2008-11-26 01:32:15 UTC (rev 13409)
@@ -0,0 +1,13 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = density vs vp
+ value-units = none none none
+ num-locs = 1
+ data-dim = 0
+ space-dim = 2
+ cs-data = cartesian {
+ space-dim = 2
+ }
+}
+0.0 0.0 1.0 1.0 1.8
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_dislocation.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_dislocation.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_dislocation.cfg 2008-11-26 01:32:15 UTC (rev 13409)
@@ -0,0 +1,119 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshiocubit = 1
+implicitelasticity = 1
+#quadrature2d = 1
+#fiatlagrange = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+importer = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.importer]
+filename = slice_quad4.exo
+coordsys.space_dim = 2
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 2
+bc = [bc_east,bc_west]
+interfaces = [fault]
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [elastic]
+materials.elastic = pylith.materials.ElasticPlaneStrain
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+
+db = spatialdata.spatialdb.SimpleDB
+db.iohandler.filename = matprops.spatialdb
+
+quadrature = pylith.feassemble.quadrature.Quadrature2D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.bc_east]
+fixed_dof = [0, 1]
+label = 21
+db.label = Dirichlet BC east surface
+
+[pylithapp.timedependent.bc.bc_west]
+fixed_dof = [0, 1]
+label = 22
+db.label = Dirichlet BC west surface
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+fault = pylith.faults.FaultCohesiveKin
+
+[pylithapp.timedependent.interfaces.fault]
+id = 2
+label = 10
+quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+mat_db = spatialdata.spatialdb.SimpleDB
+mat_db.iohandler.filename = matprops.spatialdb
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+slip.iohandler.filename = dislocation_slip.spatialdb
+slip_time.iohandler.filename = dislocation_sliptime.spatialdb
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_rtol = 1.0e-8
+pc_type = asm
+ksp_max_it = 800
+ksp_gmres_restart = 100
+ksp_monitor = true
+ksp_view = true
+log_summary = true
+# start_in_debugger = true
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output.writer]
+filename = slice-dislocation.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = slice-dislocation-statevars-elastic.vtk
+
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = slice-dislocation-fault.vtk
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_dislocation_norm.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_dislocation_norm.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_dislocation_norm.cfg 2008-11-26 01:32:15 UTC (rev 13409)
@@ -0,0 +1,119 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshiocubit = 1
+implicitelasticity = 1
+#quadrature2d = 1
+#fiatlagrange = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+importer = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.importer]
+filename = slice_quad4_norm.exo
+coordsys.space_dim = 2
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 2
+bc = [bc_east,bc_west]
+interfaces = [fault]
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [elastic]
+materials.elastic = pylith.materials.ElasticPlaneStrain
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+
+db = spatialdata.spatialdb.SimpleDB
+db.iohandler.filename = matprops_norm.spatialdb
+
+quadrature = pylith.feassemble.quadrature.Quadrature2D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.bc_east]
+fixed_dof = [0, 1]
+label = 21
+db.label = Dirichlet BC east surface
+
+[pylithapp.timedependent.bc.bc_west]
+fixed_dof = [0, 1]
+label = 22
+db.label = Dirichlet BC west surface
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+fault = pylith.faults.FaultCohesiveKin
+
+[pylithapp.timedependent.interfaces.fault]
+id = 2
+label = 10
+quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+mat_db = spatialdata.spatialdb.SimpleDB
+mat_db.iohandler.filename = matprops_norm.spatialdb
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+slip.iohandler.filename = dislocation_slip_norm.spatialdb
+slip_time.iohandler.filename = dislocation_sliptime.spatialdb
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_rtol = 1.0e-8
+pc_type = asm
+ksp_max_it = 800
+ksp_gmres_restart = 100
+ksp_monitor = true
+ksp_view = true
+log_summary = true
+# start_in_debugger = true
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output.writer]
+filename = slice-dislocation-norm.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = slice-dislocation-norm-statevars-elastic.vtk
+
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = slice-dislocation-norm-fault.vtk
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_geometry.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_geometry.jou 2008-11-26 01:32:15 UTC (rev 13409)
@@ -0,0 +1,23 @@
+# ----------------------------------------------------------------------
+# Create surface using vertices
+# ----------------------------------------------------------------------
+
+# Surface is 100km x 100km
+# -50 km <= x <= 50 km
+# -50 km <= y <= 50 km
+# Geometry is created in units of km (scale to meters at end).
+reset
+create vertex -50.0 -50.0 0.0
+create vertex -50.0 +50.0 0.0
+create vertex +50.0 +50.0 0.0
+create vertex +50.0 -50.0 0.0
+create surface vertex 1 2 3 4
+delete vertex all
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create vertex 0.0 -50.0 0.0
+create vertex 0.0 +50.0 0.0
+split surface 1 across location vertex 9 location vertex 10
+delete vertex all
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_quad4.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_quad4.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_quad4.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_quad4.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_quad4.jou 2008-11-26 01:32:15 UTC (rev 13409)
@@ -0,0 +1,58 @@
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'slice_geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+surface all size 1.0
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+surface all scheme submap
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Generate blocks and nodesets
+# ----------------------------------------------------------------------
+block 1 surface 2 3
+block 1 name "elastic"
+
+# ----------------------------------------------------------------------
+# Create nodeset for fault
+# ----------------------------------------------------------------------
+group "fault" add node in curve 5
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for BC
+# ----------------------------------------------------------------------
+group "surf_east" add node in curve 3
+nodeset 21 group surf_east
+nodeset 21 name "surf_east"
+
+group "surf_west" add node in curve 1
+nodeset 22 group surf_west
+nodeset 22 name "surf_west"
+
+group "surf_north" add node in curve 7
+group "surf_north" add node in curve 8
+group "surf_north" remove node in fault
+nodeset 23 group surf_north
+nodeset 23 name "surf_north"
+
+group "surf_south" add node in curve 6
+group "surf_south" add node in curve 9
+group "surf_south" remove node in fault
+nodeset 24 group surf_south
+nodeset 24 name "surf_south"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "slice_quad4_normalized.exo" dimension 2 overwrite
+volume all scale 1.0e+3
+export mesh "slice_quad4.exo" dimension 2 overwrite
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_quad4_norm.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/static/nondimensional/slice_quad4_norm.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/volume_dislocation.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/nondimensional/volume_dislocation.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/nondimensional/volume_dislocation.cfg 2008-11-26 01:32:15 UTC (rev 13409)
@@ -0,0 +1,145 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshiocubit = 1
+implicitelasticity = 1
+#quadrature3d = 1
+#fiatlagrange = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+importer = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.importer]
+filename = volume_hex8.exo
+
+coordsys = spatialdata.geocoords.CSGeoProj
+coordsys.projector.projection = tmerc
+coordsys.projector.proj-options = +lon_0=-118.2467 +lat_0=34.0550
+coordsys.datum_vert = mean sea level
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 3
+bc = [bc_east,bc_west,bc_bottom]
+interfaces = [fault]
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [elastic]
+materials.elastic = pylith.materials.ElasticIsotropic3D
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+
+db = spatialdata.spatialdb.SCECCVMH
+db.data_dir = /home/brad/data/sceccvm-h/vx53/bin
+db.squash = True
+db.squash_limit = -8.0*km
+#db = spatialdata.spatialdb.SimpleDB
+#db.iohandler.filename = mat_homogeneous.spatialdb
+
+quadrature = pylith.feassemble.quadrature.Quadrature3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.bc_east]
+fixed_dof = [0, 1]
+label = 21
+db = spatialdata.spatialdb.UniformDB
+db.label = Dirichlet BC east surface
+db.values = [dof-0,dof-1]
+db.data = [0.0,-2.0]
+
+[pylithapp.timedependent.bc.bc_west]
+fixed_dof = [0, 1]
+label = 22
+db = spatialdata.spatialdb.UniformDB
+db.label = Dirichlet BC west surface
+db.values = [dof-0,dof-1]
+db.data = [0.0,+2.0]
+
+[pylithapp.timedependent.bc.bc_bottom]
+fixed_dof = [2]
+label = 26
+db = spatialdata.spatialdb.UniformDB
+db.label = Dirichlet BC bottom surface
+db.values = [dof-2]
+db.data = [0.0]
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+fault = pylith.faults.FaultCohesiveKin
+
+[pylithapp.timedependent.interfaces.fault]
+id = 2
+label = 10
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+#mat_db = spatialdata.spatialdb.SCECCVMH
+#mat_db.data_dir = /home/brad/data/sceccvm-h/vx53/bin
+#mat_db.squash = True
+#mat_db.squash_limit = -8.0*km
+mat_db = spatialdata.spatialdb.SimpleDB
+mat_db.iohandler.filename = mat_homogeneous.spatialdb
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+slip.iohandler.filename = dislocation_slip.spatialdb
+slip_time.iohandler.filename = dislocation_sliptime.spatialdb
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_rtol = 1.0e-8
+pc_type = asm
+ksp_max_it = 200
+ksp_gmres_restart = 50
+ksp_monitor = true
+ksp_view = true
+log_summary = true
+# start_in_debugger = true
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output.writer]
+filename = volume-dislocation.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = volume-dislocation-statevars-elastic.vtk
+
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = volume-dislocation-fault.vtk
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/volume_geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/nondimensional/volume_geometry.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/nondimensional/volume_geometry.jou 2008-11-26 01:32:15 UTC (rev 13409)
@@ -0,0 +1,28 @@
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+reset
+brick x 48000 y 48000 z 16000
+volume 1 move x -1000 y -30000 z -8000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create vertex -8000.0 -6000.0 0.0
+create vertex 6000.0 -54000.0 0.0
+create vertex 6000.0 -54000.0 -16000.0
+create planar surface with plane vertex 9 vertex 10 vertex 11
+delete vertex all
+surface 7 name "fault_surface"
+
+# ----------------------------------------------------------------------
+# Divide volumes using interfaces surfaces
+# ----------------------------------------------------------------------
+webcut volume 1 with plane surface fault_surface
+
+# ----------------------------------------------------------------------
+# Imprint curves on volumes
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+delete body 2
Added: short/3D/PyLith/benchmarks/trunk/static/nondimensional/volume_hex8.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/nondimensional/volume_hex8.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/nondimensional/volume_hex8.jou 2008-11-26 01:32:15 UTC (rev 13409)
@@ -0,0 +1,67 @@
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'volume_geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+surface 9 11 15 16 scheme pave
+surface all size 1000
+volume all size 1000
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Generate blocks and nodesets
+# ----------------------------------------------------------------------
+block 1 volume 1 3
+block 1 name "elastic"
+
+# ----------------------------------------------------------------------
+# Create nodeset for fault
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for BC
+# ----------------------------------------------------------------------
+group "surf_east" add node in surface 6
+nodeset 21 group surf_east
+nodeset 21 name "surf_east"
+
+group "surf_west" add node in surface 4
+nodeset 22 group surf_west
+nodeset 22 name "surf_west"
+
+group "surf_north" add node in surface 14
+group "surf_north" add node in surface 12
+group "surf_north" remove node in fault
+nodeset 23 group surf_north
+nodeset 23 name "surf_north"
+
+group "surf_south" add node in surface 17
+group "surf_south" add node in surface 10
+group "surf_south" remove node in fault
+nodeset 24 group surf_south
+nodeset 24 name "surf_south"
+
+group "surf_ground" add node in surface 16
+group "surf_ground" add node in surface 9
+nodeset 25 group surf_ground
+nodeset 25 name "surf_ground"
+
+group "surf_bottom" add node in surface 15
+group "surf_bottom" add node in surface 11
+nodeset 26 group surf_bottom
+nodeset 26 name "surf_bottom"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "volume_hex8.exo" dimension 3 overwrite
More information about the CIG-COMMITS
mailing list