!-------------------------------------- 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 fuvacart  - compute the winds in the cartesian space from
!                  the tangential wind vector components.
!                  Fast version of uvacart special to Var analysis.
!

      subroutine fuvacart(pcart,pu,pv,psinlon,pcoslon,psinlat,pcoslat,knk) 1
!
      IMPLICIT NONE
#include "taglam4d.cdk"
      integer knk
      real*8 psinlon,pcoslon,psinlat,pcoslat
      real*8 pcart(3,knk),pu(knk),pv(knk)
!
!author Luc Fillion -ARMA/MSC 7 Feb 05. 
!            From ez_uvacart michel roch - april 90
!
!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.
!
!
      integer ji,jj,jk
      real*8  zfac1,zfac2
! 
      zfac1=pcoslon*psinlat
      zfac2=psinlon*psinlat
!
      do jk=1,knk
        pcart(1,jk) = -(pu(jk)*psinlon) - (pv(jk)*zfac1)
        pcart(2,jk) =  (pu(jk)*pcoslon) - (pv(jk)*zfac2)
        pcart(3,jk) =   pv(jk)*pcoslat
      enddo
!
      return
      end