source: LMDZ6/branches/Amaury_dev/libf/dyn3d/lmdz_dteta1.f90 @ 5220

MODULE lmdz_dteta1
  IMPLICIT NONE; PRIVATE
  PUBLIC dteta1

CONTAINS


  SUBROUTINE dteta1(teta, pbaru, pbarv, dteta)
    USE lmdz_filtreg, ONLY: filtreg
    USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
    USE lmdz_paramet
    IMPLICIT NONE

    !=======================================================================

    !   Auteur:  P. Le Van
    !   -------
    ! Modif F.Forget 03/94 (on retire q et dq  pour construire dteta1)

    !   ********************************************************************
    !   ... calcul du terme de convergence horizontale du flux d'enthalpie
    !    potentielle   ......
    !   ********************************************************************
    !  .. teta,pbaru et pbarv sont des arguments d'entree  pour le s-pg ....
    ! dteta               sont des arguments de sortie pour le s-pg ....

    !=======================================================================

    REAL :: teta(ip1jmp1, llm), pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
    REAL :: dteta(ip1jmp1, llm)
    INTEGER :: l, ij

    REAL :: hbyv(ip1jm, llm), hbxu(ip1jmp1, llm)

    !

    DO l = 1, llm

      DO ij = iip2, ip1jm - 1
        hbxu(ij, l) = pbaru(ij, l) * 0.5 * (teta(ij, l) + teta(ij + 1, l))
      END DO

      !    .... correction pour  hbxu(iip1,j,l)  .....
      !    ....   hbxu(iip1,j,l)= hbxu(1,j,l) ....

      !DIR$ IVDEP
      DO ij = iip1 + iip1, ip1jm, iip1
        hbxu(ij, l) = hbxu(ij - iim, l)
      END DO

      DO ij = 1, ip1jm
        hbyv(ij, l) = pbarv(ij, l) * 0.5 * (teta(ij, l) + teta(ij + iip1, l))
      END DO

    END DO

    CALL  convflu (hbxu, hbyv, llm, dteta)


    !    stockage dans  dh de la convergence horizont. filtree' du  flux
    ! ....                           ...........
    ! d'enthalpie potentielle .

    CALL filtreg(dteta, jjp1, llm, 2, 2, .TRUE., 1)

    !

  END SUBROUTINE dteta1


END MODULE lmdz_dteta1