dyn3d/Dissipation/divgrad2.f

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_gam_m, only: laplacien_gam
    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	module divgrad2_m
2
3	IMPLICIT NONE
4
5	contains
6
7	SUBROUTINE divgrad2(klevel, h, deltapres, lh, divgra, cdivh)
8
9	! 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_gam_m, only: laplacien_gam
17	USE laplacien_m, ONLY: laplacien
18	USE paramet_m, ONLY: ip1jmp1
19
20	INTEGER, intent(in):: klevel
21	REAL, intent(in):: h(ip1jmp1, klevel), deltapres(ip1jmp1, klevel)
22	integer, intent(in):: lh
23	REAL, intent(out):: divgra(ip1jmp1, klevel)
24	real, intent(in):: cdivh
25
26	! Variables locales
27	REAL sqrtps(ip1jmp1, klevel)
28	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