The shallow-water equations are solved using the "ocean" version of Gerris. The tidal elevations for the lunar (M2) component obtained from a larger-area tidal model are imposed as conditions on the boundaries of the domain.
The comments in the tides.sh script describe how to generate the appropriate GTS files from the tidal elevation and bathymetry data.
After an initial transient (t < ≈ 1 day) due to relaxation of the model toward a state consistent with the mathematical model and with the imposed boundary conditions, the model reaches a periodic regime (Figure 29).
Online harmonic decomposition can then be used to extract the amplitudes and phases of the computed M2 tidal components. The simulation stops automatically when convergence of the harmonic decomposition is reached (Figure 30).
The final tidal amplitudes and phases are illustrated in Figures 31 and 32 respectively. The harmonic decomposition is also applied to the velocity field. The results can be represented as tidal ellipses (Figure 33) and residual currents (Figure 34).
Note that the results for this simulation will not be as good as these described in Rym Msadek’s technical report because iterative Flather conditions have not been applied. See the report for details.