trunk/dyn3d/laplacien.f

module laplacien_m

  IMPLICIT NONE

contains

  SUBROUTINE laplacien(klevel, teta)

    ! From LMDZ4/libf/dyn3d/laplacien.F, version 1.1.1.1 2004/05/19 12:53:06
    ! P. Le Van
    ! Calcul de div(grad) de teta.

    use dimens_m, only: iim, jjm
    use divergf_m, only: divergf
    use filtreg_m, only: filtreg
    use grad_m, only: grad
    USE paramet_m, ONLY: ip1jm, ip1jmp1, jjp1

    INTEGER, intent(in):: klevel
    REAL, intent(inout):: teta(iim + 1, jjm + 1, klevel)

    ! Variables locales:
    REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel)

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

    CALL filtreg(teta, jjp1, klevel, 2, 1, .TRUE.)
    CALL grad(klevel, teta, ghx, ghy)
    CALL divergf(klevel, ghx, ghy, teta)

  END SUBROUTINE laplacien

end module laplacien_m
1	module laplacien_m
2
3	IMPLICIT NONE
4
5	contains
6
7	SUBROUTINE laplacien(klevel, teta)
8
9	! From LMDZ4/libf/dyn3d/laplacien.F, version 1.1.1.1 2004/05/19 12:53:06
10	! P. Le Van
11	! Calcul de div(grad) de teta.
12
13	use dimens_m, only: iim, jjm
14	use divergf_m, only: divergf
15	use filtreg_m, only: filtreg
16	use grad_m, only: grad
17	USE paramet_m, ONLY: ip1jm, ip1jmp1, jjp1
18
19	INTEGER, intent(in):: klevel
20	REAL, intent(inout):: teta(iim + 1, jjm + 1, klevel)
21
22	! Variables locales:
23	REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel)
24
25	!-----------------------------------------------------------------
26
27	CALL filtreg(teta, jjp1, klevel, 2, 1, .TRUE.)
28	CALL grad(klevel, teta, ghx, ghy)
29	CALL divergf(klevel, ghx, ghy, teta)
30
31	END SUBROUTINE laplacien
32
33	end module laplacien_m