--- trunk/Sources/filtrez/filtreg_hemisph.f 2015/04/30 18:35:49 136 +++ trunk/Sources/filtrez/filtreg_hemisph.f 2015/06/05 19:21:08 141 @@ -4,12 +4,12 @@ contains - subroutine filtreg_hemisph(champ, sdd1, sdd2, matri) + subroutine filtreg_hemisph(champ, sdd, matri) USE dimens_m, ONLY: iim REAL, intent(inout):: champ(:, :, :) ! (iim + 1, :, :) - REAL, intent(in):: sdd1(:), sdd2(:) ! (iim) + REAL, intent(in):: sdd(:) ! (iim) real, intent(in), dimension(:, :, :):: matri ! (iim, iim, :) ! Local: @@ -19,9 +19,9 @@ DO l = 1, size(champ, 3) DO j = 1, size(champ, 2) - champ(:iim, j, l) = champ(:iim, j, l) * sdd1 + champ(:iim, j, l) = champ(:iim, j, l) * sdd champ(:iim, j, l) = (champ(:iim, j, l) & - + matmul(matri(:, :, j), champ(:iim, j, l))) * sdd2 + + matmul(matri(:, :, j), champ(:iim, j, l))) / sdd champ(iim + 1, j, l) = champ(1, j, l) END DO END DO