Sources/dyn3d/vitvert.f

module vitvert_m

  IMPLICIT NONE

contains

  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 orientée de haut en bas. Au sol, au
    ! niveau sigma(1), w(i, j, 1) = 0. Au sommet, au niveau
    ! sigma(llm+1), la vitesse verticale est aussi égale à 0 et n'est
    ! pas stockée dans le tableau w.

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

    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

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