! $Id: pres2lev.F 1179 2009-06-11 14:18:47Z jghattas $

MODULE lmdz_pres2lev
  IMPLICIT NONE; PRIVATE
  PUBLIC pres2lev

CONTAINS

  !******************************************************
  SUBROUTINE pres2lev(varo, varn, lmo, lmn, po, pn, ni, nj, ok_invertp)

    ! interpolation lineaire pour passer
    ! a une nouvelle discretisation verticale pour
    ! les variables de GCM
    ! Francois Forget (01/1995)
    ! MOdif remy roca 12/97 pour passer de pres2sig
    ! Modif F.Codron 07/08 po en 3D
    !**********************************************************

    IMPLICIT NONE

    !   Declarations:
    ! ==============

    !  ARGUMENTS
    !  """""""""
    LOGICAL, INTENT(IN) :: ok_invertp
    INTEGER, INTENT(IN) :: lmo ! dimensions ancienne couches
    INTEGER, INTENT(IN) :: lmn ! dimensions nouvelle couches

    INTEGER, INTENT(IN) :: ni, nj ! nombre de point horizontal
    REAL, INTENT(IN) :: po(ni * nj, lmo) ! niveau de pression ancienne grille
    REAL, INTENT(IN) :: pn(ni * nj, lmn) ! niveau de pression nouvelle grille

    REAL, INTENT(IN) :: varo(ni * nj, lmo) ! var dans l'ancienne grille
    REAL, INTENT(OUT) :: varn(ni * nj, lmn) ! var dans la nouvelle grille

    REAL :: zvaro(ni * nj, lmo), zpo(ni * nj, lmo)

    ! Autres variables
    ! """"""""""""""""
    INTEGER :: ln, lo, k
    REAL :: coef


    ! Inversion de l'ordre des niveaux verticaux
    IF (ok_invertp) THEN
      DO lo = 1, lmo
        DO k = 1, ni * nj
          zpo(k, lo) = po(k, lmo + 1 - lo)
          zvaro(k, lo) = varo(k, lmo + 1 - lo)
        ENDDO
      ENDDO
    ELSE
      DO lo = 1, lmo
        DO k = 1, ni * nj
          zpo(k, lo) = po(k, lo)
          zvaro(k, lo) = varo(k, lo)
        ENDDO
      ENDDO
    ENDIF

    DO ln = 1, lmn
      DO lo = 1, lmo - 1
        DO k = 1, ni * nj
          IF (pn(k, ln) >= zpo(k, 1)) THEN
            varn(k, ln) = zvaro(k, 1)
          ELSE IF (pn(k, ln) <= zpo(k, lmo)) THEN
            varn(k, ln) = zvaro(k, lmo)
          ELSE IF (pn(k, ln) <= zpo(k, lo) .AND. pn(k, ln) > zpo(k, lo + 1)) THEN
            coef = (pn(k, ln) - zpo(k, lo)) / (zpo(k, lo + 1) - zpo(k, lo))
            varn(k, ln) = zvaro(k, lo) + coef * (zvaro(k, lo + 1) - zvaro(k, lo))
          ENDIF

        ENDDO
      ENDDO
    ENDDO

  END SUBROUTINE pres2lev

END MODULE lmdz_pres2lev