1 |
C $Header: /u/gcmpack/models/MITgcmUV/pkg/zonal_filt/ZONAL_FILT.h,v 1.3 2001/02/04 14:38:50 cnh Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
#ifdef ALLOW_ZONAL_FILT |
5 |
|
6 |
C- Package flag and logical parameters : |
7 |
C zonal_filt_uvStar :: filter applied to u*,v* (before SOLVE_FOR_P) |
8 |
C zonal_filt_TrStagg :: if using a Stager time-step, filter T,S before |
9 |
C computing PhiHyd ; |
10 |
C has no effect if syncr. time step is used |
11 |
|
12 |
LOGICAL zonal_filt_uvStar, zonal_filt_TrStagg |
13 |
COMMON /ZONAL_FILT_PARM_L/ |
14 |
& zonal_filt_uvStar, zonal_filt_TrStagg |
15 |
|
16 |
C- Zonal Filter integer parameters : |
17 |
C zonal_filt_cospow :: Latitude dependance of the damping function |
18 |
C = ( cos Lat / cos zonal_filt_lat )**cospow |
19 |
C zonal_filt_sinpow :: zonal mode dependance of the damping function |
20 |
C = 1 / ( sin pi.kx/Nx )**sinpow |
21 |
C zonal_filt_mode2dx :: to specify how to treat the 2.dx mode : |
22 |
C = 0 : damped like other modes. |
23 |
C = 1 : removed in regions where Zonal_filt apply |
24 |
C = 2 : removed every where. |
25 |
|
26 |
INTEGER zonal_filt_cospow, zonal_filt_sinpow, zonal_filt_mode2dx |
27 |
COMMON /ZONAL_FILT_PARM_I/ |
28 |
& zonal_filt_cospow, zonal_filt_sinpow, zonal_filt_mode2dx |
29 |
|
30 |
C- Zonal Filter (real) parameters : |
31 |
C zonal_filt_lat :: Low latitude for FFT filtering of latitude |
32 |
C circles |
33 |
_RL zonal_filt_lat |
34 |
COMMON /ZONAL_FILT_PARAMS/ |
35 |
& zonal_filt_lat |
36 |
|
37 |
C Amplitude factor as function of mode number and latitude (U,T points) |
38 |
COMMON /ZONAL_FFT/ ampFactor,ampFactorV |
39 |
C _RL ampFactor( Nx, 1-Oly:sNy+Oly, nSx, nSy ) |
40 |
C _RL ampFactorV( Nx, 1-Oly:sNy+Oly, nSx, nSy ) |
41 |
_RL ampFactor( 1-Olx:sNx+Olx, 1-Oly:sNy+Oly, nSx, nSy ) |
42 |
_RL ampFactorV( 1-Olx:sNx+Olx, 1-Oly:sNy+Oly, nSx, nSy ) |
43 |
|
44 |
#endif /* ALLOW_ZONAL_FILT */ |