pkg/timeave/timeave_statvars.F

C $Header: /u/gcmpack/MITgcm/pkg/timeave/timeave_statvars.F,v 1.19 2005/02/20 11:46:25 dimitri Exp $
C $Name:  $
#include "TIMEAVE_OPTIONS.h"

      SUBROUTINE TIMEAVE_STATVARS(
     I     myTime, myIter, bi, bj, myThid)
C     /==========================================================\
C     | SUBROUTINE TIMEAVE_STATVARS                              |
C     | o Time averaging routine for eta, U, V, W, T, S, UT, VT  |
C     |   in model main time-stepping                            |
C     \==========================================================/
      IMPLICIT NONE

C     == Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DYNVARS.h"
#include "GRID.h"
#include "TIMEAVE_STATV.h"

      LOGICAL  DIFF_BASE_MULTIPLE
      EXTERNAL DIFF_BASE_MULTIPLE

C     == Routine arguments ==
C     myThid - Thread number for this instance of the routine.
C     myIter - Iteration number
C     myTime - Current time of simulation ( s )
      INTEGER myThid
      INTEGER myIter, bi, bj
      _RL     myTime

#ifdef ALLOW_TIMEAVE

C     == Local variables ==
      INTEGER I, J, K
      _RL DDTT
      _RL tempArray (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      INTEGER thisdate(4), prevdate(4)
 
C-    Initialize fields for the first call ever
      IF ( myIter .EQ. nIter0 ) THEN
        CALL TIMEAVE_RESET(uFluxtave, 1,  bi, bj, myThid)
        CALL TIMEAVE_RESET(vFluxtave, 1,  bi, bj, myThid)
        CALL TIMEAVE_RESET(tFluxtave, 1,  bi, bj, myThid)
        CALL TIMEAVE_RESET(sFluxtave, 1,  bi, bj, myThid)
        CALL TIMEAVE_RESET(etatave,   1,  bi, bj, myThid)
        CALL TIMEAVE_RESET(thetatave, Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(salttave,  Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(uVeltave,  Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(vVeltave,  Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(wVeltave,  Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(phiHydLowtave,1, bi, bj, myThid)
#ifndef MINIMAL_TAVE_OUTPUT
        CALL TIMEAVE_RESET(UTtave,    Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(VTtave,    Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(WTtave,    Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(UStave,    Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(VStave,    Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(WStave,    Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(Eta2tave,  1,  bi, bj, myThid)
        CALL TIMEAVE_RESET(TTtave,    Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(UUtave,    Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(VVtave,    Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(UVtave,    Nr, bi, bj, myThid)
c       CALL TIMEAVE_RESET(KEtave,    Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(TdiffRtave,Nr, bi, bj, myThid)
#ifdef ALLOW_MOM_VECINV
        CALL TIMEAVE_RESET(uZetatave, Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(vZetatave, Nr, bi, bj, myThid)
#endif
        CALL TIMEAVE_RESET(phiHydtave,Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(phiHydLow2Tave,1, bi, bj, myThid)
        CALL TIMEAVE_RESET(ConvectCountTave,Nr,bi,bj,myThid)
#ifdef NONLIN_FRSURF
        CALL TIMEAVE_RESET(hUtave,    Nr, bi, bj, myThid)
        CALL TIMEAVE_RESET(hVtave,    Nr, bi, bj, myThid)
c       CALL TIMEAVE_RESET(hFacCtave, Nr, bi, bj, myThid)
c       CALL TIMEAVE_RESET(hFacWtave, Nr, bi, bj, myThid)
c       CALL TIMEAVE_RESET(hFacStave, Nr, bi, bj, myThid)
#endif /* NONLIN_FRSURF */
#endif /* ndef MINIMAL_TAVE_OUTPUT */
        DO K=1,Nr
         TimeAve_half(k,bi,bj)=0.
         TimeAve_full(k,bi,bj)=0.
        ENDDO
      ENDIF
      
C--   Cumulate state-variables with Half or Full time step : 
      DDTT = deltaTclock
      IF ( myIter .EQ. nIter0 ) THEN
        DDTT = deltaTclock*(1. _d 0 - tave_lastIter)
#ifdef ALLOW_CAL
      ELSEIF ( calendarDumps .AND. (
     &       ( taveFreq.GE. 2592000 .AND. taveFreq.LE. 2678400 ) .OR.
     &       ( taveFreq.GE.31104000 .AND. taveFreq.LE.31968000 ))) THEN
C--   Convert approximate months (30-31 days) and years (360-372 days)
C     to exact calendar months and years.
C-    First determine calendar dates for this and previous time step.
       call cal_GetDate( myiter  ,mytime            ,thisdate,mythid )
       call cal_GetDate( myiter-1,mytime-deltaTClock,prevdate,mythid )
C-    Monthly taveFreq:
       IF(taveFreq.GE. 2592000 .AND. taveFreq.LE. 2678400 .AND.
     &  (thisdate(1)-prevdate(1)).GT.50  )DDTT=deltaTclock*tave_lastIter
C-    Yearly  taveFreq:
       IF(taveFreq.GE.31104000 .AND. taveFreq.LE.31968000 .AND.
     &  (thisdate(1)-prevdate(1)).GT.5000)DDTT=deltaTclock*tave_lastIter
#endif
      ELSEIF ( 
     &  DIFF_BASE_MULTIPLE(baseTime,taveFreq,myTime,deltaTClock) ) THEN
        DDTT = deltaTclock*tave_lastIter
      ENDIF

      IF ( DDTT .NE. 0. _d 0) THEN

C-    Time Averages of single fields (no hFactor)
      CALL TIMEAVE_CUMULATE(etatave,  etaN,  1 , DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMULATE(thetatave,theta, Nr, DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMULATE(salttave, salt,  Nr, DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMULATE(uVeltave, uVel,  Nr, DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMULATE(vVeltave, vVel,  Nr, DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMULATE(wVeltave, wVel,  Nr, DDTT, bi, bj, myThid)

#ifndef MINIMAL_TAVE_OUTPUT
C-    Time Averages of "double" fields (no hFactor)
      CALL TIMEAVE_CUMUL_2V(Eta2tave, etaN,etaN, 1,  0,
     &     DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMUL_2V(TTtave, theta,theta, Nr, 0,
     &     DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMUL_2V(UUtave, uVel,  uVel, Nr, 0,
     &     DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMUL_2V(VVtave, vVel,  vVel, Nr, 0,
     &     DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMUL_2V(UVtave, uVel,  vVel, Nr, 12,
     &     DDTT, bi, bj, myThid)
c     CALL TIMEAVE_CUMUL_KE(KEtave, uVel,  vVel, Nr,
c    &     DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMUL_2V(WTtave, theta, wVel, Nr, 3,
     &     DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMUL_2V(WStave, salt,  wVel, Nr, 3,
     &     DDTT, bi, bj, myThid)

#ifdef NONLIN_FRSURF

c     CALL TIMEAVE_CUMUL_FC(hFacCtave,hFacC, Nr, DDTT, bi, bj, myThid)
c     CALL TIMEAVE_CUMUL_FC(hFacWtave,hFacW, Nr, DDTT, bi, bj, myThid)
c     CALL TIMEAVE_CUMUL_FC(hFacStave,hFacS, Nr, DDTT, bi, bj, myThid)

C-    Time Averages of single fields (* hFactor)
      CALL TIMEAVE_CUMUL_1VFC(hUtave,  uVel,  hFacW, Nr,
     &     DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMUL_1VFC(hVtave,  vVel,  hFacS, Nr,
     &     DDTT, bi, bj, myThid)

#endif /* NONLIN_FRSURF */

C-    Time Averages of "double" fields (* hFactor)
      CALL TIMEAVE_CUMUL_2VFC(UTtave, theta, uVel,  hFacW, Nr, 1,
     &     DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMUL_2VFC(VTtave, theta, vVel,  hFacS, Nr, 2,
     &     DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMUL_2VFC(UStave, salt, uVel,  hFacW, Nr, 1,
     &     DDTT, bi, bj, myThid)
      CALL TIMEAVE_CUMUL_2VFC(VStave, salt, vVel,  hFacS, Nr, 2,
     &     DDTT, bi, bj, myThid)

C-    Time Averages of "double" fields (no hFactor)
c     CALL TIMEAVE_CUMUL_2V(UTtave, theta, uVel, Nr, 1,
c    &     DDTT, bi, bj, myThid)
c     CALL TIMEAVE_CUMUL_2V(VTtave, theta, vVel, Nr, 2,
c    &     DDTT, bi, bj, myThid)
c     CALL TIMEAVE_CUMUL_2V(UStave, salt, uVel, Nr, 1,
c    &     DDTT, bi, bj, myThid)
c     CALL TIMEAVE_CUMUL_2V(VStave, salt, vVel, Nr, 2,
c    &     DDTT, bi, bj, myThid)
#endif /* ndef MINIMAL_TAVE_OUTPUT */

C-    Keep record of how much time has been integrated over
      DO K=1,Nr
        TimeAve_half(k,bi,bj)=TimeAve_half(k,bi,bj)+DDTT
      ENDDO

C-- end if DDTT ...
      ENDIF

C-    Time Averages of "intermediate" fields (no hFactor)
      IF ( myIter .NE. nIter0 ) THEN
 
C-    Time Averages of surface fluxes
C     <- moved to external_forcing_surf

        CALL TIMEAVE_CUMULATE(phiHydLowtave, phiHydLow, 1,
     &                                   deltaTclock, bi, bj, myThid)
#ifndef MINIMAL_TAVE_OUTPUT
        CALL TIMEAVE_CUMULATE(phiHydtave, totPhihyd,  Nr, 
     &                                   deltaTclock, bi, bj, myThid)
        CALL TIMEAVE_CUMUL_2V(phiHydLow2Tave, 
     &       phiHydLow,phiHydLow, 1,  0, deltaTclock, bi, bj, myThid)
#endif /* ndef MINIMAL_TAVE_OUTPUT */
        DO K=1,Nr
         TimeAve_full(k,bi,bj)=TimeAve_full(k,bi,bj)+deltaTclock
        ENDDO
      ENDIF

#endif /* ALLOW_TIMEAVE */ 

      RETURN
      END
1	jmc	1.20	C $Header: /u/gcmpack/MITgcm/pkg/timeave/timeave_statvars.F,v 1.19 2005/02/20 11:46:25 dimitri Exp $
2	jmc	1.1	C $Name: $
3	edhill	1.9	#include "TIMEAVE_OPTIONS.h"
4	jmc	1.1
5			SUBROUTINE TIMEAVE_STATVARS(
6			I myTime, myIter, bi, bj, myThid)
7			C /==========================================================\
8			C \| SUBROUTINE TIMEAVE_STATVARS \|
9			C \| o Time averaging routine for eta, U, V, W, T, S, UT, VT \|
10			C \| in model main time-stepping \|
11			C \==========================================================/
12			IMPLICIT NONE
13
14			C == Global variables ===
15			#include "SIZE.h"
16			#include "EEPARAMS.h"
17			#include "PARAMS.h"
18			#include "DYNVARS.h"
19	jmc	1.3	#include "GRID.h"
20	jmc	1.1	#include "TIMEAVE_STATV.h"
21
22	jmc	1.20	LOGICAL DIFF_BASE_MULTIPLE
23			EXTERNAL DIFF_BASE_MULTIPLE
24	jmc	1.1
25			C == Routine arguments ==
26			C myThid - Thread number for this instance of the routine.
27			C myIter - Iteration number
28			C myTime - Current time of simulation ( s )
29			INTEGER myThid
30			INTEGER myIter, bi, bj
31			_RL myTime
32
33	jmc	1.3	#ifdef ALLOW_TIMEAVE
34
35	jmc	1.1	C == Local variables ==
36	dimitri	1.13	INTEGER I, J, K
37	jmc	1.1	_RL DDTT
38	dimitri	1.13	_RL tempArray (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
39	dimitri	1.19	INTEGER thisdate(4), prevdate(4)
40	jmc	1.1
41			C- Initialize fields for the first call ever
42			IF ( myIter .EQ. nIter0 ) THEN
43	dimitri	1.13	CALL TIMEAVE_RESET(uFluxtave, 1, bi, bj, myThid)
44			CALL TIMEAVE_RESET(vFluxtave, 1, bi, bj, myThid)
45			CALL TIMEAVE_RESET(tFluxtave, 1, bi, bj, myThid)
46			CALL TIMEAVE_RESET(sFluxtave, 1, bi, bj, myThid)
47	jmc	1.3	CALL TIMEAVE_RESET(etatave, 1, bi, bj, myThid)
48	jmc	1.1	CALL TIMEAVE_RESET(thetatave, Nr, bi, bj, myThid)
49			CALL TIMEAVE_RESET(salttave, Nr, bi, bj, myThid)
50			CALL TIMEAVE_RESET(uVeltave, Nr, bi, bj, myThid)
51			CALL TIMEAVE_RESET(vVeltave, Nr, bi, bj, myThid)
52			CALL TIMEAVE_RESET(wVeltave, Nr, bi, bj, myThid)
53	dimitri	1.18	CALL TIMEAVE_RESET(phiHydLowtave,1, bi, bj, myThid)
54			#ifndef MINIMAL_TAVE_OUTPUT
55	dimitri	1.15	CALL TIMEAVE_RESET(UTtave, Nr, bi, bj, myThid)
56			CALL TIMEAVE_RESET(VTtave, Nr, bi, bj, myThid)
57			CALL TIMEAVE_RESET(WTtave, Nr, bi, bj, myThid)
58	adcroft	1.16	CALL TIMEAVE_RESET(UStave, Nr, bi, bj, myThid)
59			CALL TIMEAVE_RESET(VStave, Nr, bi, bj, myThid)
60	adcroft	1.17	CALL TIMEAVE_RESET(WStave, Nr, bi, bj, myThid)
61	jmc	1.3	CALL TIMEAVE_RESET(Eta2tave, 1, bi, bj, myThid)
62	adcroft	1.2	CALL TIMEAVE_RESET(TTtave, Nr, bi, bj, myThid)
63			CALL TIMEAVE_RESET(UUtave, Nr, bi, bj, myThid)
64			CALL TIMEAVE_RESET(VVtave, Nr, bi, bj, myThid)
65	jmc	1.7	CALL TIMEAVE_RESET(UVtave, Nr, bi, bj, myThid)
66	jmc	1.3	c CALL TIMEAVE_RESET(KEtave, Nr, bi, bj, myThid)
67	jmc	1.7	CALL TIMEAVE_RESET(TdiffRtave,Nr, bi, bj, myThid)
68	edhill	1.12	#ifdef ALLOW_MOM_VECINV
69	jmc	1.8	CALL TIMEAVE_RESET(uZetatave, Nr, bi, bj, myThid)
70			CALL TIMEAVE_RESET(vZetatave, Nr, bi, bj, myThid)
71			#endif
72	jmc	1.3	CALL TIMEAVE_RESET(phiHydtave,Nr, bi, bj, myThid)
73	mlosch	1.4	CALL TIMEAVE_RESET(phiHydLow2Tave,1, bi, bj, myThid)
74	jmc	1.1	CALL TIMEAVE_RESET(ConvectCountTave,Nr,bi,bj,myThid)
75	jmc	1.3	#ifdef NONLIN_FRSURF
76			CALL TIMEAVE_RESET(hUtave, Nr, bi, bj, myThid)
77			CALL TIMEAVE_RESET(hVtave, Nr, bi, bj, myThid)
78			c CALL TIMEAVE_RESET(hFacCtave, Nr, bi, bj, myThid)
79			c CALL TIMEAVE_RESET(hFacWtave, Nr, bi, bj, myThid)
80			c CALL TIMEAVE_RESET(hFacStave, Nr, bi, bj, myThid)
81			#endif /* NONLIN_FRSURF */
82	dimitri	1.18	#endif /* ndef MINIMAL_TAVE_OUTPUT */
83	jmc	1.1	DO K=1,Nr
84			TimeAve_half(k,bi,bj)=0.
85			TimeAve_full(k,bi,bj)=0.
86			ENDDO
87			ENDIF
88
89			C-- Cumulate state-variables with Half or Full time step :
90	dimitri	1.19	DDTT = deltaTclock
91	jmc	1.5	IF ( myIter .EQ. nIter0 ) THEN
92			DDTT = deltaTclock*(1. _d 0 - tave_lastIter)
93	dimitri	1.19	#ifdef ALLOW_CAL
94			ELSEIF ( calendarDumps .AND. (
95			& ( taveFreq.GE. 2592000 .AND. taveFreq.LE. 2678400 ) .OR.
96			& ( taveFreq.GE.31104000 .AND. taveFreq.LE.31968000 ))) THEN
97			C-- Convert approximate months (30-31 days) and years (360-372 days)
98			C to exact calendar months and years.
99			C- First determine calendar dates for this and previous time step.
100			call cal_GetDate( myiter ,mytime ,thisdate,mythid )
101			call cal_GetDate( myiter-1,mytime-deltaTClock,prevdate,mythid )
102			C- Monthly taveFreq:
103			IF(taveFreq.GE. 2592000 .AND. taveFreq.LE. 2678400 .AND.
104			& (thisdate(1)-prevdate(1)).GT.50 )DDTT=deltaTclock*tave_lastIter
105			C- Yearly taveFreq:
106			IF(taveFreq.GE.31104000 .AND. taveFreq.LE.31968000 .AND.
107			& (thisdate(1)-prevdate(1)).GT.5000)DDTT=deltaTclock*tave_lastIter
108			#endif
109	jmc	1.5	ELSEIF (
110	jmc	1.20	& DIFF_BASE_MULTIPLE(baseTime,taveFreq,myTime,deltaTClock) ) THEN
111	jmc	1.5	DDTT = deltaTclock*tave_lastIter
112	jmc	1.1	ENDIF
113
114	jmc	1.5	IF ( DDTT .NE. 0. _d 0) THEN
115
116	jmc	1.3	C- Time Averages of single fields (no hFactor)
117			CALL TIMEAVE_CUMULATE(etatave, etaN, 1 , DDTT, bi, bj, myThid)
118	jmc	1.1	CALL TIMEAVE_CUMULATE(thetatave,theta, Nr, DDTT, bi, bj, myThid)
119			CALL TIMEAVE_CUMULATE(salttave, salt, Nr, DDTT, bi, bj, myThid)
120			CALL TIMEAVE_CUMULATE(uVeltave, uVel, Nr, DDTT, bi, bj, myThid)
121			CALL TIMEAVE_CUMULATE(vVeltave, vVel, Nr, DDTT, bi, bj, myThid)
122	jmc	1.7	CALL TIMEAVE_CUMULATE(wVeltave, wVel, Nr, DDTT, bi, bj, myThid)
123	jmc	1.1
124	dimitri	1.18	#ifndef MINIMAL_TAVE_OUTPUT
125	jmc	1.3	C- Time Averages of "double" fields (no hFactor)
126			CALL TIMEAVE_CUMUL_2V(Eta2tave, etaN,etaN, 1, 0,
127			& DDTT, bi, bj, myThid)
128			CALL TIMEAVE_CUMUL_2V(TTtave, theta,theta, Nr, 0,
129			& DDTT, bi, bj, myThid)
130			CALL TIMEAVE_CUMUL_2V(UUtave, uVel, uVel, Nr, 0,
131			& DDTT, bi, bj, myThid)
132			CALL TIMEAVE_CUMUL_2V(VVtave, vVel, vVel, Nr, 0,
133	jmc	1.7	& DDTT, bi, bj, myThid)
134			CALL TIMEAVE_CUMUL_2V(UVtave, uVel, vVel, Nr, 12,
135	jmc	1.3	& DDTT, bi, bj, myThid)
136			c CALL TIMEAVE_CUMUL_KE(KEtave, uVel, vVel, Nr,
137			c & DDTT, bi, bj, myThid)
138	jmc	1.5	CALL TIMEAVE_CUMUL_2V(WTtave, theta, wVel, Nr, 3,
139	adcroft	1.17	& DDTT, bi, bj, myThid)
140			CALL TIMEAVE_CUMUL_2V(WStave, salt, wVel, Nr, 3,
141	jmc	1.5	& DDTT, bi, bj, myThid)
142	jmc	1.3
143			#ifdef NONLIN_FRSURF
144
145			c CALL TIMEAVE_CUMUL_FC(hFacCtave,hFacC, Nr, DDTT, bi, bj, myThid)
146			c CALL TIMEAVE_CUMUL_FC(hFacWtave,hFacW, Nr, DDTT, bi, bj, myThid)
147			c CALL TIMEAVE_CUMUL_FC(hFacStave,hFacS, Nr, DDTT, bi, bj, myThid)
148
149			C- Time Averages of single fields (* hFactor)
150			CALL TIMEAVE_CUMUL_1VFC(hUtave, uVel, hFacW, Nr,
151			& DDTT, bi, bj, myThid)
152			CALL TIMEAVE_CUMUL_1VFC(hVtave, vVel, hFacS, Nr,
153			& DDTT, bi, bj, myThid)
154
155			#endif /* NONLIN_FRSURF */
156
157			C- Time Averages of "double" fields (* hFactor)
158			CALL TIMEAVE_CUMUL_2VFC(UTtave, theta, uVel, hFacW, Nr, 1,
159			& DDTT, bi, bj, myThid)
160			CALL TIMEAVE_CUMUL_2VFC(VTtave, theta, vVel, hFacS, Nr, 2,
161			& DDTT, bi, bj, myThid)
162	adcroft	1.16	CALL TIMEAVE_CUMUL_2VFC(UStave, salt, uVel, hFacW, Nr, 1,
163			& DDTT, bi, bj, myThid)
164			CALL TIMEAVE_CUMUL_2VFC(VStave, salt, vVel, hFacS, Nr, 2,
165			& DDTT, bi, bj, myThid)
166	jmc	1.3
167			C- Time Averages of "double" fields (no hFactor)
168			c CALL TIMEAVE_CUMUL_2V(UTtave, theta, uVel, Nr, 1,
169			c & DDTT, bi, bj, myThid)
170			c CALL TIMEAVE_CUMUL_2V(VTtave, theta, vVel, Nr, 2,
171	adcroft	1.16	c & DDTT, bi, bj, myThid)
172			c CALL TIMEAVE_CUMUL_2V(UStave, salt, uVel, Nr, 1,
173			c & DDTT, bi, bj, myThid)
174			c CALL TIMEAVE_CUMUL_2V(VStave, salt, vVel, Nr, 2,
175	jmc	1.3	c & DDTT, bi, bj, myThid)
176	dimitri	1.18	#endif /* ndef MINIMAL_TAVE_OUTPUT */
177	jmc	1.1
178			C- Keep record of how much time has been integrated over
179			DO K=1,Nr
180			TimeAve_half(k,bi,bj)=TimeAve_half(k,bi,bj)+DDTT
181			ENDDO
182
183	jmc	1.5	C-- end if DDTT ...
184			ENDIF
185	jmc	1.3
186			C- Time Averages of "intermediate" fields (no hFactor)
187	jmc	1.1	IF ( myIter .NE. nIter0 ) THEN
188	dimitri	1.13
189			C- Time Averages of surface fluxes
190	jmc	1.14	C <- moved to external_forcing_surf
191	dimitri	1.13
192	dimitri	1.18	CALL TIMEAVE_CUMULATE(phiHydLowtave, phiHydLow, 1,
193			& deltaTclock, bi, bj, myThid)
194			#ifndef MINIMAL_TAVE_OUTPUT
195	jmc	1.6	CALL TIMEAVE_CUMULATE(phiHydtave, totPhihyd, Nr,
196			& deltaTclock, bi, bj, myThid)
197	mlosch	1.4	CALL TIMEAVE_CUMUL_2V(phiHydLow2Tave,
198			& phiHydLow,phiHydLow, 1, 0, deltaTclock, bi, bj, myThid)
199	dimitri	1.18	#endif /* ndef MINIMAL_TAVE_OUTPUT */
200	jmc	1.1	DO K=1,Nr
201			TimeAve_full(k,bi,bj)=TimeAve_full(k,bi,bj)+deltaTclock
202			ENDDO
203			ENDIF
204	jmc	1.3
205			#endif /* ALLOW_TIMEAVE */
206	jmc	1.1
207			RETURN
208			END