“…As outlined previously [11,12] in this numerical model, the rock profile, y (z, t + T) represented as a column of horizontally aligned layers of height Δz = 0.05 m, moves a quantity δy(z,T) in the points z affected by the erosive processes after one tidal period T. This quantity can be defined as the erosion of each element of the cliff front and depends on H b, which is the breaking wave height; T b, which is the wave period; K, which is a calibration term representing hydrodynamic constants (100 m 13/4 s 7/2 /kg for the application site; see Refs. [11,12]); σ c (z), which is the uniaxial compressive strength (USC) of the rock mass; ∂ z y(z,t) −1 , which is the slope of each rock element z (which changes through simulation time in response to the calculated erosion and consequently to profiles evolution); p w (z,t), which is the shape function (which include from experimental data- [41,42]-the erosion of different type waves) and w t (t), which is the tidal expression introduced as a sinusoid that oscillates about mean sealevel with an amplitude A m and period T (both can be obtained through governmental marine agencies). The rock profile evolution can be described as seen in Eq.…”