The atmospheric boundary layer is characterized by a logarithmic increase in wind speed within the surface layer and a turning of the wind vector aloft. The simplest canonical configuration featuring both a logarithmic layer and wind turning is Ekman flow. A formulation of the mean velocity profile was derived from scaling considerations and calibrated using direct numerical simulation (DNS) in Part I of this work. Here, we explore the extrapolation of this formulation to atmospheric Reynolds numbers (Re) using large-eddy simulation (LES). Theoretical profiles of the wind vector are compared to LES results at intermediate Re, which require consideration of viscous effects along with modifications to the standard bottom boundary condition. A grid convergence analysis shows that LES data converge towards the theoretical profiles for intermediate and high Re. The LES thus confirms that (i) the spanwise velocity scales as Reτ −1, where Reτ is the friction Reynolds number, and (ii) one third of the wind veer is confined to the surface layer. Such agreement between the theoretical formulation and LES data increases confidence in the underlying scaling assumptions, reinforcing the utility of the theoretical profiles as a reference for intermediate and a quasi-reference for both idealized simulation and field observation.