10 |
! P. Le Van |
! P. Le Van |
11 |
! Calcul du gradient de la divergence du vecteur v. |
! Calcul du gradient de la divergence du vecteur v. |
12 |
|
|
13 |
USE dimens_m, ONLY : iim, jjm, llm |
USE comgeom, ONLY: cuvscvgam1, cvuscugam1, unsair_gam1, unsapolnga1, & |
|
use divergf_m, only: divergf |
|
|
USE comgeom, ONLY : cuvscvgam1, cvuscugam1, unsair_gam1, unsapolnga1, & |
|
14 |
unsapolsga1 |
unsapolsga1 |
15 |
USE filtreg_m, ONLY : filtreg |
USE dimens_m, ONLY: iim, jjm, llm |
16 |
|
use divergf_m, only: divergf |
17 |
|
USE filtreg_m, ONLY: filtreg |
18 |
use grad_m, only: grad |
use grad_m, only: grad |
19 |
|
use laplacien_m, only: laplacien |
20 |
use nr_util, only: assert_eq, assert |
use nr_util, only: assert_eq, assert |
21 |
|
|
22 |
! Composantes covariantes de v : |
! Composantes covariantes de v : |
44 |
CALL divergf(klevel, xcov, ycov, div) |
CALL divergf(klevel, xcov, ycov, div) |
45 |
|
|
46 |
IF (ld > 1) THEN |
IF (ld > 1) THEN |
47 |
CALL laplacien(klevel, div, div) |
CALL laplacien(klevel, div) |
48 |
|
|
49 |
! Itération de l'opérateur laplacien_gam |
! Itération de l'opérateur laplacien_gam |
50 |
DO iter = 1, ld -2 |
DO iter = 1, ld -2 |
53 |
END DO |
END DO |
54 |
ENDIF |
ENDIF |
55 |
|
|
56 |
CALL filtreg(div, jjm + 1, klevel, 2, 1, .TRUE., 1) |
CALL filtreg(div, jjm + 1, klevel, 2, 1, .TRUE.) |
57 |
CALL grad(klevel, div, gdx, gdy) |
CALL grad(klevel, div, gdx, gdy) |
58 |
nugrads = (-1.)**ld * cdivu |
nugrads = (-1.)**ld * cdivu |
59 |
|
|