trunk/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	guez	91	module vitvert_m
2	guez	3
3	guez	91	IMPLICIT NONE
4	guez	3
5	guez	91	contains
6	guez	3
7	guez	91	SUBROUTINE vitvert(convm, w)
8	guez	3
9	guez	91	! From libf/dyn3d/vitvert.F, version 1.1.1.1, 2004/05/19 12:53:05
10			! Authors: P. Le Van, F. Hourdin
11	guez	3
12	guez	91	! Objet : calcul de la vitesse verticale aux niveaux sigma
13	guez	3
14	guez	91	! 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	guez	3
19	guez	91	USE dimens_m, ONLY : llm
20			USE disvert_m, ONLY : bp
21			USE paramet_m, ONLY : ip1jmp1
22	guez	3
23	guez	91	real, intent(in):: convm(ip1jmp1, llm)
24			REAL, intent(out):: w(ip1jmp1, llm)
25	guez	3
26	guez	91	! Local:
27			INTEGER l
28	guez	3
29	guez	91	!------------------------------------------------------
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