trunk/dyn3d/dteta1.f

module dteta1_m

  IMPLICIT NONE

contains

  SUBROUTINE dteta1(teta, pbaru, pbarv, dteta)

    ! From LMDZ4/libf/dyn3d/dteta1.F, version 1.1.1.1, 2004/05/19 12:53:06
    ! Authors: P. Le Van, F. Forget

    ! Calcul du terme de convergence horizontale du flux d'enthalpie
    ! potentielle.

    USE dimens_m, ONLY: iim, jjm, llm
    USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmp1, jjp1
    USE filtreg_m, ONLY: filtreg

    REAL, intent(in):: teta(ip1jmp1, llm)
    REAL, intent(in):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
    REAL, intent(out):: dteta(iim + 1, jjm + 1, llm)

    ! Local:
    INTEGER l, ij
    REAL hbyv(ip1jm, llm), hbxu(ip1jmp1, llm)

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

    DO l = 1, llm
       DO ij = iip2, ip1jm - 1
          hbxu(ij, l) = pbaru(ij, l) * 0.5 * (teta(ij, l) + teta(ij + 1, l))
       end DO

       DO ij = iip1+ iip1, ip1jm, iip1
          hbxu(ij, l) = hbxu(ij - iim, l)
       end DO

       DO ij = 1, ip1jm
          hbyv(ij, l)= pbarv(ij, l) * 0.5 * (teta(ij, l) + teta(ij + iip1, l))
       end DO
    end DO

    CALL convflu(hbxu, hbyv, llm, dteta)

    ! stockage dans dh de la convergence horizontale filtrée du flux
    ! d'enthalpie potentielle
    CALL filtreg(dteta, direct = .true., intensive = .false.)

  END SUBROUTINE dteta1

end module dteta1_m
1	guez	88	module dteta1_m
2	guez	27
3	guez	88	IMPLICIT NONE
4	guez	3
5	guez	88	contains
6	guez	3
7	guez	88	SUBROUTINE dteta1(teta, pbaru, pbarv, dteta)
8	guez	3
9	guez	88	! From LMDZ4/libf/dyn3d/dteta1.F, version 1.1.1.1, 2004/05/19 12:53:06
10			! Authors: P. Le Van, F. Forget
11	guez	3
12	guez	88	! Calcul du terme de convergence horizontale du flux d'enthalpie
13			! potentielle.
14	guez	3
15	guez	104	USE dimens_m, ONLY: iim, jjm, llm
16	guez	88	USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmp1, jjp1
17			USE filtreg_m, ONLY: filtreg
18	guez	3
19	guez	88	REAL, intent(in):: teta(ip1jmp1, llm)
20			REAL, intent(in):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
21	guez	104	REAL, intent(out):: dteta(iim + 1, jjm + 1, llm)
22	guez	3
23	guez	88	! Local:
24			INTEGER l, ij
25			REAL hbyv(ip1jm, llm), hbxu(ip1jmp1, llm)
26	guez	3
27	guez	88	!----------------------------------------------------------------
28	guez	3
29	guez	88	DO l = 1, llm
30			DO ij = iip2, ip1jm - 1
31			hbxu(ij, l) = pbaru(ij, l) * 0.5 * (teta(ij, l) + teta(ij + 1, l))
32			end DO
33	guez	3
34	guez	88	DO ij = iip1+ iip1, ip1jm, iip1
35			hbxu(ij, l) = hbxu(ij - iim, l)
36			end DO
37	guez	3
38	guez	88	DO ij = 1, ip1jm
39			hbyv(ij, l)= pbarv(ij, l) * 0.5 * (teta(ij, l) + teta(ij + iip1, l))
40			end DO
41			end DO
42	guez	3
43	guez	88	CALL convflu(hbxu, hbyv, llm, dteta)
44	guez	3
45	guez	88	! stockage dans dh de la convergence horizontale filtrée du flux
46			! d'enthalpie potentielle
47	guez	107	CALL filtreg(dteta, direct = .true., intensive = .false.)
48	guez	88
49			END SUBROUTINE dteta1
50
51			end module dteta1_m