/[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.4 by mlosch, Wed Sep 18 16:38:02 2002 UTC revision 1.5 by jmc, Mon Oct 7 16:29:04 2002 UTC
# Line 71  c       CALL TIMEAVE_RESET(hFacStave, Nr Line 71  c       CALL TIMEAVE_RESET(hFacStave, Nr
71        ENDIF        ENDIF
72                
73  C--   Cumulate state-variables with Half or Full time step :  C--   Cumulate state-variables with Half or Full time step :
74        IF ( myIter .EQ. nIter0 .OR.        IF ( myIter .EQ. nIter0 ) THEN
75            DDTT = deltaTclock*(1. _d 0 - tave_lastIter)
76          ELSEIF (
77       &  DIFFERENT_MULTIPLE(taveFreq, myTime, myTime-deltaTClock)) THEN       &  DIFFERENT_MULTIPLE(taveFreq, myTime, myTime-deltaTClock)) THEN
78         DDTT=0.5*deltaTclock          DDTT = deltaTclock*tave_lastIter
79        ELSE        ELSE
80         DDTT=deltaTclock          DDTT = deltaTclock
81        ENDIF        ENDIF
82    
83          IF ( DDTT .NE. 0. _d 0) THEN
84    
85  C-    Time Averages of single fields (no hFactor)  C-    Time Averages of single fields (no hFactor)
86        CALL TIMEAVE_CUMULATE(etatave,  etaN,  1 , DDTT, bi, bj, myThid)        CALL TIMEAVE_CUMULATE(etatave,  etaN,  1 , DDTT, bi, bj, myThid)
87        CALL TIMEAVE_CUMULATE(thetatave,theta, Nr, DDTT, bi, bj, myThid)        CALL TIMEAVE_CUMULATE(thetatave,theta, Nr, DDTT, bi, bj, myThid)
88        CALL TIMEAVE_CUMULATE(salttave, salt,  Nr, DDTT, bi, bj, myThid)        CALL TIMEAVE_CUMULATE(salttave, salt,  Nr, DDTT, bi, bj, myThid)
89        CALL TIMEAVE_CUMULATE(uVeltave, uVel,  Nr, DDTT, bi, bj, myThid)        CALL TIMEAVE_CUMULATE(uVeltave, uVel,  Nr, DDTT, bi, bj, myThid)
90        CALL TIMEAVE_CUMULATE(vVeltave, vVel,  Nr, DDTT, bi, bj, myThid)        CALL TIMEAVE_CUMULATE(vVeltave, vVel,  Nr, DDTT, bi, bj, myThid)
91          CALL TIMEAVE_CUMULATE(wVeltave, wVel, Nr,  DDTT, bi, bj, myThid)
92    
93  C-    Time Averages of "double" fields (no hFactor)  C-    Time Averages of "double" fields (no hFactor)
94        CALL TIMEAVE_CUMUL_2V(Eta2tave, etaN,etaN, 1,  0,        CALL TIMEAVE_CUMUL_2V(Eta2tave, etaN,etaN, 1,  0,
# Line 96  C-    Time Averages of "double" fields ( Line 101  C-    Time Averages of "double" fields (
101       &     DDTT, bi, bj, myThid)       &     DDTT, bi, bj, myThid)
102  c     CALL TIMEAVE_CUMUL_KE(KEtave, uVel,  vVel, Nr,  c     CALL TIMEAVE_CUMUL_KE(KEtave, uVel,  vVel, Nr,
103  c    &     DDTT, bi, bj, myThid)  c    &     DDTT, bi, bj, myThid)
104          CALL TIMEAVE_CUMUL_2V(WTtave, theta, wVel, Nr, 3,
105         &     DDTT, bi, bj, myThid)
106    
107  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
108    
# Line 129  C-    Keep record of how much time has b Line 136  C-    Keep record of how much time has b
136          TimeAve_half(k,bi,bj)=TimeAve_half(k,bi,bj)+DDTT          TimeAve_half(k,bi,bj)=TimeAve_half(k,bi,bj)+DDTT
137        ENDDO        ENDDO
138    
139  C- CAUTIOUS : wVel and theta are not synchronously updated during the model  C-- end if DDTT ...
140  C   time-stepping, so that the time-average diagnostic of the cross product        ENDIF
 C   W*Theta is not perfectly accurate.  
141    
142  C-    Time Averages of "intermediate" fields (no hFactor)  C-    Time Averages of "intermediate" fields (no hFactor)
143        IF ( myIter .NE. nIter0 ) THEN        IF ( myIter .NE. nIter0 ) THEN
         CALL TIMEAVE_CUMULATE(wVeltave, wVel, Nr, deltaTclock,  
      &                          bi, bj, myThid)  
         CALL TIMEAVE_CUMUL_2V(WTtave, theta, wVel, Nr, 3,  
      &                        deltaTclock, bi, bj, myThid)  
144          CALL TIMEAVE_CUMULATE(phiHydLowtave, phiHydLow, 1, deltaTclock,          CALL TIMEAVE_CUMULATE(phiHydLowtave, phiHydLow, 1, deltaTclock,
145       &                          bi, bj, myThid)       &                          bi, bj, myThid)
146          CALL TIMEAVE_CUMUL_2V(phiHydLow2Tave,          CALL TIMEAVE_CUMUL_2V(phiHydLow2Tave,
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22