[cig-commits] r16666 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d
surendra at geodynamics.org
surendra at geodynamics.org
Fri May 7 13:52:54 PDT 2010
Author: surendra
Date: 2010-05-07 13:52:54 -0700 (Fri, 07 May 2010)
New Revision: 16666
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m.jou
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/friction.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_offfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_onfaultdata.py
Log:
Added files for tri3 tests at 100m and 200m resolution. Moved reference cell and quadrature information from pylithapp.cfg to so as to avoid clash between tri3 and quad4 meshes
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/createbc.jou 2010-05-06 23:47:13 UTC (rev 16665)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/createbc.jou 2010-05-07 20:52:54 UTC (rev 16666)
@@ -2,12 +2,15 @@
# Create blocks for materials
# ----------------------------------------------------------------------
block 1 surface 8 9
+#block 1 surface 11 12 13 14
block 1 name "elastic"
# ----------------------------------------------------------------------
# Create nodeset for fault
# ----------------------------------------------------------------------
group "fault" add node in fault_surface
+#group "fault" add node in fault_surface at A
+#group "fault" add node in fault_surface at B
group "fault" remove node with y_coord > +15.0001e+3
group "fault" remove node with y_coord < -15.0001e+3
nodeset 10 group fault
@@ -42,3 +45,13 @@
group "face_yneg" add node in curve 20
nodeset 14 group face_yneg
nodeset 14 name "face yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for Off-Fault stations
+# ----------------------------------------------------------------------
+#group "off_fault_stat" add node in vertex 15
+#group "off_fault_stat" add node in vertex 16
+#nodeset 15 group off_fault_stat
+#nodeset 15 name "off_fault_stat"
+
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/friction.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/friction.spatialdb 2010-05-06 23:47:13 UTC (rev 16665)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/friction.spatialdb 2010-05-07 20:52:54 UTC (rev 16666)
@@ -11,7 +11,7 @@
space-dim = 2
}
}
-0.0 -15.1 1000.0 0.525 0.4 0.0
-0.0 -14.9 0.677 0.525 0.4 0.0
-0.0 +14.9 0.677 0.525 0.4 0.0
-0.0 +15.1 1000.0 0.525 0.4 0.0
+0.0 -15.01 10000.0 0.525 0.4 0.0
+0.0 -14.99 0.677 0.525 0.4 0.0
+0.0 +14.99 0.677 0.525 0.4 0.0
+0.0 +15.01 10000.0 0.525 0.4 0.0
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry.jou 2010-05-06 23:47:13 UTC (rev 16665)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry.jou 2010-05-07 20:52:54 UTC (rev 16666)
@@ -7,7 +7,7 @@
#
# We create the 2-D mesh by extracting the mid-surface from a 3-D mesh.
reset
-brick x 24000 y 48000 z 24000
+brick x 74000 y 48000 z 24000
create midsurface volume 1 surface 1 2
delete volume 1
@@ -20,7 +20,20 @@
curve 17 name "fault_surface"
# ----------------------------------------------------------------------
+# Create surface for Off-Fault stations
+# ----------------------------------------------------------------------
+
+#create vertex 3000 12000 0
+#create vertex 3000 -12000 0
+#create vertex -3000 -12000 0
+#create vertex -3000 12000 0
+
+#create surface vertex 15 16 17 18
+
+# ----------------------------------------------------------------------
# Imprint all volumes, merging surfaces
# ----------------------------------------------------------------------
imprint all
merge all
+
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/pylithapp.cfg 2010-05-06 23:47:13 UTC (rev 16665)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/pylithapp.cfg 2010-05-07 20:52:54 UTC (rev 16666)
@@ -67,6 +67,10 @@
dt = 0.01*s
+# Set subdomain component to OutputSolnSubset (subset of domain).
+#output.subdomain = pylith.meshio.OutputSolnSubset
+
+
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
@@ -79,10 +83,7 @@
label = Elastic material
id = 1
db_properties.iohandler.filename = matprops.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
@@ -90,35 +91,22 @@
label = 11
db.label = Absorbing BC +x
db.iohandler.filename = matprops.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
[pylithapp.timedependent.bc.x_neg]
label = 12
db.label = Absorbing BC -x
db.iohandler.filename = matprops.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
[pylithapp.timedependent.bc.y_pos]
label = 13
db.label = Absorbing BC +y
db.iohandler.filename = matprops.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
[pylithapp.timedependent.bc.y_neg]
label = 14
db.label = Absorbing BC -y
db.iohandler.filename = matprops.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
@@ -136,12 +124,6 @@
# Nodeset id from CUBIT
label = 10
-# 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
-
db_initial_tractions = spatialdata.spatialdb.SimpleDB
db_initial_tractions.iohandler.filename = tractions.spatialdb
db_initial_tractions.label = "Initial fault tractions"
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4.cfg 2010-05-07 20:52:54 UTC (rev 16666)
@@ -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/tpv205-2d/quad4_100m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m.cfg 2010-05-06 23:47:13 UTC (rev 16665)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_100m.cfg 2010-05-07 20:52:54 UTC (rev 16666)
@@ -13,7 +13,6 @@
[pylithapp.timedependent.formulation.time_step]
dt = 0.01*s
-
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.cfg 2010-05-06 23:47:13 UTC (rev 16665)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.cfg 2010-05-07 20:52:54 UTC (rev 16666)
@@ -21,6 +21,15 @@
[pylithapp.problem.formulation.output.output]
writer.filename = output/quad4_200m.vtk
+# Give basename for VTK domain output of solution over Off-fault stations.
+#[pylithapp.problem.formulation.output.subdomain]
+#label = 15 ; nodeset for subdomain
+#vertex_data_fields=[displacement,velocity]
+#output_freq = time_step
+#time_step = 0.04999*s
+#writer.filename = output/quad4_200m_OffFault.vtk
+#writer.time_format = %05.2f
+
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
writer.filename = output/quad4_200m-fault.vtk
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_offfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_offfaultdata.py 2010-05-06 23:47:13 UTC (rev 16665)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_offfaultdata.py 2010-05-07 20:52:54 UTC (rev 16666)
@@ -10,7 +10,7 @@
# ----------------------------------------------------------------------
#
-dx = 200
+dx = 100
dt = 0.05
outputRoot = "output/quad4_%3dm" % dx
@@ -110,10 +110,10 @@
dip))
fout.write(headerB)
data = numpy.transpose((time,
- +disp[:,iloc,1],
- +vel[:,iloc,1],
- -disp[:,iloc,2],
- -vel[:,iloc,2],
+ -disp[:,iloc,1],
+ -vel[:,iloc,1],
+ +disp[:,iloc,2],
+ +vel[:,iloc,2],
-disp[:,iloc,0],
-vel[:,iloc,0]))
numpy.savetxt(fout, data, fmt='%14.6e')
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_onfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_onfaultdata.py 2010-05-06 23:47:13 UTC (rev 16665)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_onfaultdata.py 2010-05-07 20:52:54 UTC (rev 16666)
@@ -92,10 +92,10 @@
headerA = \
"# problem = TPV205\n" + \
- "# author = BradAagaard\n" + \
+ "# author = Surendra N. Somala\n" + \
"# date = %s\n" % (time.asctime()) + \
"# code = PyLith\n" + \
- "# code_version = 1.4.3\n" + \
+ "# code_version = 1.5.0a\n" + \
"# element_size = %s\n" % dx
headerB = \
"# Time series in 7 columns of E14.6:\n" + \
@@ -118,7 +118,7 @@
timeScale = 100.0
dip = 7.5
strike = y_OnFltStat[:] / lengthScale
-time = (nfileNames[:] - 1) / timeScale
+time = nfileNames[:] / timeScale
print "time", time
# Write data
@@ -137,8 +137,8 @@
-slip[:,iloc,0],
-slip_rate[:,iloc,0],
-traction[:,iloc,0]/1e+6,
- +slip[:,iloc,1],
- +slip_rate[:,iloc,1],
- +traction[:,iloc,1]/1e+6))
+ +slip[:,iloc,2],
+ +slip_rate[:,iloc,2],
+ +traction[:,iloc,2]/1e+6))
numpy.savetxt(fout, data, fmt='%14.6e')
fout.close()
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3.cfg 2010-05-07 20:52:54 UTC (rev 16666)
@@ -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
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m.cfg 2010-05-07 20:52:54 UTC (rev 16666)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_100m.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_100m.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tri3_100m-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tri3_100m-elastic.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_100m.jou 2010-05-07 20:52:54 UTC (rev 16666)
@@ -0,0 +1,32 @@
+# ----------------------------------------------------------------------
+# Create tri3 mesh at 100m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size and scheme
+# ----------------------------------------------------------------------
+surface all size 100
+surface all scheme trimesh
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "tri3_100m.exo" dimension 2 overwrite
+
+
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m.cfg 2010-05-07 20:52:54 UTC (rev 16666)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_200m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01666666666666666666666667*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/tpv205-2d/tri3_200m.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m.jou 2010-05-07 20:52:54 UTC (rev 16666)
@@ -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