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

Last change on this file since 701 was 630, checked in by Laurent Fairhead, 20 years ago

Import d'une version parallele de la dynamique YM
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.3 KB
RevLine 
[630]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
40     
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
49      DO 5 ij  = ijb,ije
50      w(ij,1)  = 0.
515     CONTINUE
52
53      RETURN
54      END
Note: See TracBrowser for help on using the repository browser.