unexpected boundary condition values when trying to create monochromatic or monodirectional waves

When I tried to create monodirectional waves by setting BOUNDARY DIRECTIONAL SPREAD 1 to a very high value, I noticed the significant wave height calculated in the model was different form the value that was imposed using BOUNDARY SIGNIFICANT WAVE HEIGHT. Apparently the analytic expression used to calculate the directional distribution (subroutines FSPRD1.f, FSPRD2.f,FSPRD3.f) can give values

We might consider correcting by numerically integrating the directional distribution (array FRA in speini.f) and using that normalize.

Something similar can also occur for the spectrum when using a very high spectral peak factor to make monochromatic waves.

Not sure it is worth the effort though.