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

[Oliver Boyd]

Excellent! Thanks for your help Brad - worked like a charm.


On 10/22/09 5:54 PM, "Brad Aagaard" <baagaard at usgs.gov> wrote:

> 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