13 |
! pour l'opérateur filtre. |
! pour l'opérateur filtre. |
14 |
|
|
15 |
USE dimens_m, ONLY : iim, jjm |
USE dimens_m, ONLY : iim, jjm |
16 |
USE parafilt, ONLY : nfilun, nfilus, nfilvn, nfilvs |
USE parafilt, ONLY: matriceun, matriceus, matricevn, matricevs, matrinvn, & |
17 |
|
matrinvs |
18 |
USE coefils, ONLY : jfiltnu, jfiltnv, jfiltsu, jfiltsv, sddu, sddv, & |
USE coefils, ONLY : jfiltnu, jfiltnv, jfiltsu, jfiltsv, sddu, sddv, & |
19 |
unsddu, unsddv |
unsddu, unsddv |
20 |
|
|
49 |
|
|
50 |
INTEGER jdfil1, jdfil2, jffil1, jffil2, jdfil, jffil |
INTEGER jdfil1, jdfil2, jffil1, jffil2, jdfil, jffil |
51 |
INTEGER i, j, l, k |
INTEGER i, j, l, k |
|
REAL matriceun, matriceus, matricevn, matricevs, matrinvn, matrinvs |
|
|
COMMON /matrfil/matriceun(iim, iim, nfilun), matriceus(iim, iim, nfilus), & |
|
|
matricevn(iim, iim, nfilvn), matricevs(iim, iim, nfilvs), & |
|
|
matrinvn(iim, iim, nfilun), matrinvs(iim, iim, nfilus) |
|
52 |
REAL eignq(iim), sdd1(iim), sdd2(iim) |
REAL eignq(iim), sdd1(iim), sdd2(iim) |
53 |
INTEGER hemisph |
INTEGER hemisph |
54 |
|
|
153 |
END DO |
END DO |
154 |
DO i = 1, iim |
DO i = 1, iim |
155 |
DO k = 1, iim |
DO k = 1, iim |
156 |
eignq(k) = eignq(k) + matriceun(k, i, j)*champ(i, j, l) |
eignq(k) = eignq(k) + matriceun(k, i, j) & |
157 |
|
* champ(i, j, l) |
158 |
END DO |
END DO |
159 |
END DO |
END DO |
160 |
ELSE |
ELSE |
163 |
END DO |
END DO |
164 |
DO i = 1, iim |
DO i = 1, iim |
165 |
DO k = 1, iim |
DO k = 1, iim |
166 |
eignq(k) = eignq(k) + matricevn(k, i, j)*champ(i, j, l) |
eignq(k) = eignq(k) + matricevn(k, i, j) & |
167 |
|
* champ(i, j, l) |
168 |
END DO |
END DO |
169 |
END DO |
END DO |
170 |
END IF |
END IF |