[CIG-SHORT] Cubit tet mesh not being read into pylith - journal files

[Brad Aagaard]

Oliver-

You never cut out the small box from your big box. You need to add the 
following 2 lines to your geometry.jou file after deleting surface 20.

subtract volume 2 4 from volume 6 keep
delete volume 6

Your revised journal file should look like what I have below.

NOTE TO ALL USERS: There are many different ways to create the solid 
model geometry in LaGriT and CUBIT. Regardless of how you create your 
geometry, you need to make sure the volumes do *not* overlap.

Brad


# ----------------------------------------------------------------------
# Create block
# ----------------------------------------------------------------------
# Block is 1000 km x 1000 km x 100 km
# -500 km <= x <= 500 km
# -500 km <= y <= 500 km
# -100 km <= z <= 0 km
reset
brick x 1000000 y 1000000 z 100000
volume 1 move x 0 y 0 z -50000
# ----------------------------------------------------------------------
# Create fault surface and subdivide fault volume
# ----------------------------------------------------------------------
create vertex 0 -65e3 -20e3
create vertex 0 65e3 -20e3
create vertex 0 65e3 -40e3
create vertex 0 -65e3 -40e3

create surface vertex 9 to 12

sweep surface 7 vector 1 0 0 distance 100e3
sweep surface 12 vector -1 0 0 distance 100e3
create surface vertex 9 to 12

webcut volume 2 with plane surface 17

delete vertex 9 to 12
delete surface 17

# ----------------------------------------------------------------------
# Create material surface and subdivide material volume
# ----------------------------------------------------------------------
create planar surface with plane zplane offset -20000
webcut volume 1 with plane surface 20
delete surface 20
subtract volume 2 4 from volume 6 keep
delete volume 6

# ----------------------------------------------------------------------
# Name volumes and surfaces
# ----------------------------------------------------------------------
volume 1 name "elastic"
volume 2 4 6 name "visco"
surface 18 name "fault_surface"

# ----------------------------------------------------------------------
# Imprint all volumes, merging surfaces
# ----------------------------------------------------------------------
imprint all with volume all
merge all