/[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.7 by molod, Mon May 16 23:41:32 2005 UTC revision 1.11 by molod, Mon May 23 19:31:09 2005 UTC
# Line 28  C     !USES: Line 28  C     !USES:
28  #include "PARAMS.h"  #include "PARAMS.h"
29  #include "DIAGNOSTICS_SIZE.h"  #include "DIAGNOSTICS_SIZE.h"
30  #include "DIAGNOSTICS.h"  #include "DIAGNOSTICS.h"
 #ifdef ALLOW_FIZHI  
 #include "chronos.h"  
 #endif  
31    
32  C     !INPUT PARAMETERS:  C     !INPUT PARAMETERS:
33  C     myTime     :: current Time of simulation ( s )  C     myTime     :: current Time of simulation ( s )
# Line 47  C      j,m,n   :: loop index Line 44  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 j, m, n
47          INTEGER bi, bj, iSp
48        LOGICAL time4SnapShot        LOGICAL time4SnapShot
49        _RL     phiSec, freqSec        _RL     phiSec, freqSec
50          LOGICAL dBugFlag
51  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
52        integer mmdd,hhmmss,nsecf2,alarm2        logical alarm2,alarm2next
53        character *9 tagname        character *9 tagname
54  #endif  #endif
55    
# Line 73  C--    Select diagnostics list that uses Line 72  C--    Select diagnostics list that uses
72       &                                       myTime, deltaTclock )       &                                       myTime, deltaTclock )
73  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
74          if( useFIZHI) then          if( useFIZHI) then
          mmdd = int(freq(n))  
          hhmmss = int((freq(n) - int(freq(n)))*1.e6)  
          freqSec = nsecf2(hhmmss,mmdd,nymd)  
75           write(tagname,'(A,I2.2)')'diagtag',n           write(tagname,'(A,I2.2)')'diagtag',n
76           time4SnapShot = alarm2(tagname)           time4SnapShot = alarm2next(tagname,deltaT)
77          endif          endif
78  #endif  #endif
79    
# Line 100  c           IF (ndiag(j).NE.-1) WRITE(0, Line 96  c           IF (ndiag(j).NE.-1) WRITE(0,
96         ENDIF         ENDIF
97        ENDDO        ENDDO
98    
 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  
   
99        _END_MASTER(myThid)        _END_MASTER(myThid)
100    
101  C-jmc: do we need a "BARRIER" at this point ?  C-jmc: do we need a "BARRIER" at this point ?
102  c     _BARRIER  c     _BARRIER
103    
104    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
105    
106          DO n = 1,diagSt_nbLists
107    
108           IF ( diagSt_freq(n).LT.0. ) THEN
109    C--    Select diagnostics list that uses instantaneous output
110    
111            dBugFlag = debugLevel.GT.debLevB
112    
113            freqSec = diagSt_freq(n)
114            phiSec = diagSt_phase(n)
115            time4SnapShot = DIFF_PHASE_MULTIPLE( phiSec, freqSec,
116         &                                       myTime, deltaTclock )
117    
118    #ifdef ALLOW_FIZHI
119            if( useFIZHI) then
120             write(tagname,'(A,I2.2)')'diagStg',n
121             time4SnapShot = alarm2(tagname)
122            endif
123    #endif
124    
125            DO bj=myByLo(myThid), myByHi(myThid)
126             DO bi=myBxLo(myThid), myBxHi(myThid)
127              dBugFlag = dBugFlag.AND.(bi.EQ.1.AND.bj.EQ.1.AND.myThid.EQ.1)
128              IF ( time4SnapShot ) THEN
129    C--     switch ON diagnostics of output-stream # n
130                DO m=1,diagSt_nbActv(n)
131                  j = jSdiag(m,n)
132                  iSp = iSdiag(j)
133                  IF (dBugFlag.AND.qSdiag(0,0,iSp,bi,bj).NE. 0.)
134         &          WRITE(0,'(I8,I4,3A,1PE10.3,A,I3)') myIter,
135         &            j,' ',cdiag(j),' :',qSdiag(0,0,iSp,bi,bj),' ->', 0
136                  qSdiag(0,0,iSp,bi,bj) = 0.
137                ENDDO
138              ELSE
139    C--     switch OFF diagnostics of output-stream # n
140                DO m=1,diagSt_nbActv(n)
141                  j = jSdiag(m,n)
142                  iSp = iSdiag(j)
143                  IF (dBugFlag.AND.qSdiag(0,0,iSp,bi,bj).NE.-1.)
144         &          WRITE(0,'(I8,I4,3A,1PE10.3,A,I3)') myIter,
145         &            j,' ',cdiag(j),' :',qSdiag(0,0,iSp,bi,bj),' ->',-1
146                  qSdiag(0,0,iSp,bi,bj) = -1.
147                ENDDO
148              ENDIF
149             ENDDO
150            ENDDO
151    
152           ENDIF
153          ENDDO
154    
155        RETURN        RETURN
156        END        END
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22