source: LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv2.F90 @ 5118

! $Header$

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

  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  ....
  !
  !=======================================================================
  !
  include "dimensions.h"
  include "paramet.h"

  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