copyright (C) 2001  MSC-RPN COMM  %%%MC2%%%
***s/r slabsym
*

      subroutine slabsym 1
      implicit none
*
*OBJECT
*
*METHOD
*
*EXTERNALS
*
*AUTHOR   Claude Girard
*
*HISTORY
*
**
#include "dynmem.cdk"
*
      integer i,j,k,n
*----------------------------------------------------------------------
*
* slab-symetry

      if(north.eq.1) then
         do k=1,gnk
            do j=1,hy+1
               do i=1-hx,ldni+hx
                  ppp(i,ldnj-1+j,k)     =    ppp(i,ldnj-1,k)
                  uup(i,ldnj-1+j,k)     =    uup(i,ldnj-1,k)
                  vvp(i,ldnj-1+j,k)     =    vvp(i,ldnj-1,k)
                  wwp(i,ldnj-1+j,k)     =    wwp(i,ldnj-1,k)
                  bbp(i,ldnj-1+j,k)     =    bbp(i,ldnj-1,k)
                  hup(i,ldnj-1+j,k)     =    hup(i,ldnj-1,k)
                  do n = 1, ntr
                    trp(i,ldnj-1+j,k,n) =    trp(i,ldnj-1,k,n)
                  enddo
               enddo
            enddo
         enddo
         do j=1,hy+1
            do i=1-hx,ldni+hx
               ppp(i,ldnj-1+j,0)    =    ppp(i,ldnj-1,0)
            enddo
         enddo
      endif
*
      if(south.eq.1) then
         do k=1,gnk
            do i=1-hx,ldni+hx
              vvp(i,1,k)   =  vvp(i,2,k)
            enddo
            do j=1,hy
               do i=1-hx,ldni+hx
                  ppp(i,1-j,k)   =   ppp(i,1,k)
                  uup(i,1-j,k)   =   uup(i,1,k)
                  vvp(i,1-j,k)   =   vvp(i,1,k)
                  wwp(i,1-j,k)   =   wwp(i,1,k)
                  bbp(i,1-j,k)   =   bbp(i,1,k)
                  hup(i,1-j,k)   =   hup(i,1,k)
                  do n = 1, ntr
                    trp(i,1-j,k,n) = trp(i,1,k,n)
                  enddo
               enddo
            enddo
         enddo
         do j=1,hy
            do i=1-hx,ldni+hx
               ppp(i,1-j,0)   =   ppp(i,1,0)
            enddo
         enddo
      endif
*
*----------------------------------------------------------------------
      return
      end