## 4 Met-Ocean Strategy for Assessing the Feasibility of Marine Renewable Energy near Ports
### Post-Processing  
The system was included within a web-based graphical user interface (GUI), capable of providing the following visualisations: i. hourly agitation prediction maps for the next 72 h for Hs; ii. hourly time series of Hs, Tp, and _Pw_ at specific points of interest (pre-selected and agreed upon with APV for each study port); iii. Ocean-meteorological prediction reports featuring agitation maps across the port domain and time series of Hs, Tp, and _Pw_ at points of interest; iv. directional wave disaggregation; and v. the generation of hourly energy potential maps covering a minimum of one year (see Figure 7).  
For the calculation of directional energy potential (_Pw_), the method employs a sea state represented hourly by its directional energy spectrum \(S(\omega,\theta)\), which represents the  
\begin{table}
\begin{tabular}{c c c c c c c} \hline \hline
**Port** & **\# of Hours** & **BIAS** & **RMSE** & **SI** & **CORR** & \(\mathbf{R^{2}}\) \\ \hline
**Valencia** & 26,962 & 0.035 & 0.122 & 0.555 & 0.81 & 0.637 \\ \hline
**Sagunto** & 26,879 & \(-\)0.039 & 0.126 & 0.465 & 0.82 & 0.667 \\ \hline
**Gandia** & 26,745 & \(-\)0.010 & 0.065 & 0.400 & 0.88 & 0.778 \\ \hline \hline \end{tabular}
\end{table}
Table 2: Post-validation statistical values for Hs validation for Valencia, Sagunto y Gandía.  
Figure 7: Valencia Port GUI example for (**a**) regional wave downscaling, (**b**) local wave propagation, (**c**) high-resolution agitation, and (**d**) wave forcing, wave agitation, and (total, incident, and reflected) wave potential energy forecast.  
temporal energy flux given by the following expression, integrating the energy fluxes from all components of the spectrum:  
\[Pw=\rho g\int_{0}^{2\pi}\int_{0}^{\infty}c_{g}(\omega,z)\ S(\omega,\theta)d \omega d\theta \tag{1}\]  
\[c_{g}=\frac{c}{2}\bigg{(}1+\frac{2 kz}{\sinh\ (2 kz)}\bigg{)} \tag{2}\]  
where  
\(\rho\) is the density of seawater, \(\rho=1025\) kg/m\({}^{3}\);  
\(g\) is the acceleration due to gravity, \(g=9.81\) m/s\({}^{2}\);  
\(S(\omega,\theta)\) is the directional energy spectrum representing the energy density (energy per unit area) assigned to each frequency \(\omega\) and direction \(\theta\) of the sea state.  
\(Cg(\omega,z)\) is the group celerity representing the propagation velocity of wave energy, calculated by the following:  
\(c\) is the wave celerity, \(c=L/T=\omega/k\);  
\(\omega\) is the angular frequency of the wave, \(\omega=2\pi/T\);  
\(k\) is the wave number, \(k=2\pi/L\);  
\(T\) is the wave period;  
\(L\) is the wavelength;  
\(z\) is the water depth (positive distance from the surface to the seabed).  
Angular frequency, wave number, and depth are related by the dispersion relation:  
\[\omega^{2}=gk\tanh\ (kz) \tag{3}\]  
In cases where the directional spectrum of the sea state is unavailable, wave power can be determined using spectral parameters representative of the sea state (obtained directly through the spectral reconstruction method developed by IHCantabria and implemented in Puertos del Estado's operations), as follows:  
\[Pw=\frac{1}{8}\rho g\ Hm0^{2}c_{g}(\omega,z) \tag{4}\]  
where the moment of order \(n\) of the energy spectrum is as follows:  
\[m\ n=\int_{0}^{2\pi}\int_{0}^{\infty}\omega^{n}S(\omega,\theta)\ d\omega d \theta;\ \ n=0,\ 1,\ 2,\dots \tag{5}\]