[cig-commits] r16695 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d
surendra at geodynamics.org
surendra at geodynamics.org
Wed May 12 18:46:04 PDT 2010
Author: surendra
Date: 2010-05-12 18:46:03 -0700 (Wed, 12 May 2010)
New Revision: 16695
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_200m.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_200m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_200m.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_200m.jou
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_100m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tabulate_onfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_100m.cfg
Log:
Added files to run test tpv210-2d at 200m resolution. Moved reference cell and quadrature information to separate files one each for quad4 and tri3. Need to rename this directory to tpv210-2d.
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4.cfg 2010-05-13 01:46:03 UTC (rev 16695)
@@ -0,0 +1,42 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.elastic]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_neg]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_100m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_100m.cfg 2010-05-12 19:29:28 UTC (rev 16694)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_100m.cfg 2010-05-13 01:46:03 UTC (rev 16695)
@@ -14,45 +14,6 @@
dt = 0.01*s
# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials.elastic]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.x_neg]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.y_neg]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-
-# Provide information on the fault (interface).
-[pylithapp.timedependent.interfaces.fault]
-# Specify the basis functions and quadrature:
-# 1-D Lagrange cell in 2-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_200m.cfg 2010-05-13 01:46:03 UTC (rev 16695)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = quad4_200m.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.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/quad4_200m-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/quad4_200m-elastic.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_200m.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_200m.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_200m.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/quad4_200m.jou 2010-05-13 01:46:03 UTC (rev 16695)
@@ -0,0 +1,34 @@
+# ----------------------------------------------------------------------
+# Create quad4 mesh at 200m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set meshing scheme
+# ----------------------------------------------------------------------
+surface 18 19 scheme pave
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+surface all size 200
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "quad4_200m.exo" dimension 2 overwrite
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tabulate_onfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tabulate_onfaultdata.py 2010-05-12 19:29:28 UTC (rev 16694)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tabulate_onfaultdata.py 2010-05-13 01:46:03 UTC (rev 16695)
@@ -10,7 +10,7 @@
# ----------------------------------------------------------------------
#
-cell = "quad4"
+cell = "tri3"
dx = 100
dt = 0.05
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3.cfg 2010-05-13 01:46:03 UTC (rev 16695)
@@ -0,0 +1,36 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.elastic]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+
+[pylithapp.timedependent.bc.x_neg]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+
+[pylithapp.timedependent.bc.y_pos]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+
+[pylithapp.timedependent.bc.y_neg]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_100m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_100m.cfg 2010-05-12 19:29:28 UTC (rev 16694)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_100m.cfg 2010-05-13 01:46:03 UTC (rev 16695)
@@ -8,42 +8,6 @@
filename = tri3_100m.exo
# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.formulation.time_step]
-dt = 0.01*s
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials.elastic]
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
-
-[pylithapp.timedependent.bc.x_neg]
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
-
-[pylithapp.timedependent.bc.y_neg]
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault]
-# Specify the basis functions and quadrature:
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
-
-# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_200m.cfg 2010-05-13 01:46:03 UTC (rev 16695)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_200m.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_200m.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tri3_200m-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tri3_200m-elastic.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_200m.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_200m.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_200m.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv10-2d/tri3_200m.jou 2010-05-13 01:46:03 UTC (rev 16695)
@@ -0,0 +1,32 @@
+# ----------------------------------------------------------------------
+# Create tri3 mesh at 200m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size and scheme
+# ----------------------------------------------------------------------
+surface all size 200
+surface all scheme trimesh
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "tri3_200m.exo" dimension 2 overwrite
+
+
+
More information about the CIG-COMMITS
mailing list