[cig-commits] r19566 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d
brad at geodynamics.org
brad at geodynamics.org
Wed Feb 1 15:34:04 PST 2012
Author: brad
Date: 2012-02-01 15:34:04 -0800 (Wed, 01 Feb 2012)
New Revision: 19566
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/gradient.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_050m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_050m.jou
Removed:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/refinecells.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.jou
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m.jou
Log:
Reworked tpv205-2d cubit journal files.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/createbc.jou 2012-02-01 23:28:22 UTC (rev 19565)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/createbc.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -1,44 +1,44 @@
# ----------------------------------------------------------------------
# Create blocks for materials
# ----------------------------------------------------------------------
-block 1 surface 18 19
+block 1 surface 8 9
block 1 name "elastic"
# ----------------------------------------------------------------------
# Create nodeset for fault
# ----------------------------------------------------------------------
-group "fault" add node in curve 33
-group "fault" remove node with x_coord < -7.5001e+3
-nodeset 10 group fault
-nodeset 10 name "fault"
+group "fault" add node in curve fault_surface
+nodeset 20 group fault
+nodeset 20 name "fault"
# ----------------------------------------------------------------------
# Create nodeset for +x face
# ----------------------------------------------------------------------
-group "face_xpos" add node in curve 39
-nodeset 11 group face_xpos
-nodeset 11 name "face_xpos"
+group "face_xpos" add node in curve 14
+nodeset 10 group face_xpos
+nodeset 10 name "face_xpos"
# ----------------------------------------------------------------------
# Create nodeset for -x face
# ----------------------------------------------------------------------
-group "face_xneg" add node in curve 35
-nodeset 12 group face_xneg
-nodeset 12 name "face_xneg"
+group "face_xneg" add node in curve 16
+nodeset 11 group face_xneg
+nodeset 11 name "face_xneg"
# ----------------------------------------------------------------------
# Create nodeset for +y face
# ----------------------------------------------------------------------
-group "face_ypos" add node in curve 36
-group "face_ypos" add node in curve 40
-nodeset 13 group face_ypos
-nodeset 13 name "face_ypos"
+group "face_ypos" add node in curve 19
+group "face_ypos" add node in curve 20
+nodeset 12 group face_ypos
+nodeset 12 name "face_ypos"
# ----------------------------------------------------------------------
# Create nodeset for -y face
# ----------------------------------------------------------------------
-group "face_yneg" add node in curve 34
-group "face_yneg" add node in curve 38
-nodeset 14 group face_yneg
-nodeset 14 name "face_yneg"
+group "face_yneg" add node in curve 18
+group "face_yneg" add node in curve 21
+nodeset 13 group face_yneg
+nodeset 13 name "face_yneg"
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry.jou 2012-02-01 23:28:22 UTC (rev 19565)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -1,44 +1,60 @@
+# -*- Python -*- (syntax highlighting)
# ----------------------------------------------------------------------
-# Create block
+#
+# Brad T. Aagaard, U.S. Geological Survey
+#
# ----------------------------------------------------------------------
-# Block is 36.0 km x 18.0*sqrt(3) km
-# -21.75 km <= x <= 14.25 km
-# -18*Sqrt(3) km <= y <= 0 km
#
-# We create the 2-D mesh by extracting the mid-surface from a 3-D mesh.
-reset
-brick x 36000 y 31176.914 z 48000
-
+# CUBIT journal file with geometry for TPV210-2D.
+#
# ----------------------------------------------------------------------
-# Set z=0 at free surface
+# Set units to SI.
# ----------------------------------------------------------------------
-move volume 1 x -3750 y -15588.457
-
+# {Units('si')}
+#
# ----------------------------------------------------------------------
-# Create interface surfaces
+# Reset geometry.
# ----------------------------------------------------------------------
-create surface extended from surface 5
-body 2 rotate +60 about z
-webcut volume 1 with plane surface 7
+reset
-delete surface 7
+# Make sure undo is off to prevent errors in stitching volumes.
+undo off
# ----------------------------------------------------------------------
-# Extract the mid surfaces for 2-D
+# Create block
# ----------------------------------------------------------------------
-create midsurface volume 3 surface 15 16
-create midsurface volume 1 surface 9 11
+# Block is 48.0 km x 36.0 km
+# -24.0 km <= x <= 24.0 km
+# -36.0 km <= y <= 0.0 km
+#
+# Fault dip is 60 degrees and width is 15.0 km.
+#
+#{blockWidth=48.0*km}
+#{blockHeight=36.0*km}
+#
+#{faultDipAngle=60.0}
+#{faultWidth=15.0*km}
+brick x {blockWidth} y {blockHeight} z {blockWidth}
+# Translate block so the top is at z=0
+volume 1 move x {-0.5*faultWidth*cosd(faultDipAngle)} y {-0.5*blockHeight} z 0
+
+create midsurface volume 1 surface 1 2
delete volume 1
-delete volume 3
# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+
+# Create fault (yz plane) at x = 0.0
+split surface 7 across location position 0 0 0 location position {-blockHeight/tand(faultDipAngle)} {-blockHeight} 0
+split curve 17 at position {-faultWidth*cosd(faultDipAngle)} {-faultWidth*sind(faultDipAngle)} 0
+curve 23 name "fault_surface"
+
+# ----------------------------------------------------------------------
# Imprint all volumes, merging surfaces
# ----------------------------------------------------------------------
imprint all
merge all
-
-
-
-
+# End of file
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry_refine.jou 2012-02-01 23:28:22 UTC (rev 19565)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry_refine.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -1,44 +0,0 @@
-# ----------------------------------------------------------------------
-# Create block
-# ----------------------------------------------------------------------
-# Block is 144.0 km x 72 km
-# -72 km <= x <= 72 km
-# -72 km <= y <= 0 km
-#
-# We create the 2-D mesh by extracting the mid-surface from a 3-D mesh.
-reset
-brick x 144000 y 72000 z 72000
-
-# ----------------------------------------------------------------------
-# Set z=0 at free surface
-# ----------------------------------------------------------------------
-move volume 1 x -3750 y -36000
-
-# ----------------------------------------------------------------------
-# Create interface surfaces
-# ----------------------------------------------------------------------
-create surface extended from surface 5
-body 2 rotate +60 about z
-webcut volume 1 with plane surface 7
-
-delete surface 7
-
-# ----------------------------------------------------------------------
-# Extract the mid surfaces for 2-D
-# ----------------------------------------------------------------------
-create midsurface volume 3 surface 15 16
-create midsurface volume 1 surface 9 11
-
-delete volume 1
-delete volume 3
-
-# ----------------------------------------------------------------------
-# Imprint all volumes, merging surfaces
-# ----------------------------------------------------------------------
-imprint all
-merge all
-
-
-
-
-
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/gradient.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/gradient.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/gradient.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -0,0 +1,58 @@
+# ----------------------------------------------------------------------
+# Set vertex spacing with increasing spacing away from fault
+# ----------------------------------------------------------------------
+#{bias_factor=1.015}
+
+# Function for computing discretization size at end of curve with bias
+#{fun_dx='dxStart*bias_factor**ceil( ln(1-curveL/dxStart*(1-bias_factor))/ln(bias_factor))'}
+
+# ----------------------------------------------------------------------
+# Compute sizes at curve endpoints
+# ----------------------------------------------------------------------
+
+# dxA - size at vertex 11
+#{dxStart=dx}{curveL=Length(19)}{execute("dxA="//fun_dx)}
+
+# dxB - size at vertex 12
+#{dxStart=dx}{curveL=Length(20)}{execute("dxB="//fun_dx)}
+
+# dxC - size at vertex 13
+#{dxStart=dx}{curveL=Length(22)}{execute("dxC="//fun_dx)}
+
+# dxD - size at vertex 9
+#{dxStart=dxA}{curveL=Length(14)}{execute("dxD="//fun_dx)}
+
+# dxE - size at vertex 10
+#{dxStart=dxC}{curveL=Length(21)}{execute("dxE="//fun_dx)}
+
+# ----------------------------------------------------------------------
+# Reset sizes
+# ----------------------------------------------------------------------
+curve all scheme default
+surface all sizing function none
+
+# Set size on faults
+curve fault_surface size {dx}
+
+# Set bias on curves extending from faults
+curve 19 scheme bias fine size {dx} factor {bias_factor} start vertex 14
+curve 20 scheme bias fine size {dx} factor {bias_factor} start vertex 14
+curve 22 scheme bias fine size {dx} factor {bias_factor} start vertex 15
+
+# A to D
+curve 14 scheme bias fine size {dxA} factor {bias_factor} start vertex 11
+
+# C to E
+curve 21 scheme bias fine size {dxC} factor {bias_factor} start vertex 13
+
+# C to D
+curve 18 scheme bias fine size {dxC} coarse size {dxD} start vertex 13
+
+# B to E
+curve 16 scheme bias fine size {dxB} coarse size {dxE} start vertex 12
+
+surface 8 9 sizing function type bias start curve fault_surface factor {bias_factor}
+
+
+# End of file
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_050m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_050m.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_050m.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -0,0 +1,40 @@
+# ----------------------------------------------------------------------
+# Create quad4 mesh at 50m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size and scheme
+# ----------------------------------------------------------------------
+#{dx=50*m}
+surface all scheme pave
+playback 'gradient.jou'
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "quad4_050m.exo" dimension 2 overwrite
+
+
+# End of file
+
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.exo.gz
===================================================================
(Binary files differ)
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.jou 2012-02-01 23:28:22 UTC (rev 19565)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -8,14 +8,14 @@
playback 'geometry.jou'
# ----------------------------------------------------------------------
-# Set meshing scheme
+# Set discretization size and scheme
# ----------------------------------------------------------------------
-surface 18 19 scheme pave
+#{dx=100*m}
+surface all scheme pave
+playback 'gradient.jou'
-# ----------------------------------------------------------------------
-# Set discretization size
-# ----------------------------------------------------------------------
-surface all size 100
+# A to D
+curve 14 scheme bias fine size {dxA*1.01} factor {bias_factor} start vertex 11
# ----------------------------------------------------------------------
# Generate the mesh
@@ -23,6 +23,12 @@
mesh surface all
# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
+
+# ----------------------------------------------------------------------
# Setup boundary conditions.
# ----------------------------------------------------------------------
playback 'createbc.jou'
@@ -32,3 +38,6 @@
# ----------------------------------------------------------------------
export mesh "quad4_100m.exo" dimension 2 overwrite
+
+# End of file
+
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.exo.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.jou 2012-02-01 23:28:22 UTC (rev 19565)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m_refine.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -1,46 +0,0 @@
-# ----------------------------------------------------------------------
-# Create quad4 mesh at 100m resolution.
-# ----------------------------------------------------------------------
-
-# ----------------------------------------------------------------------
-# Generate geometry
-# ----------------------------------------------------------------------
-playback 'geometry_refine.jou'
-
-# ----------------------------------------------------------------------
-# Set meshing scheme
-# ----------------------------------------------------------------------
-surface all scheme pave
-
-# ----------------------------------------------------------------------
-# Set discretization size
-# ----------------------------------------------------------------------
-surface all size 400
-
-# ----------------------------------------------------------------------
-# Generate the mesh
-# ----------------------------------------------------------------------
-mesh surface all
-
-# ----------------------------------------------------------------------
-# Refine elements within 1/2 of rupture width
-# ----------------------------------------------------------------------
-playback 'refinecells.jou'
-
-# ----------------------------------------------------------------------
-# Smooth mesh to improve quality
-# ----------------------------------------------------------------------
-surface all smooth scheme condition number beta 1.5 cpu 10
-smooth surface all
-
-# ----------------------------------------------------------------------
-# Setup boundary conditions.
-# ----------------------------------------------------------------------
-playback 'createbc.jou'
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "quad4_100m_refine.exo" dimension 2 overwrite
-
-
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.exo.gz
===================================================================
(Binary files differ)
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.jou 2012-02-01 23:28:22 UTC (rev 19565)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -8,19 +8,22 @@
playback 'geometry.jou'
# ----------------------------------------------------------------------
-# Set meshing scheme
+# Set discretization size and scheme
# ----------------------------------------------------------------------
-surface 18 19 scheme pave
+#{dx=200*m}
+surface all scheme pave
+playback 'gradient.jou'
# ----------------------------------------------------------------------
-# Set discretization size
+# Generate the mesh
# ----------------------------------------------------------------------
-surface all size 200
+mesh surface all
# ----------------------------------------------------------------------
-# Generate the mesh
+# Smooth mesh to improve quality
# ----------------------------------------------------------------------
-mesh surface all
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
# ----------------------------------------------------------------------
# Setup boundary conditions.
@@ -32,3 +35,5 @@
# ----------------------------------------------------------------------
export mesh "quad4_200m.exo" dimension 2 overwrite
+
+# End of file
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.exo.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.jou 2012-02-01 23:28:22 UTC (rev 19565)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m_refine.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -1,45 +0,0 @@
-# ----------------------------------------------------------------------
-# Create quad4 mesh at 200m resolution.
-# ----------------------------------------------------------------------
-
-# ----------------------------------------------------------------------
-# Generate geometry
-# ----------------------------------------------------------------------
-playback 'geometry_refine.jou'
-
-# ----------------------------------------------------------------------
-# Set meshing scheme
-# ----------------------------------------------------------------------
-surface all scheme pave
-
-# ----------------------------------------------------------------------
-# Set discretization size
-# ----------------------------------------------------------------------
-surface all size 800
-
-# ----------------------------------------------------------------------
-# Generate the mesh
-# ----------------------------------------------------------------------
-mesh surface all
-
-# ----------------------------------------------------------------------
-# Refine elements within 1/2 of rupture width
-# ----------------------------------------------------------------------
-playback 'refinecells.jou'
-
-# ----------------------------------------------------------------------
-# Smooth mesh to improve quality
-# ----------------------------------------------------------------------
-surface all smooth scheme condition number beta 1.3 cpu 10
-smooth surface all
-
-# ----------------------------------------------------------------------
-# Setup boundary conditions.
-# ----------------------------------------------------------------------
-playback 'createbc.jou'
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "quad4_200m_refine.exo" dimension 2 overwrite
-
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/refinecells.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/refinecells.jou 2012-02-01 23:28:22 UTC (rev 19565)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/refinecells.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -1,8 +0,0 @@
-# ----------------------------------------------------------------------
-# Refine elements within 1/2 of rupture width
-# ----------------------------------------------------------------------
-group "refinecells1" add node with x_min >= -39.751e+3 and x_max <= 32.251e+3 and y_min >= -36.001e+3
-refine node in refinecells1 depth 1
-group "refinecells2" add node with x_min >= -21.751e+3 and x_max <= 14.251e+3 and y_min >= -24.001e+3
-refine node in refinecells2 depth 1
-
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_050m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_050m.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_050m.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -0,0 +1,40 @@
+# ----------------------------------------------------------------------
+# Create tri3 mesh at 50m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size and scheme
+# ----------------------------------------------------------------------
+#{dx=50*m}
+surface all scheme trimesh
+playback 'gradient.jou'
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "tri3_050m.exo" dimension 2 overwrite
+
+
+# End of file
+
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m.exo.gz
===================================================================
(Binary files differ)
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m.jou 2012-02-01 23:28:22 UTC (rev 19565)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -10,8 +10,9 @@
# ----------------------------------------------------------------------
# Set discretization size and scheme
# ----------------------------------------------------------------------
-surface all size 100
-surface all scheme trimesh
+#{dx=100*m}
+surface all scheme trimesh
+playback 'gradient.jou'
# ----------------------------------------------------------------------
# Generate the mesh
@@ -19,6 +20,12 @@
mesh surface all
# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
+
+# ----------------------------------------------------------------------
# Setup boundary conditions.
# ----------------------------------------------------------------------
playback 'createbc.jou'
@@ -29,4 +36,5 @@
export mesh "tri3_100m.exo" dimension 2 overwrite
+# End of file
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.exo.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.jou 2012-02-01 23:28:22 UTC (rev 19565)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m_refine.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -1,46 +0,0 @@
-# ----------------------------------------------------------------------
-# Create tri3 mesh at 100m resolution.
-# ----------------------------------------------------------------------
-
-# ----------------------------------------------------------------------
-# Generate geometry
-# ----------------------------------------------------------------------
-playback 'geometry_refine.jou'
-
-# ----------------------------------------------------------------------
-# Set discretization size and scheme
-# ----------------------------------------------------------------------
-surface all size 400
-surface all scheme trimesh
-
-# ----------------------------------------------------------------------
-# Generate the mesh
-# ----------------------------------------------------------------------
-mesh surface all
-
-# ----------------------------------------------------------------------
-# Refine elements within 1/2 of rupture width
-# ----------------------------------------------------------------------
-playback 'refinecells.jou'
-
-# ----------------------------------------------------------------------
-# Smooth mesh to improve quality
-# ----------------------------------------------------------------------
-surface all smooth scheme laplacian
-smooth surface all
-surface all smooth scheme condition number beta 1.3 cpu 10
-smooth surface all
-
-# ----------------------------------------------------------------------
-# Setup boundary conditions.
-# ----------------------------------------------------------------------
-playback 'createbc.jou'
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "tri3_100m_refine.exo" dimension 2 overwrite
-
-
-
-
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m.exo.gz
===================================================================
(Binary files differ)
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m.jou 2012-02-01 23:28:22 UTC (rev 19565)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -10,8 +10,9 @@
# ----------------------------------------------------------------------
# Set discretization size and scheme
# ----------------------------------------------------------------------
-surface all size 200
-surface all scheme trimesh
+#{dx=200*m}
+surface all scheme trimesh
+playback 'gradient.jou'
# ----------------------------------------------------------------------
# Generate the mesh
@@ -19,6 +20,12 @@
mesh surface all
# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
+
+# ----------------------------------------------------------------------
# Setup boundary conditions.
# ----------------------------------------------------------------------
playback 'createbc.jou'
@@ -29,4 +36,4 @@
export mesh "tri3_200m.exo" dimension 2 overwrite
-
+# End of file
\ No newline at end of file
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.exo.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.jou 2012-02-01 23:28:22 UTC (rev 19565)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_200m_refine.jou 2012-02-01 23:34:04 UTC (rev 19566)
@@ -1,46 +0,0 @@
-# ----------------------------------------------------------------------
-# Create tri3 mesh at 200m resolution.
-# ----------------------------------------------------------------------
-
-# ----------------------------------------------------------------------
-# Generate geometry
-# ----------------------------------------------------------------------
-playback 'geometry_refine.jou'
-
-# ----------------------------------------------------------------------
-# Set discretization size and scheme
-# ----------------------------------------------------------------------
-surface all size 800
-surface all scheme trimesh
-
-# ----------------------------------------------------------------------
-# Generate the mesh
-# ----------------------------------------------------------------------
-mesh surface all
-
-# ----------------------------------------------------------------------
-# Refine elements within 1/2 of rupture width
-# ----------------------------------------------------------------------
-playback 'refinecells.jou'
-
-# ----------------------------------------------------------------------
-# Smooth mesh to improve quality
-# ----------------------------------------------------------------------
-surface all smooth scheme laplacian
-smooth surface all
-surface all smooth scheme condition number beta 1.3 cpu 10
-smooth surface all
-
-# ----------------------------------------------------------------------
-# Setup boundary conditions.
-# ----------------------------------------------------------------------
-playback 'createbc.jou'
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "tri3_200m_refine.exo" dimension 2 overwrite
-
-
-
-
More information about the CIG-COMMITS
mailing list