12 |
! Calcul de nxgrad(rot) du vecteur v |
! Calcul de nxgrad(rot) du vecteur v |
13 |
|
|
14 |
USE dimens_m, ONLY: iim, jjm |
USE dimens_m, ONLY: iim, jjm |
15 |
USE filtreg_m, ONLY: filtreg |
USE filtreg_v_m, ONLY: filtreg_v |
16 |
use nr_util, only: assert, assert_eq |
use nr_util, only: assert, assert_eq |
17 |
|
use rotatf_m, only: rotatf |
18 |
|
|
19 |
! Composantes covariantes de v : |
! Composantes covariantes de v : |
20 |
REAL, intent(in):: xcov(:, :, :) ! (iim + 1, jjm + 1, :) |
REAL, intent(in):: xcov(:, :, :) ! (iim + 1, jjm + 1, :) |
50 |
CALL laplacien_rotgam(klevel, rot, rot) |
CALL laplacien_rotgam(klevel, rot, rot) |
51 |
ENDDO |
ENDDO |
52 |
|
|
53 |
CALL filtreg(rot, jjm, klevel, 2, 1, .FALSE.) |
CALL filtreg_v(rot, intensive = .true.) |
54 |
CALL nxgrad(klevel, rot, grx, gry) |
CALL nxgrad(klevel, rot, grx, gry) |
55 |
|
|
56 |
nugradrs = (-1.)**lr * crot |
nugradrs = (-1.)**lr * crot |