11 |
|
|
12 |
! Calcul des composantes covariantes en x et y du gradient de g. |
! Calcul des composantes covariantes en x et y du gradient de g. |
13 |
|
|
14 |
USE dimens_m, ONLY: iim |
USE dimens_m, ONLY: iim, jjm |
|
USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1 |
|
15 |
|
|
16 |
INTEGER, intent(in):: klevel |
INTEGER, intent(in):: klevel |
17 |
REAL, intent(in):: g(ip1jmp1, klevel) |
REAL, intent(in):: g(iim + 1, jjm + 1, klevel) |
18 |
REAL, intent(out):: gx(ip1jmp1, klevel) , gy(ip1jm, klevel) |
REAL, intent(out):: gx(iim + 1, jjm + 1, klevel) , gy(iim + 1, jjm, klevel) |
19 |
|
|
20 |
! Local: |
! Local: |
21 |
INTEGER l, ij |
INTEGER i, j |
22 |
|
|
23 |
!---------------------------------------------------------------- |
!---------------------------------------------------------------- |
24 |
|
|
25 |
DO l = 1, klevel |
forall (i = 1:iim) gx(i, :, :) = g(i + 1, :, :) - g(i, :, :) |
26 |
DO ij = 1, ip1jmp1 - 1 |
gx(iim + 1, :, :)= gx(1, :, :) |
27 |
gx(ij, l) = g(ij +1, l) - g(ij, l) |
|
28 |
end DO |
forall (j = 1:jjm) gy(:, j, :) = g(:, j, :) - g(:, j + 1, :) |
|
|
|
|
! 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 |
|
29 |
|
|
30 |
END SUBROUTINE grad |
END SUBROUTINE grad |
31 |
|
|