The time discretization used in TELEMAC2D is semi- implicit, with a Crank-Nicolson factor of 0.6. The equations are solved using Method of Characteristics, PSI distributive scheme and N distributive scheme for advection for flow velocities, water depth and tracers respectively. As for the PSI scheme, the Courant–Friedrichs–Lewy condition (CFL condition) must be considered for convergence, and commonly the Courant number should be less than 1. Therefore, in this study, the time step is set