diags/src/write_time_averages.F

#include "CPP_EEOPTIONS.h"

      SUBROUTINE WRITE_TIME_AVERAGES(myCurrentTime, myIter, myThid)
C     /==========================================================\
C     | SUBROUTINE WRITE_TIME_AVERAGES                           |
C     |                                                          |
C     \==========================================================/

C     == Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "AVER.h"
#include "DYNVARS.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     myCurrentTime - Current time of simulation ( s )
      INTEGER myThid
      INTEGER myIter
      REAL    myCurrentTime
      INTEGER prevPrec

C     == Local variables ==
C     suff - Hold suffix part of a filename
C     TimeAve - total time over average
      CHARACTER*(MAX_LEN_FNAM) suff
      REAL DDTT
      INTEGER bi,bj,k
      INTEGER beginIOErrCount
      INTEGER endIOErrCount
      CHARACTER*(MAX_LEN_MBUF) msgBuf
 
C     Final Time Averages and Dump Files if needed
      IF ( myIter.NE.nIter0 .AND. 
     & DIFFERENT_MULTIPLE(taveFreq,myCurrentTime,myCurrentTime-deltaTClock)
     &   ) THEN

       WRITE(suff,'(I10.10)') myIter

C      Normalize by integrated time
       DO bj = myByLo(myThid), myByHi(myThid)
        DO bi = myBxLo(myThid), myBxHi(myThid)
         DO k=1,Nz
          CALL TIMEAVER_FINAL_XYZ(thetatave,timeave,bi,bj,K,myThid)
          CALL TIMEAVER_FINAL_XYZ(salttave,timeave,bi,bj,K,myThid)
          CALL TIMEAVER_FINAL_XYZ(uVeltave,timeave,bi,bj,K,myThid)
          CALL TIMEAVER_FINAL_XYZ(vVeltave,timeave,bi,bj,K,myThid)
          CALL TIMEAVER_FINAL_XYZ(uttave,timeave,bi,bj,K,myThid)
          CALL TIMEAVER_FINAL_XYZ(vttave,timeave,bi,bj,K,myThid)
          CALL TIMEAVER_FINAL_XYZ(k13tave,intertimeave,bi,bj,K,myThid)
          CALL TIMEAVER_FINAL_XYZ(k23tave,intertimeave,bi,bj,K,myThid)
          CALL TIMEAVER_FINAL_XYZ(wVeltave,intertimeave,bi,bj,K,myThid)
         ENDDO
         CALL TIMEAVER_FINAL_XY(KapGMtave,intertimeave,bi,bj,myThid)
        ENDDO
       ENDDO

C      Write to files
       _BARRIER
       _BEGIN_MASTER( myThid )

C--     Set IO "context" for writing state
        CALL DFILE_SET_RW
        CALL DFILE_SET_CONT_ON_ERROR
C--     Read IO error counter
        beginIOErrCount = IO_ERRCOUNT(myThid)
        prevPrec = writeBinaryPrec
        writeBinaryPrec = precFloat32
        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('UTtave.',suff,uttave,myIter,myThid)
        CALL WRITE_FLD_XYZ_RL('VTtave.',suff,vttave,myIter,myThid)
        CALL WRITE_FLD_XYZ_RL('K13tave.',suff,k13tave,myIter,myThid)
        CALL WRITE_FLD_XYZ_RL('K23tave.',suff,k23tave,myIter,myThid)
        CALL WRITE_FLD_XYZ_RL('WTtave.',suff,wVeltave,myIter,myThid)
        CALL WRITE_FLD_XY_RL('KapGMtave.',suff,KapGMtave,myIter,myThid)
        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)
         DO k=1,Nz

C      Reset the averages to zero
          CALL TIMEAVER_INI_XYZ(thetatave, bi,bj,K,myThid)
          CALL TIMEAVER_INI_XYZ(salttave, bi,bj,K,myThid)
          CALL TIMEAVER_INI_XYZ(uVeltave, bi,bj,K,myThid)
          CALL TIMEAVER_INI_XYZ(vVeltave, bi,bj,K,myThid)
          CALL TIMEAVER_INI_XYZ(uttave, bi,bj,K,myThid)
          CALL TIMEAVER_INI_XYZ(vttave, bi,bj,K,myThid)
          CALL TIMEAVER_INI_XYZ(k13tave, bi,bj,K,myThid)
          CALL TIMEAVER_INI_XYZ(k23tave, bi,bj,K,myThid)
          CALL TIMEAVER_INI_XYZ(wVeltave, bi,bj,K,myThid)
          IF (K.EQ.1) CALL TIMEAVER_INI_XY(KapGMtave, bi,bj,myThid)
          TimeAve(k,bi,bj)=0.
          InterTimeAve(k,bi,bj)=0.

C      Time Averages of single fields
          DDTT=0.5*deltaTclock
          CALL TIMEAVER_1FLD_XYZ(theta, thetatave, DDTT, bi,bj,K,myThid)
          CALL TIMEAVER_1FLD_XYZ(salt,  salttave,  DDTT, bi,bj,K,myThid)
          CALL TIMEAVER_1FLD_XYZ(uVel,  uVeltave,  DDTT, bi,bj,K,myThid)
          CALL TIMEAVER_1FLD_XYZ(vVel,  vVeltave,  DDTT, bi,bj,K,myThid)

C      Time Averages of "double" fields
          CALL TIMEAVER_2FLD_XYZ(theta, uVel, uttave, 1, DDTT, bi,bj,K,myThid)
          CALL TIMEAVER_2FLD_XYZ(theta, vVel, vttave, 2, DDTT, bi,bj,K,myThid)

C      Keep record of how much time has been integrated over
          TimeAve(k,bi,bj)=TimeAve(k,bi,bj)+DDTT

         ENDDO
        ENDDO
       ENDDO

      ENDIF
 
      RETURN
      END
1	adcroft	1.1	#include "CPP_EEOPTIONS.h"
2
3			SUBROUTINE WRITE_TIME_AVERAGES(myCurrentTime, myIter, myThid)
4			C /==========================================================\
5			C \| SUBROUTINE WRITE_TIME_AVERAGES \|
6			C \| \|
7			C \==========================================================/
8
9			C == Global variables ===
10			#include "SIZE.h"
11			#include "EEPARAMS.h"
12			#include "PARAMS.h"
13			#include "AVER.h"
14			#include "DYNVARS.h"
15
16			LOGICAL DIFFERENT_MULTIPLE
17			EXTERNAL DIFFERENT_MULTIPLE
18	adcroft	1.2	INTEGER IO_ERRCOUNT
19			EXTERNAL IO_ERRCOUNT
20	adcroft	1.1
21			C == Routine arguments ==
22			C myThid - Thread number for this instance of the routine.
23			C myIter - Iteration number
24			C myCurrentTime - Current time of simulation ( s )
25			INTEGER myThid
26			INTEGER myIter
27			REAL myCurrentTime
28			INTEGER prevPrec
29
30			C == Local variables ==
31			C suff - Hold suffix part of a filename
32			C TimeAve - total time over average
33			CHARACTER*(MAX_LEN_FNAM) suff
34			REAL DDTT
35			INTEGER bi,bj,k
36	adcroft	1.2	INTEGER beginIOErrCount
37			INTEGER endIOErrCount
38			CHARACTER*(MAX_LEN_MBUF) msgBuf
39	adcroft	1.1
40			C Final Time Averages and Dump Files if needed
41			IF ( myIter.NE.nIter0 .AND.
42			& DIFFERENT_MULTIPLE(taveFreq,myCurrentTime,myCurrentTime-deltaTClock)
43			& ) THEN
44
45			WRITE(suff,'(I10.10)') myIter
46
47			C Normalize by integrated time
48			DO bj = myByLo(myThid), myByHi(myThid)
49			DO bi = myBxLo(myThid), myBxHi(myThid)
50			DO k=1,Nz
51			CALL TIMEAVER_FINAL_XYZ(thetatave,timeave,bi,bj,K,myThid)
52			CALL TIMEAVER_FINAL_XYZ(salttave,timeave,bi,bj,K,myThid)
53			CALL TIMEAVER_FINAL_XYZ(uVeltave,timeave,bi,bj,K,myThid)
54			CALL TIMEAVER_FINAL_XYZ(vVeltave,timeave,bi,bj,K,myThid)
55			CALL TIMEAVER_FINAL_XYZ(uttave,timeave,bi,bj,K,myThid)
56			CALL TIMEAVER_FINAL_XYZ(vttave,timeave,bi,bj,K,myThid)
57			CALL TIMEAVER_FINAL_XYZ(k13tave,intertimeave,bi,bj,K,myThid)
58			CALL TIMEAVER_FINAL_XYZ(k23tave,intertimeave,bi,bj,K,myThid)
59			CALL TIMEAVER_FINAL_XYZ(wVeltave,intertimeave,bi,bj,K,myThid)
60			ENDDO
61	adcroft	1.3	CALL TIMEAVER_FINAL_XY(KapGMtave,intertimeave,bi,bj,myThid)
62	adcroft	1.1	ENDDO
63			ENDDO
64
65			C Write to files
66			_BARRIER
67			_BEGIN_MASTER( myThid )
68	adcroft	1.2
69			C-- Set IO "context" for writing state
70			CALL DFILE_SET_RW
71			CALL DFILE_SET_CONT_ON_ERROR
72			C-- Read IO error counter
73			beginIOErrCount = IO_ERRCOUNT(myThid)
74	adcroft	1.1	prevPrec = writeBinaryPrec
75			writeBinaryPrec = precFloat32
76			CALL WRITE_FLD_XYZ_RL('Ttave.',suff,thetatave,myIter,myThid)
77			CALL WRITE_FLD_XYZ_RL('Stave.',suff,salttave,myIter,myThid)
78			CALL WRITE_FLD_XYZ_RL('uVeltave.',suff,uVeltave,myIter,myThid)
79			CALL WRITE_FLD_XYZ_RL('vVeltave.',suff,vVeltave,myIter,myThid)
80			CALL WRITE_FLD_XYZ_RL('UTtave.',suff,uttave,myIter,myThid)
81			CALL WRITE_FLD_XYZ_RL('VTtave.',suff,vttave,myIter,myThid)
82			CALL WRITE_FLD_XYZ_RL('K13tave.',suff,k13tave,myIter,myThid)
83			CALL WRITE_FLD_XYZ_RL('K23tave.',suff,k23tave,myIter,myThid)
84			CALL WRITE_FLD_XYZ_RL('WTtave.',suff,wVeltave,myIter,myThid)
85	adcroft	1.3	CALL WRITE_FLD_XY_RL('KapGMtave.',suff,KapGMtave,myIter,myThid)
86	adcroft	1.1	writeBinaryPrec = prevPrec
87	adcroft	1.2	C-- Reread IO error counter
88			endIOErrCount = IO_ERRCOUNT(myThid)
89
90			C-- Check for IO errors
91			IF ( endIOErrCount .NE. beginIOErrCount ) THEN
92			WRITE(msgBuf,'(A)') 'S/R WRITE_TIME_AVERAGES'
93			CALL PRINT_ERROR( msgBuf, 1 )
94			WRITE(msgBuf,'(A)') 'Error writing out data'
95			CALL PRINT_ERROR( msgBuf, 1 )
96			WRITE(msgBuf,'(A,I10)') 'Timestep ',myIter
97			CALL PRINT_ERROR( msgBuf, 1 )
98			ELSE
99			WRITE(msgBuf,'(A,I10)') '// Time-average data written, t-step', myIter
100			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT, 1 )
101			WRITE(msgBuf,'(A)') ' '
102			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT, 1 )
103			ENDIF
104
105	adcroft	1.1	_END_MASTER( myThid )
106			_BARRIER
107
108			DO bj = myByLo(myThid), myByHi(myThid)
109			DO bi = myBxLo(myThid), myBxHi(myThid)
110			DO k=1,Nz
111
112			C Reset the averages to zero
113			CALL TIMEAVER_INI_XYZ(thetatave, bi,bj,K,myThid)
114			CALL TIMEAVER_INI_XYZ(salttave, bi,bj,K,myThid)
115			CALL TIMEAVER_INI_XYZ(uVeltave, bi,bj,K,myThid)
116			CALL TIMEAVER_INI_XYZ(vVeltave, bi,bj,K,myThid)
117			CALL TIMEAVER_INI_XYZ(uttave, bi,bj,K,myThid)
118			CALL TIMEAVER_INI_XYZ(vttave, bi,bj,K,myThid)
119			CALL TIMEAVER_INI_XYZ(k13tave, bi,bj,K,myThid)
120			CALL TIMEAVER_INI_XYZ(k23tave, bi,bj,K,myThid)
121			CALL TIMEAVER_INI_XYZ(wVeltave, bi,bj,K,myThid)
122	adcroft	1.3	IF (K.EQ.1) CALL TIMEAVER_INI_XY(KapGMtave, bi,bj,myThid)
123	adcroft	1.1	TimeAve(k,bi,bj)=0.
124			InterTimeAve(k,bi,bj)=0.
125
126			C Time Averages of single fields
127			DDTT=0.5*deltaTclock
128			CALL TIMEAVER_1FLD_XYZ(theta, thetatave, DDTT, bi,bj,K,myThid)
129			CALL TIMEAVER_1FLD_XYZ(salt, salttave, DDTT, bi,bj,K,myThid)
130			CALL TIMEAVER_1FLD_XYZ(uVel, uVeltave, DDTT, bi,bj,K,myThid)
131			CALL TIMEAVER_1FLD_XYZ(vVel, vVeltave, DDTT, bi,bj,K,myThid)
132
133			C Time Averages of "double" fields
134			CALL TIMEAVER_2FLD_XYZ(theta, uVel, uttave, 1, DDTT, bi,bj,K,myThid)
135			CALL TIMEAVER_2FLD_XYZ(theta, vVel, vttave, 2, DDTT, bi,bj,K,myThid)
136
137			C Keep record of how much time has been integrated over
138			TimeAve(k,bi,bj)=TimeAve(k,bi,bj)+DDTT
139
140			ENDDO
141			ENDDO
142			ENDDO
143
144			ENDIF
145
146			RETURN
147			END