! $Header$

  SUBROUTINE interpost(q,qppm)

   implicit none


  include "dimensions.h"
  include "paramet.h"
  include "comgeom2.h"

  ! Arguments
  real :: q(iip1,jjp1,llm)
  real :: qppm(iim,jjp1,llm)
  ! Local
  integer :: l,i,j

  ! RE-INVERSION DES NIVEAUX
  ! le programme ppm3d travaille avec une 3ème coordonnée inversée par rapport
  ! de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
  ! On passe donc des niveaux de Lin à ceux du LMDZ

    do l=1,llm
      do j=1,jjp1
         do i=1,iim
             q(i,j,l)=qppm(i,j,llm-l+1)
         enddo
      enddo
     enddo

  ! BOUCLAGE EN LONGITUDE PAS EFFECTUE DANS PPM3D

     do l=1,llm
       do j=1,jjp1
        q(iip1,j,l)=q(1,j,l)
       enddo
     enddo


   return

end subroutine interpost