[cig-commits] r16792 - in short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205: . scecfiles
surendra at geodynamics.org
surendra at geodynamics.org
Wed May 26 00:39:40 PDT 2010
Author: surendra
Date: 2010-05-26 00:39:40 -0700 (Wed, 26 May 2010)
New Revision: 16792
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/geometry_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m_refine.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m_refine.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/output/
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/scecfiles/
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/scecfiles/hex8_100m/
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/scecfiles/hex8_200m/
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/scecfiles/tet4_100m/
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/scecfiles/tet4_200m/
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m_refine.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m_refine.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m_refine.jou
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_onfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m.cfg
Log:
Added refined meshes for tpv205 and post-processing script to extract off-fault data
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou 2010-05-26 06:04:39 UTC (rev 16791)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou 2010-05-26 07:39:40 UTC (rev 16792)
@@ -7,10 +7,11 @@
# ----------------------------------------------------------------------
# Create nodeset for fault
# ----------------------------------------------------------------------
-group "fault" add node in fault_surface
-group "fault" remove node with x_coord > +15.0001e+3
-group "fault" remove node with x_coord < -15.0001e+3
-group "fault" remove node with z_coord < -15.0001e+3
+group "fault" add node with x_min >= -15.001e+3 and x_max <= 15.001e+3 and z_min >= -15.001e+3
+#group "fault" add node in fault_surface
+#group "fault" remove node with x_coord > +15.0001e+3
+#group "fault" remove node with x_coord < -15.0001e+3
+#group "fault" remove node with z_coord < -15.0001e+3
nodeset 10 group fault
nodeset 10 name "fault"
@@ -69,4 +70,11 @@
group "face_zpos_subset" remove node with y_coord > +5.001e+3
group "face_zpos_subset" remove node with y_coord < -5.001e+3
nodeset 20 group face_zpos_subset
-nodeset 20 name "face_zpos_subset"
+nodeset 20 name "face zpos_subset"
+
+# ----------------------------------------------------------------------
+# Create nodeset for off-fault stations
+# ----------------------------------------------------------------------
+group "off_fault" add node with x_min >= -12.001e+3 and x_max <= 12.001e+3 and z_min >= -7.501e+3 and z_coord = 3.000e+3
+nodeset 25 group off_fault
+nodeset 25 name "off_fault"
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/geometry_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/geometry_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/geometry_refine.jou 2010-05-26 07:39:40 UTC (rev 16792)
@@ -0,0 +1,34 @@
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+# Block is 128.0 km x 128.0 km
+# -64.0 km <= x <= 64.0 km
+# -64.0 km <= y <= 64.0 km
+# -64.0 km <= z <= 0.0 km
+#
+# We create the 2-D mesh by extracting the mid-surface from a 3-D mesh.
+reset
+brick x 128000 y 128000 z 64000
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -32000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane yplane offset 0
+surface 7 name "fault_surface"
+
+# ----------------------------------------------------------------------
+# Divide volumes using interface surfaces
+# ----------------------------------------------------------------------
+webcut volume 1 with plane surface fault_surface
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all
+merge all
+delete body 2
+
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m.cfg 2010-05-26 06:04:39 UTC (rev 16791)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m.cfg 2010-05-26 07:39:40 UTC (rev 16792)
@@ -17,9 +17,13 @@
# output
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
+[pylithapp.problem.formulation.output.domain]
writer.filename = output/hex8_100m.vtk
+# Give basename for VTK domain output of solution over subdomain.
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_100m-off_fault.vtk
+
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
writer.filename = output/hex8_100m-fault.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m_refine.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m_refine.cfg 2010-05-26 07:39:40 UTC (rev 16792)
@@ -0,0 +1,33 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8_100m_refine.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8_100m_refine.vtk
+
+# Give basename for VTK domain output of solution over subdomain.
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_100m_refine-off_fault.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/hex8_100m_refine-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8_100m_refine-elastic.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m_refine.jou 2010-05-26 07:39:40 UTC (rev 16792)
@@ -0,0 +1,40 @@
+# ----------------------------------------------------------------------
+# Create hex8 mesh at 100m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry_refine.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size 400
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine elements within 1/2 of rupture width
+# ----------------------------------------------------------------------
+playback 'refinecells.jou'
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+volume all smooth scheme condition number beta 1.3 cpu 10
+smooth volume all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "hex8_100m_refine.exo" dimension 3 overwrite
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m.cfg 2010-05-26 06:04:39 UTC (rev 16791)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m.cfg 2010-05-26 07:39:40 UTC (rev 16792)
@@ -18,9 +18,13 @@
# output
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
+[pylithapp.problem.formulation.output.domain]
writer.filename = output/hex8_200m.vtk
+# Give basename for VTK domain output of solution over subdomain.
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_200m-off_fault.vtk
+
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
writer.filename = output/hex8_200m-fault.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m_refine.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m_refine.cfg 2010-05-26 07:39:40 UTC (rev 16792)
@@ -0,0 +1,34 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8_200m_refine.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01666666666666666666666667*s
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8_200m_refine.vtk
+
+# Give basename for VTK domain output of solution over subdomain.
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_200m_refine-off_fault.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/hex8_200m_refine-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8_200m_refine-elastic.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m_refine.jou 2010-05-26 07:39:40 UTC (rev 16792)
@@ -0,0 +1,41 @@
+# ----------------------------------------------------------------------
+# Create hex8 mesh at 200m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry_refine.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size 800
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine elements within 1/2 of rupture width
+# ----------------------------------------------------------------------
+playback 'refinecells.jou'
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+volume all smooth scheme condition number beta 1.3 cpu 10
+smooth volume all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "hex8_200m_refine.exo" dimension 3 overwrite
+
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/pylithapp.cfg 2010-05-26 06:04:39 UTC (rev 16791)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/pylithapp.cfg 2010-05-26 07:39:40 UTC (rev 16792)
@@ -68,8 +68,9 @@
dt = 0.01*s
-# Set subdomain component to OutputSolnSubset (subset of domain).
-#output.subdomain = pylith.meshio.OutputSolnSubset
+[pylithapp.timedependent.implicit]
+output = [domain,subdomain]
+output.subdomain = pylith.meshio.OutputSolnSubset
# ----------------------------------------------------------------------
@@ -141,12 +142,19 @@
# output
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
+[pylithapp.problem.formulation.output.domain]
vertex_data_fields=[displacement,velocity]
output_freq = time_step
time_step = 12.0*s
writer.time_format = %06.3f
+[pylithapp.problem.formulation.output.subdomain]
+label = off_fault ; nodeset for subdomain
+vertex_data_fields=[displacement,velocity]
+output_freq = time_step
+time_step = 0.04999*s
+writer.time_format = %06.3f
+
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
vertex_info_fields = [strike_dir,normal_dir,initial_traction]
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells.jou 2010-05-26 07:39:40 UTC (rev 16792)
@@ -0,0 +1,7 @@
+# ----------------------------------------------------------------------
+# Refine elements within 1/2 of rupture width
+# ----------------------------------------------------------------------
+group "refinecells1" add node with x_min >= -24.001e+3 and x_max <= 24.001e+3 and y_min >= -36.001e+3 and y_max <= 36.001e+3
+refine node in refinecells1 depth 1
+group "refinecells2" add node with x_min >= -12.001e+3 and x_max <= 12.001e+3 and y_min >= -24.001e+3 and y_max <= 24.001e+3
+refine node in refinecells2 depth 1
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_onfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_onfaultdata.py 2010-05-26 06:04:39 UTC (rev 16791)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_onfaultdata.py 2010-05-26 07:39:40 UTC (rev 16792)
@@ -10,7 +10,7 @@
# ----------------------------------------------------------------------
#
-cell = "hex8"
+cell = "tet4"
dx = 200
dt = 0.05
@@ -25,22 +25,41 @@
# ----------------------------------------------------------------------
timestamps = numpy.arange(50,12001,50)
-targets = numpy.array([[-12000.0, 0.0, 0.0],
- [ -7400.0, 0.0, 0.0],
- [ -4600.0, 0.0, 0.0],
- [ 0.0, 0.0, 0.0],
- [ 4600.0, 0.0, 0.0],
- [ 7400.0, 0.0, 0.0],
- [ 12000.0, 0.0, 0.0],
- [ 0.0, 0.0, -3000.0],
- [-12000.0, 0.0, -7600.0],
- [ -7400.0, 0.0, -7600.0],
- [ -4600.0, 0.0, -7600.0],
- [ 0.0, 0.0, -7600.0],
- [ 4600.0, 0.0, -7600.0],
- [ 7400.0, 0.0, -7600.0],
- [ 12000.0, 0.0, -7600.0],
- [ 0.0, 0.0, -12000.0]])
+if dx == 200:
+ targets = numpy.array([[-12000.0, 0.0, 0.0],
+ [ -7400.0, 0.0, 0.0],
+ [ -4600.0, 0.0, 0.0],
+ [ 0.0, 0.0, 0.0],
+ [ 4600.0, 0.0, 0.0],
+ [ 7400.0, 0.0, 0.0],
+ [ 12000.0, 0.0, 0.0],
+ [ 0.0, 0.0, -3000.0],
+ [-12000.0, 0.0, -7600.0],
+ [ -7400.0, 0.0, -7600.0],
+ [ -4600.0, 0.0, -7600.0],
+ [ 0.0, 0.0, -7600.0],
+ [ 4600.0, 0.0, -7600.0],
+ [ 7400.0, 0.0, -7600.0],
+ [ 12000.0, 0.0, -7600.0],
+ [ 0.0, 0.0, -12000.0]])
+elif dx == 100:
+ targets = numpy.array([[-12000.0, 0.0, 0.0],
+ [ -7500.0, 0.0, 0.0],
+ [ -4500.0, 0.0, 0.0],
+ [ 0.0, 0.0, 0.0],
+ [ 4500.0, 0.0, 0.0],
+ [ 7500.0, 0.0, 0.0],
+ [ 12000.0, 0.0, 0.0],
+ [ 0.0, 0.0, -3000.0],
+ [-12000.0, 0.0, -7500.0],
+ [ -7500.0, 0.0, -7500.0],
+ [ -4500.0, 0.0, -7500.0],
+ [ 0.0, 0.0, -7500.0],
+ [ 4500.0, 0.0, -7500.0],
+ [ 7500.0, 0.0, -7500.0],
+ [ 12000.0, 0.0, -7500.0],
+ [ 0.0, 0.0, -12000.0]])
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py 2010-05-26 06:04:39 UTC (rev 16791)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py 2010-05-26 07:39:40 UTC (rev 16792)
@@ -92,12 +92,9 @@
time = timestamps / timeScale
print "time", time
-filename = "%s_ruptime.dat" % (outdir)
+filename = "%sruptime.dat" % (outdir)
fout = open(filename, "w")
fout.write(headerA)
-print "distStrike.shape", distStrike.shape
-print "distDip.shape", distDip.shape
-print "ruptimeIndices.shape", ruptimeIndices.shape
-data = numpy.transpose((distStrike, distDip, ruptimeIndices))
+data = numpy.transpose((distStrike, distDip, rupTime))
numpy.savetxt(fout, data, fmt='%14.6e')
fout.close()
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m.cfg 2010-05-26 06:04:39 UTC (rev 16791)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m.cfg 2010-05-26 07:39:40 UTC (rev 16792)
@@ -17,9 +17,13 @@
# output
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
+[pylithapp.problem.formulation.output.domain]
writer.filename = output/tet4_100m.vtk
+# Give basename for VTK domain output of solution over subdomain.
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tet4_100m-off_fault.vtk
+
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
writer.filename = output/tet4_100m-fault.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m_refine.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m_refine.cfg 2010-05-26 07:39:40 UTC (rev 16792)
@@ -0,0 +1,33 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_100m_refine.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tet4_100m_refine.vtk
+
+# Give basename for VTK domain output of solution over subdomain.
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tet4_100m_refine-off_fault.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tet4_100m_refine-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tet4_100m_refine-elastic.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m_refine.jou 2010-05-26 07:39:40 UTC (rev 16792)
@@ -0,0 +1,46 @@
+# ----------------------------------------------------------------------
+# Create tet4 mesh at 100m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry_refine.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size and scheme
+# ----------------------------------------------------------------------
+volume all size 400
+volume all scheme tetmesh
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine elements within 1/2 of rupture width
+# ----------------------------------------------------------------------
+playback 'refinecells.jou'
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+volume all smooth scheme laplacian
+smooth volume all
+volume all smooth scheme condition number beta 1.3 cpu 10
+smooth volume all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "tet4_100m_refine.exo" dimension 3 overwrite
+
+
+
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m.cfg 2010-05-26 06:04:39 UTC (rev 16791)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m.cfg 2010-05-26 07:39:40 UTC (rev 16792)
@@ -17,9 +17,13 @@
# output
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
+[pylithapp.problem.formulation.output.domain]
writer.filename = output/tet4_200m.vtk
+# Give basename for VTK domain output of solution over subdomain.
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tet4_200m-off_fault.vtk
+
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
writer.filename = output/tet4_200m-fault.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m_refine.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m_refine.cfg 2010-05-26 07:39:40 UTC (rev 16792)
@@ -0,0 +1,33 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_200m_refine.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01666666666666666666666667*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tet4_200m_refine.vtk
+
+# Give basename for VTK domain output of solution over subdomain.
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tet4_200m_refine-off_fault.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tet4_200m_refine-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tet4_200m_refine-elastic.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m_refine.jou 2010-05-26 07:39:40 UTC (rev 16792)
@@ -0,0 +1,46 @@
+# ----------------------------------------------------------------------
+# Create tet4 mesh at 200m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry_refine.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size and scheme
+# ----------------------------------------------------------------------
+volume all size 800
+volume all scheme tetmesh
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine elements within 1/2 of rupture width
+# ----------------------------------------------------------------------
+playback 'refinecells.jou'
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+volume all smooth scheme laplacian
+smooth volume all
+volume all smooth scheme condition number beta 1.3 cpu 10
+smooth volume all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "tet4_200m_refine.exo" dimension 3 overwrite
+
+
+
+
More information about the CIG-COMMITS
mailing list