###
#Note: This mesh is REALLY coarse and has some funny-shaped cells.  Just for demo purposes
#--------------------------------- 
# Make volumes 
reset 
brick x 60000.000000 y 60000.000000 z 30000.000000 
move volume 1 x 0 y 0 z -15000.000000 
brick x 300000.000000 y 300000.000000 z 200000.000000 
move volume 2 x 0 y 0 z -100000.000000 
chop volume 2 with volume 1 
create vertex x 0.000000 y 0.000000 z -5000.000000 
create vertex x 8660.254038 y 5000.000000 z 0.000000 
create vertex x 8660.254038 y -5000.000000 z 0.000000 
rotate -45 about X
rotate -45 about Y
rotate 45 about Z
pause
webcut volume all with plane vertex 29 to 31 

pause

delete vertex 29 to 31 
# Make inner fault 
create vertex 2078.460969 -4900.000000 -3800.000000 
create vertex 2078.460969 4900.000000 -3800.000000 
create vertex -2078.460969 4900.000000 -6200.000000 
create vertex -2078.460969 -4900.000000 -6200.000000 
zoom volume 3
pause
create surface vertex all except in volume all 
delete vertex all except in volume all 
pause
# Make fault transition zone 
create vertex 2251.666050 -4900.000000 -3700.000000 
create vertex 2251.666050 4900.000000 -3700.000000 
create vertex 2200.935456 5041.421356 -3729.289322 
create vertex 2078.460969 5100.000000 -3800.000000 
create vertex -2078.460969 5100.000000 -6200.000000 
create vertex -2200.935456 5041.421356 -6270.710678 
create vertex -2251.666050 4900.000000 -6300.000000 
create vertex -2251.666050 -4900.000000 -6300.000000 
create vertex -2200.935456 -5041.421356 -6270.710678 
create vertex -2078.460969 -5100.000000 -6200.000000 
create vertex 2078.460969 -5100.000000 -3800.000000 
create vertex 2200.935456 -5041.421356 -3729.289322 
create surface vertex all except in volume all 
delete vertex all except in volume all 
pause
imprint volume 3 5 with volume 7 8 
imprint volume 3 5 with curve all except in volume all 
delete volume 7 8 
merge volume all 
# Rename Volumes 
volume 3 name "foot_inner" 
volume 4 name "foot_outer" 
volume 5 name "hang_inner" 
volume 6 name "hang_outer" 
surface 49 name "fault_trans" 
surface 50 name "fault_inner" 
# Meshing 
curve all size 30000.000000 
curve all in foot_inner hang_inner size 2000.000000 
curve all in surface fault_trans fault_inner size 200.000000 
curve 64 scheme bias fine size 2000.000000 coarse size 60000.000000 start vertex 34 
curve 60 scheme bias fine size 2000.000000 coarse size 60000.000000 start vertex 32 
mesh curve all
preview mesh curve all
pause 
surface all scheme tridelaunay 
mesh surface fault_inner
surface all sizing function type skeleton
mesh surface all in foot_inner
volume hang_inner hang_outer visibility off
display
zoom volume foot_inner
pause

mesh surface all in foot_outer
zoom foot_outer
pause

mesh surface all
block 1 volume hang_inner hang_outer 
block 1 name "hanging_walls" 

group "fault" add node in fault_inner 
nodeset 10 group fault 
nodeset 10 name "fault" 
group "face_xpos" add node in surface 12 
nodeset 11 group face_xpos 
nodeset 11 name "face xpos"