Metasurface Modeling for the Manipulation of Goos–Hänchen and Imbert–Fedorov Shifts