A new ordinary differential equation for the evaluation of the frequency-domain Green function

Clément (2013) derived a second order ordinary differential equation (ODE) satisfied by the free-surface Green function in the frequency domain. Since then, similar ODEs for the gradient of the Green function have been developed. Unfortunately, all these ODEs degenerate at zero frequency. Therefore, it is not possible to initialize the numerical solution of these ODEs from this zero frequency. Alternative methods based on the shifting of the initial condition to frequencies strictly greater than zero have then been developed.