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 dimens_m, 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 |