! 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.-

　　　　 ! 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.