pkg/timeave/timeave_cumulate.F

30e883b651 Jean*0001 #include "CPP_EEOPTIONS.h"
d092edd79c Jean*0002
                0003 CStartofinterface
                0004       SUBROUTINE TIMEAVE_CUMULATE(
                0005      O              fldtave,
873fffa31d Patr*0006      I              fld, Ksize, deltaTloc,
d092edd79c Jean*0007      I              bi, bj, myThid )
                0008 C     /==========================================================\
                0009 C     | SUBROUTINE TIMEAVE_CUMULATE                              |
                0010 C     | o Sum over time a single array                           |
                0011 C     \==========================================================/
                0012       IMPLICIT NONE
                0013
                0014 C     == Global variables ===
                0015 #include "SIZE.h"
                0016 #include "EEPARAMS.h"
                0017
                0018 C     == Routine arguments ==
                0019 C     myThid - Thread number for this instance of the routine.
                0020 C     fldtave - time averaged Field
                0021 C     fld  - Input Field
                0022 C     Ksize - 3rd dimension of local arrays (Input and Output fields)
                0023       INTEGER Ksize
                0024       _RL fld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
                0025       _RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
873fffa31d Patr*0026       _RL deltaTloc
d092edd79c Jean*0027       INTEGER bi, bj, myThid
                0028
                0029 CEndofinterface
                0030
                0031 C     == Local variables ==
                0032 C     i,j,k,bi,bj  - Loop counters
                0033       INTEGER i, j, k
                0034
                0035 C     DO bj = myByLo(myThid), myByHi(myThid)
                0036 C      DO bi = myBxLo(myThid), myBxHi(myThid)
                0037         DO k=1,Ksize
                0038          DO j=1,sNy
                0039           DO i=1,sNx
                0040              fldtave(i,j,k,bi,bj)=
873fffa31d Patr*0041      &       fldtave(i,j,k,bi,bj)+fld(i,j,k,bi,bj)*deltaTloc
d092edd79c Jean*0042            ENDDO
                0043           ENDDO
                0044          ENDDO
                0045 C       ENDDO
                0046 C      ENDDO
                0047
                0048       RETURN
                0049       END