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 |
30 |
use filtreg_m, only: filtreg |
use filtreg_scal_m, only: filtreg_scal |
31 |
|
use filtreg_v_m, only: filtreg_v |
32 |
use gradiv2_m, only: gradiv2 |
use gradiv2_m, only: gradiv2 |
33 |
use jumble, only: new_unit |
use jumble, only: new_unit |
34 |
use nxgraro2_m, only: nxgraro2 |
use nxgraro2_m, only: nxgraro2 |
|
USE paramet_m, ONLY: jjp1 |
|
35 |
|
|
36 |
! Variables local to the procedure: |
! Variables local to the procedure: |
37 |
REAL zvert(llm), max_zvert ! no dimension |
REAL zvert(llm), max_zvert ! no dimension |
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_scal(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_scal(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_v(zv, 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_scal(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_v(zv, 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.) |