pkg/zonal_filt/ZONAL_FILT.h

C $Header: /u/gcmpack/MITgcm/pkg/zonal_filt/ZONAL_FILT.h,v 1.4 2001/12/11 14:47:32 jmc Exp $
C $Name:  $

#ifdef ALLOW_ZONAL_FILT

C-    Package flag and logical parameters :
C     zonal_filt_uvStar  :: filter applied to u*,v* (before SOLVE_FOR_P)
C     zonal_filt_TrStagg :: if using a Stager time-step, filter T,S before
C                           computing PhiHyd ; 
C                           has no effect if syncr. time step is used

      LOGICAL zonal_filt_uvStar, zonal_filt_TrStagg
      COMMON /ZONAL_FILT_PARM_L/
     &        zonal_filt_uvStar, zonal_filt_TrStagg

C-    Zonal Filter integer parameters :
C     zonal_filt_cospow  :: Latitude dependance of the damping function
C                           = ( cos Lat / cos zonal_filt_lat )**cospow 
C     zonal_filt_sinpow  :: zonal mode dependance of the damping function
C                           = 1 / ( sin pi.kx/Nx )**sinpow
C     zonal_filt_mode2dx :: to specify how to treat the 2.dx mode :
C                         = 0 : damped like other modes.
C                         = 1 : removed in regions where Zonal_filt apply
C                         = 2 : removed every where.

      INTEGER zonal_filt_cospow, zonal_filt_sinpow, zonal_filt_mode2dx
      COMMON /ZONAL_FILT_PARM_I/
     &        zonal_filt_cospow, zonal_filt_sinpow, zonal_filt_mode2dx
 
C-    Zonal Filter (real) parameters :
C     zonal_filt_lat :: Low latitude for FFT filtering of latitude
C                          circles
      _RL zonal_filt_lat  
      COMMON /ZONAL_FILT_PARAMS/
     & zonal_filt_lat

C Amplitude factor as function of mode number and latitude (U,T points)
      COMMON /ZONAL_FFT/ ampFactor,ampFactorV
C     _RL ampFactor( Nx, 1-Oly:sNy+Oly, nSx, nSy )
C     _RL ampFactorV( Nx, 1-Oly:sNy+Oly, nSx, nSy )
      _RL ampFactor( 1-Olx:sNx+Olx, 1-Oly:sNy+Oly, nSx, nSy )
      _RL ampFactorV( 1-Olx:sNx+Olx, 1-Oly:sNy+Oly, nSx, nSy )

#endif /* ALLOW_ZONAL_FILT */
1	C $Header: /u/gcmpack/MITgcm/pkg/zonal_filt/ZONAL_FILT.h,v 1.4 2001/12/11 14:47:32 jmc 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 */