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

Diff of /trunk/filtrez/filtreg_v.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 137 by guez, Wed May 6 15:51:03 2015 UTC revision 141 by guez, Fri Jun 5 19:21:08 2015 UTC
# Line 9  contains Line 9  contains
9      ! From filtrez/filtreg.F, version 1.1.1.1, 2004/05/19 12:53:09      ! From filtrez/filtreg.F, version 1.1.1.1, 2004/05/19 12:53:09
10      ! Author: P. Le Van      ! Author: P. Le Van
11    
12      ! Filtre matriciel longitudinal, avec les matrices      ! Matrix filter on longitudes. Matrices have already been
13      ! pr\'ecalcul\'ees pour l'op\'erateur filtre. On v grid, direct      ! computed. On v grid, direct filter only.
     ! only.  
14    
     USE coefils, ONLY: sddu, unsddu  
15      USE dimens_m, ONLY: iim, jjm      USE dimens_m, ONLY: iim, jjm
16      use filtreg_hemisph_m, only: filtreg_hemisph      use filtreg_hemisph_m, only: filtreg_hemisph
17        USE inifgn_m, ONLY: sddu, unsddu
18      use inifilr_m, only: jfiltnv, jfiltsv, matricevn, matricevs      use inifilr_m, only: jfiltnv, jfiltsv, matricevn, matricevs
19      use nr_util, only: assert      use nr_util, only: assert
20    
# Line 23  contains Line 22  contains
22      ! en entr\'ee : champ \`a filtrer, en sortie : champ filtr\'e      ! en entr\'ee : champ \`a filtrer, en sortie : champ filtr\'e
23    
24      logical, intent(in):: intensive      logical, intent(in):: intensive
25      ! champ intensif ou extensif (pond\'er\'e par les aires)      ! false means the field is weighted by the area of the mesh
26    
27      ! Local:      ! Local:
28      REAL sdd1(iim), sdd2(iim)      REAL sdd(iim)
29    
30      !-----------------------------------------------------------      !-----------------------------------------------------------
31    
# Line 34  contains Line 33  contains
33      call assert(size(champ, 2) == jjm, "filtreg_v jjm")      call assert(size(champ, 2) == jjm, "filtreg_v jjm")
34    
35      IF (intensive) THEN      IF (intensive) THEN
36         sdd1 = sddu         sdd = sddu
        sdd2 = unsddu  
37      ELSE      ELSE
38         sdd1 = unsddu         sdd = unsddu
        sdd2 = sddu  
39      END IF      END IF
40    
41      call filtreg_hemisph(champ(:, :jfiltnv, :), sdd1, sdd2, matricevn)      call filtreg_hemisph(champ(:, :jfiltnv, :), sdd, matricevn)
42      call filtreg_hemisph(champ(:, jfiltsv:jjm, :), sdd1, sdd2, matricevs)      call filtreg_hemisph(champ(:, jfiltsv:jjm, :), sdd, matricevs)
43    
44    END SUBROUTINE filtreg_v    END SUBROUTINE filtreg_v
45    
Legend:
Removed from v.137  
changed lines
  Added in v.141
  ViewVC Help
Powered by ViewVC 1.1.21