Fotran code for calculation internal waves for horizontal disk