source: LMDZ5/trunk/libf/dyn3dpar/vitvert_p.F @ 4735

Last change on this file since 4735 was 2600, checked in by Ehouarn Millour, 8 years ago

Cleanup in the dynamics: turn comvert.h into module comvert_mod.F90
EM

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.4 KB
Line 
1      SUBROUTINE vitvert_p ( convm , w )
2c
3      USE parallel_lmdz
4      USE comvert_mod, ONLY: bp
5     
6      IMPLICIT NONE
7
8c=======================================================================
9c
10c   Auteurs:  P. Le Van , F. Hourdin .
11c   -------
12c
13c   Objet:
14c   ------
15c
16c    *******************************************************************
17c  .... calcul de la vitesse verticale aux niveaux sigma  ....
18c    *******************************************************************
19c     convm   est un argument  d'entree pour le s-pg  ......
20c       w     est un argument de sortie pour le s-pg  ......
21c
22c    la vitesse verticale est orientee de  haut en bas .
23c    au sol, au niveau sigma(1),   w(i,j,1) = 0.
24c    au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi
25c    egale a 0. et n'est pas stockee dans le tableau w  .
26c
27c
28c=======================================================================
29
30#include "dimensions.h"
31#include "paramet.h"
32
33      REAL w(ip1jmp1,llm),convm(ip1jmp1,llm)
34      INTEGER   l, ij,ijb,ije
35
36
37      ijb=ij_begin
38      ije=ij_end+iip1
39     
40      if (pole_sud) ije=ij_end
41c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)     
42      DO 2  l = 1,llmm1
43
44      DO 1 ij = ijb,ije
45      w( ij, l+1 ) = convm( ij, l+1 ) - bp(l+1) * convm( ij, 1 )
46   1  CONTINUE
47
48   2  CONTINUE
49c$OMP END DO
50c$OMP MASTER
51      DO 5 ij  = ijb,ije
52      w(ij,1)  = 0.
535     CONTINUE
54c$OMP END MASTER
55c$OMP BARRIER
56      RETURN
57      END
Note: See TracBrowser for help on using the repository browser.