[cig-commits] r19829 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210
brad at geodynamics.org
brad at geodynamics.org
Tue Mar 20 16:04:39 PDT 2012
Author: brad
Date: 2012-03-20 16:04:38 -0700 (Tue, 20 Mar 2012)
New Revision: 19829
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/gradient.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_friction.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tet4_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tractions.spatialdb
Removed:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/friction.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tractions.spatialdb
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/geometry.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4_200m.jou
Log:
Updated mesh generation for tet4.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/createbc.jou 2012-03-20 16:41:07 UTC (rev 19828)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/createbc.jou 2012-03-20 23:04:38 UTC (rev 19829)
@@ -1,68 +1,57 @@
# ----------------------------------------------------------------------
# Create blocks for materials
# ----------------------------------------------------------------------
-block 1 volume 1 3
+block 1 volume 1 4 9 10
block 1 name "elastic"
# ----------------------------------------------------------------------
# Create nodeset for fault
# ----------------------------------------------------------------------
-group "fault" add node in fault_surface
-group "fault" remove node with x_coord < -7.5001e+3
-group "fault" remove node with y_coord < -15.0001e+3
-group "fault" remove node with y_coord > +15.0001e+3
+group "fault" add node in fault_surface fault_surface at A fault_surface at B
nodeset 10 group fault
nodeset 10 name "fault"
# ----------------------------------------------------------------------
# Create nodeset for +x face
# ----------------------------------------------------------------------
-group "face_xpos" add node in surface 6
-nodeset 11 group face_xpos
-nodeset 11 name "face_xpos"
+group "face_xpos" add node in surface 35 47
+nodeset 20 group face_xpos
+nodeset 20 name "face_xpos"
# ----------------------------------------------------------------------
# Create nodeset for -x face
# ----------------------------------------------------------------------
-group "face_xneg" add node in surface 4
-nodeset 12 group face_xneg
-nodeset 12 name "face_xneg"
+group "face_xneg" add node in surface 49 57
+nodeset 22 group face_xneg
+nodeset 22 name "face_xneg"
# ----------------------------------------------------------------------
# Create nodeset for +y face
# ----------------------------------------------------------------------
-group "face_ypos" add node in surface 10
-group "face_ypos" add node in surface 17
-nodeset 13 group face_ypos
-nodeset 13 name "face_ypos"
+group "face_ypos" add node in surface 40 42 50 56
+nodeset 23 group face_ypos
+nodeset 23 name "face_ypos"
# ----------------------------------------------------------------------
# Create nodeset for -y face
# ----------------------------------------------------------------------
-group "face_yneg" add node in surface 12
-group "face_yneg" add node in surface 14
-nodeset 14 group face_yneg
-nodeset 14 name "face_yneg"
+group "face_yneg" add node in surface 36 46 52 54
+nodeset 24group face_yneg
+nodeset 24name "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"
+group "face_zneg" add node in surface 12 16
+nodeset 25 group face_zneg
+nodeset 25 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"
+group "face_zpos" add node in surface 10 17
+nodeset 26 group face_zpos
+nodeset 26 name "face_zpos"
-# ----------------------------------------------------------------------
-# Create nodeset for off-fault stations
-# ----------------------------------------------------------------------
-group "off_fault" add node with x_min >= -4.001e+3 and x_max <= 4.001e+3 and z_min >= -1.001e+3 and y_coord = 0.0
-nodeset 25 group off_fault
-nodeset 25 name "off_fault"
\ No newline at end of file
+
+# End of file
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/friction.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/friction.spatialdb 2012-03-20 16:41:07 UTC (rev 19828)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/friction.spatialdb 2012-03-20 23:04:38 UTC (rev 19829)
@@ -1,22 +0,0 @@
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 4
- value-names = static-coefficient dynamic-coefficient slip-weakening-parameter cohesion
- value-units = none none m MPa
- num-locs = 8
- data-dim = 2
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0e+3 // x, y, z are in km
- space-dim = 3
- }
-}
--7.4999 -15.001 0.0 1000.0 0.448 0.5 0.2
--7.4999 -15.0 0.0 0.760 0.448 0.5 0.2
--7.4999 +15.0 0.0 0.760 0.448 0.5 0.2
--7.4999 +15.001 0.0 1000.0 0.448 0.5 0.2
-
--7.5001 -15.001 0.0 10000.0 0.448 0.5 1000.0
--7.5001 -15.0 0.0 10000.0 0.448 0.5 1000.0
--7.5001 +15.0 0.0 10000.0 0.448 0.5 1000.0
--7.5001 +15.001 0.0 10000.0 0.448 0.5 1000.0
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/geometry.jou 2012-03-20 16:41:07 UTC (rev 19828)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/geometry.jou 2012-03-20 23:04:38 UTC (rev 19829)
@@ -1,31 +1,103 @@
+# -*- Python -*- (syntax highlighting)
# ----------------------------------------------------------------------
-# Create block
+#
+# Brad T. Aagaard, U.S. Geological Survey
+#
# ----------------------------------------------------------------------
-# Block is 36.0 km x 56000 km x 18.0*sqrt(3) km
-# -21.75 km <= x <= 14.25 km
-# -28 km <= y <= 28km
-# -18*Sqrt(3) km <= z <= 0 km
#
+# CUBIT journal file with geometry for TPV210.
+#
+# ----------------------------------------------------------------------
+# Set units to SI.
+# ----------------------------------------------------------------------
+#{Units('si')}
+#
+# ----------------------------------------------------------------------
+# Reset geometry.
+# ----------------------------------------------------------------------
reset
-brick x 36000 y 56000 z 31176.914
+# Make sure undo is off to prevent errors in stitching volumes.
+undo off
+
# ----------------------------------------------------------------------
-# Set z=0 at free surface
+# Create block
# ----------------------------------------------------------------------
-move volume 1 x -3750 z -15588.457
+# Block is 48.0 km x 64.0 km x 36.0 km
+# Not including offset (top of fault is at x=0)
+# -24 km <= x <= 24 km
+# -32 km <= y <= 32 km
+# -36 km <= z <= 0 km
+#
+# Fault dip is 60 degrees and width is 15.0 km.
+#
+#{blockWidth=48.0*km}
+#{blockLength=64.0*km}
+#{blockHeight=36.0*km}
+#{faultDipAngle=60.0}
+#{faultWidth=15.0*km}
+#{faultLength=30.0*km}
+#
+#{xoffset=0.5*faultWidth*cosd(faultDipAngle)}
+brick x {blockWidth} y {blockLength} z {blockHeight}
+
+# Translate block so the top is at z=0
+volume 1 move x {-xoffset} y 0 z {-0.5*blockHeight}
+
# ----------------------------------------------------------------------
-# Create interface surfaces
+# Create fault surface
# ----------------------------------------------------------------------
-create surface extended from surface 1
-body 2 rotate -60 about y
-surface 7 name "fault_surface"
-webcut volume 1 with plane surface fault_surface
+create planar surface with plane xplane offset 0
+#{fI=Id("surface")}
+rotate surface {fI} angle {faultDipAngle-90} about vertex 12 9 include_merged
+create surface extended from surface {fI}
+delete surface {fI}
+#{fI=Id("surface")}
+webcut volume 1 with plane surface {fI}
+surface 9 name "fault_body"
+# Split fault surface
+create vertex 0.0 {-0.5*faultLength} 0
+#{vAn=Id("vertex")}
+create vertex 0.0 {+0.5*faultLength} 0
+#{vAp=Id("vertex")}
+create vertex {-blockHeight/tand(faultDipAngle)} {-0.5*faultLength} {-blockHeight}
+#{vBn=Id("vertex")}
+create vertex {-blockHeight/tand(faultDipAngle)} {+0.5*faultLength} {-blockHeight}
+#{vBp=Id("vertex")}
+create vertex {-faultWidth*cosd(faultDipAngle)} {-0.5*faultLength} {-faultWidth*sind(faultDipAngle)}
+#{vCn=Id("vertex")}
+create vertex {-faultWidth*cosd(faultDipAngle)} {+0.5*faultLength} {-faultWidth*sind(faultDipAngle)}
+#{vCp=Id("vertex")}
+split surface fault_body across location vertex {vAn} location vertex {vBn}
+split surface fault_body at A across location vertex {vAp} location vertex {vBp}
+split surface fault_body at A across location vertex {vCn} location vertex {vCp}
+surface fault_body at C rename "fault_surface"
+delete vertex all
+
+# asperity (to ensure exact alignment with mesh)
+brick x {3.0*km} y {3.0*km} z {3.0*km}
+{idChopper=Id("volume")}
+rotate volume {idChopper} angle {-faultDipAngle} about y
+volume {idChopper} move x {-12.0*km*cosd(faultDipAngle)} y 0 z {-12.0*km*sind(faultDipAngle)}
+chop volume 3 with volume {idChopper} keep
+delete volume {idChopper}
+
# ----------------------------------------------------------------------
+# Create interface surface to honor stress discontinuity in TPV13
+# ----------------------------------------------------------------------
+create planar surface with plane zplane offset {sind(faultDipAngle)*-13.8*km}
+#{sI=Id("surface")}
+webcut volume 1 4 with plane surface {sI}
+
+# ----------------------------------------------------------------------
# Imprint all volumes, merging surfaces
# ----------------------------------------------------------------------
imprint all
merge all
+delete body 3 6 7 8
-delete body 2
+# End of file
+
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/gradient.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/gradient.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/gradient.jou 2012-03-20 23:04:38 UTC (rev 19829)
@@ -0,0 +1,148 @@
+# -*- Python -*- (syntax highlighting)
+# ----------------------------------------------------------------------
+# 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 32 34
+#{dxStart=dx}{curveL=Length(39)}{execute("dxA="//fun_dx)}
+
+# dxB - size at vertex 63 66
+#{dxStart=dx}{curveL=Length(85)}{execute("dxB="//fun_dx)}
+
+# dxC - size at vertex 31 33
+#{dxStart=dx}{curveL=Length(49)}{execute("dxC="//fun_dx)}
+
+# dxD - size at vertex 17 18
+#{dxStart=dxC}{curveL=Length(38)}{execute("dxD="//fun_dx)}
+
+# dxE - size at vertex 1 2
+#{dxStart=dxA}{curveL=Length(25)}{execute("dxE="//fun_dx)}
+
+# dxF - size at vertex 3 4
+#{dxStart=dxA}{curveL=Length(36)}{execute("dxF="//fun_dx)}
+
+# dxG - size at vertex 64 65
+#{dxStart=dxB}{curveL=Length(84)}{execute("dxG="//fun_dx)}
+
+# dxH - size at vertex 75 76
+#{dxStart=dxB}{curveL=Length(105)}{execute("dxH="//fun_dx)}
+
+# dxI - size at vertex 5 6
+#{dxStart=dxG}{curveL=Length(99)}{execute("dxI="//fun_dx)}
+
+# dxJ - size at vertex 7 8
+#{dxStart=dxD}{curveL=Length(35)}{execute("dxJ="//fun_dx)}
+
+# ----------------------------------------------------------------------
+# Reset sizes
+# ----------------------------------------------------------------------
+curve all scheme default
+surface all sizing function none
+surface fault_surface fault_surface at A fault_surface at B size {dx}
+
+# Set bias on curves extending from faults
+
+curve 39 scheme bias fine size {dx} factor {bias_factor} start vertex 32
+curve 45 scheme bias fine size {dx} factor {bias_factor} start vertex 34
+
+curve 85 scheme bias fine size {dx} factor {bias_factor} start vertex 62
+curve 81 scheme bias fine size {dx} factor {bias_factor} start vertex 61
+
+curve 49 scheme bias fine size {dx} factor {bias_factor} start vertex 36
+curve 48 scheme bias fine size {dx} factor {bias_factor} start vertex 35
+
+# ----------------------------------------------------------------------
+# bias starting at A
+curve 25 scheme bias fine size {dxA} factor {bias_factor} start vertex 19
+curve 36 scheme bias fine size {dxA} factor {bias_factor} start vertex 19
+
+curve 26 scheme bias fine size {dxA} factor {bias_factor} start vertex 20
+curve 34 scheme bias fine size {dxA} factor {bias_factor} start vertex 20
+
+# bias starting at B
+curve 84 scheme bias fine size {dxB} factor {bias_factor} start vertex 66
+curve 105 scheme bias fine size {dxB} factor {bias_factor} start vertex 66
+
+curve 82 scheme bias fine size {dxB} factor {bias_factor} start vertex 63
+curve 107 scheme bias fine size {dxB} factor {bias_factor} start vertex 63
+
+# ----------------------------------------------------------------------
+# C to C
+curve 43 size {dxC}
+
+# bias starting at C
+curve 38 scheme bias fine size {dxC} factor {bias_factor} start vertex 31
+curve 46 scheme bias fine size {dxC} factor {bias_factor} start vertex 33
+
+# bias starting at D
+curve 35 scheme bias fine size {dxD} factor {bias_factor} start vertex 17
+curve 33 scheme bias fine size {dxD} factor {bias_factor} start vertex 18
+
+# bias E to E
+curve 1 size {dxE}
+
+# bias F to F
+curve 3 size {dxF}
+
+# bias G to G
+curve 83 size {dxG}
+
+# bias H to H
+curve 106 size {dxH}
+
+# bias I to I
+curve 5 size {dxI}
+
+# bias J to J
+curve 7 size {dxJ}
+
+
+# bias A to B
+curve 91 scheme bias fine size {dxA} coarse size {dxB} start vertex 19
+curve 88 scheme bias fine size {dxA} coarse size {dxB} start vertex 20
+
+# bias B to D
+curve 98 scheme bias fine size {dxB} coarse size {dxD} start vertex 66
+curve 102 scheme bias fine size {dxB} coarse size {dxD} start vertex 63
+
+# bias D to I
+curve 27 scheme bias fine size {dxD} coarse size {dxI} start vertex 17
+curve 28 scheme bias fine size {dxD} coarse size {dxI} start vertex 18
+
+# bias E to G
+curve 87 scheme bias fine size {dxE} coarse size {dxG} start vertex 1
+curve 86 scheme bias fine size {dxE} coarse size {dxG} start vertex 2
+
+# bias F to H
+curve 108 scheme bias fine size {dxF} coarse size {dxH} start vertex 4
+curve 109 scheme bias fine size {dxF} coarse size {dxH} start vertex 3
+
+# bias G to I
+curve 99 scheme bias fine size {dxG} coarse size {dxI} start vertex 65
+curve 103 scheme bias fine size {dxG} coarse size {dxI} start vertex 64
+
+# bias H to J
+curve 119 scheme bias fine size {dxD} coarse size {dxI} start vertex 75
+curve 117 scheme bias fine size {dxD} coarse size {dxI} start vertex 76
+
+# Set bias on surfaces
+surface fault_body sizing function type bias start curve 90 factor {bias_factor}
+surface fault_body at B sizing function type bias start curve 89 factor {bias_factor}
+surface fault_body at A sizing function type bias start curve 47 factor {bias_factor}
+surface fault_body at C sizing function type bias start curve 100 factor {bias_factor}
+surface fault_body at D sizing function type bias start curve 101 factor {bias_factor}
+surface 34 sizing function type bias start curve 80 factor {bias_factor}
+surface 48 sizing function type bias start curve 80 factor {bias_factor}
+
+# End of file
+
+
+
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg 2012-03-20 16:41:07 UTC (rev 19828)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg 2012-03-20 23:04:38 UTC (rev 19829)
@@ -1,4 +1,3 @@
-# -*- Python -*-
[pylithapp]
# ----------------------------------------------------------------------
@@ -6,17 +5,21 @@
# ----------------------------------------------------------------------
[pylithapp.journal.info]
timedependent = 1
-explicitlumped = 1
petsc = 1
+meshimporter = 1
meshiocubit = 1
+explicitlumped = 1
+explicitlumpedtet4 = 1
elasticityexplicit = 1
-meshimporter = 1
+elasticityexplicittet4 = 1
faultcohesivedyn = 1
+fiatsimplex = 1
[pylithapp.journal.debug]
#pylithapp = 1
#problem = 1
-#explicit = 1
+#explicitlumped = 1
+#explicitlumpedtet4 = 1
# ----------------------------------------------------------------------
# mesh_generator
@@ -25,6 +28,11 @@
reader = pylith.meshio.MeshIOCubit
reorder_mesh = True
+distributor.partitioner = parmetis
+distributor.write_partition = True
+distributor.data_writer = pylith.meshio.DataWriterHDF5Mesh
+distributor.data_writer.filename = output/partition.h5
+
[pylithapp.mesh_generator.reader]
coordsys.space_dim = 3
@@ -32,59 +40,44 @@
# problem
# ----------------------------------------------------------------------
[pylithapp.timedependent]
-
dimension = 3
-# Set interfaces to an array with 1 faults: 'fault'.
-interfaces = [fault]
-
-# Set materials to an array with 1 material 'elastic'.
-materials = [elastic]
-
formulation = pylith.problems.ExplicitLumped
-
-# Nondimensionalize problem using wave propagation parameters.
+formulation.norm_viscosity = 0.1
normalizer = spatialdata.units.NondimElasticDynamic
-#normalizer.mass_density = 1000.0*kg/m**3
-#normalizer.shear_wave_speed = 1.0*m/s
-#normalizer.wave_period = 1.0*s
+normalizer.shear_wave_speed = 3333*m/s
+normalizer.mass_density = 2700.0*kg/m**3
-bc = [x_neg,x_pos,z_neg,y_neg,y_pos]
-
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_neg = pylith.bc.AbsorbingDampers
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.z_neg = pylith.bc.AbsorbingDampers
-bc.y_neg = pylith.bc.AbsorbingDampers
-bc.y_pos = pylith.bc.AbsorbingDampers
-
[pylithapp.timedependent.formulation.time_step]
-total_time = 15.01*s
+#total_time = 15.01*s
+total_time = 0.01*s
dt = 0.01*s
-[pylithapp.timedependent.formulation]
-output = [groundsurf,offfault]
-output.groundsurf = pylith.meshio.OutputSolnSubset
-output.offfault = pylith.meshio.OutputSolnSubset
-
-
# ----------------------------------------------------------------------
# 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.z_neg = pylith.bc.AbsorbingDampers
+bc.y_neg = pylith.bc.AbsorbingDampers
+bc.y_pos = pylith.bc.AbsorbingDampers
+
[pylithapp.timedependent.bc.x_pos]
label = face_xpos
db.label = Absorbing BC +x
@@ -114,8 +107,10 @@
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+interfaces = [fault]
+interfaces.fault = pylith.faults.FaultCohesiveDyn
-# Provide information on the fault (interface).
[pylithapp.timedependent.interfaces]
fault = pylith.faults.FaultCohesiveDyn
@@ -124,55 +119,48 @@
id = 10
label = fault
-db_initial_tractions = spatialdata.spatialdb.SimpleDB
-db_initial_tractions.iohandler.filename = tractions.spatialdb
-db_initial_tractions.label = "Initial fault tractions"
-db_initial_tractions.query_type = linear
-
friction = pylith.friction.SlipWeakening
+friction.label = Slip weakening
friction.db_properties = spatialdata.spatialdb.SimpleDB
-friction.db_properties.iohandler.filename = friction.spatialdb
friction.db_properties.label = Slip weakening
-
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
# Ground surface
-[pylithapp.problem.formulation.output.groundsurf]
+[pylithapp.timedependent.formulation.output]
+output = pylith.meshio.OutputSolnSubset
+
+[pylithapp.problem.formulation.output.output]
label = face_zpos
vertex_data_fields=[displacement,velocity]
output_freq = time_step
time_step = 0.04999*s
-writer.time_format = %06.3f
+writer = pylith.meshio.DataWriterHDF5SubMesh
-# Off-fault locations not on ground surface
-[pylithapp.problem.formulation.output.offfault]
-label = off_fault
-vertex_data_fields=[displacement,velocity]
-output_freq = time_step
-time_step = 0.04999*s
-writer.time_format = %06.3f
-
-# fault
+# Fault
[pylithapp.timedependent.interfaces.fault.output]
vertex_info_fields = [strike_dir,normal_dir,initial_traction]
vertex_data_fields = [slip,slip_rate,traction]
output_freq = time_step
time_step = 0.04999*s
-writer.time_format = %06.3f
+writer = pylith.meshio.DataWriterHDF5SubSubMesh
# State variables.
[pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
-output_freq = time_step
-writer.time_format = %06.3f
cell_info_fields = []
cell_data_fields = []
+cell_filter = pylith.meshio.CellFilterAvgMesh
+output_freq = time_step
+time_step = 15.0*s
+writer = pylith.meshio.DataWriterHDF5Mesh
# ----------------------------------------------------------------------
# PETSc
# ----------------------------------------------------------------------
[pylithapp.petsc]
log_summary = true
+
+
+# End of file
\ No newline at end of file
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4.cfg 2012-03-20 16:41:07 UTC (rev 19828)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4.cfg 2012-03-20 23:04:38 UTC (rev 19829)
@@ -1,4 +1,3 @@
-# -*- Python -*-
[pylithapp]
# ----------------------------------------------------------------------
@@ -58,3 +57,5 @@
cell.degree = 1
cell.quad_order = 1
+
+# End of file
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4_200m.cfg 2012-03-20 16:41:07 UTC (rev 19828)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4_200m.cfg 2012-03-20 23:04:38 UTC (rev 19829)
@@ -1,33 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator.reader]
-filename = tet4_200m.exo
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.formulation.time_step]
-dt = 0.01666666666666666666666667*s
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Ground surface
-[pylithapp.problem.formulation.output.groundsurf]
-writer.filename = output/tet4_200m-groundsurf.vtk
-
-# Off-fault
-[pylithapp.problem.formulation.output.offfault]
-writer.filename = output/tet4_200m-offfault.vtk
-
-# Fault
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/tet4_200m-fault.vtk
-
-# State variables
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/tet4_200m-elastic.vtk
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4_200m.jou 2012-03-20 16:41:07 UTC (rev 19828)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4_200m.jou 2012-03-20 23:04:38 UTC (rev 19829)
@@ -1,4 +1,9 @@
+# -*- Python -*- (syntax highlighting)
# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+#
+# ----------------------------------------------------------------------
# Create tet4 mesh at 200m resolution.
# ----------------------------------------------------------------------
@@ -10,18 +15,30 @@
# ----------------------------------------------------------------------
# Set discretization size and scheme
# ----------------------------------------------------------------------
-volume all size 200
+#{dx=200*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
# ----------------------------------------------------------------------
-# Improve element quality.
+# Smooth mesh to improve quality
# ----------------------------------------------------------------------
-playback 'tet4_smooth.jou'
+#{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.
@@ -34,6 +51,5 @@
export mesh "tet4_200m.exo" dimension 3 overwrite
+# End of file
-
-
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10.cfg 2012-03-20 23:04:38 UTC (rev 19829)
@@ -0,0 +1,23 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+db_initial_tractions = spatialdata.spatialdb.SimpleDB
+db_initial_tractions.label = Initial fault tractions
+db_initial_tractions.iohandler.filename = tpv10_tractions.spatialdb
+db_initial_tractions.query_type = linear
+
+friction.db_properties.iohandler.filename = tpv10_friction.spatialdb
+friction.db_properties.query_type = nearest
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_info_fields = [strike_dir,dip_dir,normal_dir,initial_traction]
+
+
+# End of file
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_friction.spatialdb (from rev 17829, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/friction.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_friction.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_friction.spatialdb 2012-03-20 23:04:38 UTC (rev 19829)
@@ -0,0 +1,22 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = static-coefficient dynamic-coefficient slip-weakening-parameter cohesion
+ value-units = none none m MPa
+ num-locs = 8
+ data-dim = 2
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y, z are in km
+ space-dim = 3
+ }
+}
+-7.499 -15.000 0.0 10000.0 0.448 0.5 0.2
+-7.499 -14.999 0.0 0.760 0.448 0.5 0.2
+-7.499 +14.999 0.0 0.760 0.448 0.5 0.2
+-7.499 +15.000 0.0 10000.0 0.448 0.5 0.2
+
+-7.500 -15.000 0.0 10000.0 0.448 0.5 0.2
+-7.500 -14.999 0.0 10000.0 0.448 0.5 0.2
+-7.500 +14.999 0.0 10000.0 0.448 0.5 0.2
+-7.500 +15.000 0.0 10000.0 0.448 0.5 0.2
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tet4_200m.cfg (from rev 17829, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tet4_200m.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tet4_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tet4_200m.cfg 2012-03-20 23:04:38 UTC (rev 19829)
@@ -0,0 +1,28 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_200m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.normalizer]
+wave_period = 0.6*s
+
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01666666666666666666666667*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv10_tet4_200m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv10_tet4_200m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv10_tet4_200m-elastic.h5
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tractions.spatialdb (from rev 17829, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tractions.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tractions.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tractions.spatialdb 2012-03-20 23:04:38 UTC (rev 19829)
@@ -0,0 +1,129 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = traction-shear-leftlateral traction-shear-updip traction-normal
+ value-units = MPa MPa MPa
+ num-locs = 108
+ data-dim = 2
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y, z are in km
+ space-dim = 3
+ }
+}
++0.000 -15.000 0.0 0.0 0.0 0.0
++0.000 -6.000 0.0 0.0 0.0 0.0
++0.000 -3.000 0.0 0.0 0.0 0.0
++0.000 -1.501 0.0 0.0 0.0 0.0
++0.000 -1.500 0.0 0.0 0.0 0.0
++0.000 -1.499 0.0 0.0 0.0 0.0
++0.000 +1.499 0.0 0.0 0.0 0.0
++0.000 +1.500 0.0 0.0 0.0 0.0
++0.000 +1.501 0.0 0.0 0.0 0.0
++0.000 +3.000 0.0 0.0 0.0 0.0
++0.000 +6.000 0.0 0.0 0.0 0.0
++0.000 +15.000 0.0 0.0 0.0 0.0
+
+-3.000 -15.000 0.0 0.0 24.34740000 -44.2680000
+-3.000 -6.000 0.0 0.0 24.34740000 -44.2680000
+-3.000 -3.000 0.0 0.0 24.34740000 -44.2680000
+-3.000 -1.501 0.0 0.0 24.34740000 -44.2680000
+-3.000 -1.500 0.0 0.0 24.34740000 -44.2680000
+-3.000 -1.499 0.0 0.0 24.34740000 -44.2680000
+-3.000 +1.499 0.0 0.0 24.34740000 -44.2680000
+-3.000 +1.500 0.0 0.0 24.34740000 -44.2680000
+-3.000 +1.501 0.0 0.0 24.34740000 -44.2680000
+-3.000 +3.000 0.0 0.0 24.34740000 -44.2680000
+-3.000 +6.000 0.0 0.0 24.34740000 -44.2680000
+-3.000 +15.000 0.0 0.0 24.34740000 -44.2680000
+
+-5.249 -15.000 0.0 0.0 42.59983420 -77.4542440
+-5.249 -6.000 0.0 0.0 42.59983420 -77.4542440
+-5.249 -3.000 0.0 0.0 42.59983420 -77.4542440
+-5.249 -1.501 0.0 0.0 42.59983420 -77.4542440
+-5.249 -1.500 0.0 0.0 42.59983420 -77.4542440
+-5.249 -1.499 0.0 0.0 42.59983420 -77.4542440
+-5.249 +1.499 0.0 0.0 42.59983420 -77.4542440
+-5.249 +1.500 0.0 0.0 42.59983420 -77.4542440
+-5.249 +1.501 0.0 0.0 42.59983420 -77.4542440
+-5.249 +3.000 0.0 0.0 42.59983420 -77.4542440
+-5.249 +6.000 0.0 0.0 42.59983420 -77.4542440
+-5.249 +15.000 0.0 0.0 42.59983420 -77.4542440
+
+-5.250 -15.000 0.0 0.0 42.60795000 -77.4690000
+-5.250 -6.000 0.0 0.0 42.60795000 -77.4690000
+-5.250 -3.000 0.0 0.0 42.60795000 -77.4690000
+-5.250 -1.501 0.0 0.0 42.60795000 -77.4690000
+-5.250 -1.500 0.0 0.0 51.32321250 -77.4690000
+-5.250 -1.499 0.0 0.0 60.03847500 -77.4690000
+-5.250 +1.499 0.0 0.0 60.03847500 -77.4690000
+-5.250 +1.500 0.0 0.0 51.32321250 -77.4690000
+-5.250 +1.501 0.0 0.0 42.60795000 -77.4690000
+-5.250 +3.000 0.0 0.0 42.60795000 -77.4690000
+-5.250 +6.000 0.0 0.0 42.60795000 -77.4690000
+-5.250 +15.000 0.0 0.0 42.60795000 -77.4690000
+
+-5.251 -15.000 0.0 0.0 42.61606580 -77.4837560
+-5.251 -6.000 0.0 0.0 42.61606580 -77.4837560
+-5.251 -3.000 0.0 0.0 42.61606580 -77.4837560
+-5.251 -1.501 0.0 0.0 42.61606580 -77.4837560
+-5.251 -1.500 0.0 0.0 51.07268888 -77.4837560
+-5.251 -1.499 0.0 0.0 59.52931197 -77.4837560
+-5.251 +1.499 0.0 0.0 59.52931197 -77.4837560
+-5.251 +1.500 0.0 0.0 51.07268888 -77.4837560
+-5.251 +1.501 0.0 0.0 42.61606580 -77.4837560
+-5.251 +3.000 0.0 0.0 42.61606580 -77.4837560
+-5.251 +6.000 0.0 0.0 42.61606580 -77.4837560
+-5.251 +15.000 0.0 0.0 42.61606580 -77.4837560
+
+-6.749 -15.000 0.0 0.0 54.77353420 -99.5882440
+-6.749 -6.000 0.0 0.0 54.77353420 -99.5882440
+-6.749 -3.000 0.0 0.0 54.77353420 -99.5882440
+-6.749 -1.501 0.0 0.0 54.77353420 -99.5882440
+-6.749 -1.500 0.0 0.0 65.61412632 -99.5882440
+-6.749 -1.499 0.0 0.0 76.45471840 -99.5882440
+-6.749 +1.499 0.0 0.0 76.45471840 -99.5882440
+-6.749 +1.500 0.0 0.0 65.61412632 -99.5882440
+-6.749 +1.501 0.0 0.0 54.77353420 -99.5882440
+-6.749 +3.000 0.0 0.0 54.77353420 -99.5882440
+-6.749 +6.000 0.0 0.0 54.77353420 -99.5882440
+-6.749 +15.000 0.0 0.0 54.77353420 -99.5882440
+
+-6.750 -15.000 0.0 0.0 54.78165000 -99.6030000
+-6.750 -6.000 0.0 0.0 54.78165000 -99.6030000
+-6.750 -3.000 0.0 0.0 54.78165000 -99.6030000
+-6.750 -1.501 0.0 0.0 54.78165000 -99.6030000
+-6.750 -1.500 0.0 0.0 60.20274177 -99.6030000
+-6.750 -1.499 0.0 0.0 65.62383355 -99.6030000
+-6.750 +1.499 0.0 0.0 65.62383355 -99.6030000
+-6.750 +1.500 0.0 0.0 60.20274177 -99.6030000
+-6.750 +1.501 0.0 0.0 54.78165000 -99.6030000
+-6.750 +3.000 0.0 0.0 54.78165000 -99.6030000
+-6.750 +6.000 0.0 0.0 54.78165000 -99.6030000
+-6.750 +15.000 0.0 0.0 54.78165000 -99.6030000
+
+-6.751 -15.000 0.0 0.0 54.78976580 -99.6177560
+-6.751 -6.000 0.0 0.0 54.78976580 -99.6177560
+-6.751 -3.000 0.0 0.0 54.78976580 -99.6177560
+-6.751 -1.501 0.0 0.0 54.78976580 -99.6177560
+-6.751 -1.500 0.0 0.0 54.78976580 -99.6177560
+-6.751 -1.499 0.0 0.0 54.78976580 -99.6177560
+-6.751 +1.499 0.0 0.0 54.78976580 -99.6177560
+-6.751 +1.500 0.0 0.0 54.78976580 -99.6177560
+-6.751 +1.501 0.0 0.0 54.78976580 -99.6177560
+-6.751 +3.000 0.0 0.0 54.78976580 -99.6177560
+-6.751 +6.000 0.0 0.0 54.78976580 -99.6177560
+-6.751 +15.000 0.0 0.0 54.78976580 -99.6177560
+
+-7.500 -15.000 0.0 0.0 60.86850000 -110.6700000
+-7.500 -6.000 0.0 0.0 60.86850000 -110.6700000
+-7.500 -3.000 0.0 0.0 60.86850000 -110.6700000
+-7.500 -1.501 0.0 0.0 60.86850000 -110.6700000
+-7.500 -1.500 0.0 0.0 60.86850000 -110.6700000
+-7.500 -1.499 0.0 0.0 60.86850000 -110.6700000
+-7.500 +1.499 0.0 0.0 60.86850000 -110.6700000
+-7.500 +1.500 0.0 0.0 60.86850000 -110.6700000
+-7.500 +1.501 0.0 0.0 60.86850000 -110.6700000
+-7.500 +3.000 0.0 0.0 60.86850000 -110.6700000
+-7.500 +6.000 0.0 0.0 60.86850000 -110.6700000
+-7.500 +15.000 0.0 0.0 60.86850000 -110.6700000
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tractions.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tractions.spatialdb 2012-03-20 16:41:07 UTC (rev 19828)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tractions.spatialdb 2012-03-20 23:04:38 UTC (rev 19829)
@@ -1,59 +0,0 @@
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 3
- value-names = traction-shear-leftlateral traction-shear-updip traction-normal
- value-units = MPa MPa MPa
- num-locs = 36
- data-dim = 2
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0e+3 // x, y, z are in km
- space-dim = 3
- }
-}
-0.0 -15.0 0.0 0.0 0.0 0.0
-0.0 -1.5001 0.0 0.0 0.0 0.0
-0.0 -1.5 0.0 0.0 0.0 0.0
-0.0 +1.5 0.0 0.0 0.0 0.0
-0.0 +1.5001 0.0 0.0 0.0 0.0
-0.0 +15.0 0.0 0.0 0.0 0.0
-
--5.249500000000000 -15.0 -9.092400714332822 0.0 42.603892100000 -77.461622000000
--5.249500000000000 -1.5001 -9.092400714332822 0.0 42.603892100000 -77.461622000000
--5.249500000000000 -1.5 -9.092400714332822 0.0 42.603892100000 -77.461622000000
--5.249500000000000 +1.5 -9.092400714332822 0.0 42.603892100000 -77.461622000000
--5.249500000000000 +1.5001 -9.092400714332822 0.0 42.603892100000 -77.461622000000
--5.249500000000000 +15.0 -9.092400714332822 0.0 42.603892100000 -77.461622000000
-
--5.250000000000000 -15.0 -9.093266739736606 0.0 42.407950000000 -77.469000000000
--5.250000000000000 -1.5001 -9.093266739736606 0.0 42.407950000000 -77.469000000000
--5.250000000000000 -1.5 -9.093266739736606 0.0 59.518013300000 -77.469000000000
--5.250000000000000 +1.5 -9.093266739736606 0.0 59.518013300000 -77.469000000000
--5.250000000000000 +1.5001 -9.093266739736606 0.0 42.407950000000 -77.469000000000
--5.250000000000000 +15.0 -9.093266739736606 0.0 42.407950000000 -77.469000000000
-
--6.750000000000000 -15.0 -11.691342951089920 0.0 54.781650000000 -99.603000000000
--6.750000000000000 -1.5001 -11.691342951089920 0.0 54.781650000000 -99.603000000000
--6.750000000000000 -1.5 -11.691342951089920 0.0 76.466017100000 -99.603000000000
--6.750000000000000 +1.5 -11.691342951089920 0.0 76.466017100000 -99.603000000000
--6.750000000000000 +1.5001 -11.691342951089920 0.0 54.781650000000 -99.603000000000
--6.750000000000000 +15.0 -11.691342951089920 0.0 54.781650000000 -99.603000000000
-
--6.750500000000000 -15.0 -11.692208976493704 0.0 54.785707900000 -99.610378000000
--6.750500000000000 -1.5001 -11.692208976493704 0.0 54.785707900000 -99.610378000000
--6.750500000000000 -1.5 -11.692208976493704 0.0 54.785707900000 -99.610378000000
--6.750500000000000 +1.5 -11.692208976493704 0.0 54.785707900000 -99.610378000000
--6.750500000000000 +1.5001 -11.692208976493704 0.0 54.785707900000 -99.610378000000
--6.750500000000000 +15.0 -11.692208976493704 0.0 54.785707900000 -99.610378000000
-
--7.500000000000000 -15.0 -12.990381056766578 0.0 60.868500000000 -110.670000000000
--7.500000000000000 -1.5001 -12.990381056766578 0.0 60.868500000000 -110.670000000000
--7.500000000000000 -1.5 -12.990381056766578 0.0 60.868500000000 -110.670000000000
--7.500000000000000 +1.5 -12.990381056766578 0.0 60.868500000000 -110.670000000000
--7.500000000000000 +1.5001 -12.990381056766578 0.0 60.868500000000 -110.670000000000
--7.500000000000000 +15.0 -12.990381056766578 0.0 60.868500000000 -110.670000000000
-
-
-
-
-
More information about the CIG-COMMITS
mailing list