[cig-commits] r16832 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205
surendra at geodynamics.org
surendra at geodynamics.org
Sat May 29 16:24:31 PDT 2010
Author: surendra
Date: 2010-05-29 16:24:30 -0700 (Sat, 29 May 2010)
New Revision: 16832
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells_hex8.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells_tet4.jou
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_onfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m_refine.jou
Log:
Added separate files eof refining for hex8 and tet4 meshes (tpv205)
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m_refine.jou 2010-05-29 22:02:06 UTC (rev 16831)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_100m_refine.jou 2010-05-29 23:24:30 UTC (rev 16832)
@@ -20,7 +20,7 @@
# ----------------------------------------------------------------------
# Refine elements within 1/2 of rupture width
# ----------------------------------------------------------------------
-playback 'refinecells.jou'
+playback 'refinecells_hex8.jou'
# ----------------------------------------------------------------------
# Smooth mesh to improve quality
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m_refine.jou 2010-05-29 22:02:06 UTC (rev 16831)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m_refine.jou 2010-05-29 23:24:30 UTC (rev 16832)
@@ -20,7 +20,7 @@
# ----------------------------------------------------------------------
# Refine elements within 1/2 of rupture width
# ----------------------------------------------------------------------
-playback 'refinecells.jou'
+playback 'refinecells_hex8.jou'
# ----------------------------------------------------------------------
# Smooth mesh to improve quality
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells.jou 2010-05-29 22:02:06 UTC (rev 16831)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells.jou 2010-05-29 23:24:30 UTC (rev 16832)
@@ -1,7 +1,9 @@
# ----------------------------------------------------------------------
# Refine elements within 1/2 of rupture width
# ----------------------------------------------------------------------
-group "refinecells1" add node with x_min >= -24.001e+3 and x_max <= 24.001e+3 and y_min >= -36.001e+3 and y_max <= 36.001e+3
-refine node in refinecells1 depth 1
-group "refinecells2" add node with x_min >= -12.001e+3 and x_max <= 12.001e+3 and y_min >= -24.001e+3 and y_max <= 24.001e+3
-refine node in refinecells2 depth 1
+group "refinecells1" add tet with x_min >= -36.001e+3 and x_max <= 36.001e+3 and y_min >= -24.001e+3 and y_max <= 24.001e+3 and z_min >= -36.001e+3
+refine tet in refinecells1 depth 1
+group "refinecells2" add tet with x_min >= -24.001e+3 and x_max <= 24.001e+3 and y_min >= -12.001e+3 and y_max <= 12.001e+3 and z_min >= -24.001e+3
+refine tet in refinecells2 depth 1
+
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells_hex8.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells_hex8.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells_hex8.jou 2010-05-29 23:24:30 UTC (rev 16832)
@@ -0,0 +1,8 @@
+# ----------------------------------------------------------------------
+# Refine elements within 1/2 of rupture width
+# ----------------------------------------------------------------------
+group "refinecells1" add hex with x_min >= -36.001e+3 and x_max <= 36.001e+3 and y_min >= -24.001e+3 and y_max <= 24.001e+3 and z_min >= -36.001e+3
+refine hex in refinecells1 depth 1
+group "refinecells2" add hex with x_min >= -24.001e+3 and x_max <= 24.001e+3 and y_min >= -12.001e+3 and y_max <= 12.001e+3 and z_min >= -24.001e+3
+refine hex in refinecells2 depth 1
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells_tet4.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells_tet4.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells_tet4.jou 2010-05-29 23:24:30 UTC (rev 16832)
@@ -0,0 +1,9 @@
+# ----------------------------------------------------------------------
+# Refine elements within 1/2 of rupture width
+# ----------------------------------------------------------------------
+group "refinecells1" add tet with x_min >= -36.001e+3 and x_max <= 36.001e+3 and y_min >= -24.001e+3 and y_max <= 24.001e+3 and z_min >= -36.001e+3
+refine tet in refinecells1 depth 1
+group "refinecells2" add tet with x_min >= -24.001e+3 and x_max <= 24.001e+3 and y_min >= -12.001e+3 and y_max <= 12.001e+3 and z_min >= -24.001e+3
+refine tet in refinecells2 depth 1
+
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_onfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_onfaultdata.py 2010-05-29 22:02:06 UTC (rev 16831)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_onfaultdata.py 2010-05-29 23:24:30 UTC (rev 16832)
@@ -10,7 +10,7 @@
# ----------------------------------------------------------------------
#
-cell = "hex8"
+cell = "tet4"
dx = 200
dt = 0.05
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py 2010-05-29 22:02:06 UTC (rev 16831)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py 2010-05-29 23:24:30 UTC (rev 16832)
@@ -10,7 +10,7 @@
# ----------------------------------------------------------------------
#
-cell = "hex8"
+cell = "tet4"
dx = 200
dt = 0.05
@@ -37,8 +37,9 @@
# Extract values
nsteps = timestamps.shape[0]
-# Set default rupture time to a large value (1.0e+30)
+# Set default rupture time to a large value (1.0e+9)
rupTime = 1.0e+9 * numpy.ones( (npts,), dtype=numpy.float64)
+#rupTime = numpy.zeros( (npts,), dtype=numpy.float64)
# Create buffer for current rupture time
tmpTime = numpy.zeros( (npts,), dtype=numpy.float64)
@@ -59,14 +60,30 @@
# Set rupture time at locations where threshold is exceeded
tmpTime[slipRateMag > threshold] = t
+ #print "slipRateMag \n", slipRateMag
+ #print "slipRateMag > threshold \n ", slipRateMag > threshold
+
+ #print "tmpTime \n", tmpTime[1:50]
+ #print "time", t
+
# Get indicates where current time is less than current rupture
# time (this is only the locations that just started slipping)
- indices = numpy.where(tmpTime < rupTime)[0]
+ #indices = numpy.where((tmpTime < rupTime) and (rupTime < 0.00001))[0]
+ #indices = ((tmpTime < rupTime) & (rupTime < testarr))
+ indices = numpy.where((tmpTime < rupTime) & (tmpTime > 0.00001))[0]
+ #print "indices \n", indices[1:200]
+
+
rupTime[indices] = t
+ #print "rupTime \n", rupTime[1:200]
+
itime += 1
+# print "rupTime \n", rupTime
+# rupTime[rupTime < 0.00001] = 1.0e+9
+
# Write data
headerA = \
"# problem = TPV205\n" + \
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m_refine.jou 2010-05-29 22:02:06 UTC (rev 16831)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_100m_refine.jou 2010-05-29 23:24:30 UTC (rev 16832)
@@ -26,10 +26,7 @@
# ----------------------------------------------------------------------
# Smooth mesh to improve quality
# ----------------------------------------------------------------------
-volume all smooth scheme laplacian
-smooth volume all
-volume all smooth scheme condition number beta 1.3 cpu 10
-smooth volume all
+playback 'tet4_smooth.jou'
# ----------------------------------------------------------------------
# Setup boundary conditions.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m_refine.jou 2010-05-29 22:02:06 UTC (rev 16831)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m_refine.jou 2010-05-29 23:24:30 UTC (rev 16832)
@@ -21,15 +21,12 @@
# ----------------------------------------------------------------------
# Refine elements within 1/2 of rupture width
# ----------------------------------------------------------------------
-playback 'refinecells.jou'
+playback 'refinecells_tet4.jou'
# ----------------------------------------------------------------------
# Smooth mesh to improve quality
# ----------------------------------------------------------------------
-volume all smooth scheme laplacian
-smooth volume all
-volume all smooth scheme condition number beta 1.3 cpu 10
-smooth volume all
+playback 'tet4_smooth.jou'
# ----------------------------------------------------------------------
# Setup boundary conditions.
More information about the CIG-COMMITS
mailing list