/[MITgcm]/MITgcm/pkg/exf/exf_filter_rs.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_filter_rs.F

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

revision 1.3 by dimitri, Sat Dec 28 10:11:11 2002 UTC revision 1.7 by heimbach, Thu May 25 18:32:55 2006 UTC
# Line 1  Line 1 
1  c $Header$  C $Header$
2    
3  #include "EXF_CPPOPTIONS.h"  #include "EXF_OPTIONS.h"
4    
5        subroutine exf_filter_rs(        subroutine exf_filter_rs(
6       I                          arr,       I                          arr,
# Line 15  c Line 15  c
15  c     o Read a flux record for external forcing.  c     o Read a flux record for external forcing.
16  c  c
17  c     started: Ralf.Giering@FastOpt.de 24-Mai-2000  c     started: Ralf.Giering@FastOpt.de 24-Mai-2000
18  c     changed: menemenlis@jpl.nasa.gov 20-Dec-2002  c     mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
19  c  c
20  c     ==================================================================  c     ==================================================================
21  c     SUBROUTINE exf_filter_rs  c     SUBROUTINE exf_filter_rs
# Line 53  c     == end of interface == Line 53  c     == end of interface ==
53        itlo = mybxlo(mythid)        itlo = mybxlo(mythid)
54        ithi = mybxhi(mythid)        ithi = mybxhi(mythid)
55    
56  c     Don't filter with pkg/seaice because of B/C-grid interpolation  c     Do not filter with pkg/seaice because of B/C-grid interpolation
57        IF ( .NOT. useSEAICE ) THEN        IF ( .NOT. useSEAICE ) THEN
58    
59  c     filter forcing field array  c     filter forcing field array
# Line 76  c         Set land points to zero Line 76  c         Set land points to zero
76    
77              do j = 1,sny              do j = 1,sny
78                do i = 1,snx                do i = 1,snx
79                  if ( _hFacC(i,j,1,bi,bj) .eq. 0. ) then                  if ( maskC(i,j,1,bi,bj) .eq. 0. ) then
80                    arr(i,j,bi,bj) = 0. _d 0                    arr(i,j,bi,bj) = 0. _d 0
81                  endif                  endif
82                enddo                enddo
# Line 86  c         Set land points to zero Line 86  c         Set land points to zero
86    
87              do j = 1,sny              do j = 1,sny
88                do i = 1,snx                do i = 1,snx
89                  if ( maskw(i,j,1,bi,bj) .eq. 0. ) then                  if ( maskW(i,j,1,bi,bj) .eq. 0. ) then
90                    arr(i,j,bi,bj) = 0. _d 0                    arr(i,j,bi,bj) = 0. _d 0
91                  endif                  endif
92                enddo                enddo
# Line 96  c         Set land points to zero Line 96  c         Set land points to zero
96    
97              do j = 1,sny              do j = 1,sny
98                do i = 1,snx                do i = 1,snx
99                  if ( masks(i,j,1,bi,bj) .eq. 0. ) then                  if ( maskS(i,j,1,bi,bj) .eq. 0. ) then
100                    arr(i,j,bi,bj) = 0. _d 0                    arr(i,j,bi,bj) = 0. _d 0
101                  endif                  endif
102                enddo                enddo
# Line 108  c         Set land points to zero Line 108  c         Set land points to zero
108        enddo        enddo
109    
110        ENDIF        ENDIF
111    c     END IF ( .NOT. useSEAICE )
112    
113        end        end
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22