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

Diff of /MITgcm/pkg/timeave/timeave_statvars.F

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

revision 1.5 by jmc, Mon Oct 7 16:29:04 2002 UTC revision 1.11 by edhill, Fri Oct 24 05:29:36 2003 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3  #include "CPP_OPTIONS.h"  #include "TIMEAVE_OPTIONS.h"
4    
5        SUBROUTINE TIMEAVE_STATVARS(        SUBROUTINE TIMEAVE_STATVARS(
6       I     myTime, myIter, bi, bj, myThid)       I     myTime, myIter, bi, bj, myThid)
# Line 49  C-    Initialize fields for the first ca Line 49  C-    Initialize fields for the first ca
49          CALL TIMEAVE_RESET(TTtave,    Nr, bi, bj, myThid)          CALL TIMEAVE_RESET(TTtave,    Nr, bi, bj, myThid)
50          CALL TIMEAVE_RESET(UUtave,    Nr, bi, bj, myThid)          CALL TIMEAVE_RESET(UUtave,    Nr, bi, bj, myThid)
51          CALL TIMEAVE_RESET(VVtave,    Nr, bi, bj, myThid)          CALL TIMEAVE_RESET(VVtave,    Nr, bi, bj, myThid)
52            CALL TIMEAVE_RESET(UVtave,    Nr, bi, bj, myThid)
53  c       CALL TIMEAVE_RESET(KEtave,    Nr, bi, bj, myThid)  c       CALL TIMEAVE_RESET(KEtave,    Nr, bi, bj, myThid)
54          CALL TIMEAVE_RESET(UTtave,    Nr, bi, bj, myThid)          CALL TIMEAVE_RESET(UTtave,    Nr, bi, bj, myThid)
55          CALL TIMEAVE_RESET(VTtave,    Nr, bi, bj, myThid)          CALL TIMEAVE_RESET(VTtave,    Nr, bi, bj, myThid)
56          CALL TIMEAVE_RESET(WTtave,    Nr, bi, bj, myThid)          CALL TIMEAVE_RESET(WTtave,    Nr, bi, bj, myThid)
57            CALL TIMEAVE_RESET(TdiffRtave,Nr, bi, bj, myThid)
58    #ifndef DISABLE_MOM_VECINV
59            CALL TIMEAVE_RESET(uZetatave, Nr, bi, bj, myThid)
60            CALL TIMEAVE_RESET(vZetatave, Nr, bi, bj, myThid)
61    #endif
62          CALL TIMEAVE_RESET(phiHydtave,Nr, bi, bj, myThid)          CALL TIMEAVE_RESET(phiHydtave,Nr, bi, bj, myThid)
63          CALL TIMEAVE_RESET(phiHydLowtave,1, bi, bj, myThid)          CALL TIMEAVE_RESET(phiHydLowtave,1, bi, bj, myThid)
64          CALL TIMEAVE_RESET(phiHydLow2Tave,1, bi, bj, myThid)          CALL TIMEAVE_RESET(phiHydLow2Tave,1, bi, bj, myThid)
# Line 88  C-    Time Averages of single fields (no Line 94  C-    Time Averages of single fields (no
94        CALL TIMEAVE_CUMULATE(salttave, salt,  Nr, DDTT, bi, bj, myThid)        CALL TIMEAVE_CUMULATE(salttave, salt,  Nr, DDTT, bi, bj, myThid)
95        CALL TIMEAVE_CUMULATE(uVeltave, uVel,  Nr, DDTT, bi, bj, myThid)        CALL TIMEAVE_CUMULATE(uVeltave, uVel,  Nr, DDTT, bi, bj, myThid)
96        CALL TIMEAVE_CUMULATE(vVeltave, vVel,  Nr, DDTT, bi, bj, myThid)        CALL TIMEAVE_CUMULATE(vVeltave, vVel,  Nr, DDTT, bi, bj, myThid)
97        CALL TIMEAVE_CUMULATE(wVeltave, wVel, Nr,  DDTT, bi, bj, myThid)        CALL TIMEAVE_CUMULATE(wVeltave, wVel,  Nr, DDTT, bi, bj, myThid)
98    
99  C-    Time Averages of "double" fields (no hFactor)  C-    Time Averages of "double" fields (no hFactor)
100        CALL TIMEAVE_CUMUL_2V(Eta2tave, etaN,etaN, 1,  0,        CALL TIMEAVE_CUMUL_2V(Eta2tave, etaN,etaN, 1,  0,
# Line 99  C-    Time Averages of "double" fields ( Line 105  C-    Time Averages of "double" fields (
105       &     DDTT, bi, bj, myThid)       &     DDTT, bi, bj, myThid)
106        CALL TIMEAVE_CUMUL_2V(VVtave, vVel,  vVel, Nr, 0,        CALL TIMEAVE_CUMUL_2V(VVtave, vVel,  vVel, Nr, 0,
107       &     DDTT, bi, bj, myThid)       &     DDTT, bi, bj, myThid)
108          CALL TIMEAVE_CUMUL_2V(UVtave, uVel,  vVel, Nr, 12,
109         &     DDTT, bi, bj, myThid)
110  c     CALL TIMEAVE_CUMUL_KE(KEtave, uVel,  vVel, Nr,  c     CALL TIMEAVE_CUMUL_KE(KEtave, uVel,  vVel, Nr,
111  c    &     DDTT, bi, bj, myThid)  c    &     DDTT, bi, bj, myThid)
112        CALL TIMEAVE_CUMUL_2V(WTtave, theta, wVel, Nr, 3,        CALL TIMEAVE_CUMUL_2V(WTtave, theta, wVel, Nr, 3,
# Line 141  C-- end if DDTT ... Line 149  C-- end if DDTT ...
149    
150  C-    Time Averages of "intermediate" fields (no hFactor)  C-    Time Averages of "intermediate" fields (no hFactor)
151        IF ( myIter .NE. nIter0 ) THEN        IF ( myIter .NE. nIter0 ) THEN
152          CALL TIMEAVE_CUMULATE(phiHydLowtave, phiHydLow, 1, deltaTclock,          CALL TIMEAVE_CUMULATE(phiHydtave, totPhihyd,  Nr,
153       &                          bi, bj, myThid)       &                                   deltaTclock, bi, bj, myThid)
154            CALL TIMEAVE_CUMULATE(phiHydLowtave, phiHydLow, 1,
155         &                                   deltaTclock, bi, bj, myThid)
156          CALL TIMEAVE_CUMUL_2V(phiHydLow2Tave,          CALL TIMEAVE_CUMUL_2V(phiHydLow2Tave,
157       &       phiHydLow,phiHydLow, 1,  0, deltaTclock, bi, bj, myThid)       &       phiHydLow,phiHydLow, 1,  0, deltaTclock, bi, bj, myThid)
158          DO K=1,Nr          DO K=1,Nr
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22