pkg/zonal_filt/zonal_filt_postsmooth.F

C $Header: /u/gcmpack/models/MITgcmUV/pkg/zonal_filt/Attic/zonal_filt_postsmooth.F,v 1.1.2.1 2001/01/24 16:56:07 adcroft Exp $

#include "ZONAL_FILT_OPTIONS.h"

      SUBROUTINE ZONAL_FILT_POSTSMOOTH( 
     I                           holeMask, 
     U                           field, 
     O                           avgField,
     I                           lField,
     I                           myThid )
C     /==========================================================\
C     | S/R ZONAL_FILT_POSTSMOOTH                                |
C     | o After FFT filtering processing.                        |
C     |==========================================================|
C     | Restore the mean value that was removed prior to FFT.    |
C     \==========================================================/
      IMPLICIT NONE

C     == Global data ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"

C     == Routine arguments ==
C     holeMask - Array with 0 for holes and != 0 for valid data.
C     lField   - Length of field to smooth (assumed periodic)
C     field    - Field smoothed.
C     avgField - Mean value of the field.
C     myThid   - Thread number of this instance of FFT_PRESMOOTH_IN_X
      INTEGER lField
      Real*8  holeMask(lField)
      Real*8  field(lField)
      Real*8  avgField
      INTEGER myThid

#ifdef ALLOW_ZONAL_FILT

C     == Local variables ====
C     I         - Loop counter
C 
      INTEGER I

C
C     Add fields mean value back in.
      DO I=1,lField
       IF ( holeMask(I) .NE. 0. ) THEN 
        field(I) = field(I) + avgField
       ELSE
        field(I) = 0.
       ENDIF
      ENDDO

#endif /* ALLOW_ZONAL_FILT */

      RETURN
      END
1	adcroft	1.2	C $Header: /u/gcmpack/models/MITgcmUV/pkg/zonal_filt/Attic/zonal_filt_postsmooth.F,v 1.1.2.1 2001/01/24 16:56:07 adcroft Exp $
2
3			#include "ZONAL_FILT_OPTIONS.h"
4
5			SUBROUTINE ZONAL_FILT_POSTSMOOTH(
6			I holeMask,
7			U field,
8			O avgField,
9			I lField,
10			I myThid )
11			C /==========================================================\
12			C \| S/R ZONAL_FILT_POSTSMOOTH \|
13			C \| o After FFT filtering processing. \|
14			C \|==========================================================\|
15			C \| Restore the mean value that was removed prior to FFT. \|
16			C \==========================================================/
17			IMPLICIT NONE
18
19			C == Global data ==
20			#include "SIZE.h"
21			#include "EEPARAMS.h"
22			#include "PARAMS.h"
23			#include "GRID.h"
24
25			C == Routine arguments ==
26			C holeMask - Array with 0 for holes and != 0 for valid data.
27			C lField - Length of field to smooth (assumed periodic)
28			C field - Field smoothed.
29			C avgField - Mean value of the field.
30			C myThid - Thread number of this instance of FFT_PRESMOOTH_IN_X
31			INTEGER lField
32			Real*8 holeMask(lField)
33			Real*8 field(lField)
34			Real*8 avgField
35			INTEGER myThid
36
37			#ifdef ALLOW_ZONAL_FILT
38
39			C == Local variables ====
40			C I - Loop counter
41			C
42			INTEGER I
43
44			C
45			C Add fields mean value back in.
46			DO I=1,lField
47			IF ( holeMask(I) .NE. 0. ) THEN
48			field(I) = field(I) + avgField
49			ELSE
50			field(I) = 0.
51			ENDIF
52			ENDDO
53
54			#endif /* ALLOW_ZONAL_FILT */
55
56			RETURN
57			END