The modelling of large amplitude internal waves (LAIW) propagating in the Strait of Gibraltar is carried out using a fully-nonlinear nonhydrostatic numerical model.The vertical fluid stratification, parameters of the propagating waves and bottom topography were taken close to those observed in the strait.Amplitudes of the first mode LAIWs radiating from Camarinal Sill and propagating along the strait towards the Alboran Sea were taken in a wide range, from 10 to 90m (as observed in the studied area).The focus of the modelling efforts was on three-dimensional peculiarities of LAIWs evolution, viz.cross-strait variability, interaction with lateral boundaries (including wave breaking and water mixing), radiation of secondary waves from orographic features and interaction of secondary scattered internal waves.