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 |
use rotatf_m, only: rotatf |
18 |
|
|
50 |
CALL laplacien_rotgam(klevel, rot, rot) |
CALL laplacien_rotgam(klevel, rot, rot) |
51 |
ENDDO |
ENDDO |
52 |
|
|
53 |
CALL filtreg(rot, direct = .true., intensive = .true.) |
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 |