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
    USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1

    INTEGER, intent(in):: klevel
    REAL, intent(in):: g(ip1jmp1, klevel)
    REAL, intent(out):: gx(ip1jmp1, klevel) , gy(ip1jm, klevel)

    ! Local:
    INTEGER l, ij

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

    DO l = 1, klevel
       DO ij = 1, ip1jmp1 - 1
          gx(ij, l) = g(ij +1, l) - g(ij, l)
       end DO

       ! correction pour gx(ip1, j, l)
       ! gx(iip1, j, l)= gx(1, j, l)
       DO ij = iip1, ip1jmp1, iip1
          gx(ij, l) = gx(ij -iim, l)
       end DO

       DO ij = 1, ip1jm
          gy(ij, l) = g(ij, l) - g(ij +iip1, l)
       end DO
    end DO

  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	65	USE dimens_m, ONLY: iim
15			USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1
16	guez	60
17			INTEGER, intent(in):: klevel
18	guez	65	REAL, intent(in):: g(ip1jmp1, klevel)
19			REAL, intent(out):: gx(ip1jmp1, klevel) , gy(ip1jm, klevel)
20	guez	60
21			! Local:
22			INTEGER l, ij
23
24			!----------------------------------------------------------------
25
26			DO l = 1, klevel
27			DO ij = 1, ip1jmp1 - 1
28	guez	65	gx(ij, l) = g(ij +1, l) - g(ij, l)
29	guez	60	end DO
30
31	guez	65	! correction pour gx(ip1, j, l)
32			! gx(iip1, j, l)= gx(1, j, l)
33	guez	60	DO ij = iip1, ip1jmp1, iip1
34	guez	65	gx(ij, l) = gx(ij -iim, l)
35	guez	60	end DO
36
37			DO ij = 1, ip1jm
38	guez	65	gy(ij, l) = g(ij, l) - g(ij +iip1, l)
39	guez	60	end DO
40			end DO
41
42			END SUBROUTINE grad
43
44			end module grad_m