A Numerical Implementation of a Nonlinear Mild Slope Model for Shoaling Directional Waves

We describe the numerical implementation of a phase-resolving, nonlinear spectral model for shoaling directional waves over a mild sloping beach with straight parallel isobaths. The model accounts for non-linear, quadratic (triad) wave interactions as well as shoaling and refraction. The model integrates the coupled, nonlinear hyperbolic evolution equations that describe the transformation of the complex Fourier amplitudes of the deep-water directional wave field. Because typical directional wave spectra (observed or produced by deep-water forecasting models such as WAVEWATCH IIITM) do not contain phase information, individual realizations are generated by associating a random phase to each Fourier mode. The approach provides a natural extension to the deep-water spectral wave models, and has the advantage of fully describing the shoaling wave stochastic process, i.e., the evolution of both the variance and higher order statistics (phase correlations), the latter related to the evolution of the wave shape. The numerical implementation (a Fortran 95/2003 code) includes unidirectional (shore-perpendicular) propagation as a special case. Interoperability, both with post-processing programs (e.g., MATLAB/Tecplot 360) and future model coupling (e.g., offshore wave conditions from WAVEWATCH IIITM), is promoted by using NetCDF-4/HD5 formatted output files. The capabilities of the model are demonstrated using a JONSWAP spectrum with a cos(2s) directional distribution, for shore-perpendicular and oblique propagation. The simulated wave transformation under combined shoaling, refraction and nonlinear interactions shows the expected generation of directional harmonics of the spectral peak and of infragravity (frequency < 0.05 Hz) waves. Current development efforts focus on analytic testing, development of additional physics modules essential for applications and validation with laboratory and field observations.