# Use top piece from 25m grid for top of this mesh
ln -s ../3Dgrid_octree/tet_final.inp tet_top_25m.inp
# Use 25m grid to interpolate and find elem outside 
ln -s ../3Dgrid_octree/tet_final_zones.inp tet_25m.inp

lagrit < mesh_pieces.lgi
cp outx3dgen mesh_pieces.out
cp logx3dgen mesh_pieces.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