Matlab code to produce the figures from Three-dimensional time-domain scattering of waves in the marginal ice zone

Three-dimensional scattering of ocean surface waves in the marginal ice zone (MIZ) is determined in the time-domain. The solution is found using spectral analysis of the linear operator for the Boltzmann equation. The method to calculate the scattering kernel that arises in the Boltzmann model from the single floe solution is also presented along with new identities for the far-field scattering, which can be used to validate the single floe solution. The spectrum of the operator is computed, and it is shown to have a universal structure under a special non-dimensionalization. This universal structure implies that under a scaling wave scattering in the MIZ has similar properties for a large range of ice types and wave periods. A scattering theory solution using fast Fourier transforms is given to find the solution for directional incident wave packets. A numerical solution method is also given using the split-step method and this is used to validate the spectral solution. Numerical calculations of the evolution of a typical wave field are presented.This article is part of the theme issue ‘Modelling of sea-ice phenomena’.