* Read 
read/gmv/s_r_mass1a.gmv/cmo_1
resetpts / parents
rmpoint / compress
cmo/setatt/cmo_1/zic/1 0 0/0.0
cmo/setatt/cmo_1/itetclr/1 0 0/1
dump/avs/cmo_1.inp/cmo_1
cmo / delete / cmo_1

read/gmv/s_r_mass1a.gmv/cmo_2
resetpts / parents
rmpoint / compress
cmo/setatt/cmo_2/zic/1 0 0/2.E+04
cmo/setatt/cmo_2/itetclr/1 0 0/1
dump/avs/cmo_2.inp/cmo_2
cmo / delete / cmo_2

* ----------------------------------------------------
* STACKING LAYERS
* ----------------------------------------------------

* Create stack object
cmo/create/cmo_stack

* Stack surfaces on top of each other starting from last one
stack/layers/avs/ &
cmo_1.inp 1/ &
cmo_2.inp 1 10

* Fill with prisms
stack/fill/cmo_prism/cmo_stack
cmo/setatt/cmo_prism/ipolydat/no

* Break prisms into tetrahedra
hextotet//cmo_tet/cmo_prism

mmprint
cmo status

cmo / delete / cmo_prism

* Reset
resetpts/itp

* Dump resulting tet mesh to a gmv and avs file
cmo/setatt/cmo_tet/ipolydat/no
dump/avs/cmo_tet.inp/cmo_tet
dump/gmv/cmo_tet.gmv/cmo_tet

* Boundary
* Check for holes. If mesh had holes you would have
* more than 1 boundary component.
boundary_components

define / XMAX / 200000.
define / YMAX / 200000.
define / ZMAX / 40000.
define / R1 / 1000.
define / R2 / 100000.
define / R3 / 200000.
*
* Define Geometry
*
surface / cyl_r1 / intrface / cylinder / 0. 0. 0. / 0. 0. ZMAX / R1
surface / cyl_r2 / intrface / cylinder / 0. 0. 0. / 0. 0. ZMAX / R2
surface / cyl_r3 / intrface / cylinder / 0. 0. 0. / 0. 0. ZMAX / R3

surface / b_box / reflect / box / 0. 0. 0. / XMAX YMAX ZMAX

region / r1 / le b_box and le cyl_r1
region / r2 / le b_box and ge cyl_r1 and le cyl_r2
region / r3 / le b_box and ge cyl_r2 and le cyl_r3

region / r4 / le b_box and ge cyl_r3

mregion / mr1 / le b_box and lt cyl_r1
mregion / mr2 / le b_box and gt cyl_r1 and lt cyl_r2
mregion / mr3 / le b_box and gt cyl_r2 and lt cyl_r3

mregion / mr4 / le b_box and gt cyl_r3
settets / geometry
resetpts/cell_color

dump/avs/cmo_settet.inp/cmo_tet
dump/gmv/cmo_settet.gmv/cmo_tet

finish