!-------------------------------------- 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 --------------------------------------
!

      subroutine bld_zd(pgdvort,pgddiv,pu,pv) 1,3
!
      IMPLICIT NONE
#if defined (DOC)
!
!s/r bld_zd : Shell controlling the passage from uv to zd.
!     Author L. Fillion  - ARMA/EC - 4 Feb 2010.
!Revision:
!
#endif
!
#include "taglam4d.cdk"
#include "pardim.cdk"
#include "comdim.cdk"
#include "comlun.cdk"
!
      real*8 pu(ni,nflev,nj)
      real*8 pv(ni,nflev,nj)
      real*8 pgdvort(ni,nflev,nj)
      real*8 pgddiv(ni,nflev,nj)
!
      integer ji,jj,jlev
      integer idum1,idum2,idum3,idum4,itrunc,index
      real*8  zmin,zmax
      real*8 zu2(-1:ni+2,nflev,-1:nj+2)
      real*8 zv2(-1:ni+2,nflev,-1:nj+2)
      real*8 zu9(0:ni+1,nflev,0:nj+1)
      real*8 zv9(0:ni+1,nflev,0:nj+1)
      real*8 zvort9(ni,nflev,nj)
      real*8 zdiv9(ni,nflev,nj)
!
!!
      write(nulout,*) 'Begin bld_zd'
!
!*1.  Build (VORT, DIV)
!     -----------------
!
      zu2(:,:,:) = 0.0
      zv2(:,:,:) = 0.0
!
      call symmetrize(zu2,pu,ni,nj,nflev)
      call symmetrize(zv2,pv,ni,nj,nflev)
!!
      zu9(:,:,:) = 0.0
      zv9(:,:,:) = 0.0
!
      do jlev=1,nflev
        do ji=0,ni+1
        do jj=0,nj+1
          zu9(ji,jlev,jj)=zu2(ji,jlev,jj)
          zv9(ji,jlev,jj)=zv2(ji,jlev,jj)
        enddo
        enddo
      enddo
!
      call UV2zdg_c(zvort9,zdiv9,zu9,zv9)
!
!      call maxmin(zvort9,ni,nj,nflev,zmin,zmax,
!     &            idum1,idum2,idum3,idum4,'initgdla    ',
!     &            'Z9')
!
!      call maxmin(zdiv9,ni,nj,nflev,zmin,zmax,
!     &            idum1,idum2,idum3,idum4,'initgdla    ',
!     &            'D9')
!
      do ji=1,ni
        do jj=1,nj
          do jlev=1,nflev
            pgdvort(ji,jlev,jj)=zvort9(ji,jlev,jj)
            pgddiv(ji,jlev,jj)=zdiv9(ji,jlev,jj)
          enddo
        enddo
      enddo
      write(nulout,*) 'End bld_zd'
!
      return
      end