File contents

*make a two-region mesh with a 45 degree fault
* create a 3-D tetrahedral mesh object and name it mymesh
cmo/create/mymesh/tet/

surface/bigbox/reflect/box/-50.0,-50.0,-100.0/50.0,50.0,0.0/
surface/innerbox/intrface/box/-10.0,-10.0,-20.0/10.0,10.0,0.0/
*surface/fault/intrface/plane/-10.0,-10.0,-20.0/10.0,-10.0,-20.0/-10.0,10.0,0.0/
surface/fault/intrface/plane/-15.0,-15.0,-25.0/15.0,-15.0,-25.0/-15.0,15.0,5.0/

region/outside/ le bigbox and gt innerbox /
*region/inside / le innerbox 
region/foot / lt fault and le innerbox
region/hang / ge fault and le innerbox

mregion/moutside/ le bigbox and gt innerbox
*mregion/minside / lt innerbox
mregion/mfoot / lt fault and lt innerbox
mregion/mhang / gt fault and lt innerbox

createpts/xyz/11,11,11/-50.,-50.,-100./50.,50.,0. / 1 1 1
createpts/xyz/11,11,11/-10.,-10., -20./10.,10.,0. / 1 1 1

*  eliminate coincident or nearly coincident points
*  1,0,0 means consider all points
filter/1,0,0/
rmpoint / compress
*
*dump / gmv / points.gmv / mymesh
*
* Set imt, itp based on mregion definitions
*
setpts
*
*dump / gmv / setpts.gmv / mymesh
*
*  connect the points into a Delaunay tetrahedral mesh
*  do not connect across material interfaces -
*  add points if necessary to resolve material interfaces

connect / check_interface

*dump / gmv / connect.gmv / mymesh
*
*  set element (tetrahedral) type (itetclr) based on mregion
*
settets
* now try smoothing
*smooth
smooth/position/aspect/
recon 0
massage/10./1./0.01/0.1/
*  dump mesh to some output form
dump/gmv/nestbox2s.gmv/mymesh/ascii/
dump / geofest / nestbox2s.gft / mymesh
* begin compare here
cmo/status
cmo/printatt//-all-/minmax
quality
finish