pkg/longstep/longstep_reset_3d.F

4e66ab0b67 Oliv*0001 #include "LONGSTEP_OPTIONS.h"
                0002
                0003 CBOP
                0004 C     !ROUTINE: LONGSTEP_RESET_3D
                0005 C     !INTERFACE:
                0006       SUBROUTINE LONGSTEP_RESET_3D(
                0007      O      lsCnt, lsFld,
                0008      I      myNr, myThid)
                0009
                0010 C     !DESCRIPTION: \bv
                0011 C     *==========================================================*
                0012 C     | SUBROUTINE LONGSTEP_RESET_3D
                0013 C     | o reset longstep average for one field
                0014 C     *==========================================================*
                0015 C     \ev
                0016
                0017 C     !USES:
                0018       IMPLICIT NONE
                0019 C     == Global variables ===
                0020 #include "SIZE.h"
                0021 #include "EEPARAMS.h"
                0022 #include "PARAMS.h"
                0023 #include "GRID.h"
                0024
                0025 C     !INPUT/OUTPUT PARAMETERS:
                0026 C     == Routine arguments ==
                0027 C      lsCnt :: Counter to update
                0028 C      lsFld :: average to update
                0029 C       myNr :: number of levels to update
                0030 C     myThid :: my Thread Id number
af5eb09265 Jean*0031       INTEGER myNr
4e66ab0b67 Oliv*0032       _RL lsFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
7b3defa94d Jean*0033       INTEGER lsCnt(nSx,nSy)
4e66ab0b67 Oliv*0034       INTEGER myThid
                0035
                0036 #ifdef ALLOW_LONGSTEP
                0037 C     !LOCAL VARIABLES:
                0038 C     == Local variables ==
                0039       INTEGER i,j,k,bi,bj
                0040 CEOP
                0041
                0042       DO bj=myByLo(myThid), myByHi(myThid)
                0043        DO bi=myBxLo(myThid), myBxHi(myThid)
7b3defa94d Jean*0044         lsCnt(bi,bj) = 0
4e66ab0b67 Oliv*0045         DO k = 1,myNr
                0046          DO j = 1-OLy,sNy+OLy
                0047           DO i = 1-OLx,sNx+OLx
                0048            lsFld(i,j,k,bi,bj) = 0. _d 0
                0049           ENDDO
                0050          ENDDO
                0051         ENDDO
                0052        ENDDO
                0053       ENDDO
                0054 #endif /* ALLOW_LONGSTEP */
                0055
                0056       RETURN
                0057       END
                0058