This paper will present the development of a simple subsonic boundary layer method suitable to be used coupled with panel methods in order to estimate the aerodynamic characteristics, including viscous drag and maximum lift coefficient, of 3D wings. The proposed method does not require viscous-inviscid iterations and is based on classical integral bi-dimensional boundary layer theory using Thwaites and Head ́s models with bi-dimensional empirical corrections applied to each wing strip being therefor robust and efficient to be used in the early conceptual stage of aircraft design. Presented results are compared to the Modified CS Method in an IBL scheme and experimental data and are shown to provide good results.