/[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.15 by molod, Mon May 16 23:41:32 2005 UTC revision 1.38 by jmc, Wed Aug 14 00:57:33 2013 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "DIAG_OPTIONS.h"  #include "DIAG_OPTIONS.h"
5    
6        SUBROUTINE DIAGNOSTICS_WRITE ( myTime, myIter, myThid )        SUBROUTINE DIAGNOSTICS_WRITE (
7         I                               modelEnd,
8         I                               myTime, myIter, myThid )
9  C***********************************************************************  C***********************************************************************
10  C  Purpose  C  Purpose
11  C  -------  C  -------
# Line 11  C    Output sequence for the (multiple) Line 13  C    Output sequence for the (multiple)
13  C  C
14  C  Arguments  Description  C  Arguments  Description
15  C  ----------------------  C  ----------------------
16  C     myTime :: Current time of simulation ( s )  C     modelEnd :: true if call at end of model run.
17  C     myIter :: Current Iteration Number  C     myTime   :: Current time of simulation ( s )
18  C     myThid :: my thread Id number  C     myIter   :: Current Iteration Number
19    C     myThid   :: my Thread Id number
20  C***********************************************************************  C***********************************************************************
21         IMPLICIT NONE         IMPLICIT NONE
22    #include "EEPARAMS.h"
23  #include "SIZE.h"  #include "SIZE.h"
24  #include "DIAGNOSTICS_SIZE.h"  #include "DIAGNOSTICS_SIZE.h"
 #include "DIAGNOSTICS.h"  
 #include "EEPARAMS.h"  
25  #include "PARAMS.h"  #include "PARAMS.h"
26  #ifdef ALLOW_FIZHI  #include "DIAGNOSTICS.h"
 #include "chronos.h"  
 #endif  
27    
28  C     !INPUT PARAMETERS:  C     !INPUT PARAMETERS:
29          LOGICAL modelEnd
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
44        INTEGER   myItM1, wrIter        INTEGER   myItM1, wrIter
45        LOGICAL   dump2fileNow        LOGICAL   dump2fileNow, write2file
46        _RL       phiSec, freqSec, wrTime        _RL       phiSec, freqSec, wrTime
 #ifdef ALLOW_CAL  
       INTEGER thisdate(4), prevdate(4)  
 #endif  
47  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
48        integer nsecf2,mmdd,hhmmss,alarm2        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    
54  C***********************************************************************  C***********************************************************************
55  C***   Check to see IF its time for Diagnostic Output                ***  C***   Check to see if its time for Diagnostic Output                ***
56  C***********************************************************************  C***********************************************************************
57    
58  #ifdef ALLOW_CAL          write2file = .FALSE.
         IF ( calendarDumps ) THEN  
 C-    Determine calendar dates for this and previous time step.  
            call cal_GetDate(myiter  ,mytime            ,thisdate,mythid)  
            call cal_GetDate(myiter-1,mytime-deltaTClock,prevdate,mythid)  
         ENDIF  
 #endif  
   
59          DO n = 1,nlists          DO n = 1,nlists
60            freqSec = freq(n)            freqSec = freq(n)
61            phiSec = phase(n)            phiSec = phase(n)
62    
63              IF ( freqSec.LT.0. ) THEN
64    C--     write snap-shot with suffix = myIter to be consistent with
65    C       time-average diagnostics (e.g., freq=-1 & freq=1):
66    c           wrIter = myIter
67    c           wrTime = myTime
68    C--     write snap-shot with suffix = myIter-1 to be consistent with
69    C       state-variable time-step:
70                wrIter = myItM1
71                wrTime = myTime - deltaTClock
72              ELSE
73                wrIter = myIter
74                wrTime = myTime
75              ENDIF
76              dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,
77         &                                        wrTime, deltaTClock )
78  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
79           if( useFIZHI) then            IF ( useFIZHI ) THEN
80            mmdd = int(freq(n))             WRITE(tagname,'(A,I2.2)')'diagtag',n
81            hhmmss = int((freq(n) - int(freq(n)))*1.e6)             dump2fileNow = ALARM2(tagname)
82            freqSec = nsecf2(hhmmss,mmdd,nymd)            ENDIF
          endif  
83  #endif  #endif
84    #ifdef ALLOW_CAL
85              IF ( useCAL ) THEN
86                CALL CAL_TIME2DUMP( phiSec, freqSec, deltaTClock,
87         U                          dump2fileNow,
88         I                          wrTime, myIter, myThid )
89              ENDIF
90    #endif /* ALLOW_CAL */
91              IF ( dumpAtLast .AND. modelEnd
92         &                    .AND. freqSec.GE.0. ) dump2fileNow = .TRUE.
93              IF ( dump2fileNow ) THEN
94                write2file = .TRUE.
95                CALL DIAGNOSTICS_OUT(n,wrTime,wrIter,myThid)
96              ENDIF
97            ENDDO
98    
99    C---   Check to see if its time for Statistics Diag. Output
100    
101            DO n = 1,diagSt_nbLists
102              freqSec = diagSt_freq(n)
103              phiSec = diagSt_phase(n)
104    
105            IF ( freqSec.LT.0. ) THEN            IF ( freqSec.LT.0. ) THEN
106  C--     write snap-shot with suffix = myIter to be consistent with  C--     write snap-shot with suffix = myIter to be consistent with
107  C       time-average diagnostics (e.g., freq=-1 & freq=1):  C       time-average diagnostics (e.g., freq=-1 & freq=1):
108  c           wrIter = myIter  c           wrIter = myIter
109  c           wrTime = myTime  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)')'diagtag',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  #ifdef ALLOW_CAL
127            IF ( calendarDumps .AND. (            IF ( useCAL ) THEN
128       &     ( freqSec.GE. 2592000 .AND. freqSec.LE. 2678400 ) .OR.              CALL CAL_TIME2DUMP( phiSec, freqSec, deltaTClock,
129       &     ( freqSec.GE.31104000 .AND. freqSec.LE.31968000 ))) THEN       U                          dump2fileNow,
130  C--   Convert approximate months (30-31 days) and years (360-372 days)       I                          wrTime, myIter, myThid )
 C     to exact calendar months and years.  
            dump2fileNow = .FALSE.  
 C-    Monthly freqSec:  
            IF( freqSec.GE. 2592000 .AND. freqSec.LE. 2678400 .AND.  
      &        (thisdate(1)-prevdate(1)).GT.50   ) dump2fileNow = .TRUE.  
 C-    Yearly  freqSec:  
            IF( freqSec.GE.31104000 .AND. freqSec.LE.31968000 .AND.  
      &        (thisdate(1)-prevdate(1)).GT.5000 ) dump2fileNow = .TRUE.  
131            ENDIF            ENDIF
132  #endif  #endif /* ALLOW_CAL */
133              IF ( dumpAtLast .AND. modelEnd
134         &                    .AND. freqSec.GE.0. ) dump2fileNow = .TRUE.
135            IF ( dump2fileNow ) THEN            IF ( dump2fileNow ) THEN
136              CALL DIAGNOSTICS_OUT(n,wrIter,myThid)              write2file = .TRUE.
137                CALL DIAGSTATS_OUTPUT(n,wrTime,wrIter,myThid)
138            ENDIF            ENDIF
139          ENDDO          ENDDO
140    
141  C-      wait for everyone before setting arrays to zero:  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
142          _BARRIER  
143            IF ( write2file ) THEN
144              IF ( debugLevel.GE.debLevC ) THEN
145                CALL DIAGNOSTICS_SUMMARY( myTime, myIter, myThid )
146              ENDIF
147    C-      wait for everyone before setting arrays to zero:
148              _BARRIER
149            ENDIF
150            IF ( modelEnd ) THEN
151    C-      Track diagnostics pkg activation status:
152              _BEGIN_MASTER(myThid)
153    c         IF ( diag_pkgStatus.NE.ready2fillDiags ) STOP
154              diag_pkgStatus = 99
155              _END_MASTER(myThid)
156    C       Close all Stat-diags output files
157              CALL DIAGSTATS_CLOSE_IO( myThid )
158            ENDIF
159    
160    C--     Clear storage space:
161    
162          DO n = 1,nlists          DO n = 1,nlists
163            freqSec = freq(n)            freqSec = freq(n)
164            phiSec = phase(n)            phiSec = phase(n)
165  #ifdef ALLOW_FIZHI  
166           if( useFIZHI) then            wrTime = myTime
167            mmdd = int(freq(n))            IF ( freqSec.LT.0. ) wrTime = myTime - deltaTClock
168            hhmmss = int((freq(n) - int(freq(n)))*1.e6)            dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,
169            freqSec = nsecf2(hhmmss,mmdd,nymd)       &                                        wrTime, deltaTClock )
170           endif  #ifdef ALLOW_FIZHI
171              IF ( useFIZHI ) THEN
172               WRITE(tagname,'(A,I2.2)')'diagtag',n
173               dump2fileNow = ALARM2(tagname)
174              ENDIF
175  #endif  #endif
176    #ifdef ALLOW_CAL
177              IF ( useCAL ) THEN
178                CALL CAL_TIME2DUMP( phiSec, freqSec, deltaTClock,
179         U                          dump2fileNow,
180         I                          wrTime, myIter, myThid )
181              ENDIF
182    #endif /* ALLOW_CAL */
183              IF ( dumpAtLast .AND. modelEnd
184         &                    .AND. freqSec.GE.0. ) dump2fileNow = .TRUE.
185              IF ( dump2fileNow ) CALL DIAGNOSTICS_CLEAR(n,myThid)
186            ENDDO
187    
188            DO n = 1,diagSt_nbLists
189              freqSec = diagSt_freq(n)
190              phiSec = diagSt_phase(n)
191            wrTime = myTime            wrTime = myTime
192            IF ( freqSec.LT.0. ) wrTime = myTime - deltaTclock            IF ( freqSec.LT.0. ) wrTime = myTime - deltaTClock
193            dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,            dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,
194       &                                        wrTime, deltaTclock )       &                                        wrTime, deltaTClock )
195  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
196            if( useFIZHI) then            IF ( useFIZHI ) THEN
197             write(tagname,'(A,I2.2)')'diagtag',n             WRITE(tagname,'(A,I2.2)')'diagStg',n
198             dump2fileNow = alarm2(tagname)             dump2fileNow = ALARM2(tagname)
199            endif            ENDIF
200  #endif  #endif
201            IF ( dump2fileNow ) CALL CLRINDX(n,myThid)  #ifdef ALLOW_CAL
202              IF ( useCAL ) THEN
203                CALL CAL_TIME2DUMP( phiSec, freqSec, deltaTClock,
204         U                          dump2fileNow,
205         I                          wrTime, myIter, myThid )
206              ENDIF
207    #endif /* ALLOW_CAL */
208              IF ( dumpAtLast .AND. modelEnd
209         &                    .AND. freqSec.GE.0. ) dump2fileNow = .TRUE.
210              IF ( dump2fileNow ) CALL DIAGSTATS_CLEAR( n, myThid )
211          ENDDO          ENDDO
212    
213  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
214        ENDIF        ENDIF
215    
216        RETURN        RETURN
217        END        END
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.38
  ViewVC Help
Powered by ViewVC 1.1.22