/[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.32 by edhill, Thu Sep 29 12:19:52 2005 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
       INTEGER thisdate(4), prevdate(4)  
54        LOGICAL dumpFiles        LOGICAL dumpFiles
55    #ifdef ALLOW_MNC
56        CHARACTER*(1) pf        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        dumpFiles = .FALSE.        dumpFiles = DIFFERENT_MULTIPLE(taveFreq,myTime,deltaTClock)
       IF ( DIFFERENT_MULTIPLE(taveFreq,myTime,deltaTClock) )  
      &     dumpFiles = .TRUE.  
61  #ifdef ALLOW_CAL  #ifdef ALLOW_CAL
62        IF ( calendarDumps .AND. (        IF ( useCAL ) THEN
63       &       ( taveFreq.GE. 2592000 .AND. taveFreq.LE. 2678400 ) .OR.           CALL CAL_TIME2DUMP( taveFreq, deltaTClock,
64       &       ( taveFreq.GE.31104000 .AND. taveFreq.LE.31968000 ))) THEN       U                       dumpFiles,
65  C--   Convert approximate months (30-31 days) and years (360-372 days)       I                       myTime, myIter, myThid )
 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 )  
          dumpFiles = .FALSE.  
 C-    Monthly taveFreq:  
          IF( taveFreq.GE. 2592000 .AND. taveFreq.LE. 2678400 .AND.  
      &        (thisdate(1)-prevdate(1)).GT.50   ) dumpFiles = .TRUE.  
 C-    Yearly  taveFreq:  
          IF( taveFreq.GE.31104000 .AND. taveFreq.LE.31968000 .AND.  
      &        (thisdate(1)-prevdate(1)).GT.5000 ) dumpFiles = .TRUE.  
66        ENDIF        ENDIF
67  #endif  #endif
68    
# Line 83  C-    Yearly  taveFreq: Line 72  C-    Yearly  taveFreq:
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 117  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 131  C         CALL TIMEAVE_NORMALIZ(hFacStav Line 120  C         CALL TIMEAVE_NORMALIZ(hFacStav
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          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 147  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)
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)
# Line 167  c         writeBinaryPrec = precFloat32 Line 153  c         writeBinaryPrec = precFloat32
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 197  C         CALL WRITE_FLD_XYZ_RL('KEtave. Line 183  C         CALL WRITE_FLD_XYZ_RL('KEtave.
183              CALL WRITE_FLD_XYZ_RL('Convtave.',suff,ConvectCountTave,              CALL WRITE_FLD_XYZ_RL('Convtave.',suff,ConvectCountTave,
184       &           myIter,myThid)       &           myIter,myThid)
185            ENDIF            ENDIF
186              
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              
 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_mnc) THEN          IF (useMNC .AND. timeave_mnc) THEN
224              
225            IF ( writeBinaryPrec .EQ. precFloat64 ) THEN            IF ( writeBinaryPrec .EQ. precFloat64 ) THEN
226              pf(1:1) = 'D'              pf(1:1) = 'D'
227            ELSE            ELSE
# Line 298  C         CALL MNC_CW_RL_W(pf,'tave',0,0 Line 287  C         CALL MNC_CW_RL_W(pf,'tave',0,0
287  C         CALL MNC_CW_RL_W(pf,'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(pf,'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(pf,'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 */
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.32  
changed lines
  Added in v.1.38
  ViewVC Help
Powered by ViewVC 1.1.22