A general inflow turbulence synthetic method for large eddy simulation (LES) involved in parallel calculation was proposed. Firstly, fluctuating velocities of some main grid points on the inflow boundary were generated adopting POD based spectral representation method. Secondly, velocity timeseries of all grid points on the inflow panel were spatially interpolated, adopting the finite element shape functions. The synthesized velocities were then read in commercial code Fluent by user-defined-functions (UDF) program with special consideration of parallel calculation. In order to verify the efficiency of the present method, the flow around a single square section tall building with width to height aspect ratio of 1:6 was numerically simulated by adopting the present inflow turbulence generation method for LES. The simulated wind profiles, wind spectra, statistics, self-spectra, and coherence of wind pressure coefficients, were compared with both the experimental and numerical results of existing works in detail. The present inflow turbulence generation method for LES was verified as being applicable to a certain extent in both simulating turbulent wind field and predicting wind loads on buildings.