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

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

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

Revision 1.3 - (show annotations) (download)
Thu Jan 3 16:25:44 2002 UTC (22 years, 4 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint44e_post, checkpoint46g_pre, checkpoint46f_post, checkpoint44f_post, checkpoint46b_post, checkpoint43a-release1mods, chkpt44d_post, checkpoint44e_pre, checkpoint46d_pre, release1-branch_tutorials, checkpoint45d_post, chkpt44a_post, checkpoint44h_pre, checkpoint46a_post, chkpt44c_pre, checkpoint45a_post, checkpoint44g_post, checkpoint46e_pre, checkpoint45b_post, checkpoint46b_pre, release1-branch-end, release1_final_v1, checkpoint46c_pre, checkpoint46, checkpoint44b_post, checkpoint46a_pre, checkpoint45c_post, checkpoint44h_post, chkpt44a_pre, checkpoint46c_post, checkpoint46e_post, checkpoint44b_pre, checkpoint44, checkpoint45, chkpt44c_post, checkpoint44f_pre, checkpoint46d_post, release1-branch_branchpoint
Branch point for: release1_final, release1-branch
Changes since 1.2: +69 -79 lines 
change heat transport diagnostic (UTtave & VTtave):
 hFactor is now included in UTtave and VTtave output file and
 can be directly used to compute heat transport.
add new time average diagnostics :
 * Eta**2
 * wVel*Theta
 * volume transport hFac*U, hFac*V (only with NONLIN_FRSURF).

Note: since wVel and Theta are not updated synchronously in the model,
 WTtave output must be used with some cautious.
1 C $Header: /u/gcmpack/models/MITgcmUV/pkg/timeave/timeave_statv_write.F,v 1.2 2001/05/29 14:01:40 adcroft Exp $
2 C $Name: $
3 #include "CPP_OPTIONS.h"
4
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 #include "GRID.h"
19 #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 #ifdef ALLOW_TIMEAVE
35
36 C == Local variables ==
37 C suff - Hold suffix part of a filename
38 C TimeAve - total time over average
39 INTEGER prevPrec
40 CHARACTER*(MAX_LEN_FNAM) suff
41 _RL DDTT
42 INTEGER bi,bj,k
43 INTEGER beginIOErrCount
44 INTEGER endIOErrCount
45 CHARACTER*(MAX_LEN_MBUF) msgBuf
46
47 C- Final Time Averages and Dump Files if needed
48 c IF ( myIter.NE.nIter0 .AND.
49 c & DIFFERENT_MULTIPLE(taveFreq,myTime,myTime-deltaTClock) ) THEN
50 IF(DIFFERENT_MULTIPLE(taveFreq,myTime,myTime-deltaTClock) )THEN
51
52 WRITE(suff,'(I10.10)') myIter
53
54 DO bj = myByLo(myThid), myByHi(myThid)
55 DO bi = myBxLo(myThid), myBxHi(myThid)
56 C- Normalize by integrated time
57 CALL TIMEAVE_NORMALIZ(etatave, TimeAve_half,1 ,bi,bj,myThid)
58 CALL TIMEAVE_NORMALIZ(thetatave,TimeAve_half,Nr,bi,bj,myThid)
59 CALL TIMEAVE_NORMALIZ(salttave, TimeAve_half,Nr,bi,bj,myThid)
60 CALL TIMEAVE_NORMALIZ(uVeltave, TimeAve_half,Nr,bi,bj,myThid)
61 CALL TIMEAVE_NORMALIZ(vVeltave, TimeAve_half,Nr,bi,bj,myThid)
62 CALL TIMEAVE_NORMALIZ(Eta2tave, TimeAve_half,1 ,bi,bj,myThid)
63 CALL TIMEAVE_NORMALIZ(TTtave, TimeAve_half,Nr,bi,bj,myThid)
64 CALL TIMEAVE_NORMALIZ(UUtave, TimeAve_half,Nr,bi,bj,myThid)
65 CALL TIMEAVE_NORMALIZ(VVtave, TimeAve_half,Nr,bi,bj,myThid)
66 c CALL TIMEAVE_NORMALIZ(KEtave, TimeAve_half,Nr,bi,bj,myThid)
67 #ifdef NONLIN_FRSURF
68 C- Normalize by integrated time
69 CALL TIMEAVE_NORMALIZ(hUtave, TimeAve_half,Nr,bi,bj,myThid)
70 CALL TIMEAVE_NORMALIZ(hVtave, TimeAve_half,Nr,bi,bj,myThid)
71 c CALL TIMEAVE_NORMALIZ(hFacCtave,TimeAve_half,Nr,bi,bj,myThid)
72 c CALL TIMEAVE_NORMALIZ(hFacWtave,TimeAve_half,Nr,bi,bj,myThid)
73 c CALL TIMEAVE_NORMALIZ(hFacStave,TimeAve_half,Nr,bi,bj,myThid)
74 #endif /* NONLIN_FRSURF */
75 CALL TIMEAVE_NORMALIZ(UTtave, TimeAve_half,Nr,bi,bj,myThid)
76 CALL TIMEAVE_NORMALIZ(VTtave, TimeAve_half,Nr,bi,bj,myThid)
77
78 CALL TIMEAVE_NORMALIZ(wVeltave, TimeAve_full,Nr,bi,bj,myThid)
79 CALL TIMEAVE_NORMALIZ(WTtave, TimeAve_full,Nr,bi,bj,myThid)
80 CALL TIMEAVE_NORMALIZ(phiHydtave,TimeAve_full,Nr,
81 & bi,bj,myThid)
82 CALL TIMEAVE_NORMALIZ(ConvectCountTave,TimeAve_full,Nr,
83 & bi,bj,myThid)
84 ENDDO
85 ENDDO
86
87 C Write to files
88 _BARRIER
89 _BEGIN_MASTER( myThid )
90
91 C-- Set IO "context" for writing state
92 #ifdef USE_DFILE
93 CALL DFILE_SET_RW
94 CALL DFILE_SET_CONT_ON_ERROR
95 #endif
96 C-- Read IO error counter
97 beginIOErrCount = IO_ERRCOUNT(myThid)
98 c prevPrec = writeBinaryPrec
99 c writeBinaryPrec = precFloat32
100 CALL WRITE_FLD_XY_RL('ETAtave.',suff,etatave,myIter,myThid)
101 CALL WRITE_FLD_XY_RL('Eta2tave.',suff,Eta2tave,myIter,myThid)
102
103 CALL WRITE_FLD_XYZ_RL('Ttave.',suff,thetatave,myIter,myThid)
104 CALL WRITE_FLD_XYZ_RL('Stave.',suff,salttave,myIter,myThid)
105 CALL WRITE_FLD_XYZ_RL('uVeltave.',suff,uVeltave,myIter,myThid)
106 CALL WRITE_FLD_XYZ_RL('vVeltave.',suff,vVeltave,myIter,myThid)
107 CALL WRITE_FLD_XYZ_RL('wVeltave.',suff,wVeltave,myIter,myThid)
108 CALL WRITE_FLD_XYZ_RL('TTtave.',suff,TTtave,myIter,myThid)
109 CALL WRITE_FLD_XYZ_RL('UUtave.',suff,UUtave,myIter,myThid)
110 CALL WRITE_FLD_XYZ_RL('VVtave.',suff,VVtave,myIter,myThid)
111 c CALL WRITE_FLD_XYZ_RL('KEtave.',suff,KEtave,myIter,myThid)
112 CALL WRITE_FLD_XYZ_RL('UTtave.',suff,UTtave,myIter,myThid)
113 CALL WRITE_FLD_XYZ_RL('VTtave.',suff,VTtave,myIter,myThid)
114 CALL WRITE_FLD_XYZ_RL('WTtave.',suff,WTtave,myIter,myThid)
115 CALL WRITE_FLD_XYZ_RL('PhHytave.',suff,phiHydtave,myIter,myThid)
116 CALL WRITE_FLD_XYZ_RL('Convtave.',suff,ConvectCountTave,
117 & myIter,myThid)
118 #ifdef NONLIN_FRSURF
119 CALL WRITE_FLD_XYZ_RL('hUtave.',suff,hUtave,myIter,myThid)
120 CALL WRITE_FLD_XYZ_RL('hVtave.',suff,hVtave,myIter,myThid)
121 c CALL WRITE_FLD_XYZ_RL('hFacCtave.',suff,hFacCtave,myIter,myThid)
122 c CALL WRITE_FLD_XYZ_RL('hFacWtave.',suff,hFacWtave,myIter,myThid)
123 c CALL WRITE_FLD_XYZ_RL('hFacStave.',suff,hFacStave,myIter,myThid)
124 #endif /* NONLIN_FRSURF */
125
126 c writeBinaryPrec = prevPrec
127 C-- Reread IO error counter
128 endIOErrCount = IO_ERRCOUNT(myThid)
129
130 C-- Check for IO errors
131 IF ( endIOErrCount .NE. beginIOErrCount ) THEN
132 WRITE(msgBuf,'(A)') 'S/R WRITE_TIME_AVERAGES'
133 CALL PRINT_ERROR( msgBuf, 1 )
134 WRITE(msgBuf,'(A)') 'Error writing out data'
135 CALL PRINT_ERROR( msgBuf, 1 )
136 WRITE(msgBuf,'(A,I10)') 'Timestep ',myIter
137 CALL PRINT_ERROR( msgBuf, 1 )
138 ELSE
139 WRITE(msgBuf,'(A,I10)')
140 & '// Time-average data written, t-step', myIter
141 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
142 & SQUEEZE_RIGHT, 1 )
143 WRITE(msgBuf,'(A)') ' '
144 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
145 & SQUEEZE_RIGHT, 1 )
146 ENDIF
147
148 _END_MASTER( myThid )
149 _BARRIER
150
151 DO bj = myByLo(myThid), myByHi(myThid)
152 DO bi = myBxLo(myThid), myBxHi(myThid)
153
154 C-- Like before the 1rst iteration,
155 C ==> call TIMEAVE_STATVARS with myIter=nIter0 :
156 C 1) Reset the averages to zero ;
157 C 2) Start to cumulate state-variables with Half time step.
158
159 CALL TIMEAVE_STATVARS(myTime, nIter0, bi, bj, myThid)
160
161 ENDDO
162 ENDDO
163
164 ENDIF
165
166 #endif /* ALLOW_TIMEAVE */
167
168 RETURN
169 END
  ViewVC Help
Powered by ViewVC 1.1.22