trunk/dyn3d/divergf.f

module divergf_m

  IMPLICIT NONE

contains

  SUBROUTINE divergf(klevel, x, y, div)

    ! From libf/dyn3d/divergf.F, version 1.1.1.1 2004/05/19 12:53:05

    ! P. Le Van

    ! Calcule la divergence à tous les niveaux d'un vecteur de
    ! composantes x et y. x et y sont des composantes covariantes.

    USE comgeom, ONLY: apoln, apols, cuvsurcv_2d, cvusurcu_2d, unsaire_2d
    USE dimens_m, ONLY: iim, jjm
    USE filtreg_m, ONLY: filtreg

    INTEGER, intent(in):: klevel
    REAL, intent(in):: x(iim + 1, jjm + 1, klevel), y(iim + 1, jjm, klevel)
    real, intent(out):: div(iim + 1, jjm + 1, klevel) ! in (unit of x, y) m-2

    ! Variables locales :

    INTEGER l, i, j

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

    DO l = 1, klevel
       forall (i = 2:iim + 1, j = 2:jjm) div(i, j, l) = cvusurcu_2d(i, j) &
            * x(i, j, l) - cvusurcu_2d(i - 1, j) * x(i - 1, j , l) &
            + cuvsurcv_2d(i, j - 1) * y(i, j - 1, l) - cuvsurcv_2d(i, j) &
            * y(i, j, l)

       div(1, 2:jjm, l) = div(iim + 1, 2:jjm, l)

       ! Calcul aux pôles 
       div(:, 1, l) = - SUM(cuvsurcv_2d(:iim, 1) * y(:iim, 1, l)) / apoln
       div(:, jjm + 1, l) = SUM(cuvsurcv_2d(:iim, jjm) * y(:iim, jjm, l)) &
            / apols
    end DO

    CALL filtreg(div, direct = .true., intensive = .false.)

    DO l = 1, klevel
       div(:, 2:jjm, l) = div(:, 2:jjm, l) * unsaire_2d(:, 2:jjm)
    ENDDO

  END SUBROUTINE divergf

end module divergf_m
1	guez	57	module divergf_m
2	guez	3
3	guez	57	IMPLICIT NONE
4	guez	3
5	guez	57	contains
6
7			SUBROUTINE divergf(klevel, x, y, div)
8
9	guez	61	! From libf/dyn3d/divergf.F, version 1.1.1.1 2004/05/19 12:53:05
10	guez	57
11			! P. Le Van
12	guez	106
13	guez	57	! Calcule la divergence à tous les niveaux d'un vecteur de
14			! composantes x et y. x et y sont des composantes covariantes.
15
16	guez	106	USE comgeom, ONLY: apoln, apols, cuvsurcv_2d, cvusurcu_2d, unsaire_2d
17			USE dimens_m, ONLY: iim, jjm
18	guez	65	USE filtreg_m, ONLY: filtreg
19	guez	57
20			INTEGER, intent(in):: klevel
21	guez	106	REAL, intent(in):: x(iim + 1, jjm + 1, klevel), y(iim + 1, jjm, klevel)
22			real, intent(out):: div(iim + 1, jjm + 1, klevel) ! in (unit of x, y) m-2
23	guez	57
24	guez	60	! Variables locales :
25	guez	57
26	guez	106	INTEGER l, i, j
27	guez	57
28			!------------------------------------------------------------
29
30			DO l = 1, klevel
31	guez	106	forall (i = 2:iim + 1, j = 2:jjm) div(i, j, l) = cvusurcu_2d(i, j) &
32			* x(i, j, l) - cvusurcu_2d(i - 1, j) * x(i - 1, j , l) &
33			+ cuvsurcv_2d(i, j - 1) * y(i, j - 1, l) - cuvsurcv_2d(i, j) &
34			* y(i, j, l)
35	guez	57
36	guez	106	div(1, 2:jjm, l) = div(iim + 1, 2:jjm, l)
37	guez	57
38			! Calcul aux pôles
39	guez	106	div(:, 1, l) = - SUM(cuvsurcv_2d(:iim, 1) * y(:iim, 1, l)) / apoln
40			div(:, jjm + 1, l) = SUM(cuvsurcv_2d(:iim, jjm) * y(:iim, jjm, l)) &
41			/ apols
42	guez	57	end DO
43
44	guez	107	CALL filtreg(div, direct = .true., intensive = .false.)
45	guez	57
46			DO l = 1, klevel
47	guez	106	div(:, 2:jjm, l) = div(:, 2:jjm, l) * unsaire_2d(:, 2:jjm)
48	guez	57	ENDDO
49
50			END SUBROUTINE divergf
51
52			end module divergf_m