[cig-commits] r16749 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d
surendra at geodynamics.org
surendra at geodynamics.org
Wed May 19 13:19:02 PDT 2010
Author: surendra
Date: 2010-05-19 13:19:01 -0700 (Wed, 19 May 2010)
New Revision: 16749
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m_refine.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m_refine.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/refinecells.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m_refine.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m_refine.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.cfg
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_onfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.jou
Log:
Added refined mesh files at 100m resolution for tpv205-2d
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m_refine.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m_refine.cfg 2010-05-19 20:19:01 UTC (rev 16749)
@@ -0,0 +1,30 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = quad4_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.output]
+writer.filename = output/quad4_100m_refine.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/quad4_100m_refine-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/quad4_100m_refine-elastic.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m_refine.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m_refine.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m_refine.jou 2010-05-19 20:19:01 UTC (rev 16749)
@@ -0,0 +1,39 @@
+# ----------------------------------------------------------------------
+# Create quad4 mesh at 100m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry_refine.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+surface all size 400
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Refine elements within 1/2 of rupture width
+# ----------------------------------------------------------------------
+playback 'refinecells.jou'
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "quad4_100m_refine.exo" dimension 2 overwrite
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.cfg 2010-05-19 16:53:17 UTC (rev 16748)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.cfg 2010-05-19 20:19:01 UTC (rev 16749)
@@ -21,15 +21,6 @@
[pylithapp.problem.formulation.output.output]
writer.filename = output/quad4_200m_refine.vtk
-# Give basename for VTK domain output of solution over Off-fault stations.
-#[pylithapp.problem.formulation.output.subdomain]
-#label = 15 ; nodeset for subdomain
-#vertex_data_fields=[displacement,velocity]
-#output_freq = time_step
-#time_step = 0.04999*s
-#writer.filename = output/quad4_200m_OffFault.vtk
-#writer.time_format = %05.2f
-
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
writer.filename = output/quad4_200m_refine-fault.vtk
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.jou 2010-05-19 16:53:17 UTC (rev 16748)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.jou 2010-05-19 20:19:01 UTC (rev 16749)
@@ -20,10 +20,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
+playback 'refinecells.jou'
# ----------------------------------------------------------------------
# Smooth mesh to improve quality
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/refinecells.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/refinecells.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/refinecells.jou 2010-05-19 20:19:01 UTC (rev 16749)
@@ -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-2d/tabulate_onfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_onfaultdata.py 2010-05-19 16:53:17 UTC (rev 16748)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_onfaultdata.py 2010-05-19 20:19:01 UTC (rev 16749)
@@ -11,11 +11,11 @@
#
cell = "quad4"
-dx = 100
+dx = 200
dt = 0.05
-outputRoot = "output_pacha/%s_%3dm" % (cell,dx)
-outdir = "scecfiles/%s_%3dm/" % (cell,dx)
+outputRoot = "output/%s_%3dm_%s" % (cell,dx,"refine")
+outdir = "scecfiles/%s_%3dm_%s/" % (cell,dx,"refine")
import numpy
import time
@@ -23,7 +23,7 @@
from pylith.utils.VTKDataReader import VTKDataReader
# ----------------------------------------------------------------------
-timestamps = numpy.arange(50,12001,50)
+timestamps = numpy.arange(0,12001,50)
if dx == 200:
targets = numpy.array([[0.0, -12000.0, 0.0],
[0.0, -7600.0, 0.0],
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m_refine.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m_refine.cfg 2010-05-19 20:19:01 UTC (rev 16749)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_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.output]
+writer.filename = output/tri3_100m_refine.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tri3_100m_refine-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tri3_100m_refine-elastic.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m_refine.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m_refine.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m_refine.jou 2010-05-19 20:19:01 UTC (rev 16749)
@@ -0,0 +1,46 @@
+# ----------------------------------------------------------------------
+# Create tri3 mesh at 100m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry_refine.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size and scheme
+# ----------------------------------------------------------------------
+surface all size 400
+surface all scheme trimesh
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Refine elements within 1/2 of rupture width
+# ----------------------------------------------------------------------
+playback 'refinecells.jou'
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+surface all smooth scheme laplacian
+smooth surface all
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "tri3_100m_refine.exo" dimension 2 overwrite
+
+
+
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.cfg 2010-05-19 20:19:01 UTC (rev 16749)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_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.output]
+writer.filename = output/tri3_200m_refine.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tri3_200m_refine-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tri3_200m_refine-elastic.vtk
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.jou 2010-05-19 16:53:17 UTC (rev 16748)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.jou 2010-05-19 20:19:01 UTC (rev 16749)
@@ -21,10 +21,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
+playback 'refinecells.jou'
# ----------------------------------------------------------------------
# Smooth mesh to improve quality
More information about the CIG-COMMITS
mailing list