source: trunk/LMDZ.PLUTO.old/libf/dyn3d/vitvert.F @ 3436

Last change on this file since 3436 was 3175, checked in by emillour, 11 months ago

Pluto PCM:
Add the old Pluto LMDZ for reference (required prior step to making
an LMDZ.PLUTO using the same framework as the other physics packages).
TB+EM
File size: 1.2 KB
Line 
1      SUBROUTINE vitvert ( convm , w )
2c
3      IMPLICIT NONE
4
5c=======================================================================
6c
7c   Auteurs:  P. Le Van , F. Hourdin .
8c   -------
9c
10c   Objet:
11c   ------
12c
13c    *******************************************************************
14c  .... calcul de la vitesse verticale aux niveaux sigma  ....
15c    *******************************************************************
16c     convm   est un argument  d'entree pour le s-pg  ......
17c       w     est un argument de sortie pour le s-pg  ......
18c
19c    la vitesse verticale est orientee de  haut en bas .
20c    au sol, au niveau sigma(1),   w(i,j,1) = 0.
21c    au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi
22c    egale a 0. et n'est pas stockee dans le tableau w  .
23c
24c
25c=======================================================================
26
27#include "dimensions.h"
28#include "paramet.h"
29#include "comvert.h"
30
31      REAL w(ip1jmp1,llm),convm(ip1jmp1,llm)
32      INTEGER   l, ij
33
34
35
36      DO 2  l = 1,llmm1
37
38      DO 1 ij = 1,ip1jmp1
39      w( ij, l+1 ) = convm( ij, l+1 ) - bp(l+1) * convm( ij, 1 )
40   1  CONTINUE
41
42   2  CONTINUE
43
44      DO 5 ij  = 1,ip1jmp1
45      w(ij,1)  = 0.
465     CONTINUE
47
48      RETURN
49      END
Note: See TracBrowser for help on using the repository browser.