Linear gyrokinetic simulations of zonal flows in toroidal rotating plasmas