/[MITgcm]/MITgcm/pkg/timeave/timeave_statv_write.F
ViewVC logotype

Diff of /MITgcm/pkg/timeave/timeave_statv_write.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.18 by edhill, Fri Sep 10 12:19:31 2004 UTC revision 1.38 by jmc, Sun Aug 6 23:30:33 2006 UTC
# Line 7  CBOP Line 7  CBOP
7  C     !ROUTINE: TIMEAVE_STATV_WRITE  C     !ROUTINE: TIMEAVE_STATV_WRITE
8    
9  C     !INTERFACE:  C     !INTERFACE:
10        SUBROUTINE TIMEAVE_STATV_WRITE(myTime, myIter, myThid)        SUBROUTINE TIMEAVE_STATV_WRITE( myTime, myIter, myThid )
11    
12  C     !DESCRIPTION:  C     !DESCRIPTION:
13  C     At the end of average period, write the time-average  C     At the end of average period, write the time-average
# Line 18  C     !USES: Line 18  C     !USES:
18  #include "SIZE.h"  #include "SIZE.h"
19  #include "EEPARAMS.h"  #include "EEPARAMS.h"
20  #include "PARAMS.h"  #include "PARAMS.h"
 #ifdef ALLOW_MNC          
 #include "MNC_PARAMS.h"  
 #endif  
21  #include "DYNVARS.h"  #include "DYNVARS.h"
22  #include "GRID.h"  #include "GRID.h"
23  #include "TIMEAVE_STATV.h"  #include "TIMEAVE_STATV.h"
24        LOGICAL  DIFFERENT_MULTIPLE  #ifdef ALLOW_MNC
25        EXTERNAL DIFFERENT_MULTIPLE  #include "MNC_PARAMS.h"
26        INTEGER  IO_ERRCOUNT  #endif
       EXTERNAL IO_ERRCOUNT  
27    
28  C     !INPUT PARAMETERS:  C     !INPUT PARAMETERS:
29  C     myThid - Thread number for this instance of the routine.  C     myThid - Thread number for this instance of the routine.
# Line 38  C     myTime - Current time of simulatio Line 34  C     myTime - Current time of simulatio
34        _RL     myTime        _RL     myTime
35  CEOP  CEOP
36    
37  #ifdef ALLOW_TIMEAVE  #ifdef ALLOW_TIMEAVE
38    C     Functions:
39          LOGICAL  DIFFERENT_MULTIPLE
40          EXTERNAL DIFFERENT_MULTIPLE
41          INTEGER  IO_ERRCOUNT
42          EXTERNAL IO_ERRCOUNT
43    
44  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
45  C     suff         :: Hold suffix part of a filename  C     suff         :: Hold suffix part of a filename
46  C     TimeAve      :: total time over average  C     TimeAve      :: total time over average
47  C     useVariableK :: T when vertical diffusion is not constant  C     useVariableK :: T when vertical diffusion is not constant
48        LOGICAL useVariableK        LOGICAL useVariableK
       INTEGER prevPrec  
49        CHARACTER*(MAX_LEN_FNAM) suff        CHARACTER*(MAX_LEN_FNAM) suff
50        INTEGER bi,bj,k        INTEGER bi,bj
51        INTEGER beginIOErrCount        INTEGER beginIOErrCount
52        INTEGER endIOErrCount        INTEGER endIOErrCount
53        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
54          LOGICAL dumpFiles
55    #ifdef ALLOW_MNC
56          CHARACTER*(1) pf
57    #endif
58    
59  C     Final Time Averages and Dump Files if needed  C     Final Time Averages and Dump Files if needed
60        IF (DIFFERENT_MULTIPLE(taveFreq,myTime,myTime-deltaTClock)) THEN        dumpFiles = DIFFERENT_MULTIPLE(taveFreq,myTime,deltaTClock)
61    #ifdef ALLOW_CAL
62          IF ( useCAL ) THEN
63             CALL CAL_TIME2DUMP( taveFreq, deltaTClock,
64         U                       dumpFiles,
65         I                       myTime, myIter, myThid )
66          ENDIF
67    #endif
68    
69          useVariableK = useKPP .OR. useGMredi .OR. ivdc_kappa.NE.0.        IF (dumpFiles) THEN
70            useVariableK = useKPP .OR. usePP81 .OR. useMY82 .OR. useGGL90
71         &      .OR. useGMredi .OR. ivdc_kappa.NE.0.
72          WRITE(suff,'(I10.10)') myIter          WRITE(suff,'(I10.10)') myIter
73    
74          DO bj = myByLo(myThid), myByHi(myThid)          DO bj = myByLo(myThid), myByHi(myThid)
75          DO bi = myBxLo(myThid), myBxHi(myThid)           DO bi = myBxLo(myThid), myBxHi(myThid)
76    
77  C         Normalize by integrated time  C         Normalize by integrated time
78            CALL TIMEAVE_NORMALIZ(uFluxtave,TimeAve_full,1 ,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(uFluxtave,TimeAve_full,1 ,bi,bj,myThid)
# Line 71  C         Normalize by integrated time Line 85  C         Normalize by integrated time
85            CALL TIMEAVE_NORMALIZ(uVeltave, TimeAve_half,Nr,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(uVeltave, TimeAve_half,Nr,bi,bj,myThid)
86            CALL TIMEAVE_NORMALIZ(vVeltave, TimeAve_half,Nr,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(vVeltave, TimeAve_half,Nr,bi,bj,myThid)
87            CALL TIMEAVE_NORMALIZ(wVeltave, TimeAve_half,Nr,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(wVeltave, TimeAve_half,Nr,bi,bj,myThid)
88              CALL TIMEAVE_NORMALIZ(phiHydLowtave,TimeAve_full,1,
89         &         bi,bj,myThid)
90            CALL TIMEAVE_NORMALIZ(UTtave,   TimeAve_half,Nr,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(UTtave,   TimeAve_half,Nr,bi,bj,myThid)
91            CALL TIMEAVE_NORMALIZ(VTtave,   TimeAve_half,Nr,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(VTtave,   TimeAve_half,Nr,bi,bj,myThid)
92            CALL TIMEAVE_NORMALIZ(WTtave,   TimeAve_half,Nr,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(WTtave,   TimeAve_half,Nr,bi,bj,myThid)
93            CALL TIMEAVE_NORMALIZ(UStave,   TimeAve_half,Nr,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(UStave,   TimeAve_half,Nr,bi,bj,myThid)
94            CALL TIMEAVE_NORMALIZ(VStave,   TimeAve_half,Nr,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(VStave,   TimeAve_half,Nr,bi,bj,myThid)
95            CALL TIMEAVE_NORMALIZ(WStave,   TimeAve_half,Nr,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(WStave,   TimeAve_half,Nr,bi,bj,myThid)
 #ifndef HRCUBE  
96            CALL TIMEAVE_NORMALIZ(Eta2tave, TimeAve_half,1 ,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(Eta2tave, TimeAve_half,1 ,bi,bj,myThid)
97            CALL TIMEAVE_NORMALIZ(TTtave,   TimeAve_half,Nr,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(TTtave,   TimeAve_half,Nr,bi,bj,myThid)
98            CALL TIMEAVE_NORMALIZ(UUtave,   TimeAve_half,Nr,bi,bj,myThid)            CALL TIMEAVE_NORMALIZ(UUtave,   TimeAve_half,Nr,bi,bj,myThid)
# Line 91  C         Normalize by integrated time Line 106  C         Normalize by integrated time
106  C         CALL TIMEAVE_NORMALIZ(hFacCtave,TimeAve_half,Nr,bi,bj,myThid)  C         CALL TIMEAVE_NORMALIZ(hFacCtave,TimeAve_half,Nr,bi,bj,myThid)
107  C         CALL TIMEAVE_NORMALIZ(hFacWtave,TimeAve_half,Nr,bi,bj,myThid)  C         CALL TIMEAVE_NORMALIZ(hFacWtave,TimeAve_half,Nr,bi,bj,myThid)
108  C         CALL TIMEAVE_NORMALIZ(hFacStave,TimeAve_half,Nr,bi,bj,myThid)  C         CALL TIMEAVE_NORMALIZ(hFacStave,TimeAve_half,Nr,bi,bj,myThid)
109  #endif /* NONLIN_FRSURF */  #endif /* NONLIN_FRSURF */
110    
111            CALL TIMEAVE_NORMALIZ(TdiffRtave,TimeAve_full,Nr,            CALL TIMEAVE_NORMALIZ(TdiffRtave,TimeAve_full,Nr,
112       &         bi,bj,myThid)       &         bi,bj,myThid)
# Line 101  C         CALL TIMEAVE_NORMALIZ(hFacStav Line 116  C         CALL TIMEAVE_NORMALIZ(hFacStav
116  #endif  #endif
117            CALL TIMEAVE_NORMALIZ(phiHydtave,TimeAve_full,Nr,            CALL TIMEAVE_NORMALIZ(phiHydtave,TimeAve_full,Nr,
118       &         bi,bj,myThid)       &         bi,bj,myThid)
           CALL TIMEAVE_NORMALIZ(phiHydLowtave,TimeAve_full,1,  
      &         bi,bj,myThid)  
119            CALL TIMEAVE_NORMALIZ(phiHydLow2Tave,TimeAve_full,1,            CALL TIMEAVE_NORMALIZ(phiHydLow2Tave,TimeAve_full,1,
120       &         bi,bj,myThid)       &         bi,bj,myThid)
121            CALL TIMEAVE_NORMALIZ(ConvectCountTave,TimeAve_full,Nr,            CALL TIMEAVE_NORMALIZ(ConvectCountTave,TimeAve_full,Nr,
122       &         bi,bj,myThid)       &         bi,bj,myThid)
123  #endif /* ndef HRCUBE */           ENDDO
         ENDDO  
124          ENDDO          ENDDO
125    
126  C       Write to files  C       Write to files
127          _BARRIER          _BARRIER
         _BEGIN_MASTER( myThid )  
128    
129          IF (timeave_mdsio) THEN          IF (timeave_mdsio) THEN
130              
131  C         Set IO "context" for writing state  C         Set IO "context" for writing state
132  #ifdef USE_DFILE  #ifdef USE_DFILE
133            CALL DFILE_SET_RW            CALL DFILE_SET_RW
# Line 124  C         Set IO "context" for writing s Line 135  C         Set IO "context" for writing s
135  #endif  #endif
136  C         Read IO error counter  C         Read IO error counter
137            beginIOErrCount = IO_ERRCOUNT(myThid)            beginIOErrCount = IO_ERRCOUNT(myThid)
 c         prevPrec = writeBinaryPrec  
 c         writeBinaryPrec = precFloat32  
138            CALL WRITE_FLD_XY_RL('ETAtave.' ,suff,etatave  ,myIter,myThid)            CALL WRITE_FLD_XY_RL('ETAtave.' ,suff,etatave  ,myIter,myThid)
 #ifndef HRCUBE  
           CALL WRITE_FLD_XY_RL('Eta2tave.',suff,Eta2tave ,myIter,myThid)  
 #endif /* ndef HRCUBE */  
   
139            CALL WRITE_FLD_XYZ_RL('Ttave.',suff,thetatave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('Ttave.',suff,thetatave,myIter,myThid)
140            CALL WRITE_FLD_XYZ_RL('Stave.',suff,salttave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('Stave.',suff,salttave,myIter,myThid)
141            CALL WRITE_FLD_XYZ_RL('uVeltave.',suff,uVeltave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('uVeltave.',suff,uVeltave,myIter,myThid)
142            CALL WRITE_FLD_XYZ_RL('vVeltave.',suff,vVeltave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('vVeltave.',suff,vVeltave,myIter,myThid)
143            CALL WRITE_FLD_XYZ_RL('wVeltave.',suff,wVeltave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('wVeltave.',suff,wVeltave,myIter,myThid)
144              CALL WRITE_FLD_XY_RL('Eta2tave.',suff,Eta2tave ,myIter,myThid)
145            CALL WRITE_FLD_XYZ_RL('UTtave.',suff,UTtave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('UTtave.',suff,UTtave,myIter,myThid)
146            CALL WRITE_FLD_XYZ_RL('VTtave.',suff,VTtave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('VTtave.',suff,VTtave,myIter,myThid)
147            CALL WRITE_FLD_XYZ_RL('WTtave.',suff,WTtave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('WTtave.',suff,WTtave,myIter,myThid)
148            CALL WRITE_FLD_XYZ_RL('UStave.',suff,UStave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('UStave.',suff,UStave,myIter,myThid)
149            CALL WRITE_FLD_XYZ_RL('VStave.',suff,VStave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('VStave.',suff,VStave,myIter,myThid)
150            CALL WRITE_FLD_XYZ_RL('WStave.',suff,WStave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('WStave.',suff,WStave,myIter,myThid)
 #ifndef HRCUBE  
151            CALL WRITE_FLD_XYZ_RL('TTtave.',suff,TTtave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('TTtave.',suff,TTtave,myIter,myThid)
152            CALL WRITE_FLD_XYZ_RL('UUtave.',suff,UUtave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('UUtave.',suff,UUtave,myIter,myThid)
153            CALL WRITE_FLD_XYZ_RL('VVtave.',suff,VVtave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('VVtave.',suff,VVtave,myIter,myThid)
154            CALL WRITE_FLD_XYZ_RL('UVtave.',suff,UVtave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('UVtave.',suff,UVtave,myIter,myThid)
155  C         CALL WRITE_FLD_XYZ_RL('KEtave.',suff,KEtave,myIter,myThid)  C         CALL WRITE_FLD_XYZ_RL('KEtave.',suff,KEtave,myIter,myThid)
156              
157            IF (useVariableK)            IF (useVariableK)
158       &         CALL WRITE_FLD_XYZ_RL('Tdiftave.',suff,TdiffRtave,       &         CALL WRITE_FLD_XYZ_RL('Tdiftave.',suff,TdiffRtave,
159       &         myIter,myThid)       &         myIter,myThid)
# Line 162  C         CALL WRITE_FLD_XYZ_RL('KEtave. Line 167  C         CALL WRITE_FLD_XYZ_RL('KEtave.
167  #endif /* ALLOW_MOM_VECINV */  #endif /* ALLOW_MOM_VECINV */
168            CALL WRITE_FLD_XYZ_RL('PhHytave.',suff,phiHydtave,            CALL WRITE_FLD_XYZ_RL('PhHytave.',suff,phiHydtave,
169       &         myIter,myThid)       &         myIter,myThid)
170  #endif /* ndef HRCUBE */            IF ( fluidIsWater ) THEN
           IF ( buoyancyRelation(1:7) .EQ. 'OCEANIC' ) THEN  
 #ifndef HRCUBE  
171              CALL WRITE_FLD_XY_RL('PHLtave.',suff,phiHydLowtave,              CALL WRITE_FLD_XY_RL('PHLtave.',suff,phiHydLowtave,
172       &           myIter,myThid)       &           myIter,myThid)
             CALL WRITE_FLD_XY_RL('PHL2tave.',suff,phiHydLow2tave,  
      &           myIter,myThid)  
             CALL WRITE_FLD_XYZ_RL('Convtave.',suff,ConvectCountTave,  
      &           myIter,myThid)  
 #endif /* ndef HRCUBE */  
173              CALL WRITE_FLD_XY_RL('uFluxtave.',suff,uFluxtave,              CALL WRITE_FLD_XY_RL('uFluxtave.',suff,uFluxtave,
174       &           myIter,myThid)       &           myIter,myThid)
175              CALL WRITE_FLD_XY_RL('vFluxtave.',suff,vFluxtave,              CALL WRITE_FLD_XY_RL('vFluxtave.',suff,vFluxtave,
# Line 180  C         CALL WRITE_FLD_XYZ_RL('KEtave. Line 178  C         CALL WRITE_FLD_XYZ_RL('KEtave.
178       &           myIter,myThid)       &           myIter,myThid)
179              CALL WRITE_FLD_XY_RL('sFluxtave.',suff,sFluxtave,              CALL WRITE_FLD_XY_RL('sFluxtave.',suff,sFluxtave,
180       &           myIter,myThid)       &           myIter,myThid)
181                CALL WRITE_FLD_XY_RL('PHL2tave.',suff,phiHydLow2tave,
182         &           myIter,myThid)
183                CALL WRITE_FLD_XYZ_RL('Convtave.',suff,ConvectCountTave,
184         &           myIter,myThid)
185            ENDIF            ENDIF
186              
 #ifndef HRCUBE  
187  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
188            CALL WRITE_FLD_XYZ_RL('hUtave.',suff,hUtave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('hUtave.',suff,hUtave,myIter,myThid)
189            CALL WRITE_FLD_XYZ_RL('hVtave.',suff,hVtave,myIter,myThid)            CALL WRITE_FLD_XYZ_RL('hVtave.',suff,hVtave,myIter,myThid)
190  C         CALL WRITE_FLD_XYZ_RL('hFacCtave.',suff,hFacCtave,myIter,myThid)  C         CALL WRITE_FLD_XYZ_RL('hFacCtave.',suff,hFacCtave,myIter,myThid)
191  C         CALL WRITE_FLD_XYZ_RL('hFacWtave.',suff,hFacWtave,myIter,myThid)  C         CALL WRITE_FLD_XYZ_RL('hFacWtave.',suff,hFacWtave,myIter,myThid)
192  C         CALL WRITE_FLD_XYZ_RL('hFacStave.',suff,hFacStave,myIter,myThid)  C         CALL WRITE_FLD_XYZ_RL('hFacStave.',suff,hFacStave,myIter,myThid)
193  #endif /* NONLIN_FRSURF */  #endif /* NONLIN_FRSURF */
194  #endif /* ndef HRCUBE */  
             
 C         writeBinaryPrec = prevPrec  
195  C         Reread IO error counter  C         Reread IO error counter
196            endIOErrCount = IO_ERRCOUNT(myThid)            endIOErrCount = IO_ERRCOUNT(myThid)
197    
198  C         Check for IO errors  C         Check for IO errors
199            IF ( endIOErrCount .NE. beginIOErrCount ) THEN            IF ( endIOErrCount .NE. beginIOErrCount ) THEN
200    C-          any thread that detects an error should report
201              WRITE(msgBuf,'(A)')  'S/R WRITE_TIME_AVERAGES'              WRITE(msgBuf,'(A)')  'S/R WRITE_TIME_AVERAGES'
202              CALL PRINT_ERROR( msgBuf, 1 )              CALL PRINT_ERROR( msgBuf, myThid )
203              WRITE(msgBuf,'(A)')  'Error writing out data'              WRITE(msgBuf,'(A)')  'Error writing out data'
204              CALL PRINT_ERROR( msgBuf, 1 )              CALL PRINT_ERROR( msgBuf, myThid )
205              WRITE(msgBuf,'(A,I10)') 'Timestep ',myIter              WRITE(msgBuf,'(A,I10)') 'Timestep ',myIter
206              CALL PRINT_ERROR( msgBuf, 1 )              CALL PRINT_ERROR( msgBuf, myThid )
207            ELSE            ELSE
208              WRITE(msgBuf,'(A,I10)')    C-          normal case: 1 message is enough
209                _BEGIN_MASTER( myThid )
210                WRITE(msgBuf,'(A,I10)')
211       &           '// Time-average data written, t-step', myIter       &           '// Time-average data written, t-step', myIter
212              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
213       &           SQUEEZE_RIGHT, 1 )       &           SQUEEZE_RIGHT, myThid )
214              WRITE(msgBuf,'(A)')  ' '              WRITE(msgBuf,'(A)')  ' '
215              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,              CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
216       &           SQUEEZE_RIGHT, 1 )       &           SQUEEZE_RIGHT, myThid )
217                _END_MASTER( myThid )
218            ENDIF            ENDIF
219              
220          ENDIF          ENDIF
221    
222  #ifdef ALLOW_MNC          #ifdef ALLOW_MNC
223          IF (useMNC .AND. timeave_mdsio) THEN          IF (useMNC .AND. timeave_mnc) THEN
224              
225              IF ( writeBinaryPrec .EQ. precFloat64 ) THEN
226                pf(1:1) = 'D'
227              ELSE
228                pf(1:1) = 'R'
229              ENDIF
230    
231            CALL MNC_CW_SET_UDIM('tave', -1, myThid)            CALL MNC_CW_SET_UDIM('tave', -1, myThid)
232            CALL MNC_CW_I_W('I','tave',0,0,'iter', myIter, myThid)            CALL MNC_CW_RL_W_S('D','tave',0,0,'T', myTime, myThid)
233            CALL MNC_CW_SET_UDIM('tave', 0, myThid)            CALL MNC_CW_SET_UDIM('tave', 0, myThid)
234            CALL MNC_CW_RL_W('D','tave',0,0,'model_time',myTime,myThid)            CALL MNC_CW_I_W_S('I','tave',0,0,'iter', myIter, myThid)
235    C         CALL MNC_CW_RL_W_S('D','tave',0,0,'model_time',myTime,myThid)
236            CALL MNC_CW_RL_W('D','tave',0,0,'ETAtave',etatave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'ETAtave',etatave,myThid)
237  #ifndef HRCUBE  C         EH3: something is wrong here...
238            CALL MNC_CW_RL_W('D','tave',0,0,'Eta2tave',Eta2tave,myThid)  C         CALL MNC_CW_RL_W(pf,'tave',0,0,'Eta2tave',Eta2tave,myThid)
239  #endif /* ndef HRCUBE */            CALL MNC_CW_RL_W(pf,'tave',0,0,'Ttave',thetatave,myThid)
240            CALL MNC_CW_RL_W('D','tave',0,0,'Ttave',thetatave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'Stave',salttave,myThid)
241            CALL MNC_CW_RL_W('D','tave',0,0,'Stave',salttave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'uVeltave',uVeltave,myThid)
242            CALL MNC_CW_RL_W('D','tave',0,0,'uVeltave',uVeltave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'vVeltave',vVeltave,myThid)
243            CALL MNC_CW_RL_W('D','tave',0,0,'vVeltave',vVeltave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'wVeltave',wVeltave,myThid)
244            CALL MNC_CW_RL_W('D','tave',0,0,'wVeltave',wVeltave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'UTtave',UTtave,myThid)
245            CALL MNC_CW_RL_W('D','tave',0,0,'UTtave',UTtave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'VTtave',VTtave,myThid)
246            CALL MNC_CW_RL_W('D','tave',0,0,'VTtave',VTtave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'WTtave',WTtave,myThid)
247            CALL MNC_CW_RL_W('D','tave',0,0,'WTtave',WTtave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'UStave',UStave,myThid)
248            CALL MNC_CW_RL_W('D','tave',0,0,'UStave',UStave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'VStave',VStave,myThid)
249            CALL MNC_CW_RL_W('D','tave',0,0,'VStave',VStave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'WStave',WStave,myThid)
250            CALL MNC_CW_RL_W('D','tave',0,0,'WStave',WStave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'TTtave',TTtave,myThid)
251  #ifndef HRCUBE            CALL MNC_CW_RL_W(pf,'tave',0,0,'UUtave',UUtave,myThid)
252            CALL MNC_CW_RL_W('D','tave',0,0,'TTtave',TTtave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'VVtave',VVtave,myThid)
253            CALL MNC_CW_RL_W('D','tave',0,0,'UUtave',UUtave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'UVtave',UVtave,myThid)
254            CALL MNC_CW_RL_W('D','tave',0,0,'VVtave',VVtave,myThid)  C         CALL MNC_CW_RL_W(pf,'tave',0,0,'KEtave',KEtave,myThid)
           CALL MNC_CW_RL_W('D','tave',0,0,'UVtave',UVtave,myThid)  
 C         CALL MNC_CW_RL_W('D','tave',0,0,'KEtave',KEtave,myThid)  
255            IF (useVariableK) THEN            IF (useVariableK) THEN
256              CALL MNC_CW_RL_W('D','tave',0,0,'Tdiftave',              CALL MNC_CW_RL_W(pf,'tave',0,0,'Tdiftave',
257       &           TdiffRtave,myThid)       &           TdiffRtave,myThid)
258            ENDIF            ENDIF
259  #ifdef ALLOW_MOM_VECINV  #ifdef ALLOW_MOM_VECINV
260            IF (vectorInvariantMomentum) THEN            IF (vectorInvariantMomentum) THEN
261              CALL MNC_CW_RL_W('D','tave',0,0,'uZtave',uZetatave,myThid)              CALL MNC_CW_RL_W(pf,'tave',0,0,'uZtave',uZetatave,myThid)
262              CALL MNC_CW_RL_W('D','tave',0,0,'vZtave',vZetatave,myThid)              CALL MNC_CW_RL_W(pf,'tave',0,0,'vZtave',vZetatave,myThid)
263            ENDIF            ENDIF
264  #endif /* ALLOW_MOM_VECINV */  #endif /* ALLOW_MOM_VECINV */
265            CALL MNC_CW_RL_W('D','tave',0,0,'PhHytave',            CALL MNC_CW_RL_W(pf,'tave',0,0,'PhHytave',
266       &         phiHydtave,myThid)       &         phiHydtave,myThid)
267  #endif /* ndef HRCUBE */            IF ( fluidIsWater ) THEN
268            IF ( buoyancyRelation(1:7) .EQ. 'OCEANIC' ) THEN              CALL MNC_CW_RL_W(pf,'tave',0,0,'PHLtave',
 #ifndef HRCUBE  
             CALL MNC_CW_RL_W('D','tave',0,0,'PHLtave',  
269       &           phiHydLowtave,myThid)       &           phiHydLowtave,myThid)
270              CALL MNC_CW_RL_W('D','tave',0,0,'PHL2tave',              CALL MNC_CW_RL_W(pf,'tave',0,0,'PHL2tave',
271       &           phiHydLow2tave,myThid)       &           phiHydLow2tave,myThid)
272              CALL MNC_CW_RL_W('D','tave',0,0,'Convtave',              CALL MNC_CW_RL_W(pf,'tave',0,0,'Convtave',
273       &           ConvectCountTave,myThid)       &           ConvectCountTave,myThid)
274  #endif /* ndef HRCUBE */              CALL MNC_CW_RL_W(pf,'tave',0,0,'uFluxtave',
             CALL MNC_CW_RL_W('D','tave',0,0,'uFluxtave',  
275       &           uFluxtave,myThid)       &           uFluxtave,myThid)
276              CALL MNC_CW_RL_W('D','tave',0,0,'vFluxtave',              CALL MNC_CW_RL_W(pf,'tave',0,0,'vFluxtave',
277       &           vFluxtave,myThid)       &           vFluxtave,myThid)
278              CALL MNC_CW_RL_W('D','tave',0,0,'tFluxtave',              CALL MNC_CW_RL_W(pf,'tave',0,0,'tFluxtave',
279       &           tFluxtave,myThid)       &           tFluxtave,myThid)
280              CALL MNC_CW_RL_W('D','tave',0,0,'sFluxtave',              CALL MNC_CW_RL_W(pf,'tave',0,0,'sFluxtave',
281       &           sFluxtave,myThid)       &           sFluxtave,myThid)
282            ENDIF            ENDIF
283    
 #ifndef HRCUBE  
284  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
285            CALL MNC_CW_RL_W('D','tave',0,0,'hUtave',hUtave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'hUtave',hUtave,myThid)
286            CALL MNC_CW_RL_W('D','tave',0,0,'hVtave',hVtave,myThid)            CALL MNC_CW_RL_W(pf,'tave',0,0,'hVtave',hVtave,myThid)
287  C         CALL MNC_CW_RL_W('D','tave',0,0,'hFacCtave',hFacCtave,myThid)  C         CALL MNC_CW_RL_W(pf,'tave',0,0,'hFacCtave',hFacCtave,myThid)
288  C         CALL MNC_CW_RL_W('D','tave',0,0,'hFacWtave',hFacWtave,myThid)  C         CALL MNC_CW_RL_W(pf,'tave',0,0,'hFacWtave',hFacWtave,myThid)
289  C         CALL MNC_CW_RL_W('D','tave',0,0,'hFacStave',hFacStave,myThid)  C         CALL MNC_CW_RL_W(pf,'tave',0,0,'hFacStave',hFacStave,myThid)
290  #endif /* NONLIN_FRSURF */  #endif /* NONLIN_FRSURF */
 #endif /* ndef HRCUBE */  
291    
292          ENDIF          ENDIF
293  #endif /* ALLOW_MNC */  #endif /* ALLOW_MNC */
294    
         _END_MASTER( myThid )  
295          _BARRIER          _BARRIER
296            
297          DO bj = myByLo(myThid), myByHi(myThid)          DO bj = myByLo(myThid), myByHi(myThid)
298            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
299                
300  C           Like before the 1rst iteration,  C           Like before the 1rst iteration,
301  C           ==> call TIMEAVE_STATVARS with myIter=nIter0 :  C           ==> call TIMEAVE_STATVARS with myIter=nIter0 :
302  C           1) Reset the averages to zero ;  C           1) Reset the averages to zero ;
303  C           2) Start to cumulate state-variables with Half time step.  C           2) Start to cumulate state-variables with Half time step.
304                
305              CALL TIMEAVE_STATVARS(myTime, nIter0, bi, bj, myThid)              CALL TIMEAVE_STATVARS(myTime, nIter0, bi, bj, myThid)
306                
307            ENDDO            ENDDO
308          ENDDO          ENDDO
309            
310        ENDIF        ENDIF
311          
312  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
313    
314        RETURN        RETURN
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.38
  ViewVC Help
Powered by ViewVC 1.1.22