pkg/zonal_filt/zonal_filt_apply_ts.F

42c525bfb4 Alis*0001 #include "ZONAL_FILT_OPTIONS.h"
                0002
e4b65e705c Jean*0003 CBOP 0
                0004 C     !ROUTINE: ZONAL_FILT_APPLY_TS
                0005
                0006 C     !INTERFACE:
                0007       SUBROUTINE ZONAL_FILT_APPLY_TS(
                0008      U                 tFld, sFld,
                0009      I                 myThid )
                0010
                0011 C     !DESCRIPTION:
                0012 C     *==========================================================*
                0013 C     | S/R ZONAL_FILT_APPLY_TS
                0014 C     | o Apply FFT filter to a latitude circle.
                0015 C     *==========================================================*
42c525bfb4 Alis*0016
e4b65e705c Jean*0017 C     !USES:
                0018       IMPLICIT NONE
42c525bfb4 Alis*0019 C     == Global data ==
                0020 #include "SIZE.h"
                0021 #include "EEPARAMS.h"
                0022 #include "PARAMS.h"
                0023 #include "GRID.h"
                0024
e4b65e705c Jean*0025 C     !INPUT/OUTPUT PARAMETERS:
                0026 C     myThid  :: my Thread Id number
42c525bfb4 Alis*0027       _RL tFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0028       _RL sFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0029       INTEGER myThid
e4b65e705c Jean*0030 CEOP
42c525bfb4 Alis*0031
                0032 #ifdef ALLOW_ZONAL_FILT
                0033
e4b65e705c Jean*0034 C     !LOCAL VARIABLES:
                0035       INTEGER bi,bj
42c525bfb4 Alis*0036
                0037       DO bj=myByLo(myThid),myByHi(myThid)
                0038       DO bi=myBxLo(myThid),myBxHi(myThid)
                0039
                0040         IF ( tempStepping ) CALL ZONAL_FILTER(
e4b65e705c Jean*0041      U                       tFld( 1-OLx,1-OLy,1,bi,bj),
                0042      I                       hFacC(1-OLx,1-OLy,1,bi,bj),
                0043      I                       1, sNy, Nr, bi, bj, 1, myThid )
42c525bfb4 Alis*0044         IF ( saltStepping ) CALL ZONAL_FILTER(
e4b65e705c Jean*0045      U                       sFld( 1-OLx,1-OLy,1,bi,bj),
                0046      I                       hFacC(1-OLx,1-OLy,1,bi,bj),
                0047      I                       1, sNy, Nr, bi, bj, 1, myThid )
                0048
42c525bfb4 Alis*0049
                0050       ENDDO
                0051       ENDDO
                0052
                0053 #endif /* ALLOW_ZONAL_FILT */
                0054
                0055       RETURN
                0056       END