17 |
|
|
18 |
!----------------------------------------------------------------- |
!----------------------------------------------------------------- |
19 |
|
|
20 |
DO l = 1, size(champ, 3) |
forall (j = 1:size(champ, 2), l = 1:size(champ, 3)) |
21 |
DO j = 1, size(champ, 2) |
champ(:iim, j, l) = champ(:iim, j, l) & |
22 |
champ(:iim, j, l) = champ(:iim, j, l) * sdd |
+ matmul(matri(:, :, j), champ(:iim, j, l) * sdd) / sdd |
23 |
champ(:iim, j, l) = (champ(:iim, j, l) & |
END forall |
24 |
+ matmul(matri(:, :, j), champ(:iim, j, l))) / sdd |
|
25 |
champ(iim + 1, j, l) = champ(1, j, l) |
champ(iim + 1, :, :) = champ(1, :, :) |
|
END DO |
|
|
END DO |
|
26 |
|
|
27 |
end subroutine filtreg_hemisph |
end subroutine filtreg_hemisph |
28 |
|
|