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

MODULE lmdz_dudv2
  IMPLICIT NONE; PRIVATE
  PUBLIC dudv2

CONTAINS

  SUBROUTINE dudv2(teta, pkf, bern, du, dv)

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

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

    !   Auteur:  P. Le Van
    !   -------

    !   Objet:
    !   ------

    !   *****************************************************************
    !   ..... calcul du terme de pression (gradient de p/densite )   et
    !      du terme de ( -gradient de la fonction de Bernouilli ) ...
    !   *****************************************************************
    !      Ces termes sont ajoutes a  d(ucov)/dt et a d(vcov)/dt  ..


    !    teta , pkf, bern  sont des arguments d'entree  pour le s-pg  ....
    !    du et dv          sont des arguments de sortie pour le s-pg  ....

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

    REAL :: teta(ip1jmp1, llm), pkf(ip1jmp1, llm), bern(ip1jmp1, llm), &
            du(ip1jmp1, llm), dv(ip1jm, llm)
    INTEGER :: l, ij

    DO l = 1, llm

      DO ij = iip2, ip1jm - 1
        du(ij, l) = du(ij, l) + 0.5 * (teta(ij, l) + teta(ij + 1, l)) * &
                (pkf(ij, l) - pkf(ij + 1, l)) + bern(ij, l) - bern(ij + 1, l)
      END DO


      !    .....  correction  pour du(iip1,j,l),  j=2,jjm   ......
      !    ...          du(iip1,j,l) = du(1,j,l)                 ...

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

      DO ij = 1, ip1jm
        dv(ij, l) = dv(ij, l) + 0.5 * (teta(ij, l) + teta(ij + iip1, l)) * &
                (pkf(ij + iip1, l) - pkf(ij, l)) &
                + bern(ij + iip1, l) - bern(ij, l)
      END DO

    END DO
    !

  END SUBROUTINE dudv2


END MODULE lmdz_dudv2