trunk/dyn3d/grad.f

module grad_m

  IMPLICIT NONE

contains

  SUBROUTINE grad(klevel, g, gx, gy)

    ! From LMDZ4/libf/dyn3d/grad.F, version 1.1.1.1 2004/05/19 12:53:05
    ! P. Le Van

    ! Calcul des composantes covariantes en x et y du gradient de g.

    USE dimens_m, ONLY: iim, jjm

    INTEGER, intent(in):: klevel
    REAL, intent(in):: g(iim + 1, jjm + 1, klevel)
    REAL, intent(out):: gx(iim + 1, jjm + 1, klevel) , gy(iim + 1, jjm, klevel)

    ! Local:
    INTEGER i, j

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

    forall (i = 1:iim) gx(i, :, :) = g(i + 1, :, :) - g(i, :, :)
    gx(iim + 1, :, :)= gx(1, :, :)

    forall (j = 1:jjm) gy(:, j, :) = g(:, j, :) - g(:, j + 1, :)

  END SUBROUTINE grad

end module grad_m
1	guez	60	module grad_m
2
3			IMPLICIT NONE
4
5			contains
6
7	guez	65	SUBROUTINE grad(klevel, g, gx, gy)
8	guez	60
9			! From LMDZ4/libf/dyn3d/grad.F, version 1.1.1.1 2004/05/19 12:53:05
10			! P. Le Van
11
12			! Calcul des composantes covariantes en x et y du gradient de g.
13
14	guez	103	USE dimens_m, ONLY: iim, jjm
15	guez	60
16			INTEGER, intent(in):: klevel
17	guez	103	REAL, intent(in):: g(iim + 1, jjm + 1, klevel)
18			REAL, intent(out):: gx(iim + 1, jjm + 1, klevel) , gy(iim + 1, jjm, klevel)
19	guez	60
20			! Local:
21	guez	103	INTEGER i, j
22	guez	60
23			!----------------------------------------------------------------
24
25	guez	103	forall (i = 1:iim) gx(i, :, :) = g(i + 1, :, :) - g(i, :, :)
26			gx(iim + 1, :, :)= gx(1, :, :)
27	guez	60
28	guez	103	forall (j = 1:jjm) gy(:, j, :) = g(:, j, :) - g(:, j + 1, :)
29	guez	60
30			END SUBROUTINE grad
31
32			end module grad_m