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

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

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

revision 1.2 by jmc, Mon Nov 5 18:49:21 2007 UTC revision 1.3 by jmc, Sat Nov 10 22:09:32 2007 UTC
# Line 33  C     !LOCAL VARIABLES: Line 33  C     !LOCAL VARIABLES:
33        LOGICAL  DIAGNOSTICS_IS_ON        LOGICAL  DIAGNOSTICS_IS_ON
34        EXTERNAL DIAGNOSTICS_IS_ON        EXTERNAL DIAGNOSTICS_IS_ON
35        _RL dummy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL dummy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
36        INTEGER i,j,k,bi,bj,N        INTEGER i,j,k,bi,bj,iTr
37        CHARACTER*8 diagName        CHARACTER*8 diagName
38        INTEGER km1        INTEGER km1
39    
40        diagName = '        '        diagName = '        '
41    
42  c     DO N = 1,PTRACERS_numInUse        DO iTr = 1,PTRACERS_numInUse
 C-    due to format (I2) and diagnostics name limitation,  
 C      => do not go beyong 99 tracers  
       DO N = 1,MIN(99,PTRACERS_numInUse)  
43    
44          diagName = '        '          diagName = '        '
45          WRITE(diagName,'(A4,I2.2)') 'TRAC',N          WRITE(diagName,'(A4,A2)') 'TRAC',PTRACERS_ioLabel(iTr)
46          CALL DIAGNOSTICS_FILL( pTracer(1-Olx,1-Oly,1,1,1,N), diagName,          CALL DIAGNOSTICS_FILL( pTracer(1-Olx,1-Oly,1,1,1,iTr), diagName,
47       &                         0,Nr,0,1,1,myThid )       &                         0,Nr,0,1,1,myThid )
48    
49          diagName = '        '          diagName = '        '
50          WRITE(diagName,'(A5,I2.2)') 'UTRAC',N          WRITE(diagName,'(A5,A2)') 'UTRAC',PTRACERS_ioLabel(iTr)
51          IF ( DIAGNOSTICS_IS_ON(diagName,myThid) ) THEN          IF ( DIAGNOSTICS_IS_ON(diagName,myThid) ) THEN
52           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
53            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 59  C      => do not go beyong 99 tracers Line 56  C      => do not go beyong 99 tracers
56                  DO i = 1,sNx+1                  DO i = 1,sNx+1
57                    dummy(i,j,k,bi,bj) =                    dummy(i,j,k,bi,bj) =
58       &                 uVel(i,j,k,bi,bj)*hFacW(i,j,k,bi,bj) *       &                 uVel(i,j,k,bi,bj)*hFacW(i,j,k,bi,bj) *
59       &                 0.5 _d 0*(pTracer(i,j,k,bi,bj,N)       &                 0.5 _d 0*( pTracer(i,j,k,bi,bj,iTr)
60       &                       + pTracer(i-1,j,k,bi,bj,N))       &                          + pTracer(i-1,j,k,bi,bj,iTr) )
61                  ENDDO                  ENDDO
62                ENDDO                ENDDO
63              ENDDO              ENDDO
# Line 70  C      => do not go beyong 99 tracers Line 67  C      => do not go beyong 99 tracers
67          ENDIF          ENDIF
68    
69          diagName = '        '          diagName = '        '
70          WRITE(diagName,'(A5,I2.2)') 'VTRAC',N          WRITE(diagName,'(A5,A2)') 'VTRAC',PTRACERS_ioLabel(iTr)
71          IF ( DIAGNOSTICS_IS_ON(diagName,myThid) ) THEN          IF ( DIAGNOSTICS_IS_ON(diagName,myThid) ) THEN
72           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
73            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 79  C      => do not go beyong 99 tracers Line 76  C      => do not go beyong 99 tracers
76                  DO i = 1,sNx                  DO i = 1,sNx
77                    dummy(i,j,k,bi,bj) =                    dummy(i,j,k,bi,bj) =
78       &                 vVel(i,j,k,bi,bj)*hFacS(i,j,k,bi,bj) *       &                 vVel(i,j,k,bi,bj)*hFacS(i,j,k,bi,bj) *
79       &                 0.5 _d 0*(pTracer(i,j,k,bi,bj,N)       &                 0.5 _d 0*( pTracer(i,j,k,bi,bj,iTr)
80       &                       + pTracer(i,j-1,k,bi,bj,N))       &                          + pTracer(i,j-1,k,bi,bj,iTr) )
81                  ENDDO                  ENDDO
82                ENDDO                ENDDO
83              ENDDO              ENDDO
# Line 90  C      => do not go beyong 99 tracers Line 87  C      => do not go beyong 99 tracers
87          ENDIF          ENDIF
88    
89          diagName = '        '          diagName = '        '
90          WRITE(diagName,'(A5,I2.2)') 'WTRAC',N          WRITE(diagName,'(A5,A2)') 'WTRAC',PTRACERS_ioLabel(iTr)
91          IF ( DIAGNOSTICS_IS_ON(diagName,myThid) ) THEN          IF ( DIAGNOSTICS_IS_ON(diagName,myThid) ) THEN
92           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
93            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 99  C      => do not go beyong 99 tracers Line 96  C      => do not go beyong 99 tracers
96                DO j = 1,sNy                DO j = 1,sNy
97                  DO i = 1,sNx                  DO i = 1,sNx
98                    dummy(i,j,k,bi,bj) = wVel(i,j,k,bi,bj) *                    dummy(i,j,k,bi,bj) = wVel(i,j,k,bi,bj) *
99       &                 0.5 _d 0*(pTracer(i,j,k,bi,bj,N)       &                 0.5 _d 0*( pTracer(i,j,k,bi,bj,iTr)
100       &                       + pTracer(i,j,km1,bi,bj,N))       &                          + pTracer(i,j,km1,bi,bj,iTr) )
101                  ENDDO                  ENDDO
102                ENDDO                ENDDO
103              ENDDO              ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22