Kinetic shear Alfvn modes in tokamak plasmas (with or without temperature gradient) in the full gyrokinetic limit are numerically investigated by using the integral eigenvalue equations.