/[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.9 by jmc, Sun Feb 6 23:19:02 2005 UTC revision 1.14 by jmc, Sat May 14 20:45:27 2005 UTC
# Line 21  C*************************************** Line 21  C***************************************
21  #include "DIAGNOSTICS.h"  #include "DIAGNOSTICS.h"
22  #include "EEPARAMS.h"  #include "EEPARAMS.h"
23  #include "PARAMS.h"  #include "PARAMS.h"
24    #ifdef ALLOW_FIZHI
25    #include "chronos.h"
26    #endif
27    
28    C     !INPUT PARAMETERS:
29        _RL     myTime        _RL     myTime
30        INTEGER myIter, myThid        INTEGER myIter, myThid
31    
# Line 30  c =============== Line 34  c ===============
34        INTEGER   n        INTEGER   n
35        INTEGER   myItM1, wrIter        INTEGER   myItM1, wrIter
36        LOGICAL   dump2fileNow        LOGICAL   dump2fileNow
37          _RL       phiSec, freqSec, wrTime
38    #ifdef ALLOW_CAL
39          INTEGER thisdate(4), prevdate(4)
40    #endif
41    #ifdef ALLOW_FIZHI
42          integer nsecf2,yymmdd,hhmmss
43    #endif
44    
45          LOGICAL  DIFF_PHASE_MULTIPLE
46          EXTERNAL DIFF_PHASE_MULTIPLE
47    
48        IF ( myIter.NE.nIter0 ) THEN        IF ( myIter.NE.nIter0 ) THEN
49          myItM1 = myIter - 1          myItM1 = myIter - 1
# Line 38  C*************************************** Line 52  C***************************************
52  C***   Check to see IF its time for Diagnostic Output                ***  C***   Check to see IF its time for Diagnostic Output                ***
53  C***********************************************************************  C***********************************************************************
54    
55    #ifdef ALLOW_CAL
56            IF ( calendarDumps ) THEN
57    C-    Determine calendar dates for this and previous time step.
58               call cal_GetDate(myiter  ,mytime            ,thisdate,mythid)
59               call cal_GetDate(myiter-1,mytime-deltaTClock,prevdate,mythid)
60            ENDIF
61    #endif
62    
63          DO n = 1,nlists          DO n = 1,nlists
64            IF ( freq(n).LT.0 ) THEN            freqSec = freq(n)
65              dump2fileNow = MOD(myItM1,-freq(n)) .EQ. INT(-freq(n)/2)            phiSec = phase(n)
66    #ifdef ALLOW_FIZHI
67             if( useFIZHI) then
68              yymmdd = int(freq(n))
69              hhmmss = int((freq(n) - int(freq(n)))*1.e6)
70              freqSec = nsecf2(hhmmss,yymmdd,nymd)
71              yymmdd = int(phase(n))
72              hhmmss = int((phase(n) - int(phase(n)))*1.e6)
73              phiSec = nsecf2(hhmmss,yymmdd,nymd)
74             endif
75    #endif
76    
77              IF ( freqSec.LT.0. ) THEN
78  C--     write snap-shot with suffix = myIter to be consistent with  C--     write snap-shot with suffix = myIter to be consistent with
79  C       time-average diagnostics (e.g., freq=-1 & freq=1):  C       time-average diagnostics (e.g., freq=-1 & freq=1):
80  c           wrIter = myIter  c           wrIter = myIter
81    c           wrTime = myTime
82  C--     write snap-shot with suffix = myIter-1 to be consistent with  C--     write snap-shot with suffix = myIter-1 to be consistent with
83  C       state-variable time-step:  C       state-variable time-step:
84              wrIter = myItM1              wrIter = myItM1
85                wrTime = myTime - deltaTclock
86            ELSE            ELSE
             dump2fileNow = MOD(myIter,freq(n)) .EQ. 0  
87              wrIter = myIter              wrIter = myIter
88                wrTime = myTime
89              ENDIF
90              dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,
91         &                                        wrTime, deltaTclock )
92    
93    #ifdef ALLOW_CAL
94              IF ( calendarDumps .AND. (
95         &     ( freqSec.GE. 2592000 .AND. freqSec.LE. 2678400 ) .OR.
96         &     ( freqSec.GE.31104000 .AND. freqSec.LE.31968000 ))) THEN
97    C--   Convert approximate months (30-31 days) and years (360-372 days)
98    C     to exact calendar months and years.
99               dump2fileNow = .FALSE.
100    C-    Monthly freqSec:
101               IF( freqSec.GE. 2592000 .AND. freqSec.LE. 2678400 .AND.
102         &        (thisdate(1)-prevdate(1)).GT.50   ) dump2fileNow = .TRUE.
103    C-    Yearly  freqSec:
104               IF( freqSec.GE.31104000 .AND. freqSec.LE.31968000 .AND.
105         &        (thisdate(1)-prevdate(1)).GT.5000 ) dump2fileNow = .TRUE.
106            ENDIF            ENDIF
107    #endif
108    
109            IF ( dump2fileNow ) THEN            IF ( dump2fileNow ) THEN
110              CALL DIAGNOSTICS_OUT(n,wrIter,myThid)              CALL DIAGNOSTICS_OUT(n,wrIter,myThid)
111            ENDIF            ENDIF
# Line 59  C       state-variable time-step: Line 114  C       state-variable time-step:
114  C-      wait for everyone before setting arrays to zero:  C-      wait for everyone before setting arrays to zero:
115          _BARRIER          _BARRIER
116          DO n = 1,nlists          DO n = 1,nlists
117            IF ( freq(n).LT.0 ) THEN            freqSec = freq(n)
118              dump2fileNow = MOD(myItM1,-freq(n)) .EQ. INT(-freq(n)/2)            phiSec = phase(n)
119            ELSE  #ifdef ALLOW_FIZHI
120              dump2fileNow = MOD(myIter,freq(n)) .EQ. 0           if( useFIZHI) then
121            ENDIF            yymmdd = int(freq(n))
122              hhmmss = int((freq(n) - int(freq(n)))*1.e6)
123              freqSec = nsecf2(hhmmss,yymmdd,nymd)
124              yymmdd = int(phase(n))
125              hhmmss = int((phase(n) - int(phase(n)))*1.e6)
126              phiSec = nsecf2(hhmmss,yymmdd,nymd)
127             endif
128    #endif
129              wrTime = myTime
130              IF ( freqSec.LT.0. ) wrTime = myTime - deltaTclock
131              dump2fileNow = DIFF_PHASE_MULTIPLE( phiSec, freqSec,
132         &                                        wrTime, deltaTclock )
133    
134            IF ( dump2fileNow ) CALL CLRINDX(n,myThid)            IF ( dump2fileNow ) CALL CLRINDX(n,myThid)
135          ENDDO          ENDDO
136    
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22