trunk/filtrez/filtreg_hemisph.f90

module filtreg_hemisph_m

  implicit none

contains

  subroutine filtreg_hemisph(champ, sdd, matri)

    USE dimensions, ONLY: iim

    REAL, intent(inout):: champ(:, :, :) ! (iim + 1, :, :)
    REAL, intent(in):: sdd(:) ! (iim) xprim[uv]^{\pm 1/2}

    real, intent(in):: matri(:, :, :) ! (iim, iim, :) 
    ! filtering matrix, last dimension is latitude

    ! Local:
    integer l, j

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

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

  end subroutine filtreg_hemisph

end module filtreg_hemisph_m
1	module filtreg_hemisph_m
2
3	implicit none
4
5	contains
6
7	subroutine filtreg_hemisph(champ, sdd, matri)
8
9	USE dimensions, ONLY: iim
10
11	REAL, intent(inout):: champ(:, :, :) ! (iim + 1, :, :)
12	REAL, intent(in):: sdd(:) ! (iim) xprim[uv]^{\pm 1/2}
13
14	real, intent(in):: matri(:, :, :) ! (iim, iim, :)
15	! filtering matrix, last dimension is latitude
16
17	! Local:
18	integer l, j
19
20	!-----------------------------------------------------------------
21
22	forall (j = 1:size(champ, 2), l = 1:size(champ, 3)) &
23	champ(:iim, j, l) = champ(:iim, j, l) &
24	+ matmul(matri(:, :, j), champ(:iim, j, l) * sdd) / sdd
25	champ(iim + 1, :, :) = champ(1, :, :)
26
27	end subroutine filtreg_hemisph
28
29	end module filtreg_hemisph_m