Add feature in momentumaccpetance( ), so user can define the initial vertical cooridnate which is used to search 6D closed orbit.