[cig-commits] r19397 - in short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup: . tpv16
brad at geodynamics.org
brad at geodynamics.org
Fri Jan 20 08:04:38 PST 2012
Author: brad
Date: 2012-01-20 08:04:37 -0800 (Fri, 20 Jan 2012)
New Revision: 19397
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/create_spatialdb.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/geometry.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/gradient.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/hex8.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/hex8_150m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/matprops.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tet4.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tet4_150m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_friction.spatialdb.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_hex8_075m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_hex8_150m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_tet4_075m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_tet4_150m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_tractions.spatialdb.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv17.cfg
Log:
Added tpv16.
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/create_spatialdb.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/create_spatialdb.py (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/create_spatialdb.py 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,86 @@
+#!/usr/bin/env python
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard
+# U.S. Geological Survey
+#
+# ----------------------------------------------------------------------
+#
+
+# Requires input files from SCEC website.
+# http://scecdata.usc.edu/cvws/download/tpv16/tpv16_17_input_files.zip
+
+sim = "tpv16"
+
+# ----------------------------------------------------------------------
+import numpy
+from spatialdata.spatialdb.SimpleIOAscii import SimpleIOAscii
+from spatialdata.geocoords.CSCart import CSCart
+
+data = numpy.loadtxt("%s_input_file.txt" % sim, skiprows=3)
+
+points = numpy.transpose((0*data[:,2], data[:,2]-24.0e+3, -data[:,3]))
+
+tractionNormal = -data[:,4]
+tractionShearHoriz = -data[:,5]
+tractionShearVert = data[:,6]
+
+coefStatic = data[:,9]
+coefDyn = data[:,10]
+d0 = data[:,11]
+cohesion = data[:,12]
+weakTime = data[:,13]
+
+
+cs = CSCart()
+cs._configure()
+cs.initialize()
+
+writer = SimpleIOAscii()
+writer.inventory.filename = "empty"
+writer._configure()
+
+writer.filename("%s_tractions.spatialdb" % sim)
+dataOut = {'points': points,
+ 'coordsys': cs,
+ 'data_dim': 2,
+ 'values': [{'name': 'traction-shear-leftlateral',
+ 'units': 'Pa',
+ 'data': tractionShearHoriz},
+ {'name': 'traction-shear-updip',
+ 'units': 'Pa',
+ 'data': tractionShearVert},
+ {'name': 'traction-normal',
+ 'units': 'Pa',
+ 'data': tractionNormal},
+ ],
+ }
+
+writer.write(dataOut)
+
+writer.filename("%s_friction.spatialdb" % sim)
+dataOut = {'points': points,
+ 'coordsys': cs,
+ 'data_dim': 2,
+ 'values': [{'name': 'static-coefficient',
+ 'units': 'None',
+ 'data': coefStatic},
+ {'name': 'dynamic-coefficient',
+ 'units': 'None',
+ 'data': coefDyn},
+ {'name': 'slip-weakening-parameter',
+ 'units': 'm',
+ 'data': d0},
+ {'name': 'cohesion',
+ 'units': 'Pa',
+ 'data': cohesion},
+ {'name': 'weakening-time',
+ 'units': 's',
+ 'data': weakTime},
+ ],
+ }
+writer.write(dataOut)
+
+
+# End of file
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/create_spatialdb.py
___________________________________________________________________
Name: svn:executable
+ *
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/createbc.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/createbc.jou 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,72 @@
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 3
+block 1 name "elastic"
+
+# ----------------------------------------------------------------------
+# Create nodeset for fault
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+nodeset 20 group fault
+nodeset 20 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 6
+nodeset 10 group face_xpos
+nodeset 10 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 4
+nodeset 11 group face_xneg
+nodeset 11 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 12
+group "face_ypos" add node in surface 14
+nodeset 12 group face_ypos
+nodeset 12 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 10
+group "face_yneg" add node in surface 17
+nodeset 13 group face_yneg
+nodeset 13 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 11
+group "face_zneg" add node in surface 15
+nodeset 15 group face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 9
+group "face_zpos" add node in surface 16
+nodeset 16 group face_zpos
+nodeset 16 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for subset of +z face
+# ----------------------------------------------------------------------
+group "face_zpos_subset" add node in group face_zpos
+group "face_zpos_subset" remove node with y_coord > +24.001e+3
+group "face_zpos_subset" remove node with y_coord < -24.001e+3
+group "face_zpos_subset" remove node with x_coord > +10.001e+3
+group "face_zpos_subset" remove node with x_coord < -10.001e+3
+nodeset 17 group face_zpos_subset
+nodeset 17 name "face zpos_subset"
+
+
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/geometry.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/geometry.jou 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,67 @@
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+# Block is 48.0 km x 72.0 km x 32.0 km
+# -24.0 km <= x <= 24.0 km
+# -32.0 km <= y <= 32.0 km
+# -32.0 km <= z <= 0.0 km
+#
+#{Units('si')}
+#
+#{blockWidth=48.0*km}
+#{blockLength=72.0*km}
+#{blockHeight=32.0*km}
+#
+#{faultLength=48.0*km}
+#{faultWidth=19.5*km}
+
+reset
+brick x {blockWidth} y {blockLength} z {blockHeight}
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z {-0.5*blockHeight}
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+#{sI=Id("surface")}
+
+# ----------------------------------------------------------------------
+# Divide volumes using interface surfaces
+# ----------------------------------------------------------------------
+webcut volume 1 with plane surface {sI}
+
+# ----------------------------------------------------------------------
+# Split fault surface
+# ----------------------------------------------------------------------
+create vertex 0 {-0.5*faultLength} 0
+#{vAn=Id("vertex")}
+create vertex 0 {+0.5*faultLength} 0
+#{vAp=Id("vertex")}
+create vertex 0 {-0.5*faultLength} {-blockHeight}
+#{vBn=Id("vertex")}
+create vertex 0 {+0.5*faultLength} {-blockHeight}
+#{vBp=Id("vertex")}
+create vertex 0 {-0.5*faultLength} {-faultWidth}
+#{vCn=Id("vertex")}
+create vertex 0 {+0.5*faultLength} {-faultWidth}
+#{vCp=Id("vertex")}
+split surface 7 across location vertex {vAn} location vertex {vBn}
+split surface 18 across location vertex {vAp} location vertex {vBp}
+split surface 21 across location vertex {vCn} location vertex {vCp}
+surface 22 name "fault_surface"
+delete vertex all
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all
+merge all
+delete body 2
+
+
+
+
+
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/gradient.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/gradient.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/gradient.jou 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,85 @@
+# ----------------------------------------------------------------------
+# Set vertex spacing with increasing spacing away from fault
+# ----------------------------------------------------------------------
+#{bias_factor=1.02}
+
+# Function for computing discretization size at end of curve with bias
+#{fun_dx='dxStart*bias_factor**ceil( ln(1-curveL/dxStart*(1-bias_factor))/ln(bias_factor))'}
+
+# ----------------------------------------------------------------------
+# Compute sizes at curve endpoints
+# ----------------------------------------------------------------------
+# dxA - size at vertex 9 12
+#{dxStart=dx}{curveL=Length(40)}{execute("dxA="//fun_dx)}
+
+# dxB - size at vertex 27 29
+#{dxStart=dx}{curveL=Length(46)}{execute("dxB="//fun_dx)}
+
+# dxC - size at vertex 10 11
+#{dxStart=dxB}{curveL=Length(39)}{execute("dxC="//fun_dx)}
+
+# dxD - size at vertex 1 2 3 4
+#{dxStart=dxA}{curveL=Length(22)}{execute("dxD="//fun_dx)}
+
+# dxE - size at vertex 5 6 7 8
+#{dxStart=dxC}{curveL=Length(24)}{execute("dxE="//fun_dx)}
+
+# ----------------------------------------------------------------------
+# Reset sizes
+# ----------------------------------------------------------------------
+curve all scheme default
+surface all sizing function none
+surface fault_surface size {dx}
+
+# Set bias on curves extending from faults
+
+curve 40 scheme bias fine size {dx} factor {bias_factor} start vertex 30
+curve 36 scheme bias fine size {dx} factor {bias_factor} start vertex 28
+
+curve 46 scheme bias fine size {dx} factor {bias_factor} start vertex 31
+curve 47 scheme bias fine size {dx} factor {bias_factor} start vertex 32
+
+# ----------------------------------------------------------------------
+# B to B
+curve 42 size {dxB}
+
+# bias starting at A
+curve 22 scheme bias fine size {dxA} factor {bias_factor} start vertex 12
+curve 30 scheme bias fine size {dxA} factor {bias_factor} start vertex 12
+
+curve 21 scheme bias fine size {dxA} factor {bias_factor} start vertex 9
+curve 32 scheme bias fine size {dxA} factor {bias_factor} start vertex 9
+
+# bias A to C
+curve 13 scheme bias fine size {dxA} coarse size {dxC} start vertex 9
+curve 15 scheme bias fine size {dxA} coarse size {dxC} start vertex 12
+
+# bias starting at B
+curve 39 scheme bias fine size {dxB} factor {bias_factor} start vertex 29
+curve 37 scheme bias fine size {dxB} factor {bias_factor} start vertex 27
+
+# bias starting at C
+curve 24 scheme bias fine size {dxC} factor {bias_factor} start vertex 11
+curve 29 scheme bias fine size {dxC} factor {bias_factor} start vertex 11
+
+curve 23 scheme bias fine size {dxC} factor {bias_factor} start vertex 10
+curve 31 scheme bias fine size {dxC} factor {bias_factor} start vertex 10
+
+# bias D to D
+curve 1 size {dxD}
+curve 3 size {dxD}
+
+# bias D to E
+curve 10 scheme bias fine size {dxD} coarse size {dxE} start vertex 1
+curve 12 scheme bias fine size {dxD} coarse size {dxE} start vertex 2
+curve 11 scheme bias fine size {dxD} coarse size {dxE} start vertex 3
+curve 9 scheme bias fine size {dxD} coarse size {dxE} start vertex 4
+
+# E to E
+curve 5 size {dxE}
+curve 7 size {dxE}
+
+# End of file
+
+
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/hex8.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/hex8.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/hex8.cfg 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,46 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.elastic]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_neg]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.z_neg]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/hex8_150m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/hex8_150m.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/hex8_150m.jou 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,56 @@
+# ----------------------------------------------------------------------
+# Create hex8 mesh at 200m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size and scheme
+# ----------------------------------------------------------------------
+#{dx=148.14666667.0*m}
+volume all scheme map
+volume all size {dx*3}
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine around fault
+# ----------------------------------------------------------------------
+#{refineWidth=18.001*km}
+#{refineLength=54.001*km}
+#{refineDepth=21.001*km}
+group "refinecells" add Hex with x_min >= {-0.5*refineWidth} and x_max <= {+0.5*refineWidth} and y_min >= {-0.5*refineLength} and y_max <= {+0.5*refineLength} and z_min >= {-refineDepth}
+refine hex in group refinecells depth 1
+delete group refinecells
+
+# ----------------------------------------------------------------------
+# Smooth to improve quality
+# ----------------------------------------------------------------------
+#{condnum=3.0}
+#{loop(5)}
+volume all smooth scheme condition number beta {condnum} cpu 20
+smooth volume all
+#{condnum=condnum-0.2}
+#{endloop}
+
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "hex8_150m.exo" dimension 3 overwrite
+
+
+
+
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/matprops.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/matprops.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/matprops.spatialdb 2012-01-20 16:04:37 UTC (rev 19397)
@@ -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 = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 2670.0 3464.0 6000.0
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/pylithapp.cfg 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,165 @@
+# -*- Python -*-
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+explicitlumped = 1
+petsc = 1
+meshiocubit = 1
+elasticityexplicit = 1
+elasticityexplicittet4 = 1
+meshimporter = 1
+#quadrature2d = 1
+#fiatlagrange = 1
+faultcohesivedyn = 1
+
+[pylithapp.journal.debug]
+#pylithapp = 1
+#problem = 1
+#explicit = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+reader = pylith.meshio.MeshIOCubit
+reorder_mesh = True
+
+distributor.partitioner = parmetis
+distributor.write_partition = True
+distributor.data_writer = pylith.meshio.DataWriterHDF5ExtMesh
+distributor.data_writer.filename = output/partition.h5
+
+[pylithapp.mesh_generator.reader]
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 3
+
+formulation = pylith.problems.ExplicitLumped
+formulation.norm_viscosity = 0.3
+normalizer = spatialdata.units.NondimElasticDynamic
+normalizer.shear_wave_speed = 3333*m/s
+
+[pylithapp.timedependent.formulation.time_step]
+#total_time = 15.001*s
+total_time = 0.0*s
+dt = 0.01*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [elastic]
+materials.elastic = pylith.materials.ElasticIsotropic3D
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+db_properties.label = Elastic properties
+db_properties.iohandler.filename = matprops.spatialdb
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+bc = [x_neg,x_pos,y_neg,y_pos,z_neg]
+bc.x_neg = pylith.bc.AbsorbingDampers
+bc.x_pos = pylith.bc.AbsorbingDampers
+bc.y_neg = pylith.bc.AbsorbingDampers
+bc.y_pos = pylith.bc.AbsorbingDampers
+bc.z_neg = pylith.bc.AbsorbingDampers
+
+[pylithapp.timedependent.bc.x_pos]
+label = face_xpos
+db.label = Absorbing BC +x
+db.iohandler.filename = matprops.spatialdb
+
+[pylithapp.timedependent.bc.x_neg]
+label = face_xneg
+db.label = Absorbing BC -x
+db.iohandler.filename = matprops.spatialdb
+
+[pylithapp.timedependent.bc.y_pos]
+label = face_ypos
+db.label = Absorbing BC +y
+db.iohandler.filename = matprops.spatialdb
+
+[pylithapp.timedependent.bc.y_neg]
+label = face_yneg
+db.label = Absorbing BC -y
+db.iohandler.filename = matprops.spatialdb
+
+[pylithapp.timedependent.bc.z_neg]
+label = face_zneg
+up_dir = [0, 1, 0]
+db.label = Absorbing BC -z
+db.iohandler.filename = matprops.spatialdb
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+interfaces = [fault]
+interfaces.fault = pylith.faults.FaultCohesiveDyn
+
+[pylithapp.timedependent.interfaces.fault]
+id = 10
+label = fault
+
+db_initial_tractions = spatialdata.spatialdb.SimpleDB
+db_initial_tractions.label = Initial fault tractions
+
+friction = pylith.friction.SlipWeakeningTime
+friction.label = Slip weakening with forced time weakening
+
+friction.db_properties = spatialdata.spatialdb.SimpleDB
+friction.db_properties.label = Slip weakening
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation]
+output = [domain,subdomain]
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+[pylithapp.problem.formulation.output.domain]
+vertex_data_fields=[displacement,velocity]
+output_freq = time_step
+time_step = 12.0*s
+writer = pylith.meshio.DataWriterHDF5ExtMesh
+
+[pylithapp.problem.formulation.output.subdomain]
+label = face_zpos
+vertex_data_fields=[displacement,velocity]
+output_freq = time_step
+time_step = 0.04999*s
+writer = pylith.meshio.DataWriterHDF5ExtSubMesh
+
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_info_fields = [strike_dir,normal_dir,initial_traction,static_coefficient,dynamic_coefficient,slip_weakening_parameter,cohesion,weakening_time]
+vertex_data_fields = [slip,slip_rate,traction]
+output_freq = time_step
+time_step = 0.04999*s
+writer = pylith.meshio.DataWriterHDF5ExtSubSubMesh
+
+[pylithapp.timedependent.materials.elastic.output]
+cell_data_fields = []
+cell_filter = pylith.meshio.CellFilterAvgMesh
+output_freq = time_step
+time_step = 12.0*s
+writer = pylith.meshio.DataWriterHDF5ExtMesh
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+log_summary = true
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tet4.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tet4.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tet4.cfg 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,54 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+formulation = pylith.problems.ExplicitLumpedTet4
+
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.elastic]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_neg]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.z_neg]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tet4_150m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tet4_150m.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tet4_150m.jou 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,46 @@
+# ----------------------------------------------------------------------
+# Create tet4 mesh at 100m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size and scheme
+# ----------------------------------------------------------------------
+#{dx=150*m}
+surface all scheme trimesh
+volume all scheme tetmesh
+playback 'gradient.jou'
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh surface all
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
+
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+#{condnum=2.0}
+#{loop(4)}
+cleanup volume all
+volume all smooth scheme condition number beta {condnum} cpu 10
+smooth volume all
+#{condnum=condnum-0.1}
+#{endloop}
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "tet4_150m.exo" dimension 3 overwrite
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16.cfg 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,13 @@
+# -*- Python -*-
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+db_initial_tractions.iohandler.filename = tpv16_tractions.spatialdb
+friction.db_properties.iohandler.filename = tpv16_friction.spatialdb
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_friction.spatialdb.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_friction.spatialdb.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_hex8_075m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_hex8_075m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_hex8_075m.cfg 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,30 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8_100m.exo
+refiner = pylith.topology.RefineUniform
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8_100m.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_100m-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/hex8_100m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8_100m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_hex8_150m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_hex8_150m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_hex8_150m.cfg 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,30 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8_200m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01666666666666666666666667*s
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8_200m.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_200m-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/hex8_200m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8_200m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_tet4_075m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_tet4_075m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_tet4_075m.cfg 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,33 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+reader.filename = tet4_150m.exo
+refiner = pylith.topology.RefineUniform
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.normalizer]
+wave_period = 0.3*s
+
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.005*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tpv16_tet4_075m.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tpv16_tet4_075m-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv16_tet4_075m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv16_tet4_075m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_tet4_150m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_tet4_150m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_tet4_150m.cfg 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,32 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_150m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.normalizer]
+wave_period = 0.6*s
+
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.0125*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tpv16_tet4_150m.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tpv16_tet4_150m-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv16_tet4_150m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv16_tet4_150m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_tractions.spatialdb.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv16_tractions.spatialdb.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv17.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv17.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv16/tpv17.cfg 2012-01-20 16:04:37 UTC (rev 19397)
@@ -0,0 +1,13 @@
+# -*- Python -*-
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+db_initial_tractions.iohandler.filename = tpv17_tractions.spatialdb
+friction.db_properties.iohandler.filename = tpv17_friction.spatialdb
+
+
+# End of file
More information about the CIG-COMMITS
mailing list