/[lmdze]/trunk/filtrez/filtreg_scal.f
ViewVC logotype

Annotation of /trunk/filtrez/filtreg_scal.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 137 - (hide annotations)
Wed May 6 15:51:03 2015 UTC (9 years ago) by guez
Original Path: trunk/Sources/filtrez/filtreg_scal.f
File size: 1620 byte(s) 
Removed unused argument missval in ma_fucoll_r[1-3]1, ma_fufill_r[1-3]1.

Split filtreg into two procedures: filtreg_scal and filtreg_v. I did
not like the test on the extent of the argument and there was no
common code between the two cases: jjm and jjm + 1. Also, it is
simpler now to just remove the argument "direct" from filtreg_v instead
of allowing it and then stopping the program if it is false.

Removed the computation of pkf in reanalyse2nat, was not used.

As a consequence of the split of filtreg, had to extract the
computation of pkf out of exner_hyb. This is clearer anyway because we
want to be able to call exner_hyb with any size in the first two
dimensions (as in test_disvert). But at the same time exner_hyb
required particular sizes for the computation of pkf. It was
awkward. The only computation of pkf is now in leapfrog.
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 137 USE coefils, ONLY: sddv, unsddv
15 guez 64 USE dimens_m, ONLY: iim, jjm
16 guez 136 use filtreg_hemisph_m, only: filtreg_hemisph
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 136 REAL sdd1(iim), sdd2(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     sdd1 = sddv
39     sdd2 = unsddv
40     ELSE
41     sdd1 = unsddv
42     sdd2 = sddv
43     END IF
44    
45     if (direct) then
46     call filtreg_hemisph(champ(:, 2:jfiltnu, :), sdd1, sdd2, matriceun)
47     call filtreg_hemisph(champ(:, jfiltsu:jjm, :), sdd1, sdd2, matriceus)
48 guez 136 else
49 guez 137 call filtreg_hemisph(champ(:, 2:jfiltnu, :), sdd1, sdd2, - matrinvn)
50     call filtreg_hemisph(champ(:, jfiltsu:jjm, :), sdd1, sdd2, - matrinvs)
51 guez 136 end if
52 guez 25
53 guez 137 END SUBROUTINE filtreg_scal
54 guez 27
55 guez 137 end module filtreg_scal_m
  ViewVC Help
Powered by ViewVC 1.1.21