vitvert_p.F

Go to the documentation of this file.

      SUBROUTINE vitvert_p ( convm , w )
c
      USE parallel_lmdz
      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,ijb,ije


      ijb=ij_begin
      ije=ij_end+iip1
      
      if (pole_sud) ije=ij_end
c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
      DO 2  l = 1,llmm1

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

   2  CONTINUE
c$OMP END DO
c$OMP MASTER
      DO 5 ij  = ijb,ije
      w(ij,1)  = 0.
5     CONTINUE
c$OMP END MASTER
c$OMP BARRIER
      RETURN
      END