libf/dyn3d/pression.f90

module pression_m

  IMPLICIT NONE

contains

  SUBROUTINE pression(ngrid, ap, bp, ps, p)

    ! From dyn3d/pression.F, version 1.1.1.1 2004/05/19 12:53:07

    ! Authors : P. Le Van, F. Hourdin

    ! Calcule la pression p(l) aux différents niveaux l = 1 (niveau du
    ! sol) à l = llm +1, ces niveaux correspondant aux interfaces des (llm) 
    ! couches, avec  p(ij, llm +1) = 0 et p(ij, 1) = ps(ij).

    use dimens_m, only: llm

    INTEGER, intent(in):: ngrid
    REAL, intent(in):: ap(llm + 1), bp(llm + 1)
    real, intent(in):: ps(ngrid)
    real, intent(out):: p(ngrid, llm + 1) 

    INTEGER l

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

    forall(l = 1: llm + 1) p(:, l) = ap(l) + bp(l) * ps(:)

  END SUBROUTINE pression

end module pression_m
1	guez	3	module pression_m
2
3			IMPLICIT NONE
4
5			contains
6
7			SUBROUTINE pression(ngrid, ap, bp, ps, p)
8
9			! From dyn3d/pression.F, version 1.1.1.1 2004/05/19 12:53:07
10
11			! Authors : P. Le Van, F. Hourdin
12
13			! Calcule la pression p(l) aux différents niveaux l = 1 (niveau du
14			! sol) à l = llm +1, ces niveaux correspondant aux interfaces des (llm)
15			! couches, avec p(ij, llm +1) = 0 et p(ij, 1) = ps(ij).
16
17			use dimens_m, only: llm
18
19			INTEGER, intent(in):: ngrid
20			REAL, intent(in):: ap(llm + 1), bp(llm + 1)
21			real, intent(in):: ps(ngrid)
22			real, intent(out):: p(ngrid, llm + 1)
23
24			INTEGER l
25
26			!---------------------
27
28			forall(l = 1: llm + 1) p(:, l) = ap(l) + bp(l) * ps(:)
29
30			END SUBROUTINE pression
31
32			end module pression_m