!-------------------------------------- LICENCE BEGIN ------------------------------------
!Environment Canada - Atmospheric Science and Technology License/Disclaimer,
!                     version 3; Last Modified: May 7, 2008.
!This is free but copyrighted software; you can use/redistribute/modify it under the terms
!of the Environment Canada - Atmospheric Science and Technology License/Disclaimer
!version 3 or (at your option) any later version that should be found at:
!http://collaboration.cmc.ec.gc.ca/science/rpn.comm/license.html
!
!This software 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 above mentioned License/Disclaimer for more details.
!You should have received a copy of the License/Disclaimer along with this software;
!if not, you can write to: EC-RPN COMM Group, 2121 TransCanada, suite 500, Dorval (Quebec),
!CANADA, H9P 1J3; or send e-mail to service.rpn@ec.gc.ca
!-------------------------------------- LICENCE END --------------------------------------
!
!
!**S/R afuvacart  - Adjoint of fuvacart.
!

      subroutine afuvacart(pu,pv,pcart,psinlonr,pcoslonr,psinlatr,pcoslatr,knk) 1
!
      IMPLICIT NONE
      integer knk
      real*8 psinlonr,pcoslonr,psinlatr,pcoslatr
      real*8 pcart(3,knk),pu(knk),pv(knk)
!
!author Luc Fillion -ARMA/MSC 7 Feb 05. 
!            From ez_uvacart michel roch - april 90
!
!revision: Luc Fillion -ARMA/MSC 25 Feb 05 - Fast version
!
!arguments
!    out    pvcart - Cartesian components of the tangential wind vector
!    IN     pu     - Zonal component of the wind on the tangential plane.
!           pv     - Meridional component of the wind on the tangential plane.
!           plon   - Longitudes of grid points in radians.
!           plat   - Longitudes of grid points in radians.
!           ni     - E-W DIMENSION of the grid
!           nj     - N-S dimension of the grid
!
#include "taglam4d.cdk"
      integer jk
      real*8  zfac1,zfac2
! 
      do jk=1,knk
        pu(jk)=0.0  ! dont use sub. zero.ftn here... *4 versus *8 diff...
        pv(jk)=0.0
      enddo
!
      zfac1=psinlonr*psinlatr
      zfac2=pcoslonr*psinlatr
!
      do jk=1,knk
        pu(jk) = pu(jk)+pcart(2,jk)*pcoslonr
        pv(jk) = pv(jk)-pcart(2,jk)*zfac1
        pcart(2,jk) = 0.0
        pu(jk) = pu(jk)-pcart(1,jk)*psinlonr
        pv(jk) = pv(jk)-pcart(1,jk)*zfac2
        pcart(1,jk) = 0.0
        pv(jk) = pv(jk)+pcart(3,jk)*pcoslatr
        pcart(3,jk) = 0.0
      enddo
!
      return
      end