## 2 Materials and Methods
### Hydrodynamic Model  
This research utilizdes the MIKE 21/3 coupled FM model, which includes the MIKE 21 Hydrodynamic Model (HD) and Spectral Wave (SW) and Sediment Transport (ST) modules to simulate hydrodynamics, wave characteristics, and sediment transport [24; 37]. The MIKE 21 module simulates variations in water elevation and discharge to consider the response functions of forces in rivers, lakes, estuaries, and coastal areas. It employs a numerical solution to the full, time-dependent, and nonlinear depth-averaged Navier-Stokes equations, based on the conservation of mass and momentum [25; 37]. Depth and water velocity were computed in a 2-dimensional grid using simplified equations, as outlined in Equations (2)-(4) below. \[\frac{\partial\zeta}{\partial t}+\frac{\partial p}{\partial x}+\frac{\partial q }{\partial y}=\frac{\partial d}{\partial t} \tag{2}\]  
\[\frac{\partial p}{\partial t}+\frac{\partial}{\partial x}\left(\frac{p^{2}}{h} \right)+\frac{\partial}{\partial y}\left(\frac{pq}{h}\right)+gh\frac{\partial \zeta}{\partial x}+\frac{gp\sqrt{p^{2}+q^{2}}}{C^{2}h^{2}}-\frac{1}{\rho_{w}} \left[\frac{\partial}{\partial x}(h\tau_{xx})+\frac{\partial}{\partial y} \big{(}h\tau_{xy}\big{)}\right]-\Omega q-fVV_{x}+\frac{h}{\rho_{w}}\frac{ \partial}{\partial x}(P_{a})=0 \tag{3}\]  
\[\frac{\partial q}{\partial t}+\frac{\partial}{\partial y}\left(\frac{q^{2}}{h} \right)+\frac{\partial}{\partial x}\left(\frac{pq}{h}\right)+gh\frac{\partial \zeta}{\partial y}+\frac{gq\sqrt{p^{2}+q^{2}}}{C^{2}h^{2}}-\frac{1}{\rho_{w}} \left[\frac{\partial}{\partial y}\big{(}h\tau_{yy}\big{)}+\frac{\partial}{ \partial x}\big{(}h\tau_{xy}\big{)}\right]+\Omega p-fVV_{y}+\frac{h}{\rho_{w}} \frac{\partial}{\partial y}(p_{a})=0 \tag{4}\]  
where \(h(x,y,t)\) represents the water depth (m); \(\zeta(x,y,t)\) signifies the time-varying water depth (m); \((x,y,t)\) denotes the surface elevation (m); \(p\),\(q(x,y,t)\) stand for flux density in the x, y directions (m\({}^{3}\)/s/m); \(C(x,y)\) is Chezy's resistance (m\({}^{1}/^{2}/^{2}\)/s); \(g\) is the acceleration due to gravity (m/s\({}^{2}\));\(f(V)\) represents the wind friction factor; \(V,Vx,Vy(x,y,t)\) are wind velocity and components in the x, y directions (m/s); \(\Omega(x,y,t)\) denotes the Coriolis parameter, latitude-dependent (s\({}^{-1}\)); \(pa(x,y,t)\) indicates atmospheric pressure (kg/m/s\({}^{2}\)); \(\rho w\) is the density of water (kg/m\({}^{3}\)); \(x,y\) are spatial coordinates in meters; \(t\) signifies time (s); and  
\begin{table} \begin{tabular}{c c c} \hline \hline \multirow{2}{*}{**Point Name**} & \multicolumn{2}{c}{**Coordinates**} \\ \cline{2-3}  & **Latitude** & **Longitude** \\ \hline ADCP & \(3^{\circ}53^{\prime}41.90^{\prime\prime}\) S & \(102^{\circ}16^{\prime}52.20^{\prime\prime}\) E \\ Tide & \(3^{\circ}54^{\prime}8.14^{\prime\prime}\) S & \(102^{\circ}17^{\prime\prime}28.65^{\prime\prime}\) E \\ Wind & \(3^{\circ}53^{\prime}30.50^{\prime\prime}\) S & \(102^{\circ}18^{\prime}0.39^{\prime\prime}\) E \\ \hline \hline \end{tabular}   
* The classification of tidal types is determined based on the Formzahl number divided according to the following criteria [35; 36]:
* Semi-Diurnal Tide: If the Formzahl (\(F\)) value is 0.00 \(<F\leq\) 0.25;
* Mixed Tide, Prevailing Semi-Diurnal: If the Formzahl (\(F\)) value is 0.25 \(<F\leq\) 1.5;
* Mixed Tide, Prevailing Diurnal: If the Formzahl (\(F\)) value is 1.5 \(<F\leq\) 3.0;
* Diurnal Tide: If the Formzahl (\(F\)) value is \(F>\) 3.0. This value can be calculated using the following Equation (1):  
\[F=\frac{(K1+O1)}{(M2+S2)} \tag{1}\]  
Here, \(F\) represents the Formzahl number; \(K1\) stands for the main luni-solar diurnal constituent; \(O1\) is the principal lunar diurnal constituent; \(M2\) represents the principal lunar constituent; and \(S2\) refers to the principal solar constituent [35].