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

revision 141 by guez, Fri Jun 5 19:21:08 2015 UTC revision 169 by guez, Mon Sep 14 17:13:16 2015 UTC
# Line 9  contains Line 9  contains
9      USE dimens_m, ONLY: iim      USE dimens_m, ONLY: iim
10    
11      REAL, intent(inout):: champ(:, :, :) ! (iim + 1, :, :)      REAL, intent(inout):: champ(:, :, :) ! (iim + 1, :, :)
12      REAL, intent(in):: sdd(:) ! (iim)      REAL, intent(in):: sdd(:) ! (iim) xprim[uv]^{\pm 1/2}
13      real, intent(in), dimension(:, :, :):: matri ! (iim, iim, :)  
14        real, intent(in), dimension(:, :, :):: matri ! (iim, iim, :)
15        ! filtering matrix, last dimension is latitude
16    
17      ! Local:      ! Local:
18      integer l, j      integer l, j
19    
20      !-----------------------------------------------------------------      !-----------------------------------------------------------------
21    
22      DO l = 1, size(champ, 3)      forall (j = 1:size(champ, 2), l = 1:size(champ, 3))
23         DO j = 1, size(champ, 2)         champ(:iim, j, l) = champ(:iim, j, l) &
24            champ(:iim, j, l) = champ(:iim, j, l) * sdd              + matmul(matri(:, :, j), champ(:iim, j, l) * sdd) / sdd
25            champ(:iim, j, l) = (champ(:iim, j, l) &      END forall
26                 + matmul(matri(:, :, j), champ(:iim, j, l))) / sdd  
27            champ(iim + 1, j, l) = champ(1, j, l)      champ(iim + 1, :, :) = champ(1, :, :)
        END DO  
     END DO  
28    
29    end subroutine filtreg_hemisph    end subroutine filtreg_hemisph
30    
Legend:
Removed from v.141  
changed lines
  Added in v.169
  ViewVC Help
Powered by ViewVC 1.1.21