source: LMDZ5/trunk/libf/dyn3d/vitvert.F @ 1907

Last change on this file since 1907 was 1907, checked in by lguez, 10 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.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.