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 l,ij,iter

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

    gdx = xcov
    gdy = ycov

    CALL divergf(klevel, gdx, gdy, div)

    IF(ld.GT.1) THEN
       CALL laplacien (klevel, div, div)

       ! Iteration de l'operateur laplacien_gam
       DO iter = 1, ld -2
          CALL laplacien_gam (klevel,cuvscvgam1,cvuscugam1,unsair_gam1, &
               unsapolnga1, unsapolsga1, div, div)
       ENDDO
    ENDIF

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

    DO l = 1, klevel
       DO ij = 1, ip1jmp1
          gdx(ij,l) = gdx(ij,l) * nugrads
       ENDDO
       DO ij = 1, ip1jm
          gdy(ij,l) = gdy(ij,l) * nugrads
       ENDDO
    ENDDO

  END SUBROUTINE gradiv2

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