Add domain boundary detection to Poincare event detection

When using integrate_line to find closed orbits with Poincare return
map, detect if line integration leaves flow domain.
Use standard eig function for ocean dataset demo