cmo/create/cube1/tet/
*
* Define surfaces
*
surface/outbox/reflect/box/-64.0,-64.0,-64.0/64.0,64.0,0.0/
surface/fltbot/intrface/plane/0.0,0.0,-8.0/4.0,0.0,-8.0/4.0,4.0,-8.0/
surface/inbox/intrface/box/-8.0,-4.0,-8.0/8.0,4.0,0.0/
surface/f12_plane/intrface/plane/0.0,0.0,-4.0/0.0,-4.0,-4.0/-4.0,-4.0,-8.0/
*
* Build triangulated surface that divides mb1 and mb2
* Creates surface f12
* surface/f12/intrface/sheet/CMONAME
*
infile build_plane.lgi
*
* Define regions bounded by surfaces
*
region/b1/le inbox and le f12/
region/b2/le inbox and ge f12/
region/b3a/ge inbox and le outbox and le f12_plane and ge fltbot/
region/b3b/ge inbox and le outbox and ge f12_plane and ge fltbot/
region/b4/ge inbox and le outbox and le fltbot/
*
* Define materials bounded by surfaces
*
mregion/mb1/lt inbox  and le outbox and lt f12/
mregion/mb2/lt inbox  and le outbox and gt f12/
mregion/mb3a/gt inbox and le outbox and lt f12_plane and gt fltbot/
mregion/mb3b/gt inbox and le outbox and gt f12_plane and gt fltbot/
mregion/mb4/gt inbox and le outbox and lt fltbot/
*
* Create points at a resolution of 0.5 km.
*
createpts/xyz/65,65,33/-16.0,-16.0,-16.0/16.0,16.0,0.0/1,1,1/
*
* Create points at a resolution of 1.0 km.
*
createpts/xyz/49,49,25/-24.0,-24.0,-24.0/24.0,24.0,0.0/1,1,1/
*
* Create points at a resolution of 2.0 km.
*
createpts/xyz/33,33,17/-32.0,-32.0,-32.0/32.0,32.0,0.0/1,1,1/
*
* Create points at a resolution of 4.0 km.
*
createpts/xyz/33,33,17/-64.0,-64.0,-64.0/64.0,64.0,0.0/1,1,1/
*
* With the non-planar fault, copy the points from the surface
* defining the fault.
*
copypts / cube1 / cmo_flt_surf
*
* Remove duplicate points
*
filter/1,0,0/
rmpoint/compress
cmo / setatt / -def- / imt / 1 0 0 / 0
cmo / setatt / -def- / itp / 1 0 0 / 0
*----------------------------------
* Subset the point distribution to speed up
* calculations during prototype runs.
*
* define / XMIN / -10.0
* define / XMAX /  10.0
* define / YMIN /  -5.0
* define / YMAX /   5.0
* define / ZMIN / -16.0
* define / ZMAX /   0.0
* define / XC   /   0.0
* define / XC   /   0.0
* define / ZC   /   0.0
* infile subset_cmo.mlgi
*----------------------------------
* Set node color (imt) and type (itp) based on region and mregion
*
setpts
dump / gmv / tmp_setpts.gmv / -def-
*
* Connect points into Delaunay tet mesh, without connecting across
* material interfaces, and adding points if needed to resolve
* interfaces.
*
connect/check_interface
*
* Set element colors.
*
settets

dump / gmv / tmp_settets.gmv / -def-

* Reset parent points and remove duplicate points on interfaces.
*
resetpts/parent
rmpoint/compress
*
* Define external boundaries using dump/outside_zone command.
* Generate attributes: top, bottom, left_w, right_e, back_n, front_s

dump / zone_outside / dipblock4 / -def- / keepatt
*
*
*
* Dump results to different file formats
*
dump/gmv/dipblock4.gmv/cube1/ascii
dump/avs/dipblock4.inp/cube1/
dump/lagrit/dipblock4.lagrit
quality
finish