pkg/zonal_filt/zonal_filt_apply_ts.F

C $Header: /u/gcmpack/MITgcm/pkg/zonal_filt/zonal_filt_apply_ts.F,v 1.1 2001/12/11 14:49:07 jmc Exp $
C $Name: checkpoint44 $

#include "ZONAL_FILT_OPTIONS.h"

      SUBROUTINE ZONAL_FILT_APPLY_TS( 
     U           tFld, sFld,
     I           myThid )
C     /==========================================================\
C     | S/R ZONAL_FILT_APPLY_TS                                  |
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 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
        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/MITgcm/pkg/zonal_filt/zonal_filt_apply_ts.F,v 1.1 2001/12/11 14:49:07 jmc Exp $
2	C $Name: checkpoint44 $
3
4	#include "ZONAL_FILT_OPTIONS.h"
5
6	SUBROUTINE ZONAL_FILT_APPLY_TS(
7	U tFld, sFld,
8	I myThid )
9	C /==========================================================\
10	C \| S/R ZONAL_FILT_APPLY_TS \|
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 tFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
23	_RL sFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
24	INTEGER myThid
25
26	#ifdef ALLOW_ZONAL_FILT
27
28	C == Local data ==
29	INTEGER K,bi,bj
30
31	DO bj=myByLo(myThid),myByHi(myThid)
32	DO bi=myBxLo(myThid),myBxHi(myThid)
33
34	DO k=1,Nr
35	IF ( tempStepping ) CALL ZONAL_FILTER(
36	& tFld, hFacC, 1, sNy, k, k, bi, bj, 1, myThid)
37	IF ( saltStepping ) CALL ZONAL_FILTER(
38	& sFld, hFacC, 1, sNy, k, k, bi, bj, 1, myThid)
39	ENDDO
40
41	ENDDO
42	ENDDO
43
44	#endif /* ALLOW_ZONAL_FILT */
45
46	RETURN
47	END