Hybrid mhd-gyrokinetic simulations of shear alfvén modes in tokamak plasmas