pkg/longstep/longstep_average_3d.F

C $Header$
C $Name$

#include "LONGSTEP_OPTIONS.h"

CBOP
C     !ROUTINE: LONGSTEP_AVERAGE_3D
C     !INTERFACE:
      SUBROUTINE LONGSTEP_AVERAGE_3D(
     I      lsCnt,
     U      lsFld,
     I      myNr, myThid)

C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE LONGSTEP_AVERAGE_3D
C     | o compute longstep average for one field (incl. overlap)
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     == Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C      lsCnt :: Counter to update
C      lsFld :: average to update
C       myNr :: number of levels to update
C     myThid :: my Thread Id number
      _RL lsFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
      INTEGER lsCnt
      INTEGER myNr
      INTEGER myThid

#ifdef ALLOW_LONGSTEP
C     !LOCAL VARIABLES:
C     == Local variables ==
      INTEGER i,j,k,bi,bj
      _RL factor
CEOP

      factor = FLOAT(lsCnt)
      IF (lsCnt.NE.0) factor = 1. _d 0 / factor

      DO bj=myByLo(myThid), myByHi(myThid)
       DO bi=myBxLo(myThid), myBxHi(myThid)
        DO k = 1,myNr
         DO j = 1-OLy,sNy+OLy
          DO i = 1-OLx,sNx+OLx
           lsFld(i,j,k,bi,bj) = lsFld(i,j,k,bi,bj)*factor
          ENDDO
         ENDDO
        ENDDO
       ENDDO
      ENDDO
#endif /* ALLOW_LONGSTEP */

      RETURN
      END

1	jahn	1.1	C $Header$
2			C $Name$
3
4			#include "LONGSTEP_OPTIONS.h"
5
6			CBOP
7			C !ROUTINE: LONGSTEP_AVERAGE_3D
8			C !INTERFACE:
9			SUBROUTINE LONGSTEP_AVERAGE_3D(
10			I lsCnt,
11			U lsFld,
12			I myNr, myThid)
13
14			C !DESCRIPTION: \bv
15			C ==========================================================
16			C \| SUBROUTINE LONGSTEP_AVERAGE_3D
17			C \| o compute longstep average for one field (incl. overlap)
18			C ==========================================================
19			C \ev
20
21			C !USES:
22			IMPLICIT NONE
23			C == Global variables ===
24			#include "SIZE.h"
25			#include "EEPARAMS.h"
26			#include "PARAMS.h"
27			#include "GRID.h"
28
29			C !INPUT/OUTPUT PARAMETERS:
30			C == Routine arguments ==
31			C lsCnt :: Counter to update
32			C lsFld :: average to update
33			C myNr :: number of levels to update
34			C myThid :: my Thread Id number
35			_RL lsFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
36			INTEGER lsCnt
37			INTEGER myNr
38			INTEGER myThid
39
40			#ifdef ALLOW_LONGSTEP
41			C !LOCAL VARIABLES:
42			C == Local variables ==
43			INTEGER i,j,k,bi,bj
44			_RL factor
45			CEOP
46
47			factor = FLOAT(lsCnt)
48			IF (lsCnt.NE.0) factor = 1. _d 0 / factor
49
50			DO bj=myByLo(myThid), myByHi(myThid)
51			DO bi=myBxLo(myThid), myBxHi(myThid)
52			DO k = 1,myNr
53			DO j = 1-OLy,sNy+OLy
54			DO i = 1-OLx,sNx+OLx
55			lsFld(i,j,k,bi,bj) = lsFld(i,j,k,bi,bj)*factor
56			ENDDO
57			ENDDO
58			ENDDO
59			ENDDO
60			ENDDO
61			#endif /* ALLOW_LONGSTEP */
62
63			RETURN
64			END
65