[cig-commits] r12369 - in short/2.5D/benchmarks/savageprescott: meshes parameters utils
willic3 at geodynamics.org
willic3 at geodynamics.org
Wed Jul 2 14:45:35 PDT 2008
Author: willic3
Date: 2008-07-02 14:45:34 -0700 (Wed, 02 Jul 2008)
New Revision: 12369
Added:
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_1000m.exo
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_1000m.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_3000m.exo
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_3000m.jou
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_6750m.exo
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_6750m.jou
short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.cfg
short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.cfg
Removed:
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_8000m.exo
short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_8000m.jou
Modified:
short/2.5D/benchmarks/savageprescott/parameters/finalslip_rupture.spatialdb
short/2.5D/benchmarks/savageprescott/parameters/pylithapp.cfg
short/2.5D/benchmarks/savageprescott/parameters/sliprate_creep.spatialdb
short/2.5D/benchmarks/savageprescott/utils/savpres_ss.cfg
Log:
Added graded mesh benchmarks, and modified existing uniform benchmark and
analytical parameters to use the same elastic thickness (27 km) and
locking depth (13.5 km) as the graded mesh benchmarks.
These still need to be tested and compared to the analytical solution.
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_1000m.exo
===================================================================
(Binary files differ)
Property changes on: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_1000m.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_1000m.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_1000m.jou (rev 0)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_1000m.jou 2008-07-02 21:45:34 UTC (rev 12369)
@@ -0,0 +1,159 @@
+## /Applications/Plotting/Cubit_10_2.app/Contents/MacOS/Cubit_10_2
+## Cubit Version 10.2
+## Cubit Build 8
+## Revised 10/17/2006 15:13:41 MST
+## Running 06/27/2008 11:07:35 AM
+## Command Options:
+## -warning = On
+## -information = On
+
+# ----------------------------------------------------------------------
+# Create outer block.
+# ----------------------------------------------------------------------
+
+# Block is 486 km x 486 km x 486 km
+# -243 km <= x <= 243 km
+# -243 km <= y <= 243 km
+# -486 km <= y <= 0 km
+reset
+brick x 486000 y 486000 z 486000
+
+# Translate block so the top is at z=0.
+volume 1 move x 0 y 0 z -243000
+
+# ----------------------------------------------------------------------
+# Create inner block that will be at constant resolution.
+# ----------------------------------------------------------------------
+
+# Block is 54 km x 54 km x 54 km
+# -27 km <= x <= 27 km
+# -27 km <= y <= 27 km
+# -54 km <= y <= 0 km
+brick x 54000 y 54000 z 54000
+
+# Translate block so the top is at z=0.
+volume 2 move x 0 y 0 z -27000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces.
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+surface 13 name "fault_surface"
+create planar surface with plane zplane offset -27000
+surface 14 name "material_interface"
+create planar surface with plane yplane offset 0
+surface 15 name "symmetry_plane"
+
+# ----------------------------------------------------------------------
+# Chop outer volume with inner volume and then cut both with interfaces.
+# ----------------------------------------------------------------------
+chop volume 1 with volume 2
+compress all
+webcut volume all with plane fault_surface
+webcut volume all with plane symmetry_plane
+webcut volume all with plane material_interface
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces.
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Remove surface bodies.
+# ----------------------------------------------------------------------
+delete body 1 2 3 6 7 10 11 14 17 18 21 24
+
+# ----------------------------------------------------------------------
+# Create initial mesh at constant resolution.
+# ----------------------------------------------------------------------
+volume all size 9000
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine first region inside outer box.
+# ----------------------------------------------------------------------
+group "inner1" add hex with x_coord > -81000 and with x_coord < 81000 and with y_coord > -81000 and with y_coord < 81000 and with z_coord > -81000
+refine hex in inner1 depth 0
+
+# ----------------------------------------------------------------------
+# Refine second region inside outer box.
+# ----------------------------------------------------------------------
+group "inner2" add hex with x_coord > -27000 and with x_coord < 27000 and with y_coord > -27000 and with y_coord < 27000 and with z_coord > -27000
+refine hex in inner2 depth 0
+
+# ----------------------------------------------------------------------
+# Create blocks for materials.
+# ----------------------------------------------------------------------
+block 1 volume 4 5 8 9 12 13 15 16
+block 1 name "elastic"
+block 2 volume 19 20 22 23 25 26 27 28
+block 2 name "viscoelastic"
+
+# ----------------------------------------------------------------------
+# Create fault nodeset.
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface at A
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at B
+group "fault" add node in fault_surface at C
+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 G
+group "fault" add node in fault_surface at F
+nodeset 10 fault
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face.
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 121 176 131 179
+nodeset 11 face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face.
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 150 197 153 208
+nodeset 12 face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face.
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 122 151 132 154
+nodeset 13 face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face.
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 177 198 180 207
+nodeset 14 face_yneg
+nodeset 14 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face.
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 71 93 74 103
+nodeset 15 face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face without fault.
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 16 face_zneg_nofault
+nodeset 16 name "face_zneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face.
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 59 68 83 95 61 78 90 101
+nodeset 17 face_zpos
+nodeset 17 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Export exodus file.
+# ----------------------------------------------------------------------
+export mesh "spbm_hex8_graded_1000m.exo" dimension 3 overwrite
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_3000m.exo
===================================================================
(Binary files differ)
Property changes on: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_3000m.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_3000m.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_3000m.jou (rev 0)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded_3000m.jou 2008-07-02 21:45:34 UTC (rev 12369)
@@ -0,0 +1,159 @@
+## /Applications/Plotting/Cubit_10_2.app/Contents/MacOS/Cubit_10_2
+## Cubit Version 10.2
+## Cubit Build 8
+## Revised 10/17/2006 15:13:41 MST
+## Running 06/27/2008 11:07:35 AM
+## Command Options:
+## -warning = On
+## -information = On
+
+# ----------------------------------------------------------------------
+# Create outer block.
+# ----------------------------------------------------------------------
+
+# Block is 486 km x 486 km x 486 km
+# -243 km <= x <= 243 km
+# -243 km <= y <= 243 km
+# -486 km <= y <= 0 km
+reset
+brick x 486000 y 486000 z 486000
+
+# Translate block so the top is at z=0.
+volume 1 move x 0 y 0 z -243000
+
+# ----------------------------------------------------------------------
+# Create inner block that will be at constant resolution.
+# ----------------------------------------------------------------------
+
+# Block is 54 km x 54 km x 54 km
+# -27 km <= x <= 27 km
+# -27 km <= y <= 27 km
+# -54 km <= y <= 0 km
+brick x 54000 y 54000 z 54000
+
+# Translate block so the top is at z=0.
+volume 2 move x 0 y 0 z -27000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces.
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+surface 13 name "fault_surface"
+create planar surface with plane zplane offset -27000
+surface 14 name "material_interface"
+create planar surface with plane yplane offset 0
+surface 15 name "symmetry_plane"
+
+# ----------------------------------------------------------------------
+# Chop outer volume with inner volume and then cut both with interfaces.
+# ----------------------------------------------------------------------
+chop volume 1 with volume 2
+compress all
+webcut volume all with plane fault_surface
+webcut volume all with plane symmetry_plane
+webcut volume all with plane material_interface
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces.
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Remove surface bodies.
+# ----------------------------------------------------------------------
+delete body 1 2 3 6 7 10 11 14 17 18 21 24
+
+# ----------------------------------------------------------------------
+# Create initial mesh at constant resolution.
+# ----------------------------------------------------------------------
+volume all size 27000
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine first region inside outer box.
+# ----------------------------------------------------------------------
+group "inner1" add hex with x_coord > -81000 and with x_coord < 81000 and with y_coord > -81000 and with y_coord < 81000 and with z_coord > -81000
+refine hex in inner1 depth 0
+
+# ----------------------------------------------------------------------
+# Refine second region inside outer box.
+# ----------------------------------------------------------------------
+group "inner2" add hex with x_coord > -27000 and with x_coord < 27000 and with y_coord > -27000 and with y_coord < 27000 and with z_coord > -27000
+refine hex in inner2 depth 0
+
+# ----------------------------------------------------------------------
+# Create blocks for materials.
+# ----------------------------------------------------------------------
+block 1 volume 4 5 8 9 12 13 15 16
+block 1 name "elastic"
+block 2 volume 19 20 22 23 25 26 27 28
+block 2 name "viscoelastic"
+
+# ----------------------------------------------------------------------
+# Create fault nodeset.
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface at A
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at B
+group "fault" add node in fault_surface at C
+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 G
+group "fault" add node in fault_surface at F
+nodeset 10 fault
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face.
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 121 176 131 179
+nodeset 11 face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face.
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 150 197 153 208
+nodeset 12 face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face.
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 122 151 132 154
+nodeset 13 face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face.
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 177 198 180 207
+nodeset 14 face_yneg
+nodeset 14 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face.
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 71 93 74 103
+nodeset 15 face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face without fault.
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 16 face_zneg_nofault
+nodeset 16 name "face_zneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face.
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 59 68 83 95 61 78 90 101
+nodeset 17 face_zpos
+nodeset 17 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Export exodus file.
+# ----------------------------------------------------------------------
+export mesh "spbm_hex8_graded_3000m.exo" dimension 3 overwrite
Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_6750m.exo
===================================================================
(Binary files differ)
Property changes on: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_6750m.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Copied: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_6750m.jou (from rev 12274, short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_8000m.jou)
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_6750m.jou (rev 0)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_6750m.jou 2008-07-02 21:45:34 UTC (rev 12369)
@@ -0,0 +1,133 @@
+## /Applications/Plotting/Cubit_10_2.app/Contents/MacOS/Cubit_10_2
+## Cubit Version 10.2
+## Cubit Build 8
+## Revised 10/17/2006 15:13:41 MST
+## Running 06/16/2008 02:47:47 PM
+## Command Options:
+## -warning = On
+## -information = On
+
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+
+# Block is 216 km x 216 km x 216 km
+# -108 km <= x <= 108 km
+# -108 km <= y <= 108 km
+# -216 km <= z <= 0 km
+reset
+brick x 216000 y 216000 z 216000
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -108000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+surface 7 name "fault_surface"
+create planar surface with plane zplane offset -27000
+surface 8 name "material_interface"
+
+# ----------------------------------------------------------------------
+# Divide volumes using interface surfaces
+# ----------------------------------------------------------------------
+webcut volume 1 with plane fault_surface
+webcut volume 1 with plane material_interface
+webcut volume 4 with plane material_interface
+volume 1 name "elastic_xpos"
+volume 4 name "elastic_xneg"
+volume 5 name "visco_xpos"
+volume 6 name "visco_xneg"
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Set discretization size and generate the mesh
+# ----------------------------------------------------------------------
+volume all size 6750
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 4
+block 1 name "elastic"
+block 2 volume 5 6
+block 2 name "viscoelastic"
+
+# ----------------------------------------------------------------------
+# Create fault nodeset
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at A
+nodeset 10 fault
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 20
+group "face_xpos" add node in surface 28
+nodeset 11 face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 30
+group "face_xneg" add node in surface 38
+nodeset 12 face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 21
+group "face_ypos" add node in surface 33
+group "face_ypos" add node in surface 27
+group "face_ypos" add node in surface 35
+nodeset 13 face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 23
+group "face_yneg" add node in surface 31
+group "face_yneg" add node in surface 25
+group "face_yneg" add node in surface 37
+nodeset 14 face_yneg
+nodeset 14 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 12
+group "face_zneg" add node in surface 16
+nodeset 15 face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face without fault
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 16 face_zneg_nofault
+nodeset 16 name "face_zneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 10
+group "face_zpos" add node in surface 17
+nodeset 17 face_zpos
+nodeset 17 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "spbm_hex8_unif_6750m.exo" dimension 3 overwrite
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_8000m.exo
===================================================================
(Binary files differ)
Deleted: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_8000m.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_8000m.jou 2008-07-02 01:17:36 UTC (rev 12368)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_unif_8000m.jou 2008-07-02 21:45:34 UTC (rev 12369)
@@ -1,133 +0,0 @@
-## /Applications/Plotting/Cubit_10_2.app/Contents/MacOS/Cubit_10_2
-## Cubit Version 10.2
-## Cubit Build 8
-## Revised 10/17/2006 15:13:41 MST
-## Running 06/16/2008 02:47:47 PM
-## Command Options:
-## -warning = On
-## -information = On
-
-# ----------------------------------------------------------------------
-# Create block
-# ----------------------------------------------------------------------
-
-# Block is 256 km x 256 km x 256 km
-# -128 km <= x <= 128 km
-# -128 km <= y <= 128 km
-# -256 km <= z <= 0 km
-reset
-brick x 256000 y 256000 z 256000
-
-# Translate block so the top is at z=0
-volume 1 move x 0 y 0 z -128000
-
-# ----------------------------------------------------------------------
-# Create interface surfaces
-# ----------------------------------------------------------------------
-create planar surface with plane xplane offset 0
-surface 7 name "fault_surface"
-create planar surface with plane zplane offset -32000
-surface 8 name "material_interface"
-
-# ----------------------------------------------------------------------
-# Divide volumes using interface surfaces
-# ----------------------------------------------------------------------
-webcut volume 1 with plane fault_surface
-webcut volume 1 with plane material_interface
-webcut volume 4 with plane material_interface
-volume 1 name "elastic_xpos"
-volume 4 name "elastic_xneg"
-volume 5 name "visco_xpos"
-volume 6 name "visco_xneg"
-
-# ----------------------------------------------------------------------
-# Imprint all volumes, merging surfaces
-# ----------------------------------------------------------------------
-imprint all with volume all
-merge all
-
-# ----------------------------------------------------------------------
-# Set discretization size and generate the mesh
-# ----------------------------------------------------------------------
-volume all size 8000
-mesh volume all
-
-# ----------------------------------------------------------------------
-# Create blocks for materials
-# ----------------------------------------------------------------------
-block 1 volume 1 4
-block 1 name "elastic"
-block 2 volume 5 6
-block 2 name "viscoelastic"
-
-# ----------------------------------------------------------------------
-# Create fault nodeset
-# ----------------------------------------------------------------------
-group "fault" add node in fault_surface
-group "fault" add node in fault_surface at A
-nodeset 10 fault
-
-# ----------------------------------------------------------------------
-# Create nodeset for +x face
-# ----------------------------------------------------------------------
-group "face_xpos" add node in surface 20
-group "face_xpos" add node in surface 28
-nodeset 11 face_xpos
-nodeset 11 name "face_xpos"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -x face
-# ----------------------------------------------------------------------
-group "face_xneg" add node in surface 30
-group "face_xneg" add node in surface 38
-nodeset 12 face_xneg
-nodeset 12 name "face_xneg"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face
-# ----------------------------------------------------------------------
-group "face_ypos" add node in surface 21
-group "face_ypos" add node in surface 33
-group "face_ypos" add node in surface 27
-group "face_ypos" add node in surface 35
-nodeset 13 face_ypos
-nodeset 13 name "face_ypos"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face
-# ----------------------------------------------------------------------
-group "face_yneg" add node in surface 23
-group "face_yneg" add node in surface 31
-group "face_yneg" add node in surface 25
-group "face_yneg" add node in surface 37
-nodeset 14 face_yneg
-nodeset 14 name "face_yneg"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -z face
-# ----------------------------------------------------------------------
-group "face_zneg" add node in surface 12
-group "face_zneg" add node in surface 16
-nodeset 15 face_zneg
-nodeset 15 name "face_zneg"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -z face without fault
-# ----------------------------------------------------------------------
-group "face_zneg_nofault" add node in face_zneg
-group "face_zneg_nofault" remove node in fault
-nodeset 16 face_zneg_nofault
-nodeset 16 name "face_zneg_nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +z face
-# ----------------------------------------------------------------------
-group "face_zpos" add node in surface 10
-group "face_zpos" add node in surface 17
-nodeset 17 face_zpos
-nodeset 17 name "face_zpos"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "spbm_hex8_unif_8000m.exo" dimension 3 overwrite
Modified: short/2.5D/benchmarks/savageprescott/parameters/finalslip_rupture.spatialdb
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/finalslip_rupture.spatialdb 2008-07-02 01:17:36 UTC (rev 12368)
+++ short/2.5D/benchmarks/savageprescott/parameters/finalslip_rupture.spatialdb 2008-07-02 21:45:34 UTC (rev 12369)
@@ -26,6 +26,6 @@
// (5) reverse-slip (m)
// (6) fault-opening (m)
0.0 0.0 0.00 -4.0 0.0 0.0
-0.0 0.0 -16.00 -4.0 0.0 0.0
-0.0 0.0 -16.01 0.0 0.0 0.0
-0.0 0.0 -256.00 0.0 0.0 0.0
+0.0 0.0 -13.50 -4.0 0.0 0.0
+0.0 0.0 -13.51 0.0 0.0 0.0
+0.0 0.0 -486.00 0.0 0.0 0.0
Modified: short/2.5D/benchmarks/savageprescott/parameters/pylithapp.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/pylithapp.cfg 2008-07-02 01:17:36 UTC (rev 12368)
+++ short/2.5D/benchmarks/savageprescott/parameters/pylithapp.cfg 2008-07-02 21:45:34 UTC (rev 12369)
@@ -20,6 +20,11 @@
#quadrature3d = 1
#fiatlagrange = 1
+[pylithapp.journal.debug]
+pylithapp = 1
+problem = 1
+implicit = 1
+
# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
Modified: short/2.5D/benchmarks/savageprescott/parameters/sliprate_creep.spatialdb
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/sliprate_creep.spatialdb 2008-07-02 01:17:36 UTC (rev 12368)
+++ short/2.5D/benchmarks/savageprescott/parameters/sliprate_creep.spatialdb 2008-07-02 21:45:34 UTC (rev 12369)
@@ -27,6 +27,6 @@
// (5) reverse-slip rate (m/s)
// (6) fault-opening rate (m/s)
0.0 0.0 0.00 0.00 0.0 0.0
-0.0 0.0 -16.00 0.00 0.0 0.0
-0.0 0.0 -16.01 -6.337617562806e-10 0.0 0.0
-0.0 0.0 -256.00 -6.337617562806e-10 0.0 0.0
+0.0 0.0 -13.50 0.00 0.0 0.0
+0.0 0.0 -13.51 -6.337617562806e-10 0.0 0.0
+0.0 0.0 -486.00 -6.337617562806e-10 0.0 0.0
Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.cfg (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.cfg 2008-07-02 21:45:34 UTC (rev 12369)
@@ -0,0 +1,127 @@
+# -*- Python -*-
+[pylithapp]
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to the boundary and
+# interface conditions. The general parameters are specificed in the
+# pylithapp.cfg file which PyLith reads by default. Information controlling
+# time stepping is contained in the file timedep.cfg.
+#
+# To run the simulation:
+# pylith spbm_hex8_unif_8000m.cfg timedep.cfg
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+# Change the mesh importer to the CUBIT importer.
+importer = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.importer]
+# Set filename of mesh to import.
+filename = ../meshes/spbm_hex8_graded_1000m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
+bc = [x_pos,x_neg,z_neg]
+
+# Set interfaces to an array of 1 fault: 'fault'.
+interfaces = [fault]
+
+[pylithapp.timedependent.implicit]
+# Set the output to an array of 2 output managers.
+# We will output the solution over the domain and the ground surface.
+output = [domain,subdomain]
+
+# Set subdomain component to OutputSolnSubset (subset of domain).
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+
+# We change spatial database for the rate of change from FixedDOFDB
+# (which has zero velocity) to UniformDB (which allows nonzero
+# velocities).
+
+# Set the parameters for the desired boundary conditions.
+#
+# The label corresponds to the nodeset ID in CUBIT.
+
+[pylithapp.timedependent.bc.x_pos]
+fixed_dof = [0, 1]
+label = 11
+db.label = Dirichlet BC on +x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on +x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of -1 cm/year.
+rate_db.data = [0.0,-3.1688087814e-10]
+
+[pylithapp.timedependent.bc.x_neg]
+fixed_dof = [0, 1]
+label = 12
+db.label = Dirichlet BC on -x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on -x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of 1 cm/year.
+rate_db.data = [0.0,3.1688087814e-10]
+
+[pylithapp.timedependent.bc.z_neg]
+fixed_dof = [2]
+label = 16
+db.label = Dirichlet BC on -z
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+
+# Set the parameters for the fault interface conditions.
+
+[pylithapp.timedependent.interfaces.fault]
+label = 10
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+mat_db.iohandler.filename = mat_elastic.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain.writer]
+filename = ../results/spbm_hex8_graded_1000m.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK domain output of solution over ground surface.
+[pylithapp.problem.formulation.output.subdomain]
+label = 17 ; nodeset for subdomain
+writer.filename = ../results/spbm_hex8_graded_1000m-groundsurf.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for vtk fault rupture output.
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = ../results/spbm_hex8_graded_1000m-fault.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_graded_1000m-statevars-elastic.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.materials.maxwell.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_graded_1000m-statevars-maxwell.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.cfg (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.cfg 2008-07-02 21:45:34 UTC (rev 12369)
@@ -0,0 +1,127 @@
+# -*- Python -*-
+[pylithapp]
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to the boundary and
+# interface conditions. The general parameters are specificed in the
+# pylithapp.cfg file which PyLith reads by default. Information controlling
+# time stepping is contained in the file timedep.cfg.
+#
+# To run the simulation:
+# pylith spbm_hex8_unif_8000m.cfg timedep.cfg
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+# Change the mesh importer to the CUBIT importer.
+importer = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.importer]
+# Set filename of mesh to import.
+filename = ../meshes/spbm_hex8_graded_3000m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
+bc = [x_pos,x_neg,z_neg]
+
+# Set interfaces to an array of 1 fault: 'fault'.
+interfaces = [fault]
+
+[pylithapp.timedependent.implicit]
+# Set the output to an array of 2 output managers.
+# We will output the solution over the domain and the ground surface.
+output = [domain,subdomain]
+
+# Set subdomain component to OutputSolnSubset (subset of domain).
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+
+# We change spatial database for the rate of change from FixedDOFDB
+# (which has zero velocity) to UniformDB (which allows nonzero
+# velocities).
+
+# Set the parameters for the desired boundary conditions.
+#
+# The label corresponds to the nodeset ID in CUBIT.
+
+[pylithapp.timedependent.bc.x_pos]
+fixed_dof = [0, 1]
+label = 11
+db.label = Dirichlet BC on +x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on +x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of -1 cm/year.
+rate_db.data = [0.0,-3.1688087814e-10]
+
+[pylithapp.timedependent.bc.x_neg]
+fixed_dof = [0, 1]
+label = 12
+db.label = Dirichlet BC on -x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on -x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of 1 cm/year.
+rate_db.data = [0.0,3.1688087814e-10]
+
+[pylithapp.timedependent.bc.z_neg]
+fixed_dof = [2]
+label = 16
+db.label = Dirichlet BC on -z
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+
+# Set the parameters for the fault interface conditions.
+
+[pylithapp.timedependent.interfaces.fault]
+label = 10
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+mat_db.iohandler.filename = mat_elastic.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain.writer]
+filename = ../results/spbm_hex8_graded_3000m.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK domain output of solution over ground surface.
+[pylithapp.problem.formulation.output.subdomain]
+label = 17 ; nodeset for subdomain
+writer.filename = ../results/spbm_hex8_graded_3000m-groundsurf.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for vtk fault rupture output.
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = ../results/spbm_hex8_graded_3000m-fault.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_graded_3000m-statevars-elastic.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.materials.maxwell.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_graded_3000m-statevars-maxwell.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
Modified: short/2.5D/benchmarks/savageprescott/utils/savpres_ss.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/utils/savpres_ss.cfg 2008-07-02 01:17:36 UTC (rev 12368)
+++ short/2.5D/benchmarks/savageprescott/utils/savpres_ss.cfg 2008-07-02 21:45:34 UTC (rev 12369)
@@ -2,8 +2,8 @@
[savpres_ss]
# Top-level info
-elas_thick = 32.0*km
-lock_depth = 16.0*km
+elas_thick = 27.0*km
+lock_depth = 13.5*km
recurrence_time = 200.0*year
viscosity = 4.73364e+19*Pa*s
shear_modulus = 30.0*GPa
More information about the cig-commits
mailing list