/[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

trunk/Sources/filtrez/filtreg_v.f revision 140 by guez, Fri Jun 5 18:58:06 2015 UTC trunk/filtrez/filtreg_v.f revision 265 by guez, Tue Mar 20 09:35:59 2018 UTC
# Line 12  contains Line 12  contains
12      ! Matrix filter on longitudes. Matrices have already been      ! Matrix filter on longitudes. Matrices have already been
13      ! computed. On v grid, direct filter only.      ! computed. On v grid, direct filter only.
14    
15      USE dimens_m, ONLY: iim, jjm      USE dimensions, 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      USE inifgn_m, ONLY: sddu, unsddu
18      use inifilr_m, only: jfiltnv, jfiltsv, matricevn, matricevs      use inifilr_m, only: jfiltnv, jfiltsv, matricevn, matricevs
# Line 25  contains Line 25  contains
25      ! false means the field is weighted by the area of the mesh      ! 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    
32      call assert(size(champ, 1) == iim + 1, "filtreg_v iim + 1")      call assert(size(champ, 1) == iim + 1, "filtreg_v iim + 1")
33      call assert(size(champ, 2) == jjm, "filtreg_v jjm")      call assert(size(champ, 2) == jjm, "filtreg_v jjm")
34    
35      IF (intensive) THEN      sdd = merge(sddu, unsddu, intensive)
        sdd1 = sddu  
        sdd2 = unsddu  
     ELSE  
        sdd1 = unsddu  
        sdd2 = sddu  
     END IF  
36    
37      call filtreg_hemisph(champ(:, :jfiltnv, :), sdd1, sdd2, matricevn)      call filtreg_hemisph(champ(:, :jfiltnv, :), sdd, matricevn)
38      call filtreg_hemisph(champ(:, jfiltsv:jjm, :), sdd1, sdd2, matricevs)      call filtreg_hemisph(champ(:, jfiltsv:jjm, :), sdd, matricevs)
39    
40    END SUBROUTINE filtreg_v    END SUBROUTINE filtreg_v
41    
Legend:
Removed from v.140  
changed lines
  Added in v.265
  ViewVC Help
Powered by ViewVC 1.1.21