source: LMDZ4/branches/IPSL-CM4_IPCC_patches/libf/dyn3d/vitvert.F @ 5306

!
! $Header$
!
      SUBROUTINE vitvert ( convm , w )
c
      IMPLICIT NONE

c=======================================================================
c
c   Auteurs:  P. Le Van , F. Hourdin .
c   -------
c
c   Objet:
c   ------
c
c    *******************************************************************
c  .... calcul de la vitesse verticale aux niveaux sigma  ....
c    *******************************************************************
c     convm   est un argument  d'entree pour le s-pg  ......
c       w     est un argument de sortie pour le s-pg  ......
c
c    la vitesse verticale est orientee de  haut en bas .
c    au sol, au niveau sigma(1),   w(i,j,1) = 0.
c    au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi
c    egale a 0. et n'est pas stockee dans le tableau w  .
c
c
c=======================================================================

#include "dimensions.h"
#include "paramet.h"
#include "comvert.h"

      REAL w(ip1jmp1,llm),convm(ip1jmp1,llm)
      INTEGER   l, ij



      DO 2  l = 1,llmm1

      DO 1 ij = 1,ip1jmp1
      w( ij, l+1 ) = convm( ij, l+1 ) - bp(l+1) * convm( ij, 1 )
   1  CONTINUE

   2  CONTINUE

      DO 5 ij  = 1,ip1jmp1
      w(ij,1)  = 0.
5     CONTINUE

      RETURN
      END