/[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.2 by heimbach, Tue Nov 12 20:34:41 2002 UTC revision 1.8 by jmc, Mon Apr 16 23:27:21 2007 UTC
# Line 1  Line 1 
1  c $Header$  C $Header$
2    C $Name$
3    
4  #include "EXF_CPPOPTIONS.h"  #include "EXF_OPTIONS.h"
5    
6        subroutine exf_filter_rs(        subroutine exf_filter_rs(
7       I                          arr,       I                          arr,
# Line 15  c Line 16  c
16  c     o Read a flux record for external forcing.  c     o Read a flux record for external forcing.
17  c  c
18  c     started: Ralf.Giering@FastOpt.de 24-Mai-2000  c     started: Ralf.Giering@FastOpt.de 24-Mai-2000
19    c     mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
20  c  c
21  c     ==================================================================  c     ==================================================================
22  c     SUBROUTINE exf_filter_rs  c     SUBROUTINE exf_filter_rs
# Line 27  c     == global variables == Line 29  c     == global variables ==
29  #include "EEPARAMS.h"  #include "EEPARAMS.h"
30  #include "SIZE.h"  #include "SIZE.h"
31  #include "GRID.h"  #include "GRID.h"
32  #include "exf_constants.h"  #include "PARAMS.h"
33  #include "exf.h"  #include "EXF_CONSTANTS.h"
34  #include "exf_param.h"  #include "EXF.h"
35    #include "EXF_PARAM.h"
36    
37  c     == routine arguments ==  c     == routine arguments ==
38    
# Line 51  c     == end of interface == Line 54  c     == end of interface ==
54        itlo = mybxlo(mythid)        itlo = mybxlo(mythid)
55        ithi = mybxhi(mythid)        ithi = mybxhi(mythid)
56    
57    c     Do not filter with pkg/seaice because of B/C-grid interpolation
58          IF ( .NOT. useSEAICE ) THEN
59    
60  c     filter forcing field array  c     filter forcing field array
61        do bj = jtlo,jthi        do bj = jtlo,jthi
62          do bi = itlo,ithi          do bi = itlo,ithi
# Line 71  c         Set land points to zero Line 77  c         Set land points to zero
77    
78              do j = 1,sny              do j = 1,sny
79                do i = 1,snx                do i = 1,snx
80                  if ( _hFacC(i,j,1,bi,bj) .eq. 0. ) then                  if ( maskC(i,j,1,bi,bj) .eq. 0. ) then
81                    arr(i,j,bi,bj) = 0. _d 0                    arr(i,j,bi,bj) = 0. _d 0
82                  endif                  endif
83                enddo                enddo
# Line 81  c         Set land points to zero Line 87  c         Set land points to zero
87    
88              do j = 1,sny              do j = 1,sny
89                do i = 1,snx                do i = 1,snx
90                  if ( maskw(i,j,1,bi,bj) .eq. 0. ) then                  if ( maskW(i,j,1,bi,bj) .eq. 0. ) then
91                    arr(i,j,bi,bj) = 0. _d 0                    arr(i,j,bi,bj) = 0. _d 0
92                  endif                  endif
93                enddo                enddo
# Line 91  c         Set land points to zero Line 97  c         Set land points to zero
97    
98              do j = 1,sny              do j = 1,sny
99                do i = 1,snx                do i = 1,snx
100                  if ( masks(i,j,1,bi,bj) .eq. 0. ) then                  if ( maskS(i,j,1,bi,bj) .eq. 0. ) then
101                    arr(i,j,bi,bj) = 0. _d 0                    arr(i,j,bi,bj) = 0. _d 0
102                  endif                  endif
103                enddo                enddo
# Line 102  c         Set land points to zero Line 108  c         Set land points to zero
108          enddo          enddo
109        enddo        enddo
110    
111          ENDIF
112    c     END IF ( .NOT. useSEAICE )
113    
114        end        end
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22