We present results from a high resolved Large Eddy Simulation of a freely developing Blasius profile over a concave boundary in a large spanwise domain. Due to the large initial Reynolds and Görtler numbers (Re θ ,0 = 1170, G θ ,0 = 75), we observe the onset of two dominant wavelengths, the first dominating in the linear/transition region, λ 1 , and the second dominating in the turbulent region, λ 2 . Extending previous linear stability analysis (LSA) to higher Görtler numbers and non-dimensional wavenumbers, both dominant wavelengths of the Görtler instability correspond to predictions of LSA, the latter comparable to laminar theory by replacing the kinematic viscosity with the turbulent viscosity in the definition of the Görtler number. The predicted spatial modes compare well with the computed profiles for both λ 1 and λ 2 . The skin friction coefficient C f is found heterogeneous in the spanwise direction according to the emerging wavelengths λ 1 and λ 2 of the Görtler instability. We report a smooth increase of C f from the theoretical predictions of a laminar boundary layer to those for a turbulent boundary layer over a flat plate with a slight overshoot in the domain of existence of the second dominant wavelength λ 2 .