Sources/filtrez/filtreg_hemisph.f

module filtreg_hemisph_m

  implicit none

contains

  subroutine filtreg_hemisph(champ, sdd, matri)

    USE dimens_m, ONLY: iim

    REAL, intent(inout):: champ(:, :, :) ! (iim + 1, :, :)
    REAL, intent(in):: sdd(:) ! (iim)
    real, intent(in), dimension(:, :, :):: matri ! (iim, iim, :)

    ! Local:
    integer l, j

    !-----------------------------------------------------------------

    DO l = 1, size(champ, 3)
       DO j = 1, size(champ, 2)
          champ(:iim, j, l) = champ(:iim, j, l) * sdd
          champ(:iim, j, l) = (champ(:iim, j, l) &
               + matmul(matri(:, :, j), champ(:iim, j, l))) / sdd
          champ(iim + 1, j, l) = champ(1, j, l)
       END DO
    END DO

  end subroutine filtreg_hemisph

end module filtreg_hemisph_m
1	guez	136	module filtreg_hemisph_m
2
3			implicit none
4
5			contains
6
7	guez	141	subroutine filtreg_hemisph(champ, sdd, matri)
8	guez	136
9			USE dimens_m, ONLY: iim
10
11			REAL, intent(inout):: champ(:, :, :) ! (iim + 1, :, :)
12	guez	141	REAL, intent(in):: sdd(:) ! (iim)
13	guez	136	real, intent(in), dimension(:, :, :):: matri ! (iim, iim, :)
14
15			! Local:
16			integer l, j
17
18			!-----------------------------------------------------------------
19
20			DO l = 1, size(champ, 3)
21			DO j = 1, size(champ, 2)
22	guez	141	champ(:iim, j, l) = champ(:iim, j, l) * sdd
23	guez	136	champ(:iim, j, l) = (champ(:iim, j, l) &
24	guez	141	+ matmul(matri(:, :, j), champ(:iim, j, l))) / sdd
25	guez	136	champ(iim + 1, j, l) = champ(1, j, l)
26			END DO
27			END DO
28
29			end subroutine filtreg_hemisph
30
31			end module filtreg_hemisph_m