pkg/timeave/timeave_statv_write.F

C $Header: /u/gcmpack/MITgcm/pkg/timeave/timeave_statv_write.F,v 1.13 2003/11/13 06:35:15 dimitri Exp $
C $Name:  $
#include "TIMEAVE_OPTIONS.h"

      SUBROUTINE TIMEAVE_STATV_WRITE(myTime, myIter, myThid)
C     /==========================================================\
C     | SUBROUTINE TIMEAVE_STATV_WRITE                           |
C     | o At the end of average period, write the time-average   |
C     |   state-variables on file ; then reset for next period   |
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  DIFFERENT_MULTIPLE
      EXTERNAL DIFFERENT_MULTIPLE
      INTEGER  IO_ERRCOUNT
      EXTERNAL IO_ERRCOUNT

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
      _RL     myTime

#ifdef ALLOW_TIMEAVE 

C     == Local variables ==
C     suff - Hold suffix part of a filename
C     TimeAve - total time over average
C     useVariableK   = T when vertical diffusion is not constant
      LOGICAL useVariableK
      INTEGER prevPrec
      CHARACTER*(MAX_LEN_FNAM) suff
      INTEGER bi,bj,k
      INTEGER beginIOErrCount
      INTEGER endIOErrCount
      CHARACTER*(MAX_LEN_MBUF) msgBuf
 
C-    Final Time Averages and Dump Files if needed
      IF(DIFFERENT_MULTIPLE(taveFreq,myTime,myTime-deltaTClock) )THEN

       useVariableK = useKPP .OR. useGMredi .OR. ivdc_kappa.NE.0.
       WRITE(suff,'(I10.10)') myIter

       DO bj = myByLo(myThid), myByHi(myThid)
        DO bi = myBxLo(myThid), myBxHi(myThid)
C-     Normalize by integrated time
          CALL TIMEAVE_NORMALIZ(uFluxtave,TimeAve_full,1 ,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(vFluxtave,TimeAve_full,1 ,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(tFluxtave,TimeAve_full,1 ,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(sFluxtave,TimeAve_full,1 ,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(etatave,  TimeAve_half,1 ,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(thetatave,TimeAve_half,Nr,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(salttave, TimeAve_half,Nr,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(uVeltave, TimeAve_half,Nr,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(vVeltave, TimeAve_half,Nr,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(wVeltave, TimeAve_half,Nr,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(Eta2tave, TimeAve_half,1 ,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(TTtave,   TimeAve_half,Nr,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(UUtave,   TimeAve_half,Nr,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(VVtave,   TimeAve_half,Nr,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(UVtave,   TimeAve_half,Nr,bi,bj,myThid)
c         CALL TIMEAVE_NORMALIZ(KEtave,   TimeAve_half,Nr,bi,bj,myThid)
#ifdef NONLIN_FRSURF
C-     Normalize by integrated time
          CALL TIMEAVE_NORMALIZ(hUtave,   TimeAve_half,Nr,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(hVtave,   TimeAve_half,Nr,bi,bj,myThid)
c         CALL TIMEAVE_NORMALIZ(hFacCtave,TimeAve_half,Nr,bi,bj,myThid)
c         CALL TIMEAVE_NORMALIZ(hFacWtave,TimeAve_half,Nr,bi,bj,myThid)
c         CALL TIMEAVE_NORMALIZ(hFacStave,TimeAve_half,Nr,bi,bj,myThid)
#endif /* NONLIN_FRSURF */ 
          CALL TIMEAVE_NORMALIZ(UTtave,   TimeAve_half,Nr,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(VTtave,   TimeAve_half,Nr,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(WTtave,   TimeAve_half,Nr,bi,bj,myThid)

          CALL TIMEAVE_NORMALIZ(TdiffRtave,TimeAve_full,Nr,
     &                            bi,bj,myThid)
#ifdef ALLOW_MOM_VECINV
          CALL TIMEAVE_NORMALIZ(uZetatave,TimeAve_full,Nr,bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(vZetatave,TimeAve_full,Nr,bi,bj,myThid)
#endif
          CALL TIMEAVE_NORMALIZ(phiHydtave,TimeAve_full,Nr,
     &                            bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(phiHydLowtave,TimeAve_full,1,
     &                            bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(phiHydLow2Tave,TimeAve_full,1,
     &                            bi,bj,myThid)
          CALL TIMEAVE_NORMALIZ(ConvectCountTave,TimeAve_full,Nr,
     &                            bi,bj,myThid)
        ENDDO
       ENDDO

C      Write to files
       _BARRIER
       _BEGIN_MASTER( myThid )

C--     Set IO "context" for writing state
#ifdef USE_DFILE
       CALL DFILE_SET_RW
       CALL DFILE_SET_CONT_ON_ERROR
#endif
C--    Read IO error counter
       beginIOErrCount = IO_ERRCOUNT(myThid)
c      prevPrec = writeBinaryPrec
c      writeBinaryPrec = precFloat32
       CALL WRITE_FLD_XY_RL('ETAtave.'  ,suff,etatave  ,myIter,myThid)
       CALL WRITE_FLD_XY_RL('Eta2tave.' ,suff,Eta2tave ,myIter,myThid)

       CALL WRITE_FLD_XYZ_RL('Ttave.',suff,thetatave,myIter,myThid)
       CALL WRITE_FLD_XYZ_RL('Stave.',suff,salttave,myIter,myThid)
       CALL WRITE_FLD_XYZ_RL('uVeltave.',suff,uVeltave,myIter,myThid)
       CALL WRITE_FLD_XYZ_RL('vVeltave.',suff,vVeltave,myIter,myThid)
       CALL WRITE_FLD_XYZ_RL('wVeltave.',suff,wVeltave,myIter,myThid)
       CALL WRITE_FLD_XYZ_RL('TTtave.',suff,TTtave,myIter,myThid)
       CALL WRITE_FLD_XYZ_RL('UUtave.',suff,UUtave,myIter,myThid)
       CALL WRITE_FLD_XYZ_RL('VVtave.',suff,VVtave,myIter,myThid)
       CALL WRITE_FLD_XYZ_RL('UVtave.',suff,UVtave,myIter,myThid)
c      CALL WRITE_FLD_XYZ_RL('KEtave.',suff,KEtave,myIter,myThid)
       CALL WRITE_FLD_XYZ_RL('UTtave.',suff,UTtave,myIter,myThid)
       CALL WRITE_FLD_XYZ_RL('VTtave.',suff,VTtave,myIter,myThid)
       CALL WRITE_FLD_XYZ_RL('WTtave.',suff,WTtave,myIter,myThid)

       IF (useVariableK)
     & CALL WRITE_FLD_XYZ_RL('Tdiftave.',suff,TdiffRtave,myIter,myThid)
#ifdef ALLOW_MOM_VECINV
       IF (vectorInvariantMomentum) THEN
         CALL WRITE_FLD_XYZ_RL('uZtave.',suff,uZetatave,myIter,myThid)
         CALL WRITE_FLD_XYZ_RL('vZtave.',suff,vZetatave,myIter,myThid)
       ENDIF
#endif
       CALL WRITE_FLD_XYZ_RL('PhHytave.',suff,phiHydtave,myIter,myThid)
       IF ( buoyancyRelation(1:7) .EQ. 'OCEANIC' ) THEN
        CALL WRITE_FLD_XY_RL('PHLtave.',suff,phiHydLowtave,
     &       myIter,myThid)
        CALL WRITE_FLD_XY_RL('PHL2tave.',suff,phiHydLow2tave,
     &       myIter,myThid)
        CALL WRITE_FLD_XYZ_RL('Convtave.',suff,ConvectCountTave,
     &                        myIter,myThid)
        CALL WRITE_FLD_XY_RL('uFluxtave.',suff,uFluxtave,myIter,myThid)
        CALL WRITE_FLD_XY_RL('vFluxtave.',suff,vFluxtave,myIter,myThid)
        CALL WRITE_FLD_XY_RL('tFluxtave.',suff,tFluxtave,myIter,myThid)
        CALL WRITE_FLD_XY_RL('sFluxtave.',suff,sFluxtave,myIter,myThid)
       ENDIF

#ifdef NONLIN_FRSURF
       CALL WRITE_FLD_XYZ_RL('hUtave.',suff,hUtave,myIter,myThid)
       CALL WRITE_FLD_XYZ_RL('hVtave.',suff,hVtave,myIter,myThid)
c      CALL WRITE_FLD_XYZ_RL('hFacCtave.',suff,hFacCtave,myIter,myThid)
c      CALL WRITE_FLD_XYZ_RL('hFacWtave.',suff,hFacWtave,myIter,myThid)
c      CALL WRITE_FLD_XYZ_RL('hFacStave.',suff,hFacStave,myIter,myThid)
#endif /* NONLIN_FRSURF */ 

c      writeBinaryPrec = prevPrec
C--    Reread IO error counter
       endIOErrCount = IO_ERRCOUNT(myThid)

C--    Check for IO errors
       IF ( endIOErrCount .NE. beginIOErrCount ) THEN
         WRITE(msgBuf,'(A)')  'S/R WRITE_TIME_AVERAGES'
         CALL PRINT_ERROR( msgBuf, 1 )
         WRITE(msgBuf,'(A)')  'Error writing out data'
         CALL PRINT_ERROR( msgBuf, 1 )
         WRITE(msgBuf,'(A,I10)') 'Timestep ',myIter
         CALL PRINT_ERROR( msgBuf, 1 )
       ELSE
         WRITE(msgBuf,'(A,I10)')  
     &    '// Time-average data written, t-step', myIter
         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, 
     &    SQUEEZE_RIGHT, 1 )
         WRITE(msgBuf,'(A)')  ' '
         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, 
     &    SQUEEZE_RIGHT, 1 )
       ENDIF

       _END_MASTER( myThid )
       _BARRIER

       DO bj = myByLo(myThid), myByHi(myThid)
        DO bi = myBxLo(myThid), myBxHi(myThid)

C-- Like before the 1rst iteration,
C  ==> call TIMEAVE_STATVARS with myIter=nIter0 :
C    1) Reset the averages to zero ;
C    2) Start to cumulate state-variables with Half time step.

         CALL TIMEAVE_STATVARS(myTime, nIter0, bi, bj, myThid)

        ENDDO
       ENDDO

      ENDIF
 
#endif /* ALLOW_TIMEAVE */

      RETURN
      END
1	jmc	1.14	C $Header: /u/gcmpack/MITgcm/pkg/timeave/timeave_statv_write.F,v 1.13 2003/11/13 06:35:15 dimitri Exp $
2	jmc	1.1	C $Name: $
3	edhill	1.9	#include "TIMEAVE_OPTIONS.h"
4	jmc	1.1
5			SUBROUTINE TIMEAVE_STATV_WRITE(myTime, myIter, myThid)
6			C /==========================================================\
7			C \| SUBROUTINE TIMEAVE_STATV_WRITE \|
8			C \| o At the end of average period, write the time-average \|
9			C \| state-variables on file ; then reset for next period \|
10			C \==========================================================/
11			IMPLICIT NONE
12
13			C == Global variables ===
14			#include "SIZE.h"
15			#include "EEPARAMS.h"
16			#include "PARAMS.h"
17			#include "DYNVARS.h"
18	jmc	1.3	#include "GRID.h"
19	jmc	1.1	#include "TIMEAVE_STATV.h"
20
21			LOGICAL DIFFERENT_MULTIPLE
22			EXTERNAL DIFFERENT_MULTIPLE
23			INTEGER IO_ERRCOUNT
24			EXTERNAL IO_ERRCOUNT
25
26			C == Routine arguments ==
27			C myThid - Thread number for this instance of the routine.
28			C myIter - Iteration number
29			C myTime - Current time of simulation ( s )
30			INTEGER myThid
31			INTEGER myIter
32			_RL myTime
33
34	jmc	1.3	#ifdef ALLOW_TIMEAVE
35
36	jmc	1.1	C == Local variables ==
37			C suff - Hold suffix part of a filename
38			C TimeAve - total time over average
39	jmc	1.7	C useVariableK = T when vertical diffusion is not constant
40			LOGICAL useVariableK
41	jmc	1.1	INTEGER prevPrec
42			CHARACTER*(MAX_LEN_FNAM) suff
43			INTEGER bi,bj,k
44			INTEGER beginIOErrCount
45			INTEGER endIOErrCount
46			CHARACTER*(MAX_LEN_MBUF) msgBuf
47
48			C- Final Time Averages and Dump Files if needed
49	jmc	1.3	IF(DIFFERENT_MULTIPLE(taveFreq,myTime,myTime-deltaTClock) )THEN
50	jmc	1.1
51	jmc	1.7	useVariableK = useKPP .OR. useGMredi .OR. ivdc_kappa.NE.0.
52	jmc	1.1	WRITE(suff,'(I10.10)') myIter
53
54			DO bj = myByLo(myThid), myByHi(myThid)
55			DO bi = myBxLo(myThid), myBxHi(myThid)
56	jmc	1.3	C- Normalize by integrated time
57	dimitri	1.13	CALL TIMEAVE_NORMALIZ(uFluxtave,TimeAve_full,1 ,bi,bj,myThid)
58			CALL TIMEAVE_NORMALIZ(vFluxtave,TimeAve_full,1 ,bi,bj,myThid)
59			CALL TIMEAVE_NORMALIZ(tFluxtave,TimeAve_full,1 ,bi,bj,myThid)
60			CALL TIMEAVE_NORMALIZ(sFluxtave,TimeAve_full,1 ,bi,bj,myThid)
61	jmc	1.3	CALL TIMEAVE_NORMALIZ(etatave, TimeAve_half,1 ,bi,bj,myThid)
62	jmc	1.1	CALL TIMEAVE_NORMALIZ(thetatave,TimeAve_half,Nr,bi,bj,myThid)
63			CALL TIMEAVE_NORMALIZ(salttave, TimeAve_half,Nr,bi,bj,myThid)
64			CALL TIMEAVE_NORMALIZ(uVeltave, TimeAve_half,Nr,bi,bj,myThid)
65			CALL TIMEAVE_NORMALIZ(vVeltave, TimeAve_half,Nr,bi,bj,myThid)
66	jmc	1.5	CALL TIMEAVE_NORMALIZ(wVeltave, TimeAve_half,Nr,bi,bj,myThid)
67	jmc	1.3	CALL TIMEAVE_NORMALIZ(Eta2tave, TimeAve_half,1 ,bi,bj,myThid)
68	adcroft	1.2	CALL TIMEAVE_NORMALIZ(TTtave, TimeAve_half,Nr,bi,bj,myThid)
69			CALL TIMEAVE_NORMALIZ(UUtave, TimeAve_half,Nr,bi,bj,myThid)
70			CALL TIMEAVE_NORMALIZ(VVtave, TimeAve_half,Nr,bi,bj,myThid)
71	jmc	1.6	CALL TIMEAVE_NORMALIZ(UVtave, TimeAve_half,Nr,bi,bj,myThid)
72	jmc	1.3	c CALL TIMEAVE_NORMALIZ(KEtave, TimeAve_half,Nr,bi,bj,myThid)
73			#ifdef NONLIN_FRSURF
74			C- Normalize by integrated time
75			CALL TIMEAVE_NORMALIZ(hUtave, TimeAve_half,Nr,bi,bj,myThid)
76			CALL TIMEAVE_NORMALIZ(hVtave, TimeAve_half,Nr,bi,bj,myThid)
77			c CALL TIMEAVE_NORMALIZ(hFacCtave,TimeAve_half,Nr,bi,bj,myThid)
78			c CALL TIMEAVE_NORMALIZ(hFacWtave,TimeAve_half,Nr,bi,bj,myThid)
79			c CALL TIMEAVE_NORMALIZ(hFacStave,TimeAve_half,Nr,bi,bj,myThid)
80			#endif /* NONLIN_FRSURF */
81			CALL TIMEAVE_NORMALIZ(UTtave, TimeAve_half,Nr,bi,bj,myThid)
82			CALL TIMEAVE_NORMALIZ(VTtave, TimeAve_half,Nr,bi,bj,myThid)
83	jmc	1.5	CALL TIMEAVE_NORMALIZ(WTtave, TimeAve_half,Nr,bi,bj,myThid)
84	jmc	1.1
85	jmc	1.6	CALL TIMEAVE_NORMALIZ(TdiffRtave,TimeAve_full,Nr,
86			& bi,bj,myThid)
87	edhill	1.12	#ifdef ALLOW_MOM_VECINV
88	jmc	1.8	CALL TIMEAVE_NORMALIZ(uZetatave,TimeAve_full,Nr,bi,bj,myThid)
89			CALL TIMEAVE_NORMALIZ(vZetatave,TimeAve_full,Nr,bi,bj,myThid)
90			#endif
91	jmc	1.1	CALL TIMEAVE_NORMALIZ(phiHydtave,TimeAve_full,Nr,
92			& bi,bj,myThid)
93	mlosch	1.4	CALL TIMEAVE_NORMALIZ(phiHydLowtave,TimeAve_full,1,
94			& bi,bj,myThid)
95			CALL TIMEAVE_NORMALIZ(phiHydLow2Tave,TimeAve_full,1,
96			& bi,bj,myThid)
97	jmc	1.1	CALL TIMEAVE_NORMALIZ(ConvectCountTave,TimeAve_full,Nr,
98			& bi,bj,myThid)
99			ENDDO
100			ENDDO
101
102			C Write to files
103			_BARRIER
104			_BEGIN_MASTER( myThid )
105
106			C-- Set IO "context" for writing state
107			#ifdef USE_DFILE
108	jmc	1.3	CALL DFILE_SET_RW
109			CALL DFILE_SET_CONT_ON_ERROR
110	jmc	1.1	#endif
111	jmc	1.3	C-- Read IO error counter
112			beginIOErrCount = IO_ERRCOUNT(myThid)
113			c prevPrec = writeBinaryPrec
114			c writeBinaryPrec = precFloat32
115	dimitri	1.13	CALL WRITE_FLD_XY_RL('ETAtave.' ,suff,etatave ,myIter,myThid)
116			CALL WRITE_FLD_XY_RL('Eta2tave.' ,suff,Eta2tave ,myIter,myThid)
117	jmc	1.3
118			CALL WRITE_FLD_XYZ_RL('Ttave.',suff,thetatave,myIter,myThid)
119			CALL WRITE_FLD_XYZ_RL('Stave.',suff,salttave,myIter,myThid)
120			CALL WRITE_FLD_XYZ_RL('uVeltave.',suff,uVeltave,myIter,myThid)
121			CALL WRITE_FLD_XYZ_RL('vVeltave.',suff,vVeltave,myIter,myThid)
122			CALL WRITE_FLD_XYZ_RL('wVeltave.',suff,wVeltave,myIter,myThid)
123			CALL WRITE_FLD_XYZ_RL('TTtave.',suff,TTtave,myIter,myThid)
124			CALL WRITE_FLD_XYZ_RL('UUtave.',suff,UUtave,myIter,myThid)
125			CALL WRITE_FLD_XYZ_RL('VVtave.',suff,VVtave,myIter,myThid)
126	jmc	1.6	CALL WRITE_FLD_XYZ_RL('UVtave.',suff,UVtave,myIter,myThid)
127	jmc	1.3	c CALL WRITE_FLD_XYZ_RL('KEtave.',suff,KEtave,myIter,myThid)
128			CALL WRITE_FLD_XYZ_RL('UTtave.',suff,UTtave,myIter,myThid)
129			CALL WRITE_FLD_XYZ_RL('VTtave.',suff,VTtave,myIter,myThid)
130			CALL WRITE_FLD_XYZ_RL('WTtave.',suff,WTtave,myIter,myThid)
131	jmc	1.5
132	jmc	1.7	IF (useVariableK)
133			& CALL WRITE_FLD_XYZ_RL('Tdiftave.',suff,TdiffRtave,myIter,myThid)
134	edhill	1.12	#ifdef ALLOW_MOM_VECINV
135	jmc	1.8	IF (vectorInvariantMomentum) THEN
136			CALL WRITE_FLD_XYZ_RL('uZtave.',suff,uZetatave,myIter,myThid)
137			CALL WRITE_FLD_XYZ_RL('vZtave.',suff,vZetatave,myIter,myThid)
138			ENDIF
139			#endif
140	jmc	1.3	CALL WRITE_FLD_XYZ_RL('PhHytave.',suff,phiHydtave,myIter,myThid)
141	jmc	1.14	IF ( buoyancyRelation(1:7) .EQ. 'OCEANIC' ) THEN
142			CALL WRITE_FLD_XY_RL('PHLtave.',suff,phiHydLowtave,
143			& myIter,myThid)
144			CALL WRITE_FLD_XY_RL('PHL2tave.',suff,phiHydLow2tave,
145			& myIter,myThid)
146			CALL WRITE_FLD_XYZ_RL('Convtave.',suff,ConvectCountTave,
147			& myIter,myThid)
148			CALL WRITE_FLD_XY_RL('uFluxtave.',suff,uFluxtave,myIter,myThid)
149			CALL WRITE_FLD_XY_RL('vFluxtave.',suff,vFluxtave,myIter,myThid)
150			CALL WRITE_FLD_XY_RL('tFluxtave.',suff,tFluxtave,myIter,myThid)
151			CALL WRITE_FLD_XY_RL('sFluxtave.',suff,sFluxtave,myIter,myThid)
152	jmc	1.5	ENDIF
153
154	jmc	1.3	#ifdef NONLIN_FRSURF
155			CALL WRITE_FLD_XYZ_RL('hUtave.',suff,hUtave,myIter,myThid)
156			CALL WRITE_FLD_XYZ_RL('hVtave.',suff,hVtave,myIter,myThid)
157			c CALL WRITE_FLD_XYZ_RL('hFacCtave.',suff,hFacCtave,myIter,myThid)
158			c CALL WRITE_FLD_XYZ_RL('hFacWtave.',suff,hFacWtave,myIter,myThid)
159			c CALL WRITE_FLD_XYZ_RL('hFacStave.',suff,hFacStave,myIter,myThid)
160			#endif /* NONLIN_FRSURF */
161
162			c writeBinaryPrec = prevPrec
163			C-- Reread IO error counter
164			endIOErrCount = IO_ERRCOUNT(myThid)
165	jmc	1.1
166	jmc	1.3	C-- Check for IO errors
167			IF ( endIOErrCount .NE. beginIOErrCount ) THEN
168	jmc	1.1	WRITE(msgBuf,'(A)') 'S/R WRITE_TIME_AVERAGES'
169			CALL PRINT_ERROR( msgBuf, 1 )
170			WRITE(msgBuf,'(A)') 'Error writing out data'
171			CALL PRINT_ERROR( msgBuf, 1 )
172			WRITE(msgBuf,'(A,I10)') 'Timestep ',myIter
173			CALL PRINT_ERROR( msgBuf, 1 )
174	jmc	1.3	ELSE
175	jmc	1.1	WRITE(msgBuf,'(A,I10)')
176			& '// Time-average data written, t-step', myIter
177			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
178			& SQUEEZE_RIGHT, 1 )
179			WRITE(msgBuf,'(A)') ' '
180			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
181			& SQUEEZE_RIGHT, 1 )
182	jmc	1.3	ENDIF
183	jmc	1.1
184			_END_MASTER( myThid )
185			_BARRIER
186
187			DO bj = myByLo(myThid), myByHi(myThid)
188			DO bi = myBxLo(myThid), myBxHi(myThid)
189
190	jmc	1.3	C-- Like before the 1rst iteration,
191			C ==> call TIMEAVE_STATVARS with myIter=nIter0 :
192			C 1) Reset the averages to zero ;
193			C 2) Start to cumulate state-variables with Half time step.
194
195			CALL TIMEAVE_STATVARS(myTime, nIter0, bi, bj, myThid)
196	jmc	1.1
197			ENDDO
198			ENDDO
199
200			ENDIF
201
202	jmc	1.3	#endif /* ALLOW_TIMEAVE */
203
204	jmc	1.1	RETURN
205			END