[cig-commits] r16761 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d
surendra at geodynamics.org
surendra at geodynamics.org
Sat May 22 00:19:58 PDT 2010
Author: surendra
Date: 2010-05-22 00:19:58 -0700 (Sat, 22 May 2010)
New Revision: 16761
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/refinecells.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.jou
Log:
Added REFINE meshes for TPV210-2d
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry_refine.jou 2010-05-22 07:19:58 UTC (rev 16761)
@@ -0,0 +1,44 @@
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+# Block is 144.0 km x 72 km
+# -72 km <= x <= 72 km
+# -72 km <= y <= 0 km
+#
+# We create the 2-D mesh by extracting the mid-surface from a 3-D mesh.
+reset
+brick x 144000 y 72000 z 72000
+
+# ----------------------------------------------------------------------
+# Set z=0 at free surface
+# ----------------------------------------------------------------------
+move volume 1 x -3750 y -36000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create surface extended from surface 5
+body 2 rotate +60 about z
+webcut volume 1 with plane surface 7
+
+delete surface 7
+
+# ----------------------------------------------------------------------
+# Extract the mid surfaces for 2-D
+# ----------------------------------------------------------------------
+create midsurface volume 3 surface 15 16
+create midsurface volume 1 surface 9 11
+
+delete volume 1
+delete volume 3
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all
+merge all
+
+
+
+
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.cfg 2010-05-22 07:19:58 UTC (rev 16761)
@@ -0,0 +1,29 @@
+# -*- 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/tpv210-2d/quad4_100m_refine.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.jou 2010-05-22 07:19:58 UTC (rev 16761)
@@ -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
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.cfg 2010-05-22 07:19:58 UTC (rev 16761)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = quad4_200m_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_200m_refine.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/quad4_200m_refine-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/quad4_200m_refine-elastic.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.jou 2010-05-22 07:19:58 UTC (rev 16761)
@@ -0,0 +1,39 @@
+# ----------------------------------------------------------------------
+# Create quad4 mesh at 200m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry_refine.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+surface all size 800
+
+# ----------------------------------------------------------------------
+# 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_200m_refine.exo" dimension 2 overwrite
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/refinecells.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/refinecells.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/refinecells.jou 2010-05-22 07:19:58 UTC (rev 16761)
@@ -0,0 +1,7 @@
+# ----------------------------------------------------------------------
+# Refine elements within 1/2 of rupture width
+# ----------------------------------------------------------------------
+group "refinecells1" add node with x_min >= -36.001e+3 and x_max <= 36.001e+3 and y_max <= 36.001e+3
+refine node in refinecells1 depth 1
+group "refinecells2" add node with x_min >= -18.001e+3 and x_max <= 18.001e+3 and y_max <= 24.001e+3
+refine node in refinecells2 depth 1
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.cfg 2010-05-22 07:19:58 UTC (rev 16761)
@@ -0,0 +1,23 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_100m_refine.exo
+
+# ----------------------------------------------------------------------
+# 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/tpv210-2d/tri3_100m_refine.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.jou 2010-05-22 07:19:58 UTC (rev 16761)
@@ -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/tpv210-2d/tri3_200m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.cfg 2010-05-22 07:19:58 UTC (rev 16761)
@@ -0,0 +1,23 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_200m_refine.exo
+
+# ----------------------------------------------------------------------
+# 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
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.jou 2010-05-22 07:19:58 UTC (rev 16761)
@@ -0,0 +1,46 @@
+# ----------------------------------------------------------------------
+# Create tri3 mesh at 200m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry_refine.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size and scheme
+# ----------------------------------------------------------------------
+surface all size 800
+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_200m_refine.exo" dimension 2 overwrite
+
+
+
+
More information about the CIG-COMMITS
mailing list