pkg/timeave/timeave_cumul_1k.F

C $Header: /u/gcmpack/MITgcm/pkg/timeave/timeave_cumul_1k.F,v 1.4 2005/08/19 22:50:26 heimbach Exp $
C $Name:  $
#include "CPP_EEOPTIONS.h"

CStartofinterface
      SUBROUTINE TIMEAVE_CUMUL_1K( fld, fldtave, deltaTloc, 
     &                             bi, bj, K, myThid )
C     /==========================================================\
C     | SUBROUTINE TIMEAVE_CUMUL_1K                              |
C     | o Sum over time a single array and a single level        |
C     |  (replace "timeaver_1fld_xyz" and use the same argments) |
C     \==========================================================/
      IMPLICIT NONE

C     == Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"

C     == Routine arguments ==
C     myThid - Thread number for this instance of the routine.
C     fld  - Input Field
C     fldtave - time averaged Field
      _RL fld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL deltaTloc
      INTEGER bi, bj, K, myThid

CEndofinterface

C     == Local variables ==
C     i,j,k,bi,bj  - Loop counters
      INTEGER i, j
 
C     DO bj = myByLo(myThid), myByHi(myThid)
C      DO bi = myBxLo(myThid), myBxHi(myThid)
C       DO K=1,Nr
         DO j=1,sNy
          DO i=1,sNx
             fldtave(i,j,K,bi,bj)=
     &       fldtave(i,j,K,bi,bj)+fld(i,j,K,bi,bj)*deltaTloc
           ENDDO
          ENDDO
C        ENDDO
C       ENDDO
C      ENDDO
 
      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/timeave/timeave_cumul_1k.F,v 1.4 2005/08/19 22:50:26 heimbach Exp $
2	C $Name: $
3	#include "CPP_EEOPTIONS.h"
4
5	CStartofinterface
6	SUBROUTINE TIMEAVE_CUMUL_1K( fld, fldtave, deltaTloc,
7	& bi, bj, K, myThid )
8	C /==========================================================\
9	C \| SUBROUTINE TIMEAVE_CUMUL_1K \|
10	C \| o Sum over time a single array and a single level \|
11	C \| (replace "timeaver_1fld_xyz" and use the same argments) \|
12	C \==========================================================/
13	IMPLICIT NONE
14
15	C == Global variables ===
16	#include "SIZE.h"
17	#include "EEPARAMS.h"
18
19	C == Routine arguments ==
20	C myThid - Thread number for this instance of the routine.
21	C fld - Input Field
22	C fldtave - time averaged Field
23	_RL fld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
24	_RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
25	_RL deltaTloc
26	INTEGER bi, bj, K, myThid
27
28	CEndofinterface
29
30	C == Local variables ==
31	C i,j,k,bi,bj - Loop counters
32	INTEGER i, j
33
34	C DO bj = myByLo(myThid), myByHi(myThid)
35	C DO bi = myBxLo(myThid), myBxHi(myThid)
36	C DO K=1,Nr
37	DO j=1,sNy
38	DO i=1,sNx
39	fldtave(i,j,K,bi,bj)=
40	& fldtave(i,j,K,bi,bj)+fld(i,j,K,bi,bj)*deltaTloc
41	ENDDO
42	ENDDO
43	C ENDDO
44	C ENDDO
45	C ENDDO
46
47	RETURN
48	END