[cig-commits] r19833 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210

brad at geodynamics.org brad at geodynamics.org
Wed Mar 21 09:54:04 PDT 2012


Author: brad
Date: 2012-03-21 09:54:04 -0700 (Wed, 21 Mar 2012)
New Revision: 19833

Removed:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8.cfg
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8_200m.cfg
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8_200m.jou
Log:
Removed hex8 stuff. Too many restrictions imposed by problem statement to easily generate hex8 mesh (discontinuity in stress field and apserity).

Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8.cfg	2012-03-21 16:19:41 UTC (rev 19832)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8.cfg	2012-03-21 16:54:04 UTC (rev 19833)
@@ -1,47 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials.elastic.quadrature]
-cell = pylith.feassemble.FIATLagrange
-cell.dimension = 3
-cell.quad_order = 2
-
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos.quadrature]
-cell = pylith.feassemble.FIATLagrange
-cell.dimension = 2
-cell.quad_order = 2
-
-[pylithapp.timedependent.bc.x_neg.quadrature]
-cell = pylith.feassemble.FIATLagrange
-cell.dimension = 2
-cell.quad_order = 2
-
-[pylithapp.timedependent.bc.y_pos.quadrature]
-cell = pylith.feassemble.FIATLagrange
-cell.dimension = 2
-cell.quad_order = 2
-
-[pylithapp.timedependent.bc.y_neg.quadrature]
-cell = pylith.feassemble.FIATLagrange
-cell.dimension = 2
-cell.quad_order = 2
-
-[pylithapp.timedependent.bc.z_neg.quadrature]
-cell = pylith.feassemble.FIATLagrange
-cell.dimension = 2
-cell.quad_order = 2
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault.quadrature]
-cell = pylith.feassemble.FIATLagrange
-cell.dimension = 2
-cell.quad_order = 2

Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8_200m.cfg	2012-03-21 16:19:41 UTC (rev 19832)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8_200m.cfg	2012-03-21 16:54:04 UTC (rev 19833)
@@ -1,34 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator.reader]
-filename = hex8_200m.exo
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.formulation.time_step]
-dt = 0.01666666666666666666666667*s
-
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Ground surface
-[pylithapp.problem.formulation.output.groundsurf]
-writer.filename = output/hex8_200m-groundsurf.vtk
-
-# Off-fault
-[pylithapp.problem.formulation.output.offfault]
-writer.filename = output/hex8_200m-offfault.vtk
-
-# Fault
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/hex8_200m-fault.vtk
-
-# State variables
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/hex8_200m-elastic.vtk

Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8_200m.jou	2012-03-21 16:19:41 UTC (rev 19832)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8_200m.jou	2012-03-21 16:54:04 UTC (rev 19833)
@@ -1,64 +0,0 @@
-# -*- Python -*- (syntax highlighting)
-# ----------------------------------------------------------------------
-#
-# Brad T. Aagaard, U.S. Geological Survey
-#
-# ----------------------------------------------------------------------
-# Create hex8 mesh at 200m resolution.
-# ----------------------------------------------------------------------
-
-# ----------------------------------------------------------------------
-# Generate geometry
-# ----------------------------------------------------------------------
-playback 'geometry.jou'
-
-# ----------------------------------------------------------------------
-# Set discretization size and scheme
-# ----------------------------------------------------------------------
-#{dx=200*m}
-surface all scheme map
-surface ?? scheme pave
-playback 'gradient.jou'
-
-# ----------------------------------------------------------------------
-# Generate the mesh
-# ----------------------------------------------------------------------
-mesh surface all
-surface all smooth scheme condition number beta 1.3 cpu 10
-smooth surface all
-
-mesh volume all
-
-# ----------------------------------------------------------------------
-# Refine around fault
-# ----------------------------------------------------------------------
-#{refineWidth=6.001*km}
-#{refineLength=42.001*km}
-#{refineDepth=21.001*km}
-group "refinecells" add Hex with x_min >= {-0.5*refineWidth} and x_max <= {+0.5*refineWidth} and y_min >= {-0.5*refineLength} and y_max <= {+0.5*refineLength} and z_min >= {-refineDepth}
-refine hex in group refinecells depth 1
-delete group refinecells
-
-# ----------------------------------------------------------------------
-# Smooth to improve quality
-# ----------------------------------------------------------------------
-#{condnum=3.0}
-#{loop(5)}
-volume all smooth scheme condition number beta {condnum} cpu 10
-smooth volume all
-#{condnum=condnum-0.2}
-#{endloop}
-
-
-# ----------------------------------------------------------------------
-# Setup boundary conditions.
-# ----------------------------------------------------------------------
-playback 'createbc.jou'
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "hex8_200m.exo" dimension 3 overwrite
-
-
-# End of file



More information about the CIG-COMMITS mailing list