[cig-commits] [commit] master: Updated savageprescott meshing scripts for PyLith v2.0. (2d020ad)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jul 2 15:36:54 PDT 2014
Repository : https://github.com/geodynamics/pylith_benchmarks
On branch : master
Link : https://github.com/geodynamics/pylith_benchmarks/compare/0c306e7606c6413089d7eda6f4faa6df06259470...2d020ad58bc95cac1ff40486d6b3c7259755665d
>---------------------------------------------------------------
commit 2d020ad58bc95cac1ff40486d6b3c7259755665d
Author: Brad Aagaard <baagaard at usgs.gov>
Date: Wed Jul 2 15:36:56 2014 -0700
Updated savageprescott meshing scripts for PyLith v2.0.
>---------------------------------------------------------------
2d020ad58bc95cac1ff40486d6b3c7259755665d
.../savageprescott/meshes/createbc.jou | 17 ++-
.../savageprescott/meshes/geometry.jou | 41 ++++---
.../savageprescott/meshes/gradient.jou | 132 ++++++++++-----------
.../savageprescott/meshes/hex8_20km.jou | 3 +-
.../savageprescott/meshes/hex8_6.7km.jou | 3 +-
.../savageprescott/meshes/tet4_20km.jou | 11 +-
.../savageprescott/meshes/tet4_6.7km.jou | 11 +-
7 files changed, 116 insertions(+), 102 deletions(-)
diff --git a/quasistatic/sceccrustdeform/savageprescott/meshes/createbc.jou b/quasistatic/sceccrustdeform/savageprescott/meshes/createbc.jou
index 5251ff3..84bf4c6 100644
--- a/quasistatic/sceccrustdeform/savageprescott/meshes/createbc.jou
+++ b/quasistatic/sceccrustdeform/savageprescott/meshes/createbc.jou
@@ -15,21 +15,25 @@ group "fault" add node in fault_surface at D
group "fault" add node in fault_surface at E
group "fault" add node in fault_surface at B
group "fault" add node in fault_surface at C
-group "fault" remove node with z_coord < {-4.0*lockingDepth-0.01*km}
-nodeset 10 fault
-nodeset 10 name "fault"
+nodeset 30 fault
+nodeset 30 name "fault"
+
+# Mark buried edges
+group "fault_edge" add node in curve 59 63
+nodeset 31 fault_edge
+nodeset 31 name "fault_edge"
# ----------------------------------------------------------------------
# Create nodeset for +x face
# ----------------------------------------------------------------------
-group "face_xpos" add node in surface 11 14 26 27
+group "face_xpos" add node in surface 11 14 25 26
nodeset 11 face_xpos
nodeset 11 name "face_xpos"
# ----------------------------------------------------------------------
# Create nodeset for -x face
# ----------------------------------------------------------------------
-group "face_xneg" add node in surface 19 20 32 34
+group "face_xneg" add node in surface 19 20 30 32
nodeset 12 face_xneg
nodeset 12 name "face_xneg"
@@ -50,7 +54,7 @@ nodeset 14 name "face_yneg"
# ----------------------------------------------------------------------
# Create nodeset for -z face
# ----------------------------------------------------------------------
-group "face_zneg" add node in surface 24 29 31 33
+group "face_zneg" add node in surface 24 27 29 31
nodeset 15 face_zneg
nodeset 15 name "face_zneg"
@@ -81,3 +85,4 @@ group "face_ypos_nofault" remove node in face_xneg
nodeset 19 face_ypos_nofault
nodeset 19 name "face_ypos_nofault"
+
diff --git a/quasistatic/sceccrustdeform/savageprescott/meshes/geometry.jou b/quasistatic/sceccrustdeform/savageprescott/meshes/geometry.jou
index 1d2c1b0..61b7eea 100644
--- a/quasistatic/sceccrustdeform/savageprescott/meshes/geometry.jou
+++ b/quasistatic/sceccrustdeform/savageprescott/meshes/geometry.jou
@@ -20,7 +20,7 @@
# ----------------------------------------------------------------------
# Set units to SI.
# ----------------------------------------------------------------------
-#{Units('si')}
+${Units('si')}
#
# ----------------------------------------------------------------------
# Reset geometry.
@@ -38,11 +38,11 @@ undo off
# -500 km <= y <= 500 km
# -400 km <= z <= 0 km
#
-#{blockWidth=1000.0*km}
-#{blockLength=2000.0*km}
-#{blockHeight=400.0*km}
-#{faultWidth=40.0*km}
-#{lockingDepth=20.0*km}
+${blockWidth=1000.0*km}
+${blockLength=2000.0*km}
+${blockHeight=400.0*km}
+${faultWidth=40.0*km}
+${lockingDepth=20.0*km}
brick x {blockLength} y {blockWidth} z {blockHeight}
@@ -53,35 +53,40 @@ volume 1 move x {-xoffset} y 0 z {-0.5*blockHeight}
# Create interface surfaces
# ----------------------------------------------------------------------
create planar surface with plane xplane offset 0
-#{sid=Id("surface")}
+${sid=Id("surface")}
surface {sid} name "fault_surface"
create planar surface with plane zplane offset {-faultWidth}
-#{sid=Id("surface")}
+${sid=Id("surface")}
surface {sid} name "material_interface"
create planar surface with plane yplane offset 0
-#{sid=Id("surface")}
+${sid=Id("surface")}
surface {sid} name "symmetry_plane"
-create planar surface with plane zplane offset {-lockingDepth}
-#{sid=Id("surface")}
-surface {sid} name "locking_transition"
-
# ----------------------------------------------------------------------
# Cut volume with interfaces.
# ----------------------------------------------------------------------
# Cut fault with locking transition (not needed in volume)
webcut volume 1 with plane fault_surface
-webcut volume 1 6 with plane material_interface
-webcut volume 1 6 7 8 with plane symmetry_plane
-webcut volume 2 with plane locking_transition
-delete body 5
+webcut volume 1 5 with plane material_interface
+webcut volume 1 5 6 7 with plane symmetry_plane
+
+# Inscribe locking transition on fault plane
+brick x {5.0*km} y {blockWidth+10.0*km} z {2.0*lockingDepth}
+${vid=Id("volume")}
+chop volume 2 with volume {vid}
+
+# Inscribe bottom of fault on fault plane
+brick x {5.0*km} y {blockWidth+10.0*km} z {4.0*faultWidth}
+${vid=Id("volume")}
+chop volume 14 with volume {vid}
# ----------------------------------------------------------------------
# Imprint all volumes, merging surfaces.
# ----------------------------------------------------------------------
imprint all
merge all
-delete body 2 3 4 13
+delete body 3 4 13 16 17
compress all
+
diff --git a/quasistatic/sceccrustdeform/savageprescott/meshes/gradient.jou b/quasistatic/sceccrustdeform/savageprescott/meshes/gradient.jou
index 43fd82b..22f4345 100644
--- a/quasistatic/sceccrustdeform/savageprescott/meshes/gradient.jou
+++ b/quasistatic/sceccrustdeform/savageprescott/meshes/gradient.jou
@@ -2,34 +2,34 @@
# ----------------------------------------------------------------------
# Set vertex spacing with increasing spacing away from fault
# ----------------------------------------------------------------------
-#{bias_factor=1.02}
+${bias_factor=1.02}
# 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))'}
+${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 9 28 22 12 29 21
-#{dxStart=dx}{curveL=Length(27)}{execute("dxA="//fun_dx)}
+# dxA - size at vertex 19 22 29 31 20 21 28 30
+${dxStart=dx}{curveL=Length(25)}{execute("dxA="//fun_dx)}
# dxB - size at vertex 27
-#{dxStart=dx}{curveL=Length(41)}{execute("dxB="//fun_dx)}
+${dxStart=dx}{curveL=Length(57)}{execute("dxB="//fun_dx)}
-# dxC - size at vertex 10 11
-#{dxStart=dxB}{curveL=Length(47)}{execute("dxC="//fun_dx)}
+# dxC - size at vertex 9 10
+${dxStart=dxB}{curveL=Length(44)}{execute("dxC="//fun_dx)}
-# dxD - size at vertex 28 19 23 26
-#{dxStart=dx}{curveL=Length(25)}{execute("dxD="//fun_dx)}
+# dxD - size at vertex 16 26 17 23
+${dxStart=dx}{curveL=Length(23)}{execute("dxD="//fun_dx)}
-# dxE - size at vertex 17 20
-#{dxStart=dxD}{curveL=Length(43)}{execute("dxE="//fun_dx)}
+# dxE - size at vertex 1 2 3 4
+${dxStart=dxD}{curveL=Length(24)}{execute("dxE="//fun_dx)}
-# dxF - size at vertex 1 2 3 4 13 14 15 16
-#{dxStart=dxD}{curveL=Length(26)}{execute("dxF="//fun_dx)}
+# dxF - size at vertex 15 18
+${dxStart=dxD}{curveL=Length(40)}{execute("dxF="//fun_dx)}
# dxG - size at vertex 5 6 7 8
-#{dxStart=dxF}{curveL=Length(16)}{execute("dxG="//fun_dx)}
+${dxStart=dxF}{curveL=Length(43)}{execute("dxG="//fun_dx)}
# ----------------------------------------------------------------------
# Reset sizes
@@ -42,99 +42,99 @@ surface all scheme trimesh
curve 56 57 size {dx}
# fault to A
-curve 27 scheme bias fine size {dx} coarse size {dxA} start vertex 25
-curve 31 scheme bias fine size {dx} coarse size {dxA} start vertex 25
-curve 54 scheme bias fine size {dx} coarse size {dxA} start vertex 30
-curve 58 scheme bias fine size {dx} coarse size {dxA} start vertex 30
-curve 29 scheme bias fine size {dx} coarse size {dxA} start vertex 24
-curve 30 scheme bias fine size {dx} coarse size {dxA} start vertex 24
+curve 25 scheme bias fine size {dx} coarse size {dxA} start vertex 25
+curve 29 scheme bias fine size {dx} coarse size {dxA} start vertex 25
+curve 53 scheme bias fine size {dx} coarse size {dxA} start vertex 32
+curve 61 scheme bias fine size {dx} coarse size {dxA} start vertex 32
+curve 27 scheme bias fine size {dx} coarse size {dxA} start vertex 24
+curve 28 scheme bias fine size {dx} coarse size {dxA} start vertex 24
+curve 59 scheme bias fine size {dx} coarse size {dxA} start vertex 33
+curve 63 scheme bias fine size {dx} coarse size {dxA} start vertex 33
# fault to B
-curve 41 scheme bias fine size {dx} coarse size {dxB} start vertex 24
+curve 57 scheme bias fine size {dx} coarse size {dxB} start vertex 33
# fault to D
-curve 25 scheme bias fine size {dx} coarse size {dxD} start vertex 25
-curve 35 scheme bias fine size {dx} coarse size {dxD} start vertex 25
-curve 23 scheme bias fine size {dx} coarse size {dxD} start vertex 24
-curve 34 scheme bias fine size {dx} coarse size {dxD} start vertex 24
+curve 23 scheme bias fine size {dx} coarse size {dxD} start vertex 25
+curve 33 scheme bias fine size {dx} coarse size {dxD} start vertex 25
+curve 21 scheme bias fine size {dx} coarse size {dxD} start vertex 24
+curve 32 scheme bias fine size {dx} coarse size {dxD} start vertex 24
# ----------------------------------------------------------------------
# bias starting at A
# A to C
-curve 13 scheme bias fine size {dxA} coarse size {dxC} start vertex 22
-curve 15 scheme bias fine size {dxA} coarse size {dxC} start vertex 21
+curve 51 scheme bias fine size {dxA} coarse size {dxC} start vertex 31
+curve 52 scheme bias fine size {dxA} coarse size {dxC} start vertex 30
-# A to F
-curve 1 scheme bias fine size {dxA} coarse size {dxF} start vertex 9
-curve 8 scheme bias fine size {dxA} coarse size {dxF} start vertex 9
-curve 2 scheme bias fine size {dxA} coarse size {dxF} start vertex 12
-curve 6 scheme bias fine size {dxA} coarse size {dxF} start vertex 12
+# A to E
+curve 1 scheme bias fine size {dxA} coarse size {dxE} start vertex 19
+curve 8 scheme bias fine size {dxA} coarse size {dxE} start vertex 19
+curve 2 scheme bias fine size {dxA} coarse size {dxE} start vertex 20
+curve 6 scheme bias fine size {dxA} coarse size {dxE} start vertex 20
-curve 10 scheme bias fine size {dxA} coarse size {dxF} start vertex 22
-curve 17 scheme bias fine size {dxA} coarse size {dxF} start vertex 22
-curve 9 scheme bias fine size {dxA} coarse size {dxF} start vertex 21
-curve 18 scheme bias fine size {dxA} coarse size {dxF} start vertex 21
+curve 10 scheme bias fine size {dxA} coarse size {dxE} start vertex 22
+curve 15 scheme bias fine size {dxA} coarse size {dxE} start vertex 22
+curve 9 scheme bias fine size {dxA} coarse size {dxE} start vertex 21
+curve 16 scheme bias fine size {dxA} coarse size {dxE} start vertex 21
# A to A
-curve 53 59 52 55 size {dxA}
+curve 50 62 64 49 54 60 size {dxA}
# ----------------------------------------------------------------------
# bias starting at B
# B to C
curve 44 scheme bias fine size {dxB} coarse size {dxC} start vertex 27
-curve 47 scheme bias fine size {dxB} coarse size {dxC} start vertex 27
+curve 41 scheme bias fine size {dxB} coarse size {dxC} start vertex 27
-# B to E
-curve 42 scheme bias fine size {dxB} coarse size {dxE} start vertex 27
-curve 49 scheme bias fine size {dxB} coarse size {dxE} start vertex 27
+# B to F
+curve 39 scheme bias fine size {dxB} coarse size {dxF} start vertex 27
+curve 46 scheme bias fine size {dxB} coarse size {dxF} start vertex 27
# ----------------------------------------------------------------------
# bias starting at C
# C to G
-curve 3 scheme bias fine size {dxC} coarse size {dxG} start vertex 10
-curve 7 scheme bias fine size {dxC} coarse size {dxG} start vertex 10
-curve 4 scheme bias fine size {dxC} coarse size {dxG} start vertex 11
-curve 5 scheme bias fine size {dxC} coarse size {dxG} start vertex 11
+curve 3 scheme bias fine size {dxC} coarse size {dxG} start vertex 9
+curve 7 scheme bias fine size {dxC} coarse size {dxG} start vertex 9
+curve 4 scheme bias fine size {dxC} coarse size {dxG} start vertex 10
+curve 5 scheme bias fine size {dxC} coarse size {dxG} start vertex 10
# ----------------------------------------------------------------------
# bias starting at D
-curve 26 scheme bias fine size {dxD} coarse size {dxF} start vertex 19
-curve 33 scheme bias fine size {dxD} coarse size {dxF} start vertex 19
-curve 37 scheme bias fine size {dxD} coarse size {dxF} start vertex 18
-curve 40 scheme bias fine size {dxD} coarse size {dxF} start vertex 18
+curve 31 scheme bias fine size {dxD} coarse size {dxE} start vertex 17
+curve 24 scheme bias fine size {dxD} coarse size {dxE} start vertex 17
+curve 35 scheme bias fine size {dxD} coarse size {dxE} start vertex 16
+curve 38 scheme bias fine size {dxD} coarse size {dxE} start vertex 16
-curve 32 scheme bias fine size {dxD} coarse size {dxE} start vertex 23
-curve 28 scheme bias fine size {dxD} coarse size {dxF} start vertex 23
-curve 38 scheme bias fine size {dxD} coarse size {dxF} start vertex 26
-curve 39 scheme bias fine size {dxD} coarse size {dxF} start vertex 26
+curve 26 scheme bias fine size {dxD} coarse size {dxE} start vertex 23
+curve 30 scheme bias fine size {dxD} coarse size {dxE} start vertex 23
+curve 36 scheme bias fine size {dxD} coarse size {dxE} start vertex 26
+curve 37 scheme bias fine size {dxD} coarse size {dxE} start vertex 26
-curve 43 scheme bias fine size {dxD} coarse size {dxF} start vertex 23
-curve 48 scheme bias fine size {dxD} coarse size {dxE} start vertex 26
+curve 40 scheme bias fine size {dxD} coarse size {dxF} start vertex 23
+curve 45 scheme bias fine size {dxD} coarse size {dxF} start vertex 26
# D to D
-curve 24 36 size {dxD}
+curve 22 34 size {dxD}
# ----------------------------------------------------------------------
# bias starting at E
# E to G
-curve 45 scheme bias fine size {dxE} coarse size {dxG} start vertex 20
-curve 46 scheme bias fine size {dxE} coarse size {dxG} start vertex 20
-curve 50 scheme bias fine size {dxE} coarse size {dxG} start vertex 17
-curve 51 scheme bias fine size {dxE} coarse size {dxG} start vertex 17
+curve 13 scheme bias fine size {dxE} coarse size {dxG} start vertex 13
+curve 14 scheme bias fine size {dxE} coarse size {dxG} start vertex 14
+curve 19 scheme bias fine size {dxE} coarse size {dxG} start vertex 11
+curve 20 scheme bias fine size {dxE} coarse size {dxG} start vertex 12
# ----------------------------------------------------------------------
# bias starting at F
-curve 14 scheme bias fine size {dxF} coarse size {dxG} start vertex 15
-curve 16 scheme bias fine size {dxF} coarse size {dxG} start vertex 16
-curve 22 scheme bias fine size {dxF} coarse size {dxG} start vertex 14
-curve 21 scheme bias fine size {dxF} coarse size {dxG} start vertex 13
+curve 47 scheme bias fine size {dxF} coarse size {dxG} start vertex 15
+curve 48 scheme bias fine size {dxF} coarse size {dxG} start vertex 15
+curve 42 scheme bias fine size {dxF} coarse size {dxG} start vertex 18
+curve 43 scheme bias fine size {dxF} coarse size {dxG} start vertex 18
-# F to F
-curve 12 26 19 20 size {dxF}
diff --git a/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_20km.jou b/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_20km.jou
index ecbf554..051ee6a 100644
--- a/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_20km.jou
+++ b/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_20km.jou
@@ -8,7 +8,7 @@ playback 'geometry.jou'
# ----------------------------------------------------------------------
# Set discretization size to 20 km.
# ----------------------------------------------------------------------
-#{dx=20.0*km}
+${dx=20.0*km}
volume all size {dx}
mesh volume all
@@ -25,3 +25,4 @@ export mesh "hex8_20km.exo" dimension 3 overwrite
+
diff --git a/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_6.7km.jou b/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_6.7km.jou
index d00c33c..e710ed9 100644
--- a/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_6.7km.jou
+++ b/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_6.7km.jou
@@ -10,7 +10,7 @@ playback 'geometry.jou'
# ----------------------------------------------------------------------
# Set discretization size to 20 km.
# ----------------------------------------------------------------------
-#{dx=20.0*km}
+${dx=20.0*km}
volume all size {dx}
mesh volume all
@@ -38,3 +38,4 @@ export mesh "hex8_6.7km.exo" dimension 3 overwrite
+
diff --git a/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_20km.jou b/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_20km.jou
index 021f624..dc4fc4c 100644
--- a/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_20km.jou
+++ b/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_20km.jou
@@ -9,7 +9,7 @@ playback 'geometry.jou'
# ----------------------------------------------------------------------
# Set discretization size and scheme
# ----------------------------------------------------------------------
-#{dx=20.0*km}
+${dx=20.0*km}
playback 'gradient.jou'
@@ -26,13 +26,13 @@ mesh volume all
# ----------------------------------------------------------------------
# Smooth mesh to improve quality
# ----------------------------------------------------------------------
-#{condnum=2.0}
-#{loop(4)}
+${condnum=2.0}
+${loop(4)}
cleanup volume all
volume all smooth scheme condition number beta {condnum} cpu 10
smooth volume all
-#{condnum=condnum-0.1}
-#{endloop}
+${condnum=condnum-0.1}
+${endloop}
# ----------------------------------------------------------------------
# Setup boundary conditions.
@@ -48,3 +48,4 @@ export mesh "tet4_20km.exo" dimension 3 overwrite
+
diff --git a/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_6.7km.jou b/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_6.7km.jou
index 6e49463..7a23d8f 100644
--- a/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_6.7km.jou
+++ b/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_6.7km.jou
@@ -10,7 +10,7 @@ playback 'geometry.jou'
# ----------------------------------------------------------------------
# Set discretization size and scheme
# ----------------------------------------------------------------------
-#{dx=20.0*km/3.0}
+${dx=20.0*km/3.0}
playback 'gradient.jou'
@@ -27,13 +27,13 @@ mesh volume all
# ----------------------------------------------------------------------
# Smooth mesh to improve quality
# ----------------------------------------------------------------------
-#{condnum=2.0}
-#{loop(4)}
+${condnum=2.0}
+${loop(4)}
cleanup volume all
volume all smooth scheme condition number beta {condnum} cpu 10
smooth volume all
-#{condnum=condnum-0.1}
-#{endloop}
+${condnum=condnum-0.1}
+${endloop}
# ----------------------------------------------------------------------
# Setup boundary conditions.
@@ -49,3 +49,4 @@ export mesh "tet4_6.7km.exo" dimension 3 overwrite
+
More information about the CIG-COMMITS
mailing list