# macro infile for finding fault and repo
# make sure correct mo is current

# make sure sets do not exist
# Ignore errors on first level of refinement
eltset/erefine/delete
eltset/eboth/delete
eltset/elevel/delete
eltset/ebox/delete
eltset/eflt/delete
region/rflt/delete
pset/pbox/delete
pset/pflt/delete

region / rflt / ge sleft and le sright
pset / pflt / region / rflt
eltset / eflt / inclusive / pset,get,pflt

pset/pbox/geom/xyz/1,0,0/ XMIN YMIN ZMIN / XMAX YMAX ZMAX
eltset/ebox/inclusive/pset,get,pbox
eltset/eboth/ union ebox eflt

# narrow refine area to current LEVEL

eltset/elevel/ itetlev / le / LEVEL
eltset/erefine/ inter eboth elevel

finish