dyn3d/Dissipation/divgrad2.f90

module divgrad2_m

  IMPLICIT NONE

contains

  SUBROUTINE divgrad2(klevel, h, deltapres, lh, divgra, cdivh)

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

    ! Calcul de div(grad) de (pext * h)

    USE comgeom, ONLY: cuvscvgam2, cvuscugam2, unsair_gam2, unsapolnga2, &
         unsapolsga2
    use dimensions, only: iim, jjm
    use laplacien_gam_m, only: laplacien_gam
    USE laplacien_m, ONLY: laplacien

    INTEGER, intent(in):: klevel
    REAL, intent(in), dimension(iim + 1, jjm + 1, klevel):: h, deltapres
    integer, intent(in):: lh
    REAL, intent(out):: divgra(iim + 1, jjm + 1, klevel)
    real, intent(in):: cdivh

    ! Variables locales
    REAL sqrtps(iim + 1, jjm + 1, klevel)
    INTEGER iter

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

    divgra = h
    CALL laplacien(klevel, divgra)
    sqrtps = SQRT(deltapres)
    divgra = divgra * sqrtps

    ! Itération de l'opérateur laplacien_gam
    DO iter = 1, lh - 2
       CALL laplacien_gam(klevel, cuvscvgam2, cvuscugam2, unsair_gam2, &
            unsapolnga2, unsapolsga2, divgra, divgra)
    ENDDO

    divgra = divgra * sqrtps
    CALL laplacien(klevel, divgra)
    divgra = (-1.)**lh * cdivh * divgra / deltapres

  END SUBROUTINE divgrad2

end module divgrad2_m
1	guez	65	module divgrad2_m
2	guez	3
3	guez	65	IMPLICIT NONE
4	guez	3
5	guez	65	contains
6	guez	3
7	guez	65	SUBROUTINE divgrad2(klevel, h, deltapres, lh, divgra, cdivh)
8	guez	3
9	guez	266	! From LMDZ4/libf/dyn3d/divgrad2.F, version 1.1.1.1, 2004/05/19 12:53:06
10	guez	65	! P. Le Van
11
12			! Calcul de div(grad) de (pext * h)
13
14			USE comgeom, ONLY: cuvscvgam2, cvuscugam2, unsair_gam2, unsapolnga2, &
15			unsapolsga2
16	guez	266	use dimensions, only: iim, jjm
17	guez	207	use laplacien_gam_m, only: laplacien_gam
18	guez	65	USE laplacien_m, ONLY: laplacien
19
20			INTEGER, intent(in):: klevel
21	guez	266	REAL, intent(in), dimension(iim + 1, jjm + 1, klevel):: h, deltapres
22	guez	65	integer, intent(in):: lh
23	guez	266	REAL, intent(out):: divgra(iim + 1, jjm + 1, klevel)
24	guez	65	real, intent(in):: cdivh
25
26			! Variables locales
27	guez	266	REAL sqrtps(iim + 1, jjm + 1, klevel)
28	guez	65	INTEGER iter
29
30			!-----------------------------------------------------------------
31
32			divgra = h
33			CALL laplacien(klevel, divgra)
34			sqrtps = SQRT(deltapres)
35			divgra = divgra * sqrtps
36
37			! Itération de l'opérateur laplacien_gam
38			DO iter = 1, lh - 2
39			CALL laplacien_gam(klevel, cuvscvgam2, cvuscugam2, unsair_gam2, &
40			unsapolnga2, unsapolsga2, divgra, divgra)
41			ENDDO
42
43			divgra = divgra * sqrtps
44			CALL laplacien(klevel, divgra)
45			divgra = (-1.)*lh cdivh * divgra / deltapres
46
47			END SUBROUTINE divgrad2
48
49			end module divgrad2_m