1 |
C $Header: /u/gcmpack/models/MITgcmUV/pkg/zonal_filt/zonal_filt_apply.F,v 1.3 2001/02/04 14:38:51 cnh Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
#include "ZONAL_FILT_OPTIONS.h" |
5 |
|
6 |
SUBROUTINE ZONAL_FILT_APPLY_UV( |
7 |
U uFld, vFld, |
8 |
I myThid ) |
9 |
C /==========================================================\ |
10 |
C | S/R ZONAL_FILT_APPLY_UV | |
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 |
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 |
CALL ZONAL_FILTER( |
36 |
& uFld, hFacW, 1-1, sNy+1, k, k, bi, bj, 1, myThid) |
37 |
CALL ZONAL_FILTER( |
38 |
& vFld, hFacS, 1-1, sNy+1, k, k, bi, bj, 2, myThid) |
39 |
ENDDO |
40 |
|
41 |
ENDDO |
42 |
ENDDO |
43 |
|
44 |
#endif /* ALLOW_ZONAL_FILT */ |
45 |
|
46 |
RETURN |
47 |
END |