/[MITgcm]/MITgcm/pkg/layers/layers_output.F
ViewVC logotype

Diff of /MITgcm/pkg/layers/layers_output.F

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

revision 1.8 by gforget, Wed Sep 19 18:48:18 2012 UTC revision 1.9 by rpa, Wed Oct 17 18:49:15 2012 UTC
# Line 82  cgf this needs layers_maxNum loop Line 82  cgf this needs layers_maxNum loop
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
# Line 123  c set arrays to zero if first timestep Line 123  c set arrays to zero if first timestep
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.
# Line 149  C      Normalize by integrated time Line 149  C      Normalize by integrated time
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 */
# Line 179  C      Normalize by integrated time Line 179  C      Normalize by integrated time
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    
# Line 212  C      Reset averages to zero Line 212  C      Reset averages to zero
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    
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22