source: trunk/LMDZ.COMMON/libf/dyn3d_common/vitvert.F @ 1459

Last change on this file since 1459 was 1422, checked in by milmd, 10 years ago

In GENERIC, MARS and COMMON models replace some include files by modules (usefull for decoupling physics with dynamics).
File size: 1.2 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE vitvert ( convm , w )
5c
6      USE comvert_mod, ONLY: bp
7
8      IMPLICIT NONE
9
10c=======================================================================
11c
12c   Auteurs:  P. Le Van , F. Hourdin .
13c   -------
14c
15c   Objet:
16c   ------
17c
18c    *******************************************************************
19c  .... calcul de la vitesse verticale aux niveaux sigma  ....
20c    *******************************************************************
21c     convm   est un argument  d'entree pour le s-pg  ......
22c       w     est un argument de sortie pour le s-pg  ......
23c
24c    la vitesse verticale est orientee de  haut en bas .
25c    au sol, au niveau sigma(1),   w(i,j,1) = 0.
26c    au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi
27c    egale a 0. et n'est pas stockee dans le tableau w  .
28c
29c
30c=======================================================================
31
32#include "dimensions.h"
33#include "paramet.h"
34
35      REAL w(ip1jmp1,llm),convm(ip1jmp1,llm)
36      INTEGER   l, ij
37
38
39
40      DO 2  l = 1,llmm1
41
42      DO 1 ij = 1,ip1jmp1
43      w( ij, l+1 ) = convm( ij, l+1 ) - bp(l+1) * convm( ij, 1 )
44   1  CONTINUE
45
46   2  CONTINUE
47
48      DO 5 ij  = 1,ip1jmp1
49      w(ij,1)  = 0.
505     CONTINUE
51
52      RETURN
53      END
Note: See TracBrowser for help on using the repository browser.