<sup>1</sup> Open boundary condition:

The open boundary condition is also known as the water boundary condition, where, on this boundary, either the flow rate is given, or the time series condition of the tidal level is given. For the open boundary in this work, calculation was performed in the following form of tidal harmonic analysis:

$$\eta\_i = \sum\_{i}^{m} f\_i H\_i \cos[\omega\_i t + (V\_0 + \mu)\_i - g\_i] \tag{4}$$

where *ω<sup>i</sup>* represents the angular velocity of the *i*th tidal constituent; *fi* and *u* represent the intersection factor and epoch correction of the *i*th tidal constituent, respectively; *Hi* and *gi* represent harmonic constants, which are the amplitude and epoch of each tidal constituent, respectively; *V*<sup>0</sup> represents the time angle of a tidal constituent. The time series data of the tide level at the open boundary of the model were also verified according to the global tide module.
