! $Header$

SUBROUTINE test_period(ucov, vcov, teta, q, p, phis)

  ! Auteur : P. Le Van
  !    ---------
  !  ....  Cette routine teste la periodicite en longitude des champs   ucov,
  !                       teta, q , p et phis                 ..........

  USE lmdz_infotrac, ONLY: nqtot

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


  !    ......  Arguments   ......

  REAL :: ucov(ip1jmp1, llm), vcov(ip1jm, llm), teta(ip1jmp1, llm), &
          q(ip1jmp1, llm, nqtot), p(ip1jmp1, llmp1), phis(ip1jmp1)

  !   .....  Variables  locales  .....

  INTEGER :: ij, l, nq

  DO l = 1, llm
    DO ij = 1, ip1jmp1, iip1
      IF(ucov(ij, l)/=ucov(ij + iim, l))  THEN
        PRINT *, 'STOP dans test_period car ---  UCOV  ---  n est pas', &
                ' periodique en longitude ! '
        PRINT *, ' l,  ij = ', l, ij, ij + iim
        STOP
      ENDIF
      IF(teta(ij, l)/=teta(ij + iim, l))  THEN
        PRINT *, 'STOP dans test_period car ---  TETA  ---  n est pas', &
                ' periodique en longitude ! '
        PRINT *, ' l,  ij = ', l, ij, ij + iim &
                , teta(ij, l), teta(ij + iim, l)
        STOP
      ENDIF
    ENDDO

    DO ij = 1, iim
      IF (teta(ij, l)/=teta(1, l) &
              .OR.teta(ip1jm + ij, l)/=teta(ip1jm + 1, l)) THEN
        PRINT *, 'STOP dans test_period car ---  TETA  ---  n est pas', &
                ' constant aux poles ! '
        PRINT*, 'teta(', 1, ',', l, ')=', teta(1, l)
        PRINT*, 'teta(', ij, ',', l, ')=', teta(ij, l)
        PRINT*, 'teta(', ip1jm + 1, ',', l, ')=', teta(ip1jm + 1, l)
        PRINT*, 'teta(', ip1jm + ij, ',', l, ')=', teta(ip1jm + ij, l)
        stop
      endif
    enddo
  ENDDO

  DO l = 1, llm
    DO ij = 1, ip1jm, iip1
      IF(vcov(ij, l)/=vcov(ij + iim, l))  THEN
        PRINT *, 'STOP dans test_period car ---  VCOV  ---  n est pas', &
                ' periodique en longitude !'
        PRINT *, ' l,  ij = ', l, ij, ij + iim, vcov(ij + iim, l), vcov(ij, l)
        vcov(ij + iim, l) = vcov(ij, l)
        ! STOP
      ENDIF
    ENDDO
  ENDDO

  DO nq = 1, nqtot
    DO l = 1, llm
      DO ij = 1, ip1jmp1, iip1
        IF(q(ij, l, nq)/=q(ij + iim, l, nq))  THEN
          PRINT *, 'STOP dans test_period car ---  Q  ---  n est pas ', &
                  'periodique en longitude !'
          PRINT *, ' nq , l,  ij = ', nq, l, ij, ij + iim
          STOP
        ENDIF
      ENDDO
    ENDDO
  ENDDO

  DO l = 1, llm
    DO ij = 1, ip1jmp1, iip1
      IF(p(ij, l)/=p(ij + iim, l))  THEN
        PRINT *, 'STOP dans test_period car ---  P  ---  n est pas', &
                ' periodique en longitude !'
        PRINT *, ' l ij = ', l, ij, ij + iim
        STOP
      ENDIF
      IF(phis(ij)/=phis(ij + iim))  THEN
        PRINT *, 'STOP dans test_period car ---  PHIS  ---  n est pas', &
                ' periodique en longitude !  l, IJ = ', l, ij, ij + iim
        PRINT *, ' ij = ', ij, ij + iim
        STOP
      ENDIF
    ENDDO
    DO ij = 1, iim
      IF (p(ij, l)/=p(1, l) &
              .OR.p(ip1jm + ij, l)/=p(ip1jm + 1, l)) THEN
        PRINT *, 'STOP dans test_period car ---  P     ---  n est pas', &
                ' constant aux poles ! '
        PRINT*, 'p(', 1, ',', l, ')=', p(1, l)
        PRINT*, 'p(', ij, ',', l, ')=', p(ij, l)
        PRINT*, 'p(', ip1jm + 1, ',', l, ')=', p(ip1jm + 1, l)
        PRINT*, 'p(', ip1jm + ij, ',', l, ')=', p(ip1jm + ij, l)
        stop
      endif
    enddo
  ENDDO

  RETURN
END SUBROUTINE test_period