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

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

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

Revision 1.4 - (hide 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 mlosch 1.4 C $Header: /u/gcmpack/MITgcm/pkg/timeave/timeave_statv_write.F,v 1.3 2002/01/03 16:25:44 jmc Exp $
2 jmc 1.1 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 jmc 1.3 #include "GRID.h"
19 jmc 1.1 #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 jmc 1.3 #ifdef ALLOW_TIMEAVE
35    
36 jmc 1.1 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 jmc 1.3 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 jmc 1.1
52     WRITE(suff,'(I10.10)') myIter
53    
54     DO bj = myByLo(myThid), myByHi(myThid)
55     DO bi = myBxLo(myThid), myBxHi(myThid)
56 jmc 1.3 C- Normalize by integrated time
57     CALL TIMEAVE_NORMALIZ(etatave, TimeAve_half,1 ,bi,bj,myThid)
58 jmc 1.1 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 jmc 1.3 CALL TIMEAVE_NORMALIZ(Eta2tave, TimeAve_half,1 ,bi,bj,myThid)
63 adcroft 1.2 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 jmc 1.3 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 jmc 1.1
78     CALL TIMEAVE_NORMALIZ(wVeltave, TimeAve_full,Nr,bi,bj,myThid)
79 jmc 1.3 CALL TIMEAVE_NORMALIZ(WTtave, TimeAve_full,Nr,bi,bj,myThid)
80 jmc 1.1 CALL TIMEAVE_NORMALIZ(phiHydtave,TimeAve_full,Nr,
81     & bi,bj,myThid)
82 mlosch 1.4 CALL TIMEAVE_NORMALIZ(phiHydLowtave,TimeAve_full,1,
83     & bi,bj,myThid)
84     CALL TIMEAVE_NORMALIZ(phiHydLow2Tave,TimeAve_full,1,
85     & bi,bj,myThid)
86 jmc 1.1 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 jmc 1.3 CALL DFILE_SET_RW
98     CALL DFILE_SET_CONT_ON_ERROR
99 jmc 1.1 #endif
100 jmc 1.3 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 mlosch 1.4 CALL WRITE_FLD_XY_RL('PHLtave.',suff,phiHydLowtave,
107     & myIter,myThid)
108     CALL WRITE_FLD_XY_RL('PHL2tave.',suff,phiHydLow2tave,
109     & myIter,myThid)
110 jmc 1.3
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 jmc 1.1
138 jmc 1.3 C-- Check for IO errors
139     IF ( endIOErrCount .NE. beginIOErrCount ) THEN
140 jmc 1.1 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 jmc 1.3 ELSE
147 jmc 1.1 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 jmc 1.3 ENDIF
155 jmc 1.1
156     _END_MASTER( myThid )
157     _BARRIER
158    
159     DO bj = myByLo(myThid), myByHi(myThid)
160     DO bi = myBxLo(myThid), myBxHi(myThid)
161    
162 jmc 1.3 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 jmc 1.1
169     ENDDO
170     ENDDO
171    
172     ENDIF
173    
174 jmc 1.3 #endif /* ALLOW_TIMEAVE */
175    
176 jmc 1.1 RETURN
177     END
  ViewVC Help
Powered by ViewVC 1.1.22