23 |
USE comconst, ONLY: dtvr |
USE comconst, ONLY: dtvr |
24 |
use comdissnew, only: nitergdiv, nitergrot, niterh, tetagdiv, tetagrot, & |
use comdissnew, only: nitergdiv, nitergrot, niterh, tetagdiv, tetagrot, & |
25 |
tetatemp |
tetatemp |
26 |
USE comvert, ONLY: preff, presnivs |
USE disvert_m, ONLY: preff, presnivs |
27 |
USE conf_gcm_m, ONLY: iperiod |
USE conf_gcm_m, ONLY: iperiod |
28 |
USE dimens_m, ONLY: iim, jjm |
USE dimens_m, ONLY: iim, jjm |
29 |
use divgrad2_m, only: divgrad2 |
use divgrad2_m, only: divgrad2 |
45 |
|
|
46 |
PRINT *, 'Call sequence information: inidissip' |
PRINT *, 'Call sequence information: inidissip' |
47 |
call random_seed(size=seed_size) |
call random_seed(size=seed_size) |
48 |
call random_seed(put=(/(0, ii = 1, seed_size)/)) |
call random_seed(put=(/(1, ii = 1, seed_size)/)) |
49 |
|
|
50 |
PRINT *, 'Calcul des valeurs propres de divgrad' |
PRINT *, 'Calcul des valeurs propres de divgrad' |
51 |
deltap = 1. |
deltap = 1. |
52 |
call random_number(zh) |
call random_number(zh) |
53 |
zh = zh - 0.5 |
zh = zh - 0.5 |
54 |
CALL filtreg(zh, jjp1, 1, 2, 1, .TRUE.) |
CALL filtreg(zh, direct = .true., intensive = .true.) |
55 |
|
|
56 |
DO l = 1, 50 |
DO l = 1, 50 |
57 |
CALL divgrad2(1, zh, deltap, niterh, divgra, -1.) |
CALL divgrad2(1, zh, deltap, niterh, divgra, -1.) |
65 |
PRINT *, 'Calcul des valeurs propres de gradiv' |
PRINT *, 'Calcul des valeurs propres de gradiv' |
66 |
call random_number(zu) |
call random_number(zu) |
67 |
zu = zu - 0.5 |
zu = zu - 0.5 |
68 |
CALL filtreg(zu, jjp1, 1, 2, 1, .TRUE.) |
CALL filtreg(zu, direct = .true., intensive = .true.) |
69 |
call random_number(zv) |
call random_number(zv) |
70 |
zv = zv - 0.5 |
zv = zv - 0.5 |
71 |
CALL filtreg(zv, jjm, 1, 2, 1, .FALSE.) |
CALL filtreg(zv, direct = .true., intensive = .true.) |
72 |
|
|
73 |
DO l = 1, 50 |
DO l = 1, 50 |
74 |
CALL gradiv2(zu, zv, nitergdiv, gx, gy, -1.) |
CALL gradiv2(zu, zv, nitergdiv, gx, gy, -1.) |
83 |
PRINT *, 'Calcul des valeurs propres de nxgrarot' |
PRINT *, 'Calcul des valeurs propres de nxgrarot' |
84 |
call random_number(zu) |
call random_number(zu) |
85 |
zu = zu - 0.5 |
zu = zu - 0.5 |
86 |
CALL filtreg(zu, jjp1, 1, 2, 1, .TRUE.) |
CALL filtreg(zu, direct = .true., intensive = .true.) |
87 |
call random_number(zv) |
call random_number(zv) |
88 |
zv = zv - 0.5 |
zv = zv - 0.5 |
89 |
CALL filtreg(zv, jjm, 1, 2, 1, .FALSE.) |
CALL filtreg(zv, direct = .true., intensive = .true.) |
90 |
|
|
91 |
DO l = 1, 50 |
DO l = 1, 50 |
92 |
CALL nxgraro2(zu, zv, nitergrot, gx, gy, -1.) |
CALL nxgraro2(zu, zv, nitergrot, gx, gy, -1.) |