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

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

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

revision 1.5 by jmc, Tue Jul 13 16:47:49 2004 UTC revision 1.6 by jmc, Sun Jul 18 01:12:04 2004 UTC
# Line 43  C     myTime - Current time of simulatio Line 43  C     myTime - Current time of simulatio
43  C     == Local variables ==  C     == Local variables ==
44        INTEGER iTracer, i, j, k        INTEGER iTracer, i, j, k
45        _RL DDTT        _RL DDTT
46        _RL tempArray (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  c     _RL tempArray (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
47    
48  C-    Initialize fields for the first call ever  C-    Initialize fields for the first call ever
49        IF ( myIter .EQ. nIter0 ) THEN        IF ( myIter .EQ. nIter0 ) THEN
# Line 83  C-    Time Averages of "intermediate" fi Line 83  C-    Time Averages of "intermediate" fi
83        IF ( myIter .NE. nIter0 ) THEN        IF ( myIter .NE. nIter0 ) THEN
84    
85  C-    Time Averages of surface fluxes  C-    Time Averages of surface fluxes
86         IF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN  c      IF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN
87          k=Nr  c       k=Nr
88         ELSE  c      ELSE
89          k=1  c       k=1
90         ENDIF  c      ENDIF
91    
92  C     ptracerFlux  C     ptracerFlux
93         DO iTracer=1,PTRACERS_numInUse         DO iTracer=1,PTRACERS_numInUse
94          DO j=1,sNy  c       DO j=1,sNy
95           DO i=1,sNx  c        DO i=1,sNx
96            tempArray(i,j,bi,bj)=maskC(i,j,k,bi,bj)*  c         tempArray(i,j,bi,bj)=maskC(i,j,k,bi,bj)*
97       &     surfaceTendencyPtr(i,j,bi,bj,iTracer)*  c    &     surfaceForcingPtr(i,j,bi,bj,iTracer)*
98       &     drF(k)*hFacC(i,j,k,bi,bj)  c    &     drF(k)*hFacC(i,j,k,bi,bj)
99           ENDDO  c        ENDDO
100          ENDDO  c       ENDDO
101    c       CALL TIMEAVE_CUMULATE(ptracerFluxtave(1-Olx,1-Oly,1,1,iTracer),
102    c    &                        tempArray,1,deltaTclock,bi,bj,myThid)
103          CALL TIMEAVE_CUMULATE(ptracerFluxtave(1-Olx,1-Oly,1,1,iTracer),          CALL TIMEAVE_CUMULATE(ptracerFluxtave(1-Olx,1-Oly,1,1,iTracer),
104       &                        tempArray,1,deltaTclock,bi,bj,myThid)       &                      surfaceForcingPtr(1-Olx,1-Oly,1,1,iTracer),
105         &                                  1,deltaTclock,bi,bj,myThid)
106         ENDDO         ENDDO
107    
108         DO K=1,Nr         DO K=1,Nr
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22