[cig-commits] r20019 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d
brad at geodynamics.org
brad at geodynamics.org
Tue May 1 10:20:05 PDT 2012
Author: brad
Date: 2012-05-01 10:20:05 -0700 (Tue, 01 May 2012)
New Revision: 20019
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/plot_surfvel.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/quad4_050m_gradient.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/quad4_100m_gradient.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/quad4_200m_gradient.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tabulate_offfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tabulate_onfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tri3.cfg
Log:
Small cleanup.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/plot_surfvel.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/plot_surfvel.py 2012-05-01 17:17:41 UTC (rev 20018)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/plot_surfvel.py 2012-05-01 17:20:05 UTC (rev 20019)
@@ -10,11 +10,11 @@
# ======================================================================
#
-sim = "tpv15_tri3_050m_gradient"
-tBegin = 2.0
-tEnd = 4.01
-dt = 0.1
-exaggeration = 0.05
+sim = "tpv14_tri3_100m_gradient"
+tBegin = 0.0
+tEnd = 6.01
+dt = 0.05
+exaggeration = 0.5
# ======================================================================
from enthought.mayavi.plugins.app import Mayavi
@@ -135,8 +135,8 @@
self.colorSurfTrace = (1,1,0)
self.lutReverse = False
- self.windowSize = (960, 960)
- self.aaframes = 0
+ self.windowSize = (960, 540)
+ self.aa_frames = 0
return
@@ -220,11 +220,11 @@
clipRange = numpy.array( [200, 5000] )
ptFrom = ptTo + numpy.array( [0.0, 0.0, dist])
- camera.view_up = (0,1,0)
+ camera.view_up = (-1,0,0)
camera.focal_point = ptTo
camera.position = ptFrom
camera.clipping_range = clipRange
- camera.parallel_scale = 0.5
+ camera.parallel_scale = 2.5
camera.parallel_projection = True
return
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/pylithapp.cfg 2012-05-01 17:17:41 UTC (rev 20018)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/pylithapp.cfg 2012-05-01 17:20:05 UTC (rev 20019)
@@ -27,7 +27,7 @@
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
reader = pylith.meshio.MeshIOCubit
-#reorder_mesh = True
+reorder_mesh = True
distributor.partitioner = parmetis
distributor.write_partition = True
distributor.data_writer.filename = partition.vtk
@@ -39,7 +39,6 @@
# problem
# ----------------------------------------------------------------------
[pylithapp.timedependent]
-
dimension = 2
formulation = pylith.problems.ExplicitLumped
@@ -54,23 +53,20 @@
# materials
# ----------------------------------------------------------------------
[pylithapp.timedependent]
-
materials = [elastic]
materials.elastic = pylith.materials.ElasticPlaneStrain
[pylithapp.timedependent.materials.elastic]
label = Elastic material
id = 1
+db_properties.label = Elastic properties
db_properties.iohandler.filename = matprops.spatialdb
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
[pylithapp.timedependent]
-
bc = [x_neg,x_pos,y_neg,y_pos]
-
-# Change the DirichletPoints BC to the AbsorbingDampers BC
bc.x_neg = pylith.bc.AbsorbingDampers
bc.x_pos = pylith.bc.AbsorbingDampers
bc.y_neg = pylith.bc.AbsorbingDampers
@@ -100,9 +96,7 @@
# faults
# ----------------------------------------------------------------------
[pylithapp.timedependent]
-
interfaces = [main_fault,branch_fault]
-# Change the kinematic rupture to spontaneous rupture
interfaces.main_fault = pylith.faults.FaultCohesiveDyn
interfaces.branch_fault = pylith.faults.FaultCohesiveDyn
@@ -119,6 +113,7 @@
db_initial_tractions.label = "Initial fault tractions"
friction = pylith.friction.SlipWeakening
+friction.label = Slip weakening friction on main fault
friction.db_properties = spatialdata.spatialdb.SimpleDB
friction.db_properties.iohandler.filename = friction_main.spatialdb
@@ -139,6 +134,7 @@
db_initial_tractions.data = [-70.0*MPa, -120.0*MPa]
friction = pylith.friction.SlipWeakening
+friction.label = Slip weakening friction on branch fault
friction.db_properties = spatialdata.spatialdb.SimpleDB
friction.db_properties.iohandler.filename = friction_branch.spatialdb
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/quad4_050m_gradient.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/quad4_050m_gradient.jou 2012-05-01 17:17:41 UTC (rev 20018)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/quad4_050m_gradient.jou 2012-05-01 17:20:05 UTC (rev 20019)
@@ -18,8 +18,12 @@
# ----------------------------------------------------------------------
surface all scheme pave
playback 'gradient.jou'
-curve 30 interval 26
+# Reset biases to get even number of points in loop required by paving scheme
+curve 18 scheme bias fine size {dxA} factor {0.998*bias_factor} start vertex 13
+curve 19 scheme bias fine size {dxA} factor {0.998*bias_factor} start vertex 14
+curve 20 scheme bias fine size {dxA} factor {0.998*bias_factor} start vertex 14
+
# ----------------------------------------------------------------------
# Generate the mesh
# ----------------------------------------------------------------------
@@ -41,3 +45,4 @@
# ----------------------------------------------------------------------
export mesh "quad4_050m_gradient.exo" dimension 2 overwrite
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/quad4_100m_gradient.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/quad4_100m_gradient.jou 2012-05-01 17:17:41 UTC (rev 20018)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/quad4_100m_gradient.jou 2012-05-01 17:20:05 UTC (rev 20019)
@@ -19,6 +19,10 @@
surface all scheme pave
playback 'gradient.jou'
+# Reset biases to get even number of points in loop required by paving scheme
+curve 27 scheme bias fine size {dxA} factor {0.997*bias_factor} start vertex 13
+
+
# ----------------------------------------------------------------------
# Generate the mesh
# ----------------------------------------------------------------------
@@ -41,3 +45,4 @@
export mesh "quad4_100m_gradient.exo" dimension 2 overwrite
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/quad4_200m_gradient.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/quad4_200m_gradient.jou 2012-05-01 17:17:41 UTC (rev 20018)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/quad4_200m_gradient.jou 2012-05-01 17:20:05 UTC (rev 20019)
@@ -18,8 +18,11 @@
# ----------------------------------------------------------------------
surface all scheme pave
playback 'gradient.jou'
-curve 27 interval 46
+# Reset biases to get even number of points in loop required by paving scheme
+curve 30 interval 12
+curve 14 scheme bias fine size {dxD} coarse size {0.98*dxC} start vertex 11
+
# ----------------------------------------------------------------------
# Generate the mesh
# ----------------------------------------------------------------------
@@ -43,3 +46,4 @@
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tabulate_offfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tabulate_offfaultdata.py 2012-05-01 17:17:41 UTC (rev 20018)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tabulate_offfaultdata.py 2012-05-01 17:20:05 UTC (rev 20019)
@@ -10,9 +10,9 @@
# ----------------------------------------------------------------------
#
-sim = "tpv15_alt"
+sim = "tpv15"
cell = "tri3"
-dx = 50
+dx = 200
dt = 0.05
inputRoot = "output/%s_%s_%03dm_gradient" % (sim, cell,dx)
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tabulate_onfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tabulate_onfaultdata.py 2012-05-01 17:17:41 UTC (rev 20018)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tabulate_onfaultdata.py 2012-05-01 17:20:05 UTC (rev 20019)
@@ -10,9 +10,9 @@
# ----------------------------------------------------------------------
#
-sim = "tpv15_alt"
+sim = "tpv15"
cell = "tri3"
-dx = 50
+dx = 200
dt = 0.05
inputRoot = "output/%s_%s_%03dm_gradient" % (sim, cell,dx)
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tri3.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tri3.cfg 2012-05-01 17:17:41 UTC (rev 20018)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tri3.cfg 2012-05-01 17:20:05 UTC (rev 20019)
@@ -12,35 +12,35 @@
# ----------------------------------------------------------------------
[pylithapp.timedependent.materials.elastic]
quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
+quadrature.cell.dimension = 2
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
[pylithapp.timedependent.bc.x_pos]
quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
+quadrature.cell.dimension = 1
[pylithapp.timedependent.bc.x_neg]
quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
+quadrature.cell.dimension = 1
[pylithapp.timedependent.bc.y_pos]
quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
+quadrature.cell.dimension = 1
[pylithapp.timedependent.bc.y_neg]
quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
+quadrature.cell.dimension = 1
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
[pylithapp.timedependent.interfaces.main_fault]
quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
+quadrature.cell.dimension = 1
[pylithapp.timedependent.interfaces.branch_fault]
quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
+quadrature.cell.dimension = 1
More information about the CIG-COMMITS
mailing list