! $Header$

SUBROUTINE interpost(q, qppm)
  USE lmdz_comgeom2

USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
  USE lmdz_paramet
  IMPLICIT NONE




  ! 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