!----- AGPL --------------------------------------------------------------------
!                                                                               
!  Copyright (C)  Stichting Deltares, 2017-2022.                                
!                                                                               
!  This file is part of Delft3D (D-Flow Flexible Mesh component).               
!                                                                               
!  Delft3D is free software: you can redistribute it and/or modify              
!  it under the terms of the GNU Affero General Public License as               
!  published by the Free Software Foundation version 3.                         
!                                                                               
!  Delft3D  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 Affero General Public License for more details.                          
!                                                                               
!  You should have received a copy of the GNU Affero General Public License     
!  along with Delft3D.  If not, see <http://www.gnu.org/licenses/>.             
!                                                                               
!  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",                  
!  "D-Flow Flexible Mesh" and "Deltares" are registered trademarks of Stichting 
!  Deltares, and remain the property of Stichting Deltares. All rights reserved.
!                                                                               
!-------------------------------------------------------------------------------

! $Id$
! $HeadURL$

 subroutine getucxucynoweirswrong(ku, ucxku, ucyku, ischeme)  !wrong
 use m_flow
 use m_flowgeom
 implicit none

 integer           :: ku, LL, L, Ls, ischeme
 double precision  :: ucxku, ucyku, wwx, wwy, ww, wwt, ac1, wwxt, wwyt, ux, uy

 ucxku = 0d0 ; ucyku = 0d0; wwt = 0d0; wwxt = 0d0 ; wwyt = 0d0
 do LL = 1,nd(ku)%lnx
    Ls = nd(ku)%ln(LL); L = iabs(Ls)
    if ( iadv(L) < 21 .or. iadv(L) > 25 ) then  ! no weir
       if (Ls < 0) then
          ac1 = acL(L)
       else
          ac1 = 1d0 - acL(L)
       endif
       ww    = ac1*dx(L)*wu(L)
       if (ischeme == 3) then
          wwx   = csu(L)*ww
          wwy   = snu(L)*ww
          ucxku = ucxku + wwx*u0(L)
          ucyku = ucyku + wwy*u0(L)
       else
          wwx   = abs(csu(L))*ww
          wwy   = abs(snu(L))*ww
          ux    = csu(L)*u0(L)
          uy    = snu(L)*u0(L)
          ucxku = ucxku + ux*wwx
          ucyku = ucyku + uy*wwy
       endif
       wwxt = wwxt + abs(wwx)
       wwyt = wwyt + abs(wwy)
    endif
 enddo

 if (wwxt > 0d0) ucxku = ucxku / wwxt
 if (wwyt > 0d0) ucyku = ucyku / wwyt

 end subroutine getucxucynoweirswrong