My Project
 All Classes Files Functions Variables Macros
vitvert_loc.F
Go to the documentation of this file.
1  SUBROUTINE vitvert_loc ( convm , w )
2 c
3  USE parallel
4  IMPLICIT NONE
5 
6 c=======================================================================
7 c
8 c Auteurs: P. Le Van , F. Hourdin .
9 c -------
10 c
11 c Objet:
12 c ------
13 c
14 c *******************************************************************
15 c .... calcul de la vitesse verticale aux niveaux sigma ....
16 c *******************************************************************
17 c convm est un argument d'entree pour le s-pg ......
18 c w est un argument de sortie pour le s-pg ......
19 c
20 c la vitesse verticale est orientee de haut en bas .
21 c au sol, au niveau sigma(1), w(i,j,1) = 0.
22 c au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi
23 c egale a 0. et n'est pas stockee dans le tableau w .
24 c
25 c
26 c=======================================================================
27 
28 #include "dimensions.h"
29 #include "paramet.h"
30 #include "comvert.h"
31 
32  REAL w(ijb_u:ije_u,llm),convm(ijb_u:ije_u,llm)
33  INTEGER l, ij,ijb,ije
34 
35 
36  ijb=ij_begin
37  ije=ij_end+iip1
38 
39  if (pole_sud) ije=ij_end
40 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
41  DO 2 l = 1,llmm1
42 
43  DO 1 ij = ijb,ije
44  w( ij, l+1 ) = convm( ij, l+1 ) - bp(l+1) * convm( ij, 1 )
45  1 CONTINUE
46 
47  2 CONTINUE
48 c$OMP END DO
49 c$OMP MASTER
50  DO 5 ij = ijb,ije
51  w(ij,1) = 0.
52 5 CONTINUE
53 c$OMP END MASTER
54 c$OMP BARRIER
55  RETURN
56  END