pkg/timeave/timeave_cumul_1vfc.F

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