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: iim, jjm, 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(iim + 1, jjm + 1, llm)

    ! Local:
    INTEGER l

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

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

    ! Filtrage :
    CALL filtreg(convm, direct = .true., intensive = .false.)

    ! Intégration de la convergence de masse de haut en bas :
    DO l = llmm1, 1, -1
       convm(:, :, l) = convm(:, :, l) + convm(:, :, l+1)
    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	104	USE dimens_m, ONLY: iim, jjm, llm
12	guez	91	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	guez	104	REAL, INTENT(OUT):: convm(iim + 1, jjm + 1, llm)
24	guez	81
25	guez	91	! Local:
26	guez	104	INTEGER l
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	guez	107	CALL filtreg(convm, direct = .true., intensive = .false.)
35	guez	81
36	guez	91	! Intégration de la convergence de masse de haut en bas :
37			DO l = llmm1, 1, -1
38	guez	104	convm(:, :, l) = convm(:, :, l) + convm(:, :, l+1)
39	guez	91	END DO
40	guez	81
41	guez	91	END SUBROUTINE convmas
42	guez	81
43	guez	91	end module convmas_m