Sources/dyn3d/vitvert.f

!
! $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	!
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