source: LMDZ5/branches/AI-cosp/libf/dyn3dpar/vitvert_p.F @ 5308

Last change on this file since 5308 was 1907, checked in by lguez, 11 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
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

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.

  • 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      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.