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 laplacien_m, ONLY: laplacien
    USE paramet_m, ONLY: ip1jmp1

    INTEGER, intent(in):: klevel
    REAL, intent(in):: h(ip1jmp1, klevel), deltapres(ip1jmp1, klevel)
    integer, intent(in):: lh
    REAL, intent(out):: divgra(ip1jmp1, klevel)
    real, intent(in):: cdivh

    ! Variables locales
    REAL sqrtps(ip1jmp1, 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	65	! From LMDZ4/libf/dyn3d/divgrad2.F, version 1.1.1.1 2004/05/19 12:53:06
10			! 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			USE laplacien_m, ONLY: laplacien
17			USE paramet_m, ONLY: ip1jmp1
18
19			INTEGER, intent(in):: klevel
20			REAL, intent(in):: h(ip1jmp1, klevel), deltapres(ip1jmp1, klevel)
21			integer, intent(in):: lh
22			REAL, intent(out):: divgra(ip1jmp1, klevel)
23			real, intent(in):: cdivh
24
25			! Variables locales
26			REAL sqrtps(ip1jmp1, klevel)
27			INTEGER iter
28
29			!-----------------------------------------------------------------
30
31			divgra = h
32			CALL laplacien(klevel, divgra)
33			sqrtps = SQRT(deltapres)
34			divgra = divgra * sqrtps
35
36			! Itération de l'opérateur laplacien_gam
37			DO iter = 1, lh - 2
38			CALL laplacien_gam(klevel, cuvscvgam2, cvuscugam2, unsair_gam2, &
39			unsapolnga2, unsapolsga2, divgra, divgra)
40			ENDDO
41
42			divgra = divgra * sqrtps
43			CALL laplacien(klevel, divgra)
44			divgra = (-1.)*lh cdivh * divgra / deltapres
45
46			END SUBROUTINE divgrad2
47
48			end module divgrad2_m