! $Header$

SUBROUTINE interpost(q, qppm)
  USE lmdz_comgeom2

  IMPLICIT NONE

  INCLUDE "dimensions.h"
  INCLUDE "paramet.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