[cig-commits] r17887 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14
brad at geodynamics.org
brad at geodynamics.org
Wed Feb 16 18:12:00 PST 2011
Author: brad
Date: 2011-02-16 18:11:59 -0800 (Wed, 16 Feb 2011)
New Revision: 17887
Removed:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/output/
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/friction_branch.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/friction_main.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/geometry.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/gradient.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/params.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tet4_200m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tpv14_tet4_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tpv14_tractions.spatialdb
Log:
Worked on tet4 mesh.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/createbc.jou 2011-02-17 01:13:56 UTC (rev 17886)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/createbc.jou 2011-02-17 02:11:59 UTC (rev 17887)
@@ -21,40 +21,40 @@
nodeset 101 name "branch_fault"
# ----------------------------------------------------------------------
-# Create nodeset for +x face
+# Create nodeset for +y face
# ----------------------------------------------------------------------
-group "face_xpos" add node in surface 10
-group "face_xpos" add node in surface 20
-group "face_xpos" add node in surface 25
-nodeset 11 group face_xpos
-nodeset 11 name "face_xpos"
+group "face_ypos" add node in surface 15
+group "face_ypos" add node in surface 20
+group "face_ypos" add node in surface 25
+nodeset 11 group face_ypos
+nodeset 11 name "face_ypos"
# ----------------------------------------------------------------------
-# Create nodeset for -x face
+# Create nodeset for -y face
# ----------------------------------------------------------------------
-group "face_xneg" add node in surface 12
-group "face_xneg" add node in surface 15
-nodeset 12 group face_xneg
-nodeset 12 name "face_xneg"
+group "face_yneg" add node in surface 10
+group "face_yneg" add node in surface 17
+nodeset 12 group face_yneg
+nodeset 12 name "face_yneg"
# ----------------------------------------------------------------------
-# Create nodeset for +y face
+# Create nodeset for +x face
# ----------------------------------------------------------------------
-group "face_ypos" add node in surface 5
-nodeset 13 group face_ypos
-nodeset 13 name "face_ypos"
+group "face_xpos" add node in surface 6
+nodeset 13 group face_xpos
+nodeset 13 name "face_xpos"
# ----------------------------------------------------------------------
-# Create nodeset for -y face
+# Create nodeset for -x face
# ----------------------------------------------------------------------
-group "face_yneg" add node in surface 3
-nodeset 14 group face_yneg
-nodeset 14 name "face_yneg"
+group "face_xneg" add node in surface 4
+nodeset 14 group face_xneg
+nodeset 14 name "face_xneg"
# ----------------------------------------------------------------------
# Create nodeset for -z face
# ----------------------------------------------------------------------
-group "face_zneg" add node in surface 11
+group "face_zneg" add node in surface 16
group "face_zneg" add node in surface 21
group "face_zneg" add node in surface 28
nodeset 15 group face_zneg
@@ -63,7 +63,7 @@
# ----------------------------------------------------------------------
# Create nodeset for +z face
# ----------------------------------------------------------------------
-group "face_zpos" add node in surface 9
+group "face_zpos" add node in surface 14
group "face_zpos" add node in surface 23
group "face_zpos" add node in surface 26
nodeset 16 group face_zpos
@@ -81,3 +81,4 @@
nodeset 20 name "face zpos_subset"
# End of file
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/friction_branch.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/friction_branch.spatialdb 2011-02-17 01:13:56 UTC (rev 17886)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/friction_branch.spatialdb 2011-02-17 02:11:59 UTC (rev 17887)
@@ -4,19 +4,19 @@
value-names = static-coefficient dynamic-coefficient slip-weakening-parameter cohesion
value-units = none none m MPa
num-locs = 8
- data-dim = 1
+ 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 0.1 -14.999 10000.0 0.525 0.4 0.0
-0.0 0.1001 -14.999 0.677 0.525 0.4 0.0
-0.0 5.9999 -14.999 0.677 0.525 0.4 0.0
-0.0 6.0 -14.999 10000.0 0.525 0.4 0.0
-
-0.0 0.1 -15.0 10000.0 0.525 0.4 0.0
-0.0 0.1001 -15.0 10000.0 0.525 0.4 0.0
-0.0 5.9999 -15.0 10000.0 0.525 0.4 0.0
-0.0 6.0 -15.0 10000.0 0.525 0.4 0.0
+ 0.050 0.0 -14.99 10000.0 0.525 0.4 0.0
+ 0.051 0.0 -14.99 0.677 0.525 0.4 0.0
+ 5.999 0.0 -14.99 0.677 0.525 0.4 0.0
+ 6.000 0.0 -14.99 10000.0 0.525 0.4 0.0
+
+ 0.050 0.0 -15.00 10000.0 0.525 0.4 0.0
+ 0.051 0.0 -15.00 10000.0 0.525 0.4 0.0
+ 5.999 0.0 -15.00 10000.0 0.525 0.4 0.0
+ 6.000 0.0 -15.00 10000.0 0.525 0.4 0.0
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/friction_main.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/friction_main.spatialdb 2011-02-17 01:13:56 UTC (rev 17886)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/friction_main.spatialdb 2011-02-17 02:11:59 UTC (rev 17887)
@@ -11,12 +11,12 @@
space-dim = 3
}
}
--13.99 0.0 -14.98 10000.0 0.525 0.4 0.0
--13.98 0.0 -14.98 0.677 0.525 0.4 0.0
-+13.98 0.0 -14.98 0.677 0.525 0.4 0.0
-+13.99 0.0 -14.98 10000.0 0.525 0.4 0.0
-
--13.99 0.0 -14.99 10000.0 0.525 0.4 0.0
--13.98 0.0 -14.99 10000.0 0.525 0.4 0.0
-+13.98 0.0 -14.99 10000.0 0.525 0.4 0.0
-+13.99 0.0 -14.99 10000.0 0.525 0.4 0.0
+0.0 -13.99 -14.98 10000.0 0.525 0.4 0.0
+0.0 -13.98 -14.98 0.677 0.525 0.4 0.0
+0.0 +13.98 -14.98 0.677 0.525 0.4 0.0
+0.0 +13.99 -14.98 10000.0 0.525 0.4 0.0
+
+0.0 -13.99 -14.99 10000.0 0.525 0.4 0.0
+0.0 -13.98 -14.99 10000.0 0.525 0.4 0.0
+0.0 +13.98 -14.99 10000.0 0.525 0.4 0.0
+0.0 +13.99 -14.99 10000.0 0.525 0.4 0.0
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/geometry.jou 2011-02-17 01:13:56 UTC (rev 17886)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/geometry.jou 2011-02-17 02:11:59 UTC (rev 17887)
@@ -6,7 +6,7 @@
# ----------------------------------------------------------------------
# Create and position block.
# ----------------------------------------------------------------------
-brick x {blockLength} y {blockWidth} z {blockHeight}
+brick x {blockWidth} y {blockLength} z {blockHeight}
volume 1 move z {-blockHeight/2.0}
# ----------------------------------------------------------------------
@@ -14,17 +14,17 @@
# ----------------------------------------------------------------------
# Create main fault
-create planar surface with plane yplane offset 0
+create planar surface with plane xplane offset 0
surface 7 name "temp_main"
webcut volume 1 with sheet surface temp_main
delete body 2
# Create branch fault
-create planar surface with plane xplane offset {blockLength/2.0}
+create planar surface with plane yplane offset 0.0
surface 18 name "temp_branch"
-rotate surface temp_branch about vertex 2 5 angle {branchAngle-90}
-move surface temp_branch x {branchOffsetY-blockLength/2.0} y {-blockWidth/2.0}
-webcut volume 3 with sheet surface temp_branch
+move surface temp_branch x {blockWidth/2} y {branchOffsetY} z 0
+rotate surface temp_branch about vertex 21 vertex 22 angle {90-branchAngle}
+webcut volume 1 with sheet surface temp_branch
delete body 4
# ----------------------------------------------------------------------
@@ -32,20 +32,20 @@
# ----------------------------------------------------------------------
# Split main fault surface using volume 1
-create vertex {-mainFaultLength/2} 0 0
-create vertex {+mainFaultLength/2} 0 0
-create vertex {-mainFaultLength/2} 0 {-faultWidth}
-create vertex {+mainFaultLength/2} 0 {-faultWidth}
-create vertex {-mainFaultLength/2} 0 {-blockHeight}
-create vertex {+mainFaultLength/2} 0 {-blockHeight}
-create vertex {-blockLength/2} 0 {-faultWidth}
-create vertex {+blockLength/2} 0 {-faultWidth}
+create vertex 0 {-mainFaultLength/2} 0
+create vertex 0 {+mainFaultLength/2} 0
+create vertex 0 {-mainFaultLength/2} {-faultWidth}
+create vertex 0 {+mainFaultLength/2} {-faultWidth}
+create vertex 0 {-mainFaultLength/2} {-blockHeight}
+create vertex 0 {+mainFaultLength/2} {-blockHeight}
+create vertex 0 {-blockLength/2} {-faultWidth}
+create vertex 0 {+blockLength/2} {-faultWidth}
-split surface 8 across location vertex 39 location vertex 40
-split surface 30 across location vertex 33 location vertex 35
-split surface 29 across location vertex 35 location vertex 37
-split surface 32 across location vertex 34 location vertex 36
-split surface 34 across location vertex 36 location vertex 38
+split surface 13 across location vertex 39 location vertex 40
+split surface 29 across location vertex 33 location vertex 35
+split surface 30 across location vertex 35 location vertex 37
+split surface 31 across location vertex 34 location vertex 36
+split surface 33 across location vertex 36 location vertex 38
delete vertex all
# Split branch fault surface using volume 3
@@ -66,7 +66,9 @@
imprint all
merge all
-surface 45 name "fault_surface_mainB"
-surface 48 name "fault_surface_mainA"
+surface 47 name "fault_surface_mainB"
+surface 52 name "fault_surface_mainA"
surface 41 name "fault_surface_branch"
+
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/gradient.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/gradient.jou 2011-02-17 01:13:56 UTC (rev 17886)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/gradient.jou 2011-02-17 02:11:59 UTC (rev 17887)
@@ -1,113 +1,133 @@
# ----------------------------------------------------------------------
# Set vertex spacing with increasing spacing away from fault
# ----------------------------------------------------------------------
-#{bias_factor=1.01}
+#{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 13 15 41 42
+#{dxStart=dx}{curveL=Length(68)}{execute("dxA="//fun_dx)}
+
+# dxB - size at vertex 28 45 48 58
+#{dxStart=dx}{curveL=Length(63)}{execute("dxB="//fun_dx)}
+
+# dxC - size at vertex 25 55
+#{dxStart=dx}{curveL=Length(83)}{execute("dxC="//fun_dx)}
+
+# dxD - size at vertex 14 16
+#{dxStart=dxA}{curveL=Length(57)}{execute("dxD="//fun_dx)}
+
+# dxE - size at vertex 26
+#{dxStart=dxC}{curveL=Length(76)}{execute("dxE="//fun_dx)}
+
+# dxF - size at vertex 2
+#{dxStart=dxC}{curveL=Length(49)}{execute("dxF="//fun_dx)}
+
+# dxG - size at vertex 1 3 4
+#{dxStart=dxA}{curveL=Length(30)}{execute("dxG="//fun_dx)}
+
+# dxH - size at vertex 5
+#{dxStart=dxE}{curveL=Length(50)}{execute("dxH="//fun_dx)}
+
+# dxE - size at vertex 6 7 8
+#{dxStart=dxD}{curveL=Length(32)}{execute("dxI="//fun_dx)}
+
+# ----------------------------------------------------------------------
+# Reset sizes
+# ----------------------------------------------------------------------
curve all scheme default
surface all sizing function none
surface fault_surface_mainA fault_surface_mainB fault_surface_branch size {dx}
# Set bias on curves extending from faults
-curve 70 scheme bias fine size {dx} factor {bias_factor} start vertex 46
-curve 69 scheme bias fine size {dx} factor {bias_factor} start vertex 47
-curve 71 scheme bias fine size {dx} factor {bias_factor} start vertex 47
-curve 59 scheme bias fine size {dx} factor {bias_factor} start vertex 43
-curve 60 scheme bias fine size {dx} factor {bias_factor} start vertex 44
-curve 63 scheme bias fine size {dx} factor {bias_factor} start vertex 44
+# Main fault
+curve 67 scheme bias fine size {dx} factor {bias_factor} start vertex 46
+curve 68 scheme bias fine size {dx} factor {bias_factor} start vertex 47
+curve 71 scheme bias fine size {dx} factor {bias_factor} start vertex 46
+curve 61 scheme bias fine size {dx} factor {bias_factor} start vertex 44
+curve 62 scheme bias fine size {dx} factor {bias_factor} start vertex 43
+curve 63 scheme bias fine size {dx} factor {bias_factor} start vertex 43
+
+# Branch fault
curve 83 scheme bias fine size {dx} factor {bias_factor} start vertex 56
curve 82 scheme bias fine size {dx} factor {bias_factor} start vertex 57
curve 84 scheme bias fine size {dx} factor {bias_factor} start vertex 57
+# Intersection
curve 75 scheme bias fine size {dx} factor {bias_factor} start vertex 54
# ----------------------------------------------------------------------
-# Extend bias based on 1st set of bias values
-
-# Must be set manually.
-#{dxA=452.58} # end size of curve 59
-#{dxB=440.784} # end size of curve 63
-#{dxC=511.68} # end size of curve 83
-
# A to A
-curve 56 57 size {dxA}
-curve 85 92 93 size {dxB}
+curve 54 55 size {dxA}
+curve 85 88 99 size {dxB}
# C to C
curve 77 size {dxC}
# bias starting at A
-curve 21 scheme bias fine size {dxA} factor {bias_factor} start vertex 16
-curve 22 scheme bias fine size {dxA} factor {bias_factor} start vertex 14
-curve 29 scheme bias fine size {dxA} factor {bias_factor} start vertex 14
+curve 21 scheme bias fine size {dxA} factor {bias_factor} start vertex 13
+curve 30 scheme bias fine size {dxA} factor {bias_factor} start vertex 13
+curve 29 scheme bias fine size {dxA} factor {bias_factor} start vertex 15
# bias starting at B
-curve 64 scheme bias fine size {dxB} factor {bias_factor} start vertex 45
-curve 73 scheme bias fine size {dxB} factor {bias_factor} start vertex 48
+curve 65 scheme bias fine size {dxB} factor {bias_factor} start vertex 45
+curve 72 scheme bias fine size {dxB} factor {bias_factor} start vertex 48
curve 86 scheme bias fine size {dxB} factor {bias_factor} start vertex 58
# bias starting at C
curve 49 scheme bias fine size {dxC} factor {bias_factor} start vertex 25
# bias A to C
-curve 42 scheme bias fine size {dxA} coarse size {dxC} start vertex 16
-curve 42 scheme bias fine size {dxA} coarse size {dxC} start vertex 16
+curve 42 scheme bias fine size {dxA} coarse size {dxC} start vertex 15
-# ----------------------------------------------------------------------
-# Extend bias based on 2nd set of bias values
-
-# Must be set manually.
-#{dxD=682.805} # end size of curve 64
-#{dxE=744.36} # end size of curve 86
-#{dxF=608.892} # end size of curve 49
-#{dxG=758.417} # end size of curve 22
-
# bias starting at D
-curve 24 scheme bias fine size {dxD} factor {bias_factor} start vertex 13
-curve 31 scheme bias fine size {dxD} factor {bias_factor} start vertex 13
-curve 23 scheme bias fine size {dxD} factor {bias_factor} start vertex 15
+curve 32 scheme bias fine size {dxD} factor {bias_factor} start vertex 14
+curve 23 scheme bias fine size {dxD} factor {bias_factor} start vertex 14
+curve 31 scheme bias fine size {dxD} factor {bias_factor} start vertex 16
# bias starting at E
curve 50 scheme bias fine size {dxE} factor {bias_factor} start vertex 26
# bias A to D
-curve 54 scheme bias fine size {dxA} coarse size {dxD} start vertex 41
-curve 55 scheme bias fine size {dxA} coarse size {dxD} start vertex 42
+curve 57 scheme bias fine size {dxA} coarse size {dxD} start vertex 41
+curve 56 scheme bias fine size {dxA} coarse size {dxD} start vertex 42
# bias D to E
-curve 41 scheme bias fine size {dxD} coarse size {dxE} start vertex 15
+curve 41 scheme bias fine size {dxD} coarse size {dxE} start vertex 16
# bias F to G
-curve 4 scheme bias fine size {dxF} coarse size {dxG} start vertex 1
+curve 1 scheme bias fine size {dxF} coarse size {dxG} start vertex 2
# bias C to E
curve 76 scheme bias fine size {dxC} coarse size {dxE} start vertex 55
# G to G
-curve 2 size {dxG}
+curve 3 size {dxG}
-# ----------------------------------------------------------------------
-# Extend bias based on 3rd set of bias values
-# Must be set manually.
-#{dxH=804.251} # end size of curve 50
-#{dxI=975.139} # end size of curve 24
-
# bias G to I
+curve 10 scheme bias fine size {dxG} coarse size {dxI} start vertex 1
curve 9 scheme bias fine size {dxG} coarse size {dxI} start vertex 4
curve 11 scheme bias fine size {dxG} coarse size {dxI} start vertex 3
-curve 12 scheme bias fine size {dxG} coarse size {dxI} start vertex 2
# bias F to H
-curve 10 scheme bias fine size {dxF} coarse size {dxH} start vertex 1
+curve 12 scheme bias fine size {dxF} coarse size {dxH} start vertex 2
+# bias E to H
+curve 50 scheme bias fine size {dxE} coarse size {dxH} start vertex 26
+
# bias H to I
-curve 6 scheme bias fine size {dxH} coarse size {dxI} start vertex 6
+curve 5 scheme bias fine size {dxH} coarse size {dxI} start vertex 5
# I to I
-curve 8 size {dxI}
+curve 7 size {dxI}
# End of file
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/params.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/params.jou 2011-02-17 01:13:56 UTC (rev 17886)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/params.jou 2011-02-17 02:11:59 UTC (rev 17887)
@@ -6,14 +6,14 @@
# Parameters for geometry.
# ----------------------------------------------------------------------
#
-# Block is 64.0 km x 80.0 km x 40.0 km
-# -32.0 km <= x <= 32.0 km
-# -40.0 km <= y <= 40.0 km
-# -40.0 km <= z <= 0 km
+# Block is 48.0 km x 64.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 km
#
-#{blockWidth=64*km}
-#{blockLength=80*km}
-#{blockHeight=40*km}
+#{blockWidth=48*km}
+#{blockLength=64*km}
+#{blockHeight=32*km}
#
# Faults have a common width of 15 km
#{faultWidth=15.0*km}
@@ -36,7 +36,9 @@
# ----------------------------------------------------------------------
# Parameters for off-fault stations on ground surface
# ----------------------------------------------------------------------
-#{offFaultMinX=-25.01*km}
-#{offFaultMaxX=+25.01*km}
-#{offFaultMinY=-20.01*km}
-#{offFaultMaxY=+20.01*km}
+#{offFaultMinX=-16.01*km}
+#{offFaultMaxX=+16.01*km}
+#{offFaultMinY=-24.01*km}
+#{offFaultMaxY=+24.01*km}
+
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/pylithapp.cfg 2011-02-17 01:13:56 UTC (rev 17886)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/pylithapp.cfg 2011-02-17 02:11:59 UTC (rev 17887)
@@ -27,7 +27,7 @@
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
reader = pylith.meshio.MeshIOCubit
-reorder_mesh = True
+#reorder_mesh = True
#refiner = pylith.topology.RefineUniform
[pylithapp.mesh_generator.reader]
@@ -109,6 +109,7 @@
[pylithapp.timedependent]
interfaces = [main_fault,branch_fault]
+
# Change the kinematic rupture to spontaneous rupture
interfaces.main_fault = pylith.faults.FaultCohesiveDyn
interfaces.branch_fault = pylith.faults.FaultCohesiveDyn
@@ -138,7 +139,6 @@
db_initial_tractions.values = [traction-shear-leftlateral,traction-shear-updip,traction-normal]
friction = pylith.friction.SlipWeakening
-
friction.db_properties = spatialdata.spatialdb.SimpleDB
friction.db_properties.iohandler.filename = friction_branch.spatialdb
friction.db_properties.label = Slip weakening
@@ -157,7 +157,7 @@
vertex_data_fields=[displacement,velocity]
output_freq = time_step
time_step = 0.04999*s
-#writer = pylith.meshio.DataWriterHDF5ExtSubMesh
+writer = pylith.meshio.DataWriterHDF5ExtSubMesh
[pylithapp.timedependent.interfaces.main_fault.output]
vertex_info_fields = [initial_traction,static_coefficient,dynamic_coefficient,slip_weakening_parameter]
@@ -166,7 +166,6 @@
time_step = 0.04999*s
writer = pylith.meshio.DataWriterHDF5ExtSubSubMesh
-# Give basename for VTK Branch fault output.
[pylithapp.timedependent.interfaces.branch_fault.output]
vertex_info_fields = [initial_traction,static_coefficient,dynamic_coefficient,slip_weakening_parameter]
vertex_data_fields = [slip,slip_rate,traction]
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tet4_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tet4_200m.jou 2011-02-17 01:13:56 UTC (rev 17886)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tet4_200m.jou 2011-02-17 02:11:59 UTC (rev 17887)
@@ -34,7 +34,6 @@
# ----------------------------------------------------------------------
# Smooth mesh to improve quality
# ----------------------------------------------------------------------
-cleanup volume all
volume all smooth scheme condition number beta 2.0 cpu 10
smooth volume all
volume all smooth scheme condition number beta 1.9 cpu 10
@@ -49,3 +48,4 @@
# Export exodus file
# ----------------------------------------------------------------------
export mesh "tet4_200m.exo" dimension 3 overwrite
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tpv14_tet4_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tpv14_tet4_200m.cfg 2011-02-17 01:13:56 UTC (rev 17886)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tpv14_tet4_200m.cfg 2011-02-17 02:11:59 UTC (rev 17887)
@@ -17,16 +17,37 @@
# output
# ----------------------------------------------------------------------
[pylithapp.problem.formulation.output.subdomain]
-writer.filename = output/tpv14_tet4_200m-groundsurf.h5
+writer = pylith.meshio.DataWriterVTKSubMesh
+writer.time_format = %06.2f
+writer.filename = output/tpv14_tet4_200m-groundsurf.vtk
[pylithapp.timedependent.interfaces.main_fault.output]
-writer.filename = output/tpv14_tet4_200m-main_fault.h5
+writer = pylith.meshio.DataWriterVTKSubSubMesh
+writer.time_format = %06.2f
+writer.filename = output/tpv14_tet4_200m-main_fault.vtk
[pylithapp.timedependent.interfaces.branch_fault.output]
-writer.filename = output/tpv14_tet4_200m-branch_fault.h5
+writer = pylith.meshio.DataWriterVTKSubSubMesh
+writer.time_format = %06.2f
+writer.filename = output/tpv14_tet4_200m-branch_fault.vtk
[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/tpv14_tet4_200m-elastic.h5
+writer = pylith.meshio.DataWriterVTKMesh
+writer.time_format = %06.2f
+writer.filename = output/tpv14_tet4_200m-elastic.vtk
+#[pylithapp.problem.formulation.output.subdomain]
+#writer.filename = output/tpv14_tet4_200m-groundsurf.h5
+#
+#[pylithapp.timedependent.interfaces.main_fault.output]
+#writer.filename = output/tpv14_tet4_200m-main_fault.h5
+#
+#[pylithapp.timedependent.interfaces.branch_fault.output]
+#writer.filename = output/tpv14_tet4_200m-branch_fault.h5
+#
+#[pylithapp.timedependent.materials.elastic.output]
+#writer.filename = output/tpv14_tet4_200m-elastic.h5
+
+
# End of file
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tpv14_tractions.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tpv14_tractions.spatialdb 2011-02-17 01:13:56 UTC (rev 17886)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tpv14_tractions.spatialdb 2011-02-17 02:11:59 UTC (rev 17887)
@@ -11,22 +11,22 @@
space-dim = 3
}
}
--7.501 0.0 -5.999 -70.0 0.0 -120.0
--7.5 0.0 -5.999 -70.0 0.0 -120.0
--4.5 0.0 -5.999 -70.0 0.0 -120.0
--4.499 0.0 -5.999 -70.0 0.0 -120.0
-
--7.501 0.0 -6.0 -70.0 0.0 -120.0
--7.5 0.0 -6.0 -81.6 0.0 -120.0
--4.5 0.0 -6.0 -81.6 0.0 -120.0
--4.499 0.0 -6.0 -70.0 0.0 -120.0
-
--7.501 0.0 -9.0 -70.0 0.0 -120.0
--7.5 0.0 -9.0 -81.6 0.0 -120.0
--4.5 0.0 -9.0 -81.6 0.0 -120.0
--4.499 0.0 -9.0 -70.0 0.0 -120.0
-
--7.501 0.0 -9.001 -70.0 0.0 -120.0
--7.5 0.0 -9.001 -70.0 0.0 -120.0
--4.5 0.0 -9.001 -70.0 0.0 -120.0
--4.499 0.0 -9.001 -70.0 0.0 -120.0
+0.0 -7.501 -5.999 -70.0 0.0 -120.0
+0.0 -7.5 -5.999 -70.0 0.0 -120.0
+0.0 -4.5 -5.999 -70.0 0.0 -120.0
+0.0 -4.499 -5.999 -70.0 0.0 -120.0
+
+0.0 -7.501 -6.0 -70.0 0.0 -120.0
+0.0 -7.5 -6.0 -81.6 0.0 -120.0
+0.0 -4.5 -6.0 -81.6 0.0 -120.0
+0.0 -4.499 -6.0 -70.0 0.0 -120.0
+
+0.0 -7.501 -9.0 -70.0 0.0 -120.0
+0.0 -7.5 -9.0 -81.6 0.0 -120.0
+0.0 -4.5 -9.0 -81.6 0.0 -120.0
+0.0 -4.499 -9.0 -70.0 0.0 -120.0
+
+0.0 -7.501 -9.001 -70.0 0.0 -120.0
+0.0 -7.5 -9.001 -70.0 0.0 -120.0
+0.0 -4.5 -9.001 -70.0 0.0 -120.0
+0.0 -4.499 -9.001 -70.0 0.0 -120.0
More information about the CIG-COMMITS
mailing list