Introduction - If you have any usage issues, please Google them yourself
! This FUNCTION solves the finite depth dispersion relationship:
!
! k*tanh(k*h)=(Omega^2)/g
!
! for k, the wavenumber (WaveNumber) given the frequency, Omega,
! gravitational constant, g, and water depth, h, as inputs. A
! high order initial guess is used in conjunction with a quadratic
! Newton s method for the solution with seven significant digits
! accuracy using only one iteration pass. The method is due to
! Professor J.N. Newman of M.I.T. as found in routine EIGVAL of
! the SWIM-MOTION-LINES (SML) software package in source file
! Solve.f of the SWIM module.