A numerical method for bottom interacting ocean acoustic normal modes