pkg/zonal_filt/zonal_filt_apply.F

C $Header: /u/gcmpack/models/MITgcmUV/pkg/zonal_filt/zonal_filt_apply.F,v 1.2 2001/02/02 21:36:30 adcroft Exp $
C $Name: $

#include "ZONAL_FILT_OPTIONS.h"

      SUBROUTINE ZONAL_FILT_APPLY( 
     U           uFld, vFld, tFld, sFld,
     I           myThid )
C     /==========================================================\
C     | S/R ZONAL_FILT_APPLY                                     |
C     | o Apply FFT filter to a latitude circle.                 |
C     \==========================================================/
      IMPLICIT NONE

C     == Global data ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"

C     == Routine arguments ==
      _RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL tFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL sFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      INTEGER myThid

#ifdef ALLOW_ZONAL_FILT

C     == Local data ==
      INTEGER K,bi,bj

      DO bj=myByLo(myThid),myByHi(myThid)
      DO bi=myBxLo(myThid),myBxHi(myThid)

      DO k=1,Nr
        CALL ZONAL_FILTER(
     &    uFld, hFacW, 1-1, sNy+1, k, k, bi, bj, 1, myThid)
        CALL ZONAL_FILTER(
     &    vFld, hFacS, 1-1, sNy+1, k, k, bi, bj, 2, myThid)
        IF ( tempStepping ) CALL ZONAL_FILTER(
     &    tFld, hFacC, 1, sNy, k, k, bi, bj, 1, myThid)
        IF ( saltStepping ) CALL ZONAL_FILTER(
     &    sFld, hFacC, 1, sNy, k, k, bi, bj, 1, myThid)
      ENDDO

      ENDDO
      ENDDO

#endif /* ALLOW_ZONAL_FILT */

      RETURN
      END
1	C $Header: /u/gcmpack/models/MITgcmUV/pkg/zonal_filt/zonal_filt_apply.F,v 1.2 2001/02/02 21:36:30 adcroft Exp $
2	C $Name: $
3
4	#include "ZONAL_FILT_OPTIONS.h"
5
6	SUBROUTINE ZONAL_FILT_APPLY(
7	U uFld, vFld, tFld, sFld,
8	I myThid )
9	C /==========================================================\
10	C \| S/R ZONAL_FILT_APPLY \|
11	C \| o Apply FFT filter to a latitude circle. \|
12	C \==========================================================/
13	IMPLICIT NONE
14
15	C == Global data ==
16	#include "SIZE.h"
17	#include "EEPARAMS.h"
18	#include "PARAMS.h"
19	#include "GRID.h"
20
21	C == Routine arguments ==
22	_RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
23	_RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
24	_RL tFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
25	_RL sFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
26	INTEGER myThid
27
28	#ifdef ALLOW_ZONAL_FILT
29
30	C == Local data ==
31	INTEGER K,bi,bj
32
33	DO bj=myByLo(myThid),myByHi(myThid)
34	DO bi=myBxLo(myThid),myBxHi(myThid)
35
36	DO k=1,Nr
37	CALL ZONAL_FILTER(
38	& uFld, hFacW, 1-1, sNy+1, k, k, bi, bj, 1, myThid)
39	CALL ZONAL_FILTER(
40	& vFld, hFacS, 1-1, sNy+1, k, k, bi, bj, 2, myThid)
41	IF ( tempStepping ) CALL ZONAL_FILTER(
42	& tFld, hFacC, 1, sNy, k, k, bi, bj, 1, myThid)
43	IF ( saltStepping ) CALL ZONAL_FILTER(
44	& sFld, hFacC, 1, sNy, k, k, bi, bj, 1, myThid)
45	ENDDO
46
47	ENDDO
48	ENDDO
49
50	#endif /* ALLOW_ZONAL_FILT */
51
52	RETURN
53	END