Sources/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	guez	91	module convmas_m
2	guez	3
3	guez	81	IMPLICIT NONE
4	guez	3
5	guez	91	contains
6	guez	3
7	guez	91	SUBROUTINE convmas(pbaru, pbarv, convm)
8	guez	3
9	guez	91	! From LMDZ4/libf/dyn3d/convmas.F, version 1.1.1.1, 2004/05/19 12:53:07
10	guez	3
11	guez	91	USE dimens_m, ONLY: llm
12			USE paramet_m, ONLY: ip1jm, ip1jmp1, jjp1, llmm1
13			USE filtreg_m, ONLY: filtreg
14	guez	3
15	guez	91	! Authors: P. Le Van, F. Hourdin
16			! Objet: calcul de la convergence du flux de masse aux niveaux p
17	guez	3
18	guez	91	! 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	guez	3
22	guez	91	REAL, INTENT(IN):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
23			REAL, INTENT(OUT):: convm(ip1jmp1, llm)
24	guez	81
25	guez	91	! Local:
26			INTEGER l, ij
27	guez	81
28	guez	91	!-----------------------------------------------------------------------
29	guez	81
30	guez	91	! Calcul de - (d(pbaru)/dx + d(pbarv)/dy) :
31			CALL convflu(pbaru, pbarv, llm, convm)
32	guez	81
33	guez	91	! Filtrage :
34			CALL filtreg(convm, jjp1, llm, 2, 2, .TRUE.)
35	guez	81
36	guez	91	! 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	guez	81
43	guez	91	END SUBROUTINE convmas
44	guez	81
45	guez	91	end module convmas_m