/[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.4 - (show annotations) (download)
Wed Sep 18 16:38:02 2002 UTC (21 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint46j_pre, checkpoint46h_pre, checkpoint46g_post, checkpoint46i_post, checkpoint46h_post
Changes since 1.3: +9 -1 lines 
o Include a new diagnostic variable phiHydLow for the ocean model
  - in z-coordinates, it is the bottom pressure anomaly
  - in p-coordinates, it is the sea surface elevation
  - in both cases, these variable have global drift, reflecting the mass
    drift in z-coordinates and the volume drift in p-coordinates
  - included time averaging for phiHydLow, be aware of the drift!
o depth-dependent computation of Bo_surf for pressure coordinates
  in the ocean (buoyancyRelation='OCEANICP')
  - requires a new routine (FIND_RHO_SCALAR) to compute density with only
    Theta, Salinity, and Pressure in the parameter list. This routine is
    presently contained in find_rho.F. This routine does not give the
    correct density for 'POLY3', which would be a z-dependent reference
    density.
o cleaned up find_rho
  - removed obsolete 'eqn' from the parameter list.
o added two new verification experiments: gop and goz
  (4x4 degree global ocean, 15 layers in pressure and height coordinates)
1 C $Header: /u/gcmpack/MITgcm/pkg/timeave/timeave_statv_write.F,v 1.3 2002/01/03 16:25:44 jmc 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(phiHydLowtave,TimeAve_full,1,
83 & bi,bj,myThid)
84 CALL TIMEAVE_NORMALIZ(phiHydLow2Tave,TimeAve_full,1,
85 & bi,bj,myThid)
86 CALL TIMEAVE_NORMALIZ(ConvectCountTave,TimeAve_full,Nr,
87 & bi,bj,myThid)
88 ENDDO
89 ENDDO
90
91 C Write to files
92 _BARRIER
93 _BEGIN_MASTER( myThid )
94
95 C-- Set IO "context" for writing state
96 #ifdef USE_DFILE
97 CALL DFILE_SET_RW
98 CALL DFILE_SET_CONT_ON_ERROR
99 #endif
100 C-- Read IO error counter
101 beginIOErrCount = IO_ERRCOUNT(myThid)
102 c prevPrec = writeBinaryPrec
103 c writeBinaryPrec = precFloat32
104 CALL WRITE_FLD_XY_RL('ETAtave.',suff,etatave,myIter,myThid)
105 CALL WRITE_FLD_XY_RL('Eta2tave.',suff,Eta2tave,myIter,myThid)
106 CALL WRITE_FLD_XY_RL('PHLtave.',suff,phiHydLowtave,
107 & myIter,myThid)
108 CALL WRITE_FLD_XY_RL('PHL2tave.',suff,phiHydLow2tave,
109 & myIter,myThid)
110
111 CALL WRITE_FLD_XYZ_RL('Ttave.',suff,thetatave,myIter,myThid)
112 CALL WRITE_FLD_XYZ_RL('Stave.',suff,salttave,myIter,myThid)
113 CALL WRITE_FLD_XYZ_RL('uVeltave.',suff,uVeltave,myIter,myThid)
114 CALL WRITE_FLD_XYZ_RL('vVeltave.',suff,vVeltave,myIter,myThid)
115 CALL WRITE_FLD_XYZ_RL('wVeltave.',suff,wVeltave,myIter,myThid)
116 CALL WRITE_FLD_XYZ_RL('TTtave.',suff,TTtave,myIter,myThid)
117 CALL WRITE_FLD_XYZ_RL('UUtave.',suff,UUtave,myIter,myThid)
118 CALL WRITE_FLD_XYZ_RL('VVtave.',suff,VVtave,myIter,myThid)
119 c CALL WRITE_FLD_XYZ_RL('KEtave.',suff,KEtave,myIter,myThid)
120 CALL WRITE_FLD_XYZ_RL('UTtave.',suff,UTtave,myIter,myThid)
121 CALL WRITE_FLD_XYZ_RL('VTtave.',suff,VTtave,myIter,myThid)
122 CALL WRITE_FLD_XYZ_RL('WTtave.',suff,WTtave,myIter,myThid)
123 CALL WRITE_FLD_XYZ_RL('PhHytave.',suff,phiHydtave,myIter,myThid)
124 CALL WRITE_FLD_XYZ_RL('Convtave.',suff,ConvectCountTave,
125 & myIter,myThid)
126 #ifdef NONLIN_FRSURF
127 CALL WRITE_FLD_XYZ_RL('hUtave.',suff,hUtave,myIter,myThid)
128 CALL WRITE_FLD_XYZ_RL('hVtave.',suff,hVtave,myIter,myThid)
129 c CALL WRITE_FLD_XYZ_RL('hFacCtave.',suff,hFacCtave,myIter,myThid)
130 c CALL WRITE_FLD_XYZ_RL('hFacWtave.',suff,hFacWtave,myIter,myThid)
131 c CALL WRITE_FLD_XYZ_RL('hFacStave.',suff,hFacStave,myIter,myThid)
132 #endif /* NONLIN_FRSURF */
133
134 c writeBinaryPrec = prevPrec
135 C-- Reread IO error counter
136 endIOErrCount = IO_ERRCOUNT(myThid)
137
138 C-- Check for IO errors
139 IF ( endIOErrCount .NE. beginIOErrCount ) THEN
140 WRITE(msgBuf,'(A)') 'S/R WRITE_TIME_AVERAGES'
141 CALL PRINT_ERROR( msgBuf, 1 )
142 WRITE(msgBuf,'(A)') 'Error writing out data'
143 CALL PRINT_ERROR( msgBuf, 1 )
144 WRITE(msgBuf,'(A,I10)') 'Timestep ',myIter
145 CALL PRINT_ERROR( msgBuf, 1 )
146 ELSE
147 WRITE(msgBuf,'(A,I10)')
148 & '// Time-average data written, t-step', myIter
149 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
150 & SQUEEZE_RIGHT, 1 )
151 WRITE(msgBuf,'(A)') ' '
152 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
153 & SQUEEZE_RIGHT, 1 )
154 ENDIF
155
156 _END_MASTER( myThid )
157 _BARRIER
158
159 DO bj = myByLo(myThid), myByHi(myThid)
160 DO bi = myBxLo(myThid), myBxHi(myThid)
161
162 C-- Like before the 1rst iteration,
163 C ==> call TIMEAVE_STATVARS with myIter=nIter0 :
164 C 1) Reset the averages to zero ;
165 C 2) Start to cumulate state-variables with Half time step.
166
167 CALL TIMEAVE_STATVARS(myTime, nIter0, bi, bj, myThid)
168
169 ENDDO
170 ENDDO
171
172 ENDIF
173
174 #endif /* ALLOW_TIMEAVE */
175
176 RETURN
177 END
  ViewVC Help
Powered by ViewVC 1.1.22