82 |
IF ( layers_MDSIO ) THEN |
IF ( layers_MDSIO ) THEN |
83 |
WRITE(suff,'(I10.10)') myIter |
WRITE(suff,'(I10.10)') myIter |
84 |
#ifdef LAYERS_UFLUX |
#ifdef LAYERS_UFLUX |
85 |
CALL WRITE_FLD_3D_RL( 'layers_UFlux.', suff, Nlayers, |
CALL WRITE_FLD_3D_RL( 'layers_UH.', suff, Nlayers, |
86 |
& layers_UFlux, myIter, myThid ) |
& layers_UH, myIter, myThid ) |
87 |
#ifdef LAYERS_THICKNESS |
#ifdef LAYERS_THICKNESS |
88 |
CALL WRITE_FLD_3D_RL( 'layers_HU.', suff, Nlayers, |
CALL WRITE_FLD_3D_RL( 'layers_Hw.', suff, Nlayers, |
89 |
& layers_HU, myIter, myThid ) |
& layers_Hw, myIter, myThid ) |
90 |
#endif /* LAYERS_THICKNESS */ |
#endif /* LAYERS_THICKNESS */ |
91 |
#endif /* LAYERS_UFLUX */ |
#endif /* LAYERS_UFLUX */ |
92 |
#ifdef LAYERS_VFLUX |
#ifdef LAYERS_VFLUX |
93 |
CALL WRITE_FLD_3D_RL( 'layers_VFlux.', suff, Nlayers, |
CALL WRITE_FLD_3D_RL( 'layers_VH.', suff, Nlayers, |
94 |
& layers_VFlux, myIter, myThid ) |
& layers_VH, myIter, myThid ) |
95 |
#ifdef LAYERS_THICKNESS |
#ifdef LAYERS_THICKNESS |
96 |
CALL WRITE_FLD_3D_RL( 'layers_HV.', suff, Nlayers, |
CALL WRITE_FLD_3D_RL( 'layers_Hs.', suff, Nlayers, |
97 |
& layers_HV, myIter, myThid ) |
& layers_Hs, myIter, myThid ) |
98 |
#endif /* LAYERS_THICKNESS */ |
#endif /* LAYERS_THICKNESS */ |
99 |
#endif /* LAYERS_VFLUX */ |
#endif /* LAYERS_VFLUX */ |
100 |
#ifdef LAYERS_PRHO_REF |
#ifdef LAYERS_PRHO_REF |
123 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
124 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
125 |
#ifdef LAYERS_UFLUX |
#ifdef LAYERS_UFLUX |
126 |
CALL TIMEAVE_RESET(layers_UFlux_T,Nlayers,bi,bj,myThid) |
CALL TIMEAVE_RESET(layers_UH_T,Nlayers,bi,bj,myThid) |
127 |
#ifdef LAYERS_THICKNESS |
#ifdef LAYERS_THICKNESS |
128 |
CALL TIMEAVE_RESET(layers_HU_T,Nlayers,bi,bj,myThid) |
CALL TIMEAVE_RESET(layers_Hw_T,Nlayers,bi,bj,myThid) |
129 |
#endif /* LAYERS_THICKNESS */ |
#endif /* LAYERS_THICKNESS */ |
130 |
#endif /* LAYERS_UFLUX */ |
#endif /* LAYERS_UFLUX */ |
131 |
|
|
132 |
#ifdef LAYERS_VFLUX |
#ifdef LAYERS_VFLUX |
133 |
CALL TIMEAVE_RESET(layers_VFlux_T,Nlayers,bi,bj,myThid) |
CALL TIMEAVE_RESET(layers_VH_T,Nlayers,bi,bj,myThid) |
134 |
#ifdef LAYERS_THICKNESS |
#ifdef LAYERS_THICKNESS |
135 |
CALL TIMEAVE_RESET(layers_HV_T,Nlayers,bi,bj,myThid) |
CALL TIMEAVE_RESET(layers_Hs_T,Nlayers,bi,bj,myThid) |
136 |
#endif /* LAYERS_THICKNESS */ |
#endif /* LAYERS_THICKNESS */ |
137 |
#endif /* LAYERS_VFLUX */ |
#endif /* LAYERS_VFLUX */ |
138 |
layers_TimeAve(bi,bj) = 0. |
layers_TimeAve(bi,bj) = 0. |
149 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
150 |
|
|
151 |
#ifdef LAYERS_UFLUX |
#ifdef LAYERS_UFLUX |
152 |
CALL TIMEAVE_NORMALIZE( layers_UFlux_T, layers_timeave, |
CALL TIMEAVE_NORMALIZE( layers_UH_T, layers_timeave, |
153 |
& Nlayers, bi, bj, myThid ) |
& Nlayers, bi, bj, myThid ) |
154 |
#ifdef LAYERS_THICKNESS |
#ifdef LAYERS_THICKNESS |
155 |
CALL TIMEAVE_NORMALIZE( layers_HU_T, layers_timeave, |
CALL TIMEAVE_NORMALIZE( layers_Hw_T, layers_timeave, |
156 |
& Nlayers, bi, bj, myThid ) |
& Nlayers, bi, bj, myThid ) |
157 |
#endif /* LAYERS_THICKNESS */ |
#endif /* LAYERS_THICKNESS */ |
158 |
#endif /* LAYERS_UFLUX */ |
#endif /* LAYERS_UFLUX */ |
159 |
|
|
160 |
#ifdef LAYERS_VFLUX |
#ifdef LAYERS_VFLUX |
161 |
CALL TIMEAVE_NORMALIZE( layers_VFlux_T, layers_timeave, |
CALL TIMEAVE_NORMALIZE( layers_VH_T, layers_timeave, |
162 |
& Nlayers, bi, bj, myThid ) |
& Nlayers, bi, bj, myThid ) |
163 |
#ifdef LAYERS_THICKNESS |
#ifdef LAYERS_THICKNESS |
164 |
CALL TIMEAVE_NORMALIZE( layers_HV_T, layers_timeave, |
CALL TIMEAVE_NORMALIZE( layers_Hs_T, layers_timeave, |
165 |
& Nlayers, bi, bj, myThid ) |
& Nlayers, bi, bj, myThid ) |
166 |
#endif /* LAYERS_THICKNESS */ |
#endif /* LAYERS_THICKNESS */ |
167 |
#endif /* LAYERS_VFLUX */ |
#endif /* LAYERS_VFLUX */ |
179 |
IF ( layers_MDSIO ) THEN |
IF ( layers_MDSIO ) THEN |
180 |
WRITE(suff,'(I10.10)') myIter |
WRITE(suff,'(I10.10)') myIter |
181 |
#ifdef LAYERS_UFLUX |
#ifdef LAYERS_UFLUX |
182 |
CALL WRITE_FLD_3D_RL( 'layers_UFlux-tave.', suff, Nlayers, |
CALL WRITE_FLD_3D_RL( 'layers_UH-tave.', suff, Nlayers, |
183 |
& layers_UFlux_T, myIter, myThid ) |
& layers_UH_T, myIter, myThid ) |
184 |
#ifdef LAYERS_THICKNESS |
#ifdef LAYERS_THICKNESS |
185 |
CALL WRITE_FLD_3D_RL( 'layers_HU-tave.', suff, Nlayers, |
CALL WRITE_FLD_3D_RL( 'layers_Hw-tave.', suff, Nlayers, |
186 |
& layers_HU_T, myIter, myThid ) |
& layers_Hw_T, myIter, myThid ) |
187 |
#endif /* LAYERS_THICKNESS */ |
#endif /* LAYERS_THICKNESS */ |
188 |
#endif /* LAYERS_UFLUX */ |
#endif /* LAYERS_UFLUX */ |
189 |
#ifdef LAYERS_VFLUX |
#ifdef LAYERS_VFLUX |
190 |
CALL WRITE_FLD_3D_RL( 'layers_VFlux-tave.', suff, Nlayers, |
CALL WRITE_FLD_3D_RL( 'layers_VH-tave.', suff, Nlayers, |
191 |
& layers_VFlux_T, myIter, myThid ) |
& layers_VH_T, myIter, myThid ) |
192 |
#ifdef LAYERS_THICKNESS |
#ifdef LAYERS_THICKNESS |
193 |
CALL WRITE_FLD_3D_RL( 'layers_HV-tave.', suff, Nlayers, |
CALL WRITE_FLD_3D_RL( 'layers_Hs-tave.', suff, Nlayers, |
194 |
& layers_HV_T, myIter, myThid ) |
& layers_Hs_T, myIter, myThid ) |
195 |
#endif /* LAYERS_THICKNESS */ |
#endif /* LAYERS_THICKNESS */ |
196 |
#endif /* LAYERS_VFLUX */ |
#endif /* LAYERS_VFLUX */ |
197 |
|
|
212 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
213 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
214 |
#ifdef LAYERS_UFLUX |
#ifdef LAYERS_UFLUX |
215 |
CALL TIMEAVE_RESET(layers_UFlux_T,Nlayers,bi,bj,myThid) |
CALL TIMEAVE_RESET(layers_UH_T,Nlayers,bi,bj,myThid) |
216 |
#ifdef LAYERS_THICKNESS |
#ifdef LAYERS_THICKNESS |
217 |
CALL TIMEAVE_RESET(layers_HU_T,Nlayers,bi,bj,myThid) |
CALL TIMEAVE_RESET(layers_Hw_T,Nlayers,bi,bj,myThid) |
218 |
#endif /* LAYERS_THICKNESS */ |
#endif /* LAYERS_THICKNESS */ |
219 |
#endif /* LAYERS_UFLUX */ |
#endif /* LAYERS_UFLUX */ |
220 |
|
|
221 |
#ifdef LAYERS_VFLUX |
#ifdef LAYERS_VFLUX |
222 |
CALL TIMEAVE_RESET(layers_VFlux_T,Nlayers,bi,bj,myThid) |
CALL TIMEAVE_RESET(layers_VH_T,Nlayers,bi,bj,myThid) |
223 |
#ifdef LAYERS_THICKNESS |
#ifdef LAYERS_THICKNESS |
224 |
CALL TIMEAVE_RESET(layers_HV_T,Nlayers,bi,bj,myThid) |
CALL TIMEAVE_RESET(layers_Hs_T,Nlayers,bi,bj,myThid) |
225 |
#endif /* LAYERS_THICKNESS */ |
#endif /* LAYERS_THICKNESS */ |
226 |
#endif /* LAYERS_VFLUX */ |
#endif /* LAYERS_VFLUX */ |
227 |
|
|