Computational model has been developed to simulate wind loads on a large container ship for different coastal terrains and wind incidence angles. Wind characteristics were adopted to correspond to the atmospheric boundary layer flow. Computations were performed using the software package based on solving Reynolds-Averaged-Navier-Stokes equations using the k-ε turbulence model. Prior to testing wind loading of a large container ship, wind flow through an empty computational domain without the ship was investigated in order to determine changes in the velocity, turbulent kinetic energy, and dissipation of turbulent kinetic energy profiles from inlet to the ship position. For all coastal terrain types, an increase in velocity and a decrease in turbulent kinetic energy between an inlet and the ship position are observed, while differences in dissipation of turbulent kinetic energy profiles are negligible. Differences in coastal terrain type influence drag coefficients. In particular, drag coefficients CD and CL increase with an increase in coastal urbanization, while moment coefficients CM are nearly the same for all terrain types. As the wind incidence angle increases, absolute values for CD decrease and absolute values for CL increase, while maximum values for CM are observed at 45° and 90° wind incidence angles. Computational model developed in this study has been validated against analytical expressions based on wind-tunnel experiments. Differences in computational and analytical results are possibly due to different inflow conditions in those two approaches, where wind characteristics in this study correspond to realistic atmospheric conditions.