/[MITgcm]/MITgcm/pkg/ptracers/ptracers_timeave.F
ViewVC logotype

Diff of /MITgcm/pkg/ptracers/ptracers_timeave.F

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

revision 1.1 by jmc, Sat Jan 2 23:45:46 2010 UTC revision 1.2 by jmc, Sun Jan 3 03:11:21 2010 UTC
# Line 67  C-    Initialize fields for the first ca Line 67  C-    Initialize fields for the first ca
67          ptracer_full(bi,bj) = 0.          ptracer_full(bi,bj) = 0.
68        ENDIF        ENDIF
69    
70          IF ( PTRACERS_taveFreq.GT.0. _d 0 ) THEN
71  C--   Cumulate state-variables with Half or Full time step :  C--   Cumulate state-variables with Half or Full time step :
72        IF ( myIter .EQ. nIter0 .OR.         IF ( myIter .EQ. nIter0 .OR.
73       &     DIFFERENT_MULTIPLE( PTRACERS_taveFreq, myTime, deltaTClock )       &    DIFFERENT_MULTIPLE( PTRACERS_taveFreq, myTime, deltaTClock )
74       &   ) THEN       &    ) THEN
75          DDTT = 0.5*deltaTclock          DDTT = 0.5*deltaTclock
76        ELSE         ELSE
77          DDTT = deltaTclock          DDTT = deltaTclock
78        ENDIF         ENDIF
79    
80  C-    Time Averages of state variables  C-    Time Averages of state variables
81        DO iTr=1,PTRACERS_numInUse         DO iTr=1,PTRACERS_numInUse
82          CALL TIMEAVE_CUMULATE( ptracertave(1-Olx,1-Oly,1,1,1,iTr),           CALL TIMEAVE_CUMULATE( ptracertave(1-Olx,1-Oly,1,1,1,iTr),
83       &                         pTracer(1-Olx,1-Oly,1,1,1,iTr),       &                          pTracer(1-Olx,1-Oly,1,1,1,iTr),
84       &                         Nr, DDTT, bi, bj, myThid )       &                          Nr, DDTT, bi, bj, myThid )
85        ENDDO         ENDDO
86  C-    Keep record of how much time has been integrated over  C-    Keep record of how much time has been integrated over
87        ptracer_half(bi,bj) = ptracer_half(bi,bj)+DDTT          ptracer_half(bi,bj) = ptracer_half(bi,bj)+DDTT
88    
89  C-    Time Averages of "intermediate" fields  C-    Time Averages of "intermediate" fields
90        IF ( myIter .NE. nIter0 ) THEN         IF ( myIter .NE. nIter0 ) THEN
91    
92  C-    Time Averages of surface fluxes  C-    Time Averages of surface fluxes
93         DO iTr=1,PTRACERS_numInUse          DO iTr=1,PTRACERS_numInUse
94          CALL TIMEAVE_CUMULATE( ptracerFluxtave(1-Olx,1-Oly,1,1,iTr),           CALL TIMEAVE_CUMULATE( ptracerFluxtave(1-Olx,1-Oly,1,1,iTr),
95       &                       surfaceForcingPTr(1-Olx,1-Oly,1,1,iTr),       &                        surfaceForcingPTr(1-Olx,1-Oly,1,1,iTr),
96       &                         1, deltaTclock, bi, bj, myThid )       &                          1, deltaTclock, bi, bj, myThid )
97         ENDDO          ENDDO
98         ptracer_full(bi,bj) = ptracer_full(bi,bj)+deltaTclock          ptracer_full(bi,bj) = ptracer_full(bi,bj)+deltaTclock
99    
100           ENDIF
101    
102    C-    end block if PTRACERS_taveFreq > 0
103        ENDIF        ENDIF
104    
105  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22