Wave number from wave frequency 波频率计算波数
elemental real(rk) function k01(w) result(k0) real(rk), intent(in) :: w k0 = w*w/g end function k01