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	module vitvert_m
2
3	IMPLICIT NONE
4
5	contains
6
7	pure function vitvert(convm)
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	! Purpose: Compute vertical speed at sigma levels.
13
14	! 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	USE dimens_m, ONLY : llm
19	USE disvert_m, ONLY : bp
20
21	real, intent(in):: convm(:, :, :) ! (iim + 1, jjm + 1, llm)
22	REAL vitvert(size(convm, 1), size(convm, 2), size(convm, 3))
23
24	! Local:
25	INTEGER l
26
27	!------------------------------------------------------
28
29	forall (l = 2: llm) &
30	vitvert(:, :, l) = convm(:, :, l) - bp(l) * convm(:, :, 1)
31	vitvert(:, :, 1) = 0.
32
33	END function vitvert
34
35	end module vitvert_m