[cig-commits] r17999 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d
brad at geodynamics.org
brad at geodynamics.org
Tue Mar 1 13:14:02 PST 2011
Author: brad
Date: 2011-03-01 13:14:02 -0800 (Tue, 01 Mar 2011)
New Revision: 17999
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/friction_mainbranch.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/friction_mainext.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tractions.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tri3_050m_gradient.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tri3_100m_gradient.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tri3_200m_gradient.cfg
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tabulate_offfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tabulate_onfaultdata.py
Log:
Added alternative geometry.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/createbc.jou 2011-03-01 20:59:44 UTC (rev 17998)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/createbc.jou 2011-03-01 21:14:02 UTC (rev 17999)
@@ -9,8 +9,6 @@
# ----------------------------------------------------------------------
group "main_fault" add node in fault_surface_main
group "main_fault" add node in fault_surface_main at A
-#group "main_fault" remove node with y_coord > {mainFaultLength/2.0}
-#group "main_fault" remove node with y_coord < {-mainFaultLength/2.0}
nodeset 20 group main_fault
nodeset 20 name "main_fault"
@@ -18,12 +16,27 @@
# Create nodeset for fault branch
# ----------------------------------------------------------------------
group "branch_fault" add node in fault_surface_branch
-#group "branch_fault" remove node with x_coord > {branchFaultLength*sind(branchAngle)}
group "branch_fault" remove node in main_fault
nodeset 21 group branch_fault
nodeset 21 name "branch_fault"
# ----------------------------------------------------------------------
+# Create nodeset for main-branch fault (alternative geometry)
+# ----------------------------------------------------------------------
+group "mainbranch_fault" add node in fault_surface_main at A
+group "mainbranch_fault" add node in fault_surface_branch
+nodeset 22 group mainbranch_fault
+nodeset 22 name "mainbranch_fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for main-ext fault (alternative geometry)
+# ----------------------------------------------------------------------
+group "mainext_fault" add node in fault_surface_main
+group "mainext_fault" remove node in mainbranch_fault
+nodeset 23 group mainext_fault
+nodeset 23 name "mainext_fault"
+
+# ----------------------------------------------------------------------
# Create nodeset for +x face
# ----------------------------------------------------------------------
group "face_xpos" add node in curve 14
@@ -55,3 +68,4 @@
nodeset 13 name "face_yneg"
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/friction_mainbranch.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/friction_mainbranch.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/friction_mainbranch.spatialdb 2011-03-01 21:14:02 UTC (rev 17999)
@@ -0,0 +1,17 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = static-coefficient dynamic-coefficient slip-weakening-parameter cohesion
+ value-units = none none m MPa
+ num-locs = 4
+ data-dim = 1
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y are in km
+ space-dim = 2
+ }
+}
+0.0 -13.999 10000.0 0.525 0.4 0.0
+0.0 -13.998 0.677 0.525 0.4 0.0
+0.0 +12.391 0.677 0.525 0.4 0.0
+0.0 +12.392 10000.0 0.525 0.4 0.0
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/friction_mainext.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/friction_mainext.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/friction_mainext.spatialdb 2011-03-01 21:14:02 UTC (rev 17999)
@@ -0,0 +1,17 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = static-coefficient dynamic-coefficient slip-weakening-parameter cohesion
+ value-units = none none m MPa
+ num-locs = 4
+ data-dim = 1
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y are in km
+ space-dim = 2
+ }
+}
+0.0 +2.10 10000.0 0.525 0.4 0.0
+0.0 +2.11 0.677 0.525 0.4 0.0
+0.0 +13.98 0.677 0.525 0.4 0.0
+0.0 +13.99 10000.0 0.525 0.4 0.0
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/pylithapp.cfg 2011-03-01 20:59:44 UTC (rev 17998)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/pylithapp.cfg 2011-03-01 21:14:02 UTC (rev 17999)
@@ -28,6 +28,9 @@
[pylithapp.mesh_generator]
reader = pylith.meshio.MeshIOCubit
#reorder_mesh = True
+distributor.partitioner = parmetis
+distributor.write_partition = True
+distributor.data_writer.filename = partition.vtk
[pylithapp.mesh_generator.reader]
coordsys.space_dim = 2
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 2011-03-01 20:59:44 UTC (rev 17998)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tabulate_offfaultdata.py 2011-03-01 21:14:02 UTC (rev 17999)
@@ -10,7 +10,7 @@
# ----------------------------------------------------------------------
#
-sim = "tpv14"
+sim = "tpv15_alt"
cell = "tri3"
dx = 50
dt = 0.05
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 2011-03-01 20:59:44 UTC (rev 17998)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tabulate_onfaultdata.py 2011-03-01 21:14:02 UTC (rev 17999)
@@ -10,7 +10,7 @@
# ----------------------------------------------------------------------
#
-sim = "tpv14"
+sim = "tpv15_alt"
cell = "tri3"
dx = 50
dt = 0.05
@@ -89,14 +89,31 @@
strike = (targets[:,1] - 2000) / lengthScale
elif fault == "branch_fault":
strike = 2.0 * targets[:,0] / lengthScale
+ elif fault == "mainbranch_fault":
+ maskM = targets[:,0] <= +10.0
+ strike = maskM*(targets[:,1] - 2000) / lengthScale + \
+ ~maskM*(2.0 * targets[:,0]) / lengthScale
+ elif fault == "mainext_fault":
+ strike = (targets[:,1] - 2000) / lengthScale
+ else:
+ raise ValueError("Unknown fault '%s'" % fault)
print timeStamps.shape, ntimesteps, slip.shape
zero = numpy.zeros( (ntimesteps-1, 1), dtype=numpy.float64)
for iloc in xrange(ntargets):
+ if fault == "mainbranch_fault":
+ if targets[iloc,0] <= 10.0:
+ faultName = "main_fault"
+ else:
+ faultName = "branch_fault"
+ elif fault == "mainext_fault":
+ faultName = "main_fault"
+ else:
+ faultName = fault;
pt = locName % (round(10*strike[iloc]),
round(10*dip))
- filename = "%s-%s-%s.dat" % (outputRoot, fault, pt)
+ filename = "%s-%s-%s.dat" % (outputRoot, faultName, pt)
fout = open(filename, 'w');
fout.write(headerA)
fout.write("# time_step = %14.6E\n" % dt)
@@ -116,30 +133,57 @@
return
-# ----------------------------------------------------------------------
-# MAIN FAULT
+if sim in ["tpv14", "tpv15"]:
+ # ------------------------------------------------------------------
+ # MAIN FAULT
-# Target coordinates are relative to faults intersection.
-targets = numpy.array([[0.0, -2000.0],
- [0.0, +2000.0],
- [0.0, +5000.0],
- [0.0, +9000.0]])
-# Origin of coordinate system is at center of main fault
-targets[:,1] += 2000.0
+ # Target coordinates are relative to faults intersection.
+ targets = numpy.array([[0.0, -2000.0],
+ [0.0, +2000.0],
+ [0.0, +5000.0],
+ [0.0, +9000.0]])
+ # Origin of coordinate system is at center of main fault
+ targets[:,1] += 2000.0
+
+ extract("main_fault", targets)
+
+ # ------------------------------------------------------------------
+ # BRANCH FAULT
-extract("main_fault", targets)
+ # Target coordinates are relative to faults intersection.
+ targets = numpy.array([[1000.0, 1732.0508],
+ [2500.0, 4330.1270],
+ [4500.0, 7794.2286]])
+ # Origin of coordinate system is at center of main fault
+ targets[:,1] += 2000.0
+
+ extract("branch_fault", targets)
-# ----------------------------------------------------------------------
-# BRANCH FAULT
+else:
+ # ------------------------------------------------------------------
+ # MAIN-BRANCH FAULT
-# Target coordinates are relative to faults intersection.
-targets = numpy.array([[1000.0, 1732.0508],
- [2500.0, 4330.1270],
- [4500.0, 7794.2286]])
-# Origin of coordinate system is at center of main fault
-targets[:,1] += 2000.0
+ # Target coordinates are relative to faults intersection.
+ targets = numpy.array([[0.0, -2000.0],
+ [1000.0, 1732.0508],
+ [2500.0, 4330.1270],
+ [4500.0, 7794.2286]])
+ # Origin of coordinate system is at center of main fault
+ targets[:,1] += 2000.0
-extract("branch_fault", targets)
+ extract("mainbranch_fault", targets)
+ # ------------------------------------------------------------------
+ # MAIN-EXT FAULT
+
+ # Target coordinates are relative to faults intersection.
+ targets = numpy.array([[0.0, +2000.0],
+ [0.0, +5000.0],
+ [0.0, +9000.0]])
+ # Origin of coordinate system is at center of main fault
+ targets[:,1] += 2000.0
+
+ extract("mainext_fault", targets)
+
# End of file
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt.cfg 2011-03-01 21:14:02 UTC (rev 17999)
@@ -0,0 +1,21 @@
+# Parameters specific to TPV15-2D.
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+# Main fault
+[pylithapp.timedependent.interfaces.main_fault]
+label = mainbranch_fault
+
+friction.db_properties.iohandler.filename = friction_mainbranch.spatialdb
+db_initial_tractions.iohandler.filename = tpv15_alt_tractions.spatialdb
+
+# Branch fault
+[pylithapp.timedependent.interfaces.branch_fault]
+label = mainext_fault
+
+friction.db_properties.iohandler.filename = friction_mainext.spatialdb
+db_initial_tractions.data = [+70.0*MPa, -120.0*MPa]
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tractions.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tractions.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tractions.spatialdb 2011-03-01 21:14:02 UTC (rev 17999)
@@ -0,0 +1,19 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 2
+ value-names = traction-shear traction-normal
+ value-units = MPa MPa
+ num-locs = 6
+ data-dim = 2
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y, z are in km
+ space-dim = 2
+ }
+}
+0.0 -7.501 +70.0 -120.0
+0.0 -7.500 +81.6 -120.0
+0.0 -4.500 +81.6 -120.0
+0.0 -4.499 +70.0 -120.0
+0.0 +2.000 +70.0 -120.0
+0.0 +2.001 +78.0 -120.0
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tri3_050m_gradient.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tri3_050m_gradient.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tri3_050m_gradient.cfg 2011-03-01 21:14:02 UTC (rev 17999)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_050m_gradient.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.0025*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv15_alt_tri3_050m_gradient.h5
+
+[pylithapp.timedependent.interfaces.main_fault.output]
+writer.filename = output/tpv15_alt_tri3_050m_gradient-mainbranch_fault.h5
+
+[pylithapp.timedependent.interfaces.branch_fault.output]
+writer.filename = output/tpv15_alt_tri3_050m_gradient-mainext_fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv15_alt_tri3_050m_gradient-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tri3_100m_gradient.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tri3_100m_gradient.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tri3_100m_gradient.cfg 2011-03-01 21:14:02 UTC (rev 17999)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_100m_gradient.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.005*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv15_alt_tri3_100m_gradient.h5
+
+[pylithapp.timedependent.interfaces.main_fault.output]
+writer.filename = output/tpv15_alt_tri3_100m_gradient-mainbranch_fault.h5
+
+[pylithapp.timedependent.interfaces.branch_fault.output]
+writer.filename = output/tpv15_alt_tri3_100m_gradient-mainext_fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv15_alt_tri3_100m_gradient-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tri3_200m_gradient.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tri3_200m_gradient.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14-2d/tpv15_alt_tri3_200m_gradient.cfg 2011-03-01 21:14:02 UTC (rev 17999)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_200m_gradient.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv15_alt_tri3_200m_gradient.h5
+
+[pylithapp.timedependent.interfaces.main_fault.output]
+writer.filename = output/tpv15_alt_tri3_200m_gradient-mainbranch_fault.h5
+
+[pylithapp.timedependent.interfaces.branch_fault.output]
+writer.filename = output/tpv15_alt_tri3_200m_gradient-mainext_fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv15_alt_tri3_200m_gradient-elastic.h5
More information about the CIG-COMMITS
mailing list