source: LMDZ4/trunk/libf/dyn3dpar/vitvert.F @ 1000

Last change on this file since 1000 was 1000, checked in by Laurent Fairhead, 16 years ago
  • Modifs sur le parallelisme: masquage dans la physique
  • Inclusion strato
  • mise en coherence etat0
  • le mode offline fonctionne maintenant en parallele,
  • les fichiers de la dynamiques sont correctement sortis et peuvent etre reconstruit avec rebuild
  • la version parallele de la dynamique peut s'executer sans MPI (sur 1 proc)
  • L'OPENMP fonctionne maintenant sans la parallelisation MPI.

YM
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.2 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE vitvert ( convm , w )
5c
6      IMPLICIT NONE
7
8c=======================================================================
9c
10c   Auteurs:  P. Le Van , F. Hourdin .
11c   -------
12c
13c   Objet:
14c   ------
15c
16c    *******************************************************************
17c  .... calcul de la vitesse verticale aux niveaux sigma  ....
18c    *******************************************************************
19c     convm   est un argument  d'entree pour le s-pg  ......
20c       w     est un argument de sortie pour le s-pg  ......
21c
22c    la vitesse verticale est orientee de  haut en bas .
23c    au sol, au niveau sigma(1),   w(i,j,1) = 0.
24c    au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi
25c    egale a 0. et n'est pas stockee dans le tableau w  .
26c
27c
28c=======================================================================
29
30#include "dimensions.h"
31#include "paramet.h"
32#include "comvert.h"
33
34      REAL w(ip1jmp1,llm),convm(ip1jmp1,llm)
35      INTEGER   l, ij
36
37
38
39      DO 2  l = 1,llmm1
40
41      DO 1 ij = 1,ip1jmp1
42      w( ij, l+1 ) = convm( ij, l+1 ) - bp(l+1) * convm( ij, 1 )
43   1  CONTINUE
44
45   2  CONTINUE
46
47      DO 5 ij  = 1,ip1jmp1
48      w(ij,1)  = 0.
495     CONTINUE
50
51      RETURN
52      END
Note: See TracBrowser for help on using the repository browser.