Comments:
The computational grid, named port3 for the Great Bay Estuary system is created using linear equilateral triangles fitted to the shoreline boundary. There are 8577 nodes and 14974 elements. The bathymetric depth contour reveals a center channel running south to north.

M2 tidal forcing is specified as a Dirichlet elevation boundary condition across the southern and northern inlet boundariesy. The amplitude is 1.35m at the southern boundary and 1.00m at the northern boundary. The simulation was started with fluid at rest and was terminated after six M2 tidal period (74.52hrs). The simulation setup is given below. All the simulation parameters are summarized in Table below.

Residual transport:

Comments:

Comments:
The above transient solutions are numerically well-behaved and appear to describe the realistic dynamics of the tidal flooding and dewatering process in this complex domain.