trunk/filtrez/filtreg_scal.f

module filtreg_scal_m

  IMPLICIT NONE

contains

  SUBROUTINE filtreg_scal(champ, direct, intensive)

    ! From filtrez/filtreg.F, version 1.1.1.1, 2004/05/19 12:53:09
    ! Author: P. Le Van
    ! Objet : filtre matriciel longitudinal, avec les matrices pr\'ecalcul\'ees
    ! pour l'op\'erateur filtre. 

    USE dimens_m, ONLY: iim, jjm
    use filtreg_hemisph_m, only: filtreg_hemisph
    USE inifgn_m, ONLY: sddv, unsddv
    use inifilr_m, only: jfiltnu, jfiltsu, matriceun, matriceus, matrinvn, &
         matrinvs
    use nr_util, only: assert

    REAL, intent(inout):: champ(:, :, :) ! (iim + 1, jjm + 1, :)
    ! en entr\'ee : champ \`a filtrer, en sortie : champ filtr\'e

    logical, intent(in):: direct ! filtre direct ou inverse 

    logical, intent(in):: intensive
    ! champ intensif ou extensif (pond\'er\'e par les aires)

    ! Local:
    REAL sdd(iim)

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

    call assert(size(champ, 1) == iim + 1, "filtreg_scal iim + 1")
    call assert(size(champ, 2) == jjm + 1, "filtreg_scal jjm + 1")

    IF (intensive) THEN
       sdd = sddv
    ELSE
       sdd = unsddv
    END IF

    if (direct) then
       call filtreg_hemisph(champ(:, 2:jfiltnu, :), sdd, matriceun)
       call filtreg_hemisph(champ(:, jfiltsu:jjm, :), sdd, matriceus)
    else
       call filtreg_hemisph(champ(:, 2:jfiltnu, :), sdd, - matrinvn)
       call filtreg_hemisph(champ(:, jfiltsu:jjm, :), sdd, - matrinvs)
    end if

  END SUBROUTINE filtreg_scal

end module filtreg_scal_m
1	guez	137	module filtreg_scal_m
2	guez	3
3	guez	27	IMPLICIT NONE
4	guez	3
5	guez	27	contains
6	guez	3
7	guez	137	SUBROUTINE filtreg_scal(champ, direct, intensive)
8	guez	3
9	guez	27	! From filtrez/filtreg.F, version 1.1.1.1, 2004/05/19 12:53:09
10			! Author: P. Le Van
11	guez	133	! Objet : filtre matriciel longitudinal, avec les matrices pr\'ecalcul\'ees
12			! pour l'op\'erateur filtre.
13	guez	3
14	guez	64	USE dimens_m, ONLY: iim, jjm
15	guez	136	use filtreg_hemisph_m, only: filtreg_hemisph
16	guez	140	USE inifgn_m, ONLY: sddv, unsddv
17	guez	137	use inifilr_m, only: jfiltnu, jfiltsu, matriceun, matriceus, matrinvn, &
18			matrinvs
19	guez	107	use nr_util, only: assert
20	guez	3
21	guez	137	REAL, intent(inout):: champ(:, :, :) ! (iim + 1, jjm + 1, :)
22	guez	133	! en entr\'ee : champ \`a filtrer, en sortie : champ filtr\'e
23	guez	3
24	guez	107	logical, intent(in):: direct ! filtre direct ou inverse
25	guez	3
26	guez	107	logical, intent(in):: intensive
27	guez	133	! champ intensif ou extensif (pond\'er\'e par les aires)
28	guez	3
29	guez	107	! Local:
30	guez	141	REAL sdd(iim)
31	guez	25
32	guez	27	!-----------------------------------------------------------
33	guez	25
34	guez	137	call assert(size(champ, 1) == iim + 1, "filtreg_scal iim + 1")
35			call assert(size(champ, 2) == jjm + 1, "filtreg_scal jjm + 1")
36	guez	25
37	guez	137	IF (intensive) THEN
38	guez	141	sdd = sddv
39	guez	137	ELSE
40	guez	141	sdd = unsddv
41	guez	137	END IF
42
43			if (direct) then
44	guez	141	call filtreg_hemisph(champ(:, 2:jfiltnu, :), sdd, matriceun)
45			call filtreg_hemisph(champ(:, jfiltsu:jjm, :), sdd, matriceus)
46	guez	136	else
47	guez	141	call filtreg_hemisph(champ(:, 2:jfiltnu, :), sdd, - matrinvn)
48			call filtreg_hemisph(champ(:, jfiltsu:jjm, :), sdd, - matrinvs)
49	guez	136	end if
50	guez	25
51	guez	137	END SUBROUTINE filtreg_scal
52	guez	27
53	guez	137	end module filtreg_scal_m