/[MITgcm]/MITgcm/diags/src/do_time_averages.F
ViewVC logotype

Diff of /MITgcm/diags/src/do_time_averages.F

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

revision 1.8 by adcroft, Thu Aug 26 17:47:37 1999 UTC revision 1.9 by adcroft, Wed Jun 21 19:28:40 2000 UTC
# Line 2  Line 2 
2    
3        SUBROUTINE DO_TIME_AVERAGES(        SUBROUTINE DO_TIME_AVERAGES(
4       I     myCurrentTime, myIter, bi, bj, K, kUp, kDown,       I     myCurrentTime, myIter, bi, bj, K, kUp, kDown,
5       I     K13,K23,wVel,KapGM,ConvectCount,       I     wVel,ConvectCount,
6       I     myThid)       I     myThid)
7  C     /==========================================================\  C     /==========================================================\
8  C     | SUBROUTINE DO_TIME_AVERAGES                              |  C     | SUBROUTINE DO_TIME_AVERAGES                              |
# Line 28  C     myCurrentTime - Current time of si Line 28  C     myCurrentTime - Current time of si
28        INTEGER myThid        INTEGER myThid
29        INTEGER myIter, bi, bj, K, kUp, kDown        INTEGER myIter, bi, bj, K, kUp, kDown
30        _RL     myCurrentTime        _RL     myCurrentTime
       _RL K13   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)  
       _RL K23   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)  
31        _RL wVel  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)        _RL wVel  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)
       _RL KapGM (1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
32        _RL ConvectCount (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL ConvectCount (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
33    
34  C     == Local variables ==  C     == Local variables ==
# Line 48  C     Initialize fields for the first ca Line 45  C     Initialize fields for the first ca
45           CALL TIMEAVER_INI_XYZ(vVeltave, bi, bj, K, myThid)           CALL TIMEAVER_INI_XYZ(vVeltave, bi, bj, K, myThid)
46           CALL TIMEAVER_INI_XYZ(uttave, bi, bj, K, myThid)           CALL TIMEAVER_INI_XYZ(uttave, bi, bj, K, myThid)
47           CALL TIMEAVER_INI_XYZ(vttave, bi, bj, K, myThid)           CALL TIMEAVER_INI_XYZ(vttave, bi, bj, K, myThid)
          CALL TIMEAVER_INI_XYZ(k13tave, bi, bj, K, myThid)  
          CALL TIMEAVER_INI_XYZ(k23tave, bi, bj, K, myThid)  
48           CALL TIMEAVER_INI_XYZ(wVeltave, bi, bj, K, myThid)           CALL TIMEAVER_INI_XYZ(wVeltave, bi, bj, K, myThid)
          IF (TOP_LEVEL) CALL TIMEAVER_INI_XY(KapGMtave, bi, bj, myThid)  
49           CALL TIMEAVER_INI_XYZ(ConvectCountTave, bi, bj, K, myThid)           CALL TIMEAVER_INI_XYZ(ConvectCountTave, bi, bj, K, myThid)
50           TimeAve(k,bi,bj)=0.           TimeAve(k,bi,bj)=0.
51           InterTimeAve(k,bi,bj)=0.           InterTimeAve(k,bi,bj)=0.
# Line 83  C     Keep record of how much time has b Line 77  C     Keep record of how much time has b
77    
78  C     Time Averages of "intermediate" fields  C     Time Averages of "intermediate" fields
79        DDTT=deltaTclock        DDTT=deltaTclock
       CALL TIMEAVER_INTER_XYZ(Nr, K13, k13tave, DDTT,  
      & bi, bj, K, K, myThid)  
       CALL TIMEAVER_INTER_XYZ(Nr, K23, k23tave, DDTT,  
      & bi, bj, K, K, myThid)  
80        CALL TIMEAVER_INTER_XYZ(2, wVel, wVeltave, DDTT,        CALL TIMEAVER_INTER_XYZ(2, wVel, wVeltave, DDTT,
81       & bi, bj, K, kUp, myThid)       & bi, bj, K, kUp, myThid)
       IF (TOP_LEVEL)  
      &  CALL TIMEAVER_INTER_XY(KapGM, KapGMtave, DDTT, bi, bj, myThid)  
82        InterTimeAve(k,bi,bj)=InterTimeAve(k,bi,bj)+DDTT        InterTimeAve(k,bi,bj)=InterTimeAve(k,bi,bj)+DDTT
83        CALL TIMEAVER_INTER_XYZ(Nr, ConvectCount, ConvectCountTave, DDTT,        CALL TIMEAVER_INTER_XYZ(Nr, ConvectCount, ConvectCountTave, DDTT,
84       & bi, bj, K, K, myThid)       & bi, bj, K, K, myThid)
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22