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

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

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

revision 1.33 by jmc, Mon Jun 6 15:42:58 2011 UTC revision 1.39 by jmc, Tue Jun 2 20:58:22 2015 UTC
# Line 30  C     !INPUT PARAMETERS: Line 30  C     !INPUT PARAMETERS:
30        _RL     myTime        _RL     myTime
31        INTEGER myIter, myThid        INTEGER myIter, myThid
32    
33    C     !FUNCTIONS:
34          LOGICAL  DIFF_PHASE_MULTIPLE
35          EXTERNAL DIFF_PHASE_MULTIPLE
36    #ifdef ALLOW_FIZHI
37          LOGICAL  ALARM2
38          EXTERNAL ALARM2
39    #endif
40    
41  c Local variables  c Local variables
42  c ===============  c ===============
43        INTEGER   n        INTEGER   n
# Line 37  c =============== Line 45  c ===============
45        LOGICAL   dump2fileNow, write2file        LOGICAL   dump2fileNow, write2file
46        _RL       phiSec, freqSec, wrTime        _RL       phiSec, freqSec, wrTime
47  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
       LOGICAL alarm2  
48        CHARACTER *9 tagname        CHARACTER *9 tagname
49  #endif  #endif
50    
       LOGICAL  DIFF_PHASE_MULTIPLE  
       EXTERNAL DIFF_PHASE_MULTIPLE  
   
51        IF ( myIter.NE.nIter0 ) THEN        IF ( myIter.NE.nIter0 ) THEN
52          myItM1 = myIter - 1          myItM1 = myIter - 1
53    
# Line 64  c           wrTime = myTime Line 68  c           wrTime = myTime
68  C--     write snap-shot with suffix = myIter-1 to be consistent with  C--     write snap-shot with suffix = myIter-1 to be consistent with
69  C       state-variable time-step:  C       state-variable time-step:
70              wrIter = myItM1              wrIter = myItM1
71              wrTime = myTime - deltaTclock              wrTime = myTime - deltaTClock
72            ELSE            ELSE
73              wrIter = myIter              wrIter = myIter
74              wrTime = myTime              wrTime = myTime
75            ENDIF            ENDIF
76            dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,            dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,
77       &                                        wrTime, deltaTclock )       &                                        wrTime, deltaTClock )
78  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
79            IF ( useFIZHI ) THEN            IF ( useFIZHI ) THEN
80             WRITE(tagname,'(A,I2.2)')'diagtag',n             WRITE(tagname,'(A,I2.2)')'diagtag',n
81             dump2fileNow = alarm2(tagname)             dump2fileNow = ALARM2(tagname)
82            ENDIF            ENDIF
83  #endif  #endif
84  #ifdef ALLOW_CAL  #ifdef ALLOW_CAL
85            IF ( useCAL ) THEN            IF ( useCAL ) THEN
86              CALL CAL_TIME2DUMP( freqSec, deltaTClock,              CALL CAL_TIME2DUMP( phiSec, freqSec, deltaTClock,
87       U                          dump2fileNow,       U                          dump2fileNow,
88       I                          myTime, myIter, myThid )       I                          wrTime, myIter, myThid )
89            ENDIF            ENDIF
90  #endif /* ALLOW_CAL */  #endif /* ALLOW_CAL */
91            IF ( dumpAtLast .AND. modelEnd            IF ( dumpAtLast .AND. modelEnd
92       &                    .AND. freqSec.GE.0. ) dump2fileNow = .TRUE.       &                    .AND. freqSec.GE.0. ) dump2fileNow = .TRUE.
93            IF ( dump2fileNow ) THEN            IF ( dump2fileNow ) THEN
94              write2file = .TRUE.              write2file = .TRUE.
95              CALL DIAGNOSTICS_OUT(n,wrIter,wrTime,myThid)              CALL DIAGNOSTICS_OUT(n,wrTime,wrIter,myThid)
96            ENDIF            ENDIF
97          ENDDO          ENDDO
98    
# Line 106  c           wrTime = myTime Line 110  c           wrTime = myTime
110  C--     write snap-shot with suffix = myIter-1 to be consistent with  C--     write snap-shot with suffix = myIter-1 to be consistent with
111  C       state-variable time-step:  C       state-variable time-step:
112              wrIter = myItM1              wrIter = myItM1
113              wrTime = myTime - deltaTclock              wrTime = myTime - deltaTClock
114            ELSE            ELSE
115              wrIter = myIter              wrIter = myIter
116              wrTime = myTime              wrTime = myTime
117            ENDIF            ENDIF
118            dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,            dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,
119       &                                        wrTime, deltaTclock )       &                                        wrTime, deltaTClock )
120  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
121            IF ( useFIZHI ) THEN            IF ( useFIZHI ) THEN
122             WRITE(tagname,'(A,I2.2)')'diagStg',n             WRITE(tagname,'(A,I2.2)')'diagStg',n
123             dump2fileNow = alarm2(tagname)             dump2fileNow = ALARM2(tagname)
124            ENDIF            ENDIF
125  #endif  #endif
126    #ifdef ALLOW_CAL
127              IF ( useCAL ) THEN
128                CALL CAL_TIME2DUMP( phiSec, freqSec, deltaTClock,
129         U                          dump2fileNow,
130         I                          wrTime, myIter, myThid )
131              ENDIF
132    #endif /* ALLOW_CAL */
133            IF ( dumpAtLast .AND. modelEnd            IF ( dumpAtLast .AND. modelEnd
134       &                    .AND. freqSec.GE.0. ) dump2fileNow = .TRUE.       &                    .AND. freqSec.GE.0. ) dump2fileNow = .TRUE.
135            IF ( dump2fileNow ) THEN            IF ( dump2fileNow ) THEN
# Line 137  C-      wait for everyone before setting Line 148  C-      wait for everyone before setting
148            _BARRIER            _BARRIER
149          ENDIF          ENDIF
150          IF ( modelEnd ) THEN          IF ( modelEnd ) THEN
151    C-      Track diagnostics pkg activation status:
152    c         IF ( diag_pkgStatus.NE.ready2fillDiags ) STOP
153              _BARRIER
154              _BEGIN_MASTER(myThid)
155              diag_pkgStatus = 99
156              _END_MASTER(myThid)
157              _BARRIER
158  C       Close all Stat-diags output files  C       Close all Stat-diags output files
159            CALL DIAGSTATS_CLOSE_IO( myThid )            CALL DIAGSTATS_CLOSE_IO( myThid )
160          ENDIF          ENDIF
# Line 148  C--     Clear storage space: Line 166  C--     Clear storage space:
166            phiSec = phase(n)            phiSec = phase(n)
167    
168            wrTime = myTime            wrTime = myTime
169            IF ( freqSec.LT.0. ) wrTime = myTime - deltaTclock            IF ( freqSec.LT.0. ) wrTime = myTime - deltaTClock
170            dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,            dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,
171       &                                        wrTime, deltaTclock )       &                                        wrTime, deltaTClock )
172  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
173            IF ( useFIZHI ) THEN            IF ( useFIZHI ) THEN
174             WRITE(tagname,'(A,I2.2)')'diagtag',n             WRITE(tagname,'(A,I2.2)')'diagtag',n
175             dump2fileNow = alarm2(tagname)             dump2fileNow = ALARM2(tagname)
176            ENDIF            ENDIF
177  #endif  #endif
178  #ifdef ALLOW_CAL  #ifdef ALLOW_CAL
179            IF ( useCAL ) THEN            IF ( useCAL ) THEN
180              CALL CAL_TIME2DUMP( freqSec, deltaTClock,              CALL CAL_TIME2DUMP( phiSec, freqSec, deltaTClock,
181       U                          dump2fileNow,       U                          dump2fileNow,
182       I                          myTime, myIter, myThid )       I                          wrTime, myIter, myThid )
183            ENDIF            ENDIF
184  #endif /* ALLOW_CAL */  #endif /* ALLOW_CAL */
185            IF ( dumpAtLast .AND. modelEnd            IF ( dumpAtLast .AND. modelEnd
# Line 173  C--     Clear storage space: Line 191  C--     Clear storage space:
191            freqSec = diagSt_freq(n)            freqSec = diagSt_freq(n)
192            phiSec = diagSt_phase(n)            phiSec = diagSt_phase(n)
193            wrTime = myTime            wrTime = myTime
194            IF ( freqSec.LT.0. ) wrTime = myTime - deltaTclock            IF ( freqSec.LT.0. ) wrTime = myTime - deltaTClock
195            dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,            dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,
196       &                                        wrTime, deltaTclock )       &                                        wrTime, deltaTClock )
197  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
198            IF ( useFIZHI ) THEN            IF ( useFIZHI ) THEN
199             WRITE(tagname,'(A,I2.2)')'diagStg',n             WRITE(tagname,'(A,I2.2)')'diagStg',n
200             dump2fileNow = alarm2(tagname)             dump2fileNow = ALARM2(tagname)
201            ENDIF            ENDIF
202  #endif  #endif
203    #ifdef ALLOW_CAL
204              IF ( useCAL ) THEN
205                CALL CAL_TIME2DUMP( phiSec, freqSec, deltaTClock,
206         U                          dump2fileNow,
207         I                          wrTime, myIter, myThid )
208              ENDIF
209    #endif /* ALLOW_CAL */
210            IF ( dumpAtLast .AND. modelEnd            IF ( dumpAtLast .AND. modelEnd
211       &                    .AND. freqSec.GE.0. ) dump2fileNow = .TRUE.       &                    .AND. freqSec.GE.0. ) dump2fileNow = .TRUE.
212            IF ( dump2fileNow ) CALL DIAGSTATS_CLEAR( n, myThid )            IF ( dump2fileNow ) CALL DIAGSTATS_CLEAR( n, myThid )
Legend:
Removed from v.1.33  
changed lines
  Added in v.1.39
  ViewVC Help
Powered by ViewVC 1.1.22