trunk/dyn3d/grad.f90

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 dimensions, 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	module grad_m
2
3	IMPLICIT NONE
4
5	contains
6
7	SUBROUTINE grad(klevel, g, gx, gy)
8
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	USE dimensions, ONLY: iim, jjm
15
16	INTEGER, intent(in):: klevel
17	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
20	! Local:
21	INTEGER i, j
22
23	!----------------------------------------------------------------
24
25	forall (i = 1:iim) gx(i, :, :) = g(i + 1, :, :) - g(i, :, :)
26	gx(iim + 1, :, :)= gx(1, :, :)
27
28	forall (j = 1:jjm) gy(:, j, :) = g(:, j, :) - g(:, j + 1, :)
29
30	END SUBROUTINE grad
31
32	end module grad_m