lagrit < mesh_boxes.lgi
cp outx3dgen mesh_boxes.out
cp logx3dgen mesh_boxes.log

# Create the top piece with slope
#    uses infile mass_psave.mlgi
lagrit < mesh_top_tri.lgi
cp outx3dgen mesh_top_tri.out
cp logx3dgen mesh_top_tri.log

lagrit < mesh_top.lgi
cp outx3dgen mesh_top.out
cp logx3dgen mesh_top.log

# Refine bottom piece with repo and fault
#     uses  get_eltsets_max.mlgi and get_eltsets.mlgi

lagrit < mesh_octree.lgi
cp outx3dgen mesh_octree.out
cp logx3dgen mesh_octree.log

# Join pieces and write final files
lagrit < mesh_join.lgi
cp outx3dgen mesh_join.out
cp logx3dgen mesh_join.log

# Write FEHM and Exodus files

lagrit < write_fehm.lgi
cp outx3dgen write_fehm.out
cp logx3dgen write_fehm.log

lagrit < write_exo.lgi
cp outx3dgen write_exo.out
cp logx3dgen write_exo.log

/home/tam/bin/ncdump_lin64 tet_exodus.exo > tet_exodus_ascii.exo