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

Diff of /trunk/filtrez/filtreg_hemisph.f

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

trunk/Sources/filtrez/filtreg_hemisph.f revision 169 by guez, Mon Sep 14 17:13:16 2015 UTC trunk/filtrez/filtreg_hemisph.f revision 266 by guez, Thu Apr 19 17:54:55 2018 UTC
# Line 6  contains Line 6  contains
6    
7    subroutine filtreg_hemisph(champ, sdd, matri)    subroutine filtreg_hemisph(champ, sdd, matri)
8    
9      USE dimens_m, ONLY: iim      USE dimensions, ONLY: iim
10    
11      REAL, intent(inout):: champ(:, :, :) ! (iim + 1, :, :)      REAL, intent(inout):: champ(:, :, :) ! (iim + 1, :, :)
12      REAL, intent(in):: sdd(:) ! (iim) xprim[uv]^{\pm 1/2}      REAL, intent(in):: sdd(:) ! (iim) xprim[uv]^{\pm 1/2}
13    
14      real, intent(in), dimension(:, :, :):: matri ! (iim, iim, :)      real, intent(in):: matri(:, :, :) ! (iim, iim, :)
15      ! filtering matrix, last dimension is latitude      ! filtering matrix, last dimension is latitude
16    
17      ! Local:      ! Local:
# Line 19  contains Line 19  contains
19    
20      !-----------------------------------------------------------------      !-----------------------------------------------------------------
21    
22      forall (j = 1:size(champ, 2), l = 1:size(champ, 3))      forall (j = 1:size(champ, 2), l = 1:size(champ, 3)) &
23         champ(:iim, j, l) = champ(:iim, j, l) &           champ(:iim, j, l) = champ(:iim, j, l) &
24              + matmul(matri(:, :, j), champ(:iim, j, l) * sdd) / sdd           + matmul(matri(:, :, j), champ(:iim, j, l) * sdd) / sdd
     END forall  
   
25      champ(iim + 1, :, :) = champ(1, :, :)      champ(iim + 1, :, :) = champ(1, :, :)
26    
27    end subroutine filtreg_hemisph    end subroutine filtreg_hemisph
Legend:
Removed from v.169  
changed lines
  Added in v.266
  ViewVC Help
Powered by ViewVC 1.1.21