Sources/dyn3d/vitvert.f

SUBROUTINE vitvert(convm , w)

  ! From libf/dyn3d/vitvert.F, version 1.1.1.1 2004/05/19 12:53:05
  ! Authors: P. Le Van , F. Hourdin

  ! Objet : calcul de la vitesse verticale aux niveaux sigma

  ! La vitesse verticale est orientee de haut en bas .
  ! au sol, au niveau sigma(1), w(i, j, 1) = 0.
  ! au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi
  ! egale a 0. et n'est pas stockee dans le tableau w .

  USE dimens_m, ONLY : llm
  USE paramet_m, ONLY : ip1jmp1
  USE disvert_m, ONLY : bp

  IMPLICIT NONE

  real, intent(in):: convm(ip1jmp1, llm)
  REAL, intent(out):: w(ip1jmp1, llm)

  ! Local:
  INTEGER l

  !------------------------------------------------------

  forall (l = 2: llm) w(:, l) = convm(:, l) - bp(l) * convm(:, 1)
  w(:, 1) = 0.

END SUBROUTINE vitvert
1	guez	43	SUBROUTINE vitvert(convm , w)
2	guez	3
3	guez	43	! From libf/dyn3d/vitvert.F, version 1.1.1.1 2004/05/19 12:53:05
4			! Authors: P. Le Van , F. Hourdin
5	guez	3
6	guez	43	! Objet : calcul de la vitesse verticale aux niveaux sigma
7	guez	3
8	guez	43	! La vitesse verticale est orientee de haut en bas .
9			! au sol, au niveau sigma(1), w(i, j, 1) = 0.
10			! au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi
11			! egale a 0. et n'est pas stockee dans le tableau w .
12	guez	3
13	guez	43	USE dimens_m, ONLY : llm
14			USE paramet_m, ONLY : ip1jmp1
15	guez	66	USE disvert_m, ONLY : bp
16	guez	3
17	guez	43	IMPLICIT NONE
18	guez	3
19	guez	43	real, intent(in):: convm(ip1jmp1, llm)
20			REAL, intent(out):: w(ip1jmp1, llm)
21	guez	3
22	guez	43	! Local:
23			INTEGER l
24	guez	3
25	guez	43	!------------------------------------------------------
26	guez	3
27	guez	43	forall (l = 2: llm) w(:, l) = convm(:, l) - bp(l) * convm(:, 1)
28			w(:, 1) = 0.
29	guez	3
30	guez	43	END SUBROUTINE vitvert