source: LMDZ4/trunk/libf/dyn3dpar/vitvert_p.F @ 1354

Last change on this file since 1354 was 764, checked in by Laurent Fairhead, 17 years ago

Merge entre la version V3_conv et le HEAD
YM, JG, LF

  • 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
4      IMPLICIT NONE
5
6c=======================================================================
7c
8c   Auteurs:  P. Le Van , F. Hourdin .
9c   -------
10c
11c   Objet:
12c   ------
13c
14c    *******************************************************************
15c  .... calcul de la vitesse verticale aux niveaux sigma  ....
16c    *******************************************************************
17c     convm   est un argument  d'entree pour le s-pg  ......
18c       w     est un argument de sortie pour le s-pg  ......
19c
20c    la vitesse verticale est orientee de  haut en bas .
21c    au sol, au niveau sigma(1),   w(i,j,1) = 0.
22c    au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi
23c    egale a 0. et n'est pas stockee dans le tableau w  .
24c
25c
26c=======================================================================
27
28#include "dimensions.h"
29#include "paramet.h"
30#include "comvert.h"
31
32      REAL w(ip1jmp1,llm),convm(ip1jmp1,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
40c$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
48c$OMP END DO
49c$OMP MASTER
50      DO 5 ij  = ijb,ije
51      w(ij,1)  = 0.
525     CONTINUE
53c$OMP END MASTER
54c$OMP BARRIER
55      RETURN
56      END
Note: See TracBrowser for help on using the repository browser.