pkg/ptracers/ptracers_zonal_filt_apply.F

7744ec9132 Jean*0001 #include "PTRACERS_OPTIONS.h"
                0002
                0003 CBOP 0
                0004 C     !ROUTINE: PTRACERS_ZONAL_FILT_APPLY
                0005
                0006 C     !INTERFACE:
                0007       SUBROUTINE PTRACERS_ZONAL_FILT_APPLY(
                0008      I                                      myThid )
                0009
                0010 C     !DESCRIPTION:
                0011 C     *==========================================================*
                0012 C     | S/R PTRACERS_ZONAL_FILT_APPLY
                0013 C     | o Apply FFT filter to all passive tracers
                0014 C     *==========================================================*
                0015
                0016 C     !USES:
                0017       IMPLICIT NONE
                0018 C     == Global data ==
                0019 #include "SIZE.h"
                0020 #include "EEPARAMS.h"
                0021 #include "PARAMS.h"
                0022 #include "GRID.h"
                0023 #include "PTRACERS_SIZE.h"
                0024 #include "PTRACERS_PARAMS.h"
d3a355ef62 Jean*0025 #include "PTRACERS_START.h"
7744ec9132 Jean*0026 #include "PTRACERS_FIELDS.h"
                0027
                0028 C     !INPUT/OUTPUT PARAMETERS:
                0029 C     myThid  :: my Thread Id number
                0030       INTEGER myThid
                0031 CEOP
                0032
                0033 #ifdef ALLOW_ZONAL_FILT
                0034
                0035 C     !LOCAL VARIABLES:
                0036       INTEGER iTr
                0037       INTEGER bi,bj
                0038
                0039       DO iTr = 1,PTRACERS_numInUse
d3a355ef62 Jean*0040        IF ( PTRACERS_StepFwd(iTr) ) THEN
7744ec9132 Jean*0041
d3a355ef62 Jean*0042         DO bj=myByLo(myThid),myByHi(myThid)
                0043          DO bi=myBxLo(myThid),myBxHi(myThid)
7744ec9132 Jean*0044           CALL ZONAL_FILTER(
d3a355ef62 Jean*0045      U                       pTracer(1-OLx,1-OLy,1,bi,bj,iTr),
7744ec9132 Jean*0046      I                       hFacC(1-OLx,1-OLy,1,bi,bj),
                0047      I                       1, sNy, Nr, bi, bj, 1, myThid )
d3a355ef62 Jean*0048          ENDDO
7744ec9132 Jean*0049         ENDDO
d3a355ef62 Jean*0050
                0051        ENDIF
7744ec9132 Jean*0052       ENDDO
                0053
                0054 #endif /* ALLOW_ZONAL_FILT */
                0055
                0056       RETURN
                0057       END