1 |
jmc |
1.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 |