source: LMDZ6/branches/Amaury_dev/libf/dyn3d/lmdz_dudv1.f90 @ 5218

MODULE lmdz_dudv1
  IMPLICIT NONE; PRIVATE
  PUBLIC dudv1

CONTAINS


  SUBROUTINE dudv1(vorpot, pbaru, pbarv, du, dv)
    USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
    USE lmdz_paramet
    IMPLICIT NONE

    !-----------------------------------------------------------------------

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

    !   Objet:
    !   ------
    !   calcul du terme de  rotation
    !   ce terme est ajoute a  d(ucov)/dt et a d(vcov)/dt  ..
    !   vorpot, pbaru et pbarv sont des arguments d'entree  pour le s-pg ..
    !   du  et dv              sont des arguments de sortie pour le s-pg ..

    !-----------------------------------------------------------------------

    REAL :: vorpot(ip1jm, llm), pbaru(ip1jmp1, llm), &
            pbarv(ip1jm, llm), du(ip1jmp1, llm), dv(ip1jm, llm)
    INTEGER :: l, ij

    DO l = 1, llm

      DO ij = iip2, ip1jm - 1
        du(ij, l) = 0.125 * (vorpot(ij - iip1, l) + vorpot(ij, l)) * &
                (pbarv(ij - iip1, l) + pbarv(ij - iim, l) + &
                        pbarv(ij, l) + pbarv(ij + 1, l))
      END DO

      DO ij = 1, ip1jm - 1
        dv(ij + 1, l) = - 0.125 * (vorpot(ij, l) + vorpot(ij + 1, l)) * &
                (pbaru(ij, l) + pbaru(ij + 1, l) + &
                        pbaru(ij + iip1, l) + pbaru(ij + iip2, l))
      END DO

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

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

    END DO

  END SUBROUTINE dudv1


END MODULE lmdz_dudv1