dyn3d/Dissipation/gradiv2.f

module gradiv2_m

  IMPLICIT NONE

contains

  SUBROUTINE gradiv2(klevel, xcov, ycov, ld, gdx, gdy, cdivu)

    ! From LMDZ4/libf/dyn3d/gradiv2.F, version 1.1.1.1 2004/05/19 12:53:07
    ! P. Le Van
    ! Calcul de grad div du vecteur v.

    USE dimens_m, ONLY : llm
    USE paramet_m, ONLY : ip1jm, ip1jmp1, jjp1
    USE comgeom, ONLY : cuvscvgam1, cvuscugam1, unsair_gam1, unsapolnga1, &
         unsapolsga1
    USE filtreg_m, ONLY : filtreg

    INTEGER, intent(in):: klevel

    ! composantes covariantes de v:
    REAL, intent(in):: xcov(ip1jmp1, klevel), ycov(ip1jm, klevel)

    integer, intent(in):: ld
    REAL, intent(out):: gdx(ip1jmp1, klevel), gdy(ip1jm, klevel)
    real, intent(in):: cdivu

    ! Variables locales :
    REAL div(ip1jmp1, llm)
    REAL nugrads
    INTEGER iter

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

    gdx = xcov
    gdy = ycov

    CALL divergf(klevel, gdx, gdy, div)

    IF (ld > 1) THEN
       CALL laplacien(klevel, div, div)

       ! Itération de l'opérateur laplacien_gam
       DO iter = 1, ld -2
          CALL laplacien_gam(klevel, cuvscvgam1, cvuscugam1, unsair_gam1, &
               unsapolnga1, unsapolsga1, div, div)
       END DO
    ENDIF

    CALL filtreg(div, jjp1, klevel, 2, 1, .TRUE., 1)
    CALL grad(klevel, div, gdx, gdy)
    nugrads = (-1.)**ld * cdivu

    gdx = gdx * nugrads
    gdy = gdy * nugrads

  END SUBROUTINE gradiv2

end module gradiv2_m
1	module gradiv2_m
2
3	IMPLICIT NONE
4
5	contains
6
7	SUBROUTINE gradiv2(klevel, xcov, ycov, ld, gdx, gdy, cdivu)
8
9	! From LMDZ4/libf/dyn3d/gradiv2.F, version 1.1.1.1 2004/05/19 12:53:07
10	! P. Le Van
11	! Calcul de grad div du vecteur v.
12
13	USE dimens_m, ONLY : llm
14	USE paramet_m, ONLY : ip1jm, ip1jmp1, jjp1
15	USE comgeom, ONLY : cuvscvgam1, cvuscugam1, unsair_gam1, unsapolnga1, &
16	unsapolsga1
17	USE filtreg_m, ONLY : filtreg
18
19	INTEGER, intent(in):: klevel
20
21	! composantes covariantes de v:
22	REAL, intent(in):: xcov(ip1jmp1, klevel), ycov(ip1jm, klevel)
23
24	integer, intent(in):: ld
25	REAL, intent(out):: gdx(ip1jmp1, klevel), gdy(ip1jm, klevel)
26	real, intent(in):: cdivu
27
28	! Variables locales :
29	REAL div(ip1jmp1, llm)
30	REAL nugrads
31	INTEGER iter
32
33	!--------------------------------------------------------------
34
35	gdx = xcov
36	gdy = ycov
37
38	CALL divergf(klevel, gdx, gdy, div)
39
40	IF (ld > 1) THEN
41	CALL laplacien(klevel, div, div)
42
43	! Itération de l'opérateur laplacien_gam
44	DO iter = 1, ld -2
45	CALL laplacien_gam(klevel, cuvscvgam1, cvuscugam1, unsair_gam1, &
46	unsapolnga1, unsapolsga1, div, div)
47	END DO
48	ENDIF
49
50	CALL filtreg(div, jjp1, klevel, 2, 1, .TRUE., 1)
51	CALL grad(klevel, div, gdx, gdy)
52	nugrads = (-1.)*ld cdivu
53
54	gdx = gdx * nugrads
55	gdy = gdy * nugrads
56
57	END SUBROUTINE gradiv2
58
59	end module gradiv2_m