Ekman layer equation with rough-bottom boundary condition has been solved to determine the effect of roughness on the magnitude of Ekman veering. The bottom boundary condition and the flow field were expanded in a power of roughness (h) which is always smaller than the Ekman layer thickness (${delta}_E$). By changing the magnitude of roughness parameter (h/${delta}_E$), the magnitude of the veering, which rotates counterclockwise from the interior geostrophic flow, has been computed. At a fixed depth within the Ekman layer, the magnitude of veering increases as the roughness parameter increases. However, the cross-isobar flux turns out to decrease with increasing roughness. To verify the analytic solution, laboratory experiments were carried out. Rough-bottom cylinderical container filled with homogeneous fluid was sit on a rotating table. The flow pattern during the period of steady spin-up shows that the degree of veering coincides well with the analytic results for various roughness parameters.