pkg/timeave/timeave_cumul_1vfc.F

C $Header: /u/gcmpack/MITgcm/pkg/timeave/timeave_cumul_1vfc.F,v 1.3 2003/12/05 02:27:59 jmc Exp $
C $Name:  $
#include "CPP_EEOPTIONS.h"

CStartofinterface
      SUBROUTINE TIMEAVE_CUMUL_1VFC( 
     O   fldtave,
     I   fld1, fld2, Ksize, deltaTloc,
     I   bi, bj, myThid )
C     /==========================================================\
C     | SUBROUTINE TIMEAVE_CUMUL_1VFC                            |
C     | o Sum over time a product of 1 field  (RL array)         |
C     |                          and 1 factor (RS 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     fld1  - Input Field
C     fld2  - Input Factor
C     Ksize - 3rd dimension of local arrays (Input and Output fields)
      INTEGER Ksize, dir
      _RL fld1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
      _RS fld2(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*fld1(i,j,k,bi,bj)*fld2(i,j,k,bi,bj)
           ENDDO
          ENDDO
         ENDDO
C       ENDDO
C      ENDDO

      RETURN
      END
1	heimbach	1.4	C $Header: /u/gcmpack/MITgcm/pkg/timeave/timeave_cumul_1vfc.F,v 1.3 2003/12/05 02:27:59 jmc Exp $
2	jmc	1.1	C $Name: $
3	jmc	1.3	#include "CPP_EEOPTIONS.h"
4	jmc	1.1
5			CStartofinterface
6			SUBROUTINE TIMEAVE_CUMUL_1VFC(
7			O fldtave,
8	heimbach	1.4	I fld1, fld2, Ksize, deltaTloc,
9	jmc	1.1	I bi, bj, myThid )
10			C /==========================================================\
11			C \| SUBROUTINE TIMEAVE_CUMUL_1VFC \|
12			C \| o Sum over time a product of 1 field (RL array) \|
13			C \| and 1 factor (RS array) \|
14			C \==========================================================/
15			IMPLICIT NONE
16
17			C == Global variables ===
18			#include "SIZE.h"
19			#include "EEPARAMS.h"
20
21			C == Routine arguments ==
22			C myThid - Thread number for this instance of the routine.
23			C fldtave - time averaged Field
24			C fld1 - Input Field
25			C fld2 - Input Factor
26			C Ksize - 3rd dimension of local arrays (Input and Output fields)
27			INTEGER Ksize, dir
28			_RL fld1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
29			_RS fld2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
30			_RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
31	heimbach	1.4	_RL deltaTloc
32	jmc	1.1	INTEGER bi, bj, myThid
33
34			CEndofinterface
35
36			C == Local variables ==
37			C i,j,k,bi,bj - Loop counters
38			INTEGER i, j, k
39
40			C DO bj = myByLo(myThid), myByHi(myThid)
41			C DO bi = myBxLo(myThid), myBxHi(myThid)
42			DO k=1,Ksize
43			DO j=1,sNy
44			DO i=1,sNx
45			fldtave(i,j,k,bi,bj)= fldtave(i,j,k,bi,bj)
46	heimbach	1.4	& + deltaTlocfld1(i,j,k,bi,bj)fld2(i,j,k,bi,bj)
47	jmc	1.1	ENDDO
48			ENDDO
49			ENDDO
50			C ENDDO
51			C ENDDO
52
53			RETURN
54			END