/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_switch_onoff.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_switch_onoff.F

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

revision 1.11 by molod, Mon May 23 19:31:09 2005 UTC revision 1.12 by jmc, Sun Jun 26 16:51:49 2005 UTC
# Line 43  C      newIter :: future iteration numbe Line 43  C      newIter :: future iteration numbe
43  C      j,m,n   :: loop index  C      j,m,n   :: loop index
44  c     CHARACTER*(MAX_LEN_MBUF) msgBuf  c     CHARACTER*(MAX_LEN_MBUF) msgBuf
45  c     INTEGER newIter  c     INTEGER newIter
46        INTEGER j, m, n        INTEGER m, n, nd
47        INTEGER bi, bj, iSp        INTEGER bi, bj, ip, iSp
48        LOGICAL time4SnapShot        LOGICAL time4SnapShot
49        _RL     phiSec, freqSec        _RL     phiSec, freqSec
50        LOGICAL dBugFlag        LOGICAL dBugFlag
# Line 56  c     INTEGER newIter Line 56  c     INTEGER newIter
56        LOGICAL  DIFF_PHASE_MULTIPLE        LOGICAL  DIFF_PHASE_MULTIPLE
57        EXTERNAL DIFF_PHASE_MULTIPLE        EXTERNAL DIFF_PHASE_MULTIPLE
58    
       _BEGIN_MASTER(myThid)  
   
59  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
60    
61  c     newIter = 1 + myIter  c     newIter = 1 + myIter
# Line 77  C--    Select diagnostics list that uses Line 75  C--    Select diagnostics list that uses
75          endif          endif
76  #endif  #endif
77    
78          IF ( time4SnapShot ) THEN          DO bj=myByLo(myThid), myByHi(myThid)
79  C--     switch ON diagnostics of output-stream # n           DO bi=myBxLo(myThid), myBxHi(myThid)
80            DO m=1,nActive(n)            IF ( time4SnapShot ) THEN
81              j = jdiag(m,n)  C--       switch ON diagnostics of output-stream # n
82  c           IF (ndiag(j).NE.0) WRITE(0,*) myIter,j,ndiag(j),' ->',0              DO m=1,nActive(n)
83              ndiag(j) = 0  c             nd = jdiag(m,n)
84            ENDDO  c             IF (ndiag(nd).NE.0) WRITE(0,*) myIter,nd,ndiag(nd),' ->',0
85          ELSE                ip = idiag(m,n)
86  C--     switch OFF diagnostics of output-stream # n                IF (ip.GT.0) ndiag(ip,bi,bj) = 0
87            DO m=1,nActive(n)              ENDDO
88              j = jdiag(m,n)            ELSE
89  c           IF (ndiag(j).NE.-1) WRITE(0,*) myIter,j,ndiag(j),' ->',-1  C--       switch OFF diagnostics of output-stream # n
90              ndiag(j) = -1              DO m=1,nActive(n)
91            ENDDO  c             nd = jdiag(m,n)
92          ENDIF  c             IF (ndiag(nd).NE.-1) WRITE(0,*) myIter,nd,ndiag(nd),' ->',-1
93                  ip = idiag(m,n)
94                  IF (ip.GT.0) ndiag(ip,bi,bj) = -1
95                ENDDO
96              ENDIF
97             ENDDO
98            ENDDO
99    
100         ENDIF         ENDIF
101        ENDDO        ENDDO
102    
       _END_MASTER(myThid)  
   
 C-jmc: do we need a "BARRIER" at this point ?  
 c     _BARRIER  
   
103  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
104    
105        DO n = 1,diagSt_nbLists        DO n = 1,diagSt_nbLists
# Line 128  C--    Select diagnostics list that uses Line 127  C--    Select diagnostics list that uses
127            IF ( time4SnapShot ) THEN            IF ( time4SnapShot ) THEN
128  C--     switch ON diagnostics of output-stream # n  C--     switch ON diagnostics of output-stream # n
129              DO m=1,diagSt_nbActv(n)              DO m=1,diagSt_nbActv(n)
130                j = jSdiag(m,n)               iSp = iSdiag(m,n)
131                iSp = iSdiag(j)               IF (iSp.GT.0) THEN
132                IF (dBugFlag.AND.qSdiag(0,0,iSp,bi,bj).NE. 0.)                 nd = jSdiag(m,n)
133                   IF (dBugFlag.AND.qSdiag(0,0,iSp,bi,bj).NE. 0.)
134       &          WRITE(0,'(I8,I4,3A,1PE10.3,A,I3)') myIter,       &          WRITE(0,'(I8,I4,3A,1PE10.3,A,I3)') myIter,
135       &            j,' ',cdiag(j),' :',qSdiag(0,0,iSp,bi,bj),' ->', 0       &           nd,' ',cdiag(nd),' :',qSdiag(0,0,iSp,bi,bj),' ->', 0
136                qSdiag(0,0,iSp,bi,bj) = 0.                 qSdiag(0,0,iSp,bi,bj) = 0.
137                 ENDIF
138              ENDDO              ENDDO
139            ELSE            ELSE
140  C--     switch OFF diagnostics of output-stream # n  C--     switch OFF diagnostics of output-stream # n
141              DO m=1,diagSt_nbActv(n)              DO m=1,diagSt_nbActv(n)
142                j = jSdiag(m,n)               iSp = iSdiag(m,n)
143                iSp = iSdiag(j)               IF (iSp.GT.0) THEN
144                IF (dBugFlag.AND.qSdiag(0,0,iSp,bi,bj).NE.-1.)                 nd = jSdiag(m,n)
145                   IF (dBugFlag.AND.qSdiag(0,0,iSp,bi,bj).NE.-1.)
146       &          WRITE(0,'(I8,I4,3A,1PE10.3,A,I3)') myIter,       &          WRITE(0,'(I8,I4,3A,1PE10.3,A,I3)') myIter,
147       &            j,' ',cdiag(j),' :',qSdiag(0,0,iSp,bi,bj),' ->',-1       &           nd,' ',cdiag(nd),' :',qSdiag(0,0,iSp,bi,bj),' ->',-1
148                qSdiag(0,0,iSp,bi,bj) = -1.                 qSdiag(0,0,iSp,bi,bj) = -1.
149                 ENDIF
150              ENDDO              ENDDO
151            ENDIF            ENDIF
152           ENDDO           ENDDO
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22