[cig-commits] r16854 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210
surendra at geodynamics.org
surendra at geodynamics.org
Tue Jun 1 14:37:20 PDT 2010
Author: surendra
Date: 2010-06-01 14:37:19 -0700 (Tue, 01 Jun 2010)
New Revision: 16854
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/geometry.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_offfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_onfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_ruptime.py
Log:
Adjusted journal files and post-processing files to have Z-axis up. .cfg files were already setup that way so they need noadjustment
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/createbc.jou 2010-06-01 21:36:58 UTC (rev 16853)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/createbc.jou 2010-06-01 21:37:19 UTC (rev 16854)
@@ -9,8 +9,8 @@
# ----------------------------------------------------------------------
group "fault" add node in fault_surface
group "fault" remove node with x_coord < -7.5001e+3
-group "fault" remove node with z_coord < -15.0001e+3
-group "fault" remove node with z_coord > +15.0001e+3
+group "fault" remove node with y_coord < -15.0001e+3
+group "fault" remove node with y_coord > +15.0001e+3
nodeset 10 group fault
nodeset 10 name "fault"
@@ -31,16 +31,16 @@
# ----------------------------------------------------------------------
# Create nodeset for +y face
# ----------------------------------------------------------------------
-group "face_ypos" add node in surface 12
-group "face_ypos" add node in surface 14
+group "face_ypos" add node in surface 10
+group "face_ypos" add node in surface 17
nodeset 13 group face_ypos
nodeset 13 name "face_ypos"
# ----------------------------------------------------------------------
# Create nodeset for -y face
# ----------------------------------------------------------------------
-group "face_yneg" add node in surface 10
-group "face_yneg" add node in surface 17
+group "face_yneg" add node in surface 12
+group "face_yneg" add node in surface 14
nodeset 14 group face_yneg
nodeset 14 name "face_yneg"
@@ -63,8 +63,6 @@
# ----------------------------------------------------------------------
# Create nodeset for off-fault stations
# ----------------------------------------------------------------------
-group "off_fault" add node with x_min >= -4.001e+3 and x_max <= 4.001e+3 and y_min >= -1.001e+3 and z_coord = 0.0
+group "off_fault" add node with x_min >= -4.001e+3 and x_max <= 4.001e+3 and z_min >= -1.001e+3 and y_coord = 0.0
nodeset 25 group off_fault
-nodeset 25 name "off_fault"
-
-
+nodeset 25 name "off_fault"
\ No newline at end of file
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/geometry.jou 2010-06-01 21:36:58 UTC (rev 16853)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/geometry.jou 2010-06-01 21:37:19 UTC (rev 16854)
@@ -1,24 +1,24 @@
# ----------------------------------------------------------------------
# Create block
# ----------------------------------------------------------------------
-# Block is 36.0 km x 18.0*sqrt(3) km x 56000
+# Block is 36.0 km x 56000 km x 18.0*sqrt(3) km
# -21.75 km <= x <= 14.25 km
-# -18*Sqrt(3) km <= y <= 0 km
-# -28 km <= z <= 28km
+# -28 km <= y <= 28km
+# -18*Sqrt(3) km <= z <= 0 km
#
reset
-brick x 36000 y 31176.914 z 56000
+brick x 36000 y 56000 z 31176.914
# ----------------------------------------------------------------------
# Set z=0 at free surface
# ----------------------------------------------------------------------
-move volume 1 x -3750 y -15588.457
+move volume 1 x -3750 z -15588.457
# ----------------------------------------------------------------------
# Create interface surfaces
# ----------------------------------------------------------------------
-create surface extended from surface 5
-body 2 rotate +60 about z
+create surface extended from surface 1
+body 2 rotate -60 about y
surface 7 name "fault_surface"
webcut volume 1 with plane surface fault_surface
@@ -29,8 +29,3 @@
merge all
delete body 2
-
-
-
-
-
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_offfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_offfaultdata.py 2010-06-01 21:36:58 UTC (rev 16853)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_offfaultdata.py 2010-06-01 21:37:19 UTC (rev 16854)
@@ -29,18 +29,18 @@
# ----------------------------------------------------------------------
timestamps = numpy.arange(50,15001,50)
-targets = numpy.array([[-3000.0, -0.0, 0.0],
- [-2000.0, -0.0, 0.0],
- [-1000.0, -0.0, 0.0],
- [+1000.0, -0.0, 0.0],
- [+2000.0, -0.0, 0.0],
- [+3000.0, -0.0, 0.0],
- [-1000.0,-300.0, 0.0],
- [ -500.0,-300.0, 0.0],
- [ +500.0,-300.0, 0.0],
- [+1000.0,-300.0, 0.0],
- [-3000.0, -0.0, 12000.0],
- [+3000.0, -0.0, 12000.0]])
+targets = numpy.array([[-3000.0, 0.0, -0.0],
+ [-2000.0, 0.0, -0.0],
+ [-1000.0, 0.0, -0.0],
+ [+1000.0, 0.0, -0.0],
+ [+2000.0, 0.0, -0.0],
+ [+3000.0, 0.0, -0.0],
+ [-1000.0, 0.0,-300.0],
+ [ -500.0, 0.0,-300.0],
+ [ +500.0, 0.0,-300.0],
+ [+1000.0, 0.0,-300.0],
+ [-3000.0, 12000.0, -0.0],
+ [+3000.0, 12000.0, -0.0]])
reader = VTKDataReader()
@@ -105,8 +105,8 @@
lengthScale = 1000.0
timeScale = 1000.0
-dip = -targets[:,1] / lengthScale
-strike = targets[:,2] / lengthScale
+dip = -targets[:,2] / lengthScale
+strike = targets[:,1] / lengthScale
body = -targets[:,0] / lengthScale
time = timestamps / timeScale
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_onfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_onfaultdata.py 2010-06-01 21:36:58 UTC (rev 16853)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_onfaultdata.py 2010-06-01 21:37:19 UTC (rev 16854)
@@ -28,15 +28,15 @@
# ----------------------------------------------------------------------
timestamps = numpy.arange(50,15001,50)
targets = numpy.array([[ 0.0000000, 0.0000000, 0.0000000],
- [ 0.0000000, 0.0000000, 4500.0000000],
- [ 0.0000000, 0.0000000, 12000.0000000],
- [ -750.0000000, -1299.0381057, 0.0000000],
- [ -1500.0000000, -2598.0762114, 0.0000000],
- [ -2250.0000000, -3897.1143170, 0.0000000],
- [ -3750.0000000, -6495.1905284, 0.0000000],
- [ -3750.0000000, -6495.1905284, 4500.0000000],
- [ -3750.0000000, -6495.1905284, 12000.0000000],
- [ -6000.0000000, -10392.3048454, 0.0000000]])
+ [ 0.0000000, 4500.0000000, 0.0000000],
+ [ 0.0000000, 12000.0000000, 0.0000000],
+ [ -750.0000000, 0.0000000, -1299.0381057],
+ [ -1500.0000000, 0.0000000, -2598.0762114],
+ [ -2250.0000000, 0.0000000, -3897.1143170],
+ [ -3750.0000000, 0.0000000, -6495.1905284],
+ [ -3750.0000000, 4500.0000000, -6495.1905284],
+ [ -3750.0000000, 12000.0000000, -6495.1905284],
+ [ -6000.0000000, 0.0000000, -10392.3048454]])
reader = VTKDataReader()
@@ -105,7 +105,7 @@
lengthScale = 1000.0
timeScale = 1000.0
dip = -2*targets[:,0] / lengthScale
-strike = targets[:,2] / lengthScale
+strike = targets[:,1] / lengthScale
time = timestamps / timeScale
print "time", time
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_ruptime.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_ruptime.py 2010-06-01 21:36:58 UTC (rev 16853)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_ruptime.py 2010-06-01 21:37:19 UTC (rev 16854)
@@ -105,7 +105,7 @@
lengthScale = 1000.0
timeScale = 1000.0
distDip = -2*vertices[:,0]
-distStrike = vertices[:,2]
+distStrike = vertices[:,1]
time = timestamps / timeScale
print "time", time
More information about the CIG-COMMITS
mailing list