libf/dyn3d/vitvert.f90

!
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/vitvert.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
!
      SUBROUTINE vitvert ( convm , w )
c

c=======================================================================
c
c   Auteurs:  P. Le Van , F. Hourdin .
c   -------
c
c   Objet:
c   ------
c
c    *******************************************************************
c  .... calcul de la vitesse verticale aux niveaux sigma  ....
c    *******************************************************************
c     convm   est un argument  d'entree pour le s-pg  ......
c       w     est un argument de sortie pour le s-pg  ......
c
c    la vitesse verticale est orientee de  haut en bas .
c    au sol, au niveau sigma(1),   w(i,j,1) = 0.
c    au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi
c    egale a 0. et n'est pas stockee dans le tableau w  .
c
c
c=======================================================================

      use dimens_m
      use paramet_m
      use comvert
      IMPLICIT NONE

      REAL w(ip1jmp1,llm),convm(ip1jmp1,llm)
      INTEGER   l, ij


      DO 2  l = 1,llmm1

      DO 1 ij = 1,ip1jmp1
      w( ij, l+1 ) = convm( ij, l+1 ) - bp(l+1) * convm( ij, 1 )
   1  CONTINUE

   2  CONTINUE

      DO 5 ij  = 1,ip1jmp1
      w(ij,1)  = 0.
5     CONTINUE

      RETURN
      END
1	guez	3	!
2			! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/vitvert.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
3			!
4			SUBROUTINE vitvert ( convm , w )
5			c
6
7			c=======================================================================
8			c
9			c Auteurs: P. Le Van , F. Hourdin .
10			c -------
11			c
12			c Objet:
13			c ------
14			c
15			c *******************************************************************
16			c .... calcul de la vitesse verticale aux niveaux sigma ....
17			c *******************************************************************
18			c convm est un argument d'entree pour le s-pg ......
19			c w est un argument de sortie pour le s-pg ......
20			c
21			c la vitesse verticale est orientee de haut en bas .
22			c au sol, au niveau sigma(1), w(i,j,1) = 0.
23			c au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi
24			c egale a 0. et n'est pas stockee dans le tableau w .
25			c
26			c
27			c=======================================================================
28
29			use dimens_m
30			use paramet_m
31			use comvert
32			IMPLICIT NONE
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.
49			5 CONTINUE
50
51			RETURN
52			END