pkg/timeave/timeave_cumul_fc.F

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

CStartofinterface
      SUBROUTINE TIMEAVE_CUMUL_FC( 
     O              fldtave, 
     I              fld, Ksize, deltaTloc, 
     I              bi, bj, myThid )
C     /==========================================================\
C     | SUBROUTINE TIMEAVE_CUMUL_FC                              |
C     | o Sum over time a single RS array into a RL array        |
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     fldtave - time averaged Field
C     fld  - Input Field
C     Ksize - 3rd dimension of local arrays (Input and Output fields)
      INTEGER Ksize
      _RS fld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
      _RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
      _RL deltaTloc
      INTEGER bi, bj, myThid

CEndofinterface

C     == Local variables ==
C     i,j,k,bi,bj  - Loop counters
      INTEGER i, j, k
 
C     DO bj = myByLo(myThid), myByHi(myThid)
C      DO bi = myBxLo(myThid), myBxHi(myThid)
        DO k=1,Ksize
         DO j=1,sNy
          DO i=1,sNx
           fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)
     &                       + deltaTloc*fld(i,j,k,bi,bj)
           ENDDO
          ENDDO
         ENDDO
C       ENDDO
C      ENDDO
 
      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/timeave/timeave_cumul_fc.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_FC(
7	O fldtave,
8	I fld, Ksize, deltaTloc,
9	I bi, bj, myThid )
10	C /==========================================================\
11	C \| SUBROUTINE TIMEAVE_CUMUL_FC \|
12	C \| o Sum over time a single RS array into a RL array \|
13	C \==========================================================/
14	IMPLICIT NONE
15
16	C == Global variables ===
17	#include "SIZE.h"
18	#include "EEPARAMS.h"
19
20	C == Routine arguments ==
21	C myThid - Thread number for this instance of the routine.
22	C fldtave - time averaged Field
23	C fld - Input Field
24	C Ksize - 3rd dimension of local arrays (Input and Output fields)
25	INTEGER Ksize
26	_RS fld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
27	_RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
28	_RL deltaTloc
29	INTEGER bi, bj, myThid
30
31	CEndofinterface
32
33	C == Local variables ==
34	C i,j,k,bi,bj - Loop counters
35	INTEGER i, j, k
36
37	C DO bj = myByLo(myThid), myByHi(myThid)
38	C DO bi = myBxLo(myThid), myBxHi(myThid)
39	DO k=1,Ksize
40	DO j=1,sNy
41	DO i=1,sNx
42	fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)
43	& + deltaTloc*fld(i,j,k,bi,bj)
44	ENDDO
45	ENDDO
46	ENDDO
47	C ENDDO
48	C ENDDO
49
50	RETURN
51	END