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.) |