trunk/dyn3d/convmas.f

module convmas_m

  IMPLICIT NONE

contains

  SUBROUTINE convmas(pbaru, pbarv, convm)

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

    USE dimens_m, ONLY: llm
    USE paramet_m, ONLY: ip1jm, ip1jmp1, jjp1, llmm1
    USE filtreg_m, ONLY: filtreg

    ! Authors: P. Le Van, F. Hourdin
    ! Objet: calcul de la convergence du flux de masse aux niveaux p

    ! Le calcul se fait de haut en bas, la convergence de masse au
    ! niveau p(llm+1) est égale à 0 et n'est pas stockée dans le
    ! tableau convm.

    REAL, INTENT(IN):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
    REAL, INTENT(OUT):: convm(ip1jmp1, llm)

    ! Local:
    INTEGER l, ij

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

    ! Calcul de - (d(pbaru)/dx + d(pbarv)/dy) :
    CALL convflu(pbaru, pbarv, llm, convm)

    ! Filtrage :
    CALL filtreg(convm, jjp1, llm, 2, 2, .TRUE.)

    ! Intégration de la convergence de masse de haut en bas :
    DO l = llmm1, 1, -1
       DO ij = 1, ip1jmp1
          convm(ij, l) = convm(ij, l) + convm(ij, l+1)
       END DO
    END DO

  END SUBROUTINE convmas

end module convmas_m
1	module convmas_m
2
3	IMPLICIT NONE
4
5	contains
6
7	SUBROUTINE convmas(pbaru, pbarv, convm)
8
9	! From LMDZ4/libf/dyn3d/convmas.F, version 1.1.1.1, 2004/05/19 12:53:07
10
11	USE dimens_m, ONLY: llm
12	USE paramet_m, ONLY: ip1jm, ip1jmp1, jjp1, llmm1
13	USE filtreg_m, ONLY: filtreg
14
15	! Authors: P. Le Van, F. Hourdin
16	! Objet: calcul de la convergence du flux de masse aux niveaux p
17
18	! Le calcul se fait de haut en bas, la convergence de masse au
19	! niveau p(llm+1) est égale à 0 et n'est pas stockée dans le
20	! tableau convm.
21
22	REAL, INTENT(IN):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
23	REAL, INTENT(OUT):: convm(ip1jmp1, llm)
24
25	! Local:
26	INTEGER l, ij
27
28	!-----------------------------------------------------------------------
29
30	! Calcul de - (d(pbaru)/dx + d(pbarv)/dy) :
31	CALL convflu(pbaru, pbarv, llm, convm)
32
33	! Filtrage :
34	CALL filtreg(convm, jjp1, llm, 2, 2, .TRUE.)
35
36	! Intégration de la convergence de masse de haut en bas :
37	DO l = llmm1, 1, -1
38	DO ij = 1, ip1jmp1
39	convm(ij, l) = convm(ij, l) + convm(ij, l+1)
40	END DO
41	END DO
42
43	END SUBROUTINE convmas
44
45	end module convmas_m