subroutine cur2ca(u ,v ,alfa ,k ,n , & & m ) !----- GPL --------------------------------------------------------------------- ! ! Copyright (C) Stichting Deltares, 2011-2021. ! ! This program is free software: you can redistribute it and/or modify ! it under the terms of the GNU General Public License as published by ! the Free Software Foundation version 3. ! ! This program is distributed in the hope that it will be useful, ! but WITHOUT ANY WARRANTY; without even the implied warranty of ! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the ! GNU General Public License for more details. ! ! You should have received a copy of the GNU General Public License ! along with this program. If not, see . ! ! contact: delft3d.support@deltares.nl ! Stichting Deltares ! P.O. Box 177 ! 2600 MH Delft, The Netherlands ! ! All indications and logos of, and references to, "Delft3D" and "Deltares" ! are registered trademarks of Stichting Deltares, and remain the property of ! Stichting Deltares. All rights reserved. ! !------------------------------------------------------------------------------- ! $Id$ ! $HeadURL$ !!--description----------------------------------------------------------------- ! NONE !!--pseudo code and references-------------------------------------------------- ! NONE !!--declarations---------------------------------------------------------------- use precision implicit none ! ! ! COMMON variables ! real(hp) :: facrad common /fr / facrad ! ! Global variables ! integer, intent(in) :: m integer, intent(in) :: n integer, dimension(n, m), intent(in) :: k real(hp), dimension(n, m), intent(in) :: alfa real(hp), dimension(n, m) :: u real(hp), dimension(n, m) :: v ! ! ! Local variables ! integer :: i integer :: j real(hp) :: ang real(hp) :: ux real(hp) :: uy ! ! !! executable statements ------------------------------------------------------- ! ! ! ! write (*, *) ' CUR2CA' do i = 1, m do j = 1, n if (k(j, i)==1) then ang = facrad*alfa(j, i) ux = u(j, i)*cos(ang) - v(j, i)*sin(ang) uy = u(j, i)*sin(ang) + v(j, i)*cos(ang) u(j, i) = ux v(j, i) = uy else u(j, i) = 0._hp v(j, i) = 0._hp endif enddo enddo end subroutine cur2ca