trunk/dyn3d/vitvert.f

module vitvert_m

  IMPLICIT NONE

contains

  pure function vitvert(convm)

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

    ! Purpose: Compute vertical speed at sigma levels.

    ! Vertical speed is oriented from bottom to top. At ground-level
    ! sigma(1): vitvert(i, j, 1) = 0. At top-level sigma(llm + 1), vertical
    ! speed is 0 too and is not stored in vitvert.
    
    USE dimens_m, ONLY : llm
    USE disvert_m, ONLY : bp

    real, intent(in):: convm(:, :, :) ! (iim + 1, jjm + 1, llm)
    REAL vitvert(size(convm, 1), size(convm, 2), size(convm, 3))

    ! Local:
    INTEGER l

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

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

  END function 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	252	pure function vitvert(convm)
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	252	! Purpose: Compute vertical speed at sigma levels.
13	guez	3
14	guez	252	! Vertical speed is oriented from bottom to top. At ground-level
15			! sigma(1): vitvert(i, j, 1) = 0. At top-level sigma(llm + 1), vertical
16			! speed is 0 too and is not stored in vitvert.
17
18	guez	91	USE dimens_m, ONLY : llm
19			USE disvert_m, ONLY : bp
20	guez	3
21	guez	252	real, intent(in):: convm(:, :, :) ! (iim + 1, jjm + 1, llm)
22			REAL vitvert(size(convm, 1), size(convm, 2), size(convm, 3))
23	guez	3
24	guez	91	! Local:
25			INTEGER l
26	guez	3
27	guez	91	!------------------------------------------------------
28
29	guez	252	forall (l = 2: llm) &
30			vitvert(:, :, l) = convm(:, :, l) - bp(l) * convm(:, :, 1)
31			vitvert(:, :, 1) = 0.
32	guez	91
33	guez	252	END function vitvert
34	guez	91
35			end module vitvert_m