/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/dic_diags.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/darwin/dic_diags.F

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

revision 1.6 by jahn, Wed Dec 4 21:19:12 2013 UTC revision 1.7 by jahn, Fri Dec 27 17:29:00 2013 UTC
# Line 2  Line 2 
2  #include "PTRACERS_OPTIONS.h"  #include "PTRACERS_OPTIONS.h"
3  #include "DARWIN_OPTIONS.h"  #include "DARWIN_OPTIONS.h"
4    
 #ifdef ALLOW_PTRACERS  
 #ifdef ALLOW_DARWIN  
   
 #ifdef ALLOW_CARBON  
   
5  CBOP  CBOP
 CStartOfInterface  
6        SUBROUTINE DIC_DIAGS( myTime, myIter, myThid )        SUBROUTINE DIC_DIAGS( myTime, myIter, myThid )
7  C     /==========================================================\  C     /==========================================================\
8  C     | SUBROUTINE DIC_BIOTIC_DIAGS                                  |  C     | SUBROUTINE DIC_BIOTIC_DIAGS                              |
9  C     | o Do DIC diagnostic output                           |  C     | o Do DIC diagnostic output                               |
10  C     |==========================================================|  C     |==========================================================|
11  C     \==========================================================/  C     \==========================================================/
12        IMPLICIT NONE        IMPLICIT NONE
# Line 21  C     === Global variables === Line 15  C     === Global variables ===
15  #include "SIZE.h"  #include "SIZE.h"
16  #include "EEPARAMS.h"  #include "EEPARAMS.h"
17  #include "PARAMS.h"  #include "PARAMS.h"
18    #ifdef ALLOW_DARWIN
19  #include "DARWIN_SIZE.h"  #include "DARWIN_SIZE.h"
20  #include "DARWIN_FLUX.h"  #include "DARWIN_FLUX.h"
21  #include "DARWIN_IO.h"  #include "DARWIN_IO.h"
22    #endif
23    #ifdef ALLOW_PTRACERS
24  #include "PTRACERS_SIZE.h"  #include "PTRACERS_SIZE.h"
25  #include "PTRACERS_PARAMS.h"  #include "PTRACERS_PARAMS.h"
26    #endif
27  C     == Routine arguments ==  C     == Routine arguments ==
28  C     myTime - Current time of simulation ( s )  C     myTime - Current time of simulation ( s )
29  C     myIter - Iteration number  C     myIter - Iteration number
# Line 33  C     myThid -  Number of this instance Line 31  C     myThid -  Number of this instance
31        _RL     myTime        _RL     myTime
32        INTEGER myIter        INTEGER myIter
33        INTEGER myThid        INTEGER myThid
34  CEndOfInterface  CEOP
35    
36    #ifdef ALLOW_PTRACERS
37    #ifdef ALLOW_DARWIN
38    #ifdef ALLOW_CARBON
39    #ifdef ALLOW_TIMEAVE
40    
41  C     == Local variables ==  C     == Local variables ==
42        LOGICAL  DIFFERENT_MULTIPLE        LOGICAL  DIFFERENT_MULTIPLE
# Line 50  C     == Local variables == Line 52  C     == Local variables ==
52          pf(1:1) = 'R'          pf(1:1) = 'R'
53        ENDIF        ENDIF
54    
 #ifdef ALLOW_TIMEAVE  
   
55  C     Dump files and restart average computation if needed  C     Dump files and restart average computation if needed
56        IF ( DIFFERENT_MULTIPLE(taveFreq,myTime,deltaTClock).        IF ( DIFFERENT_MULTIPLE(taveFreq,myTime,deltaTClock).
57       &                               and.myIter.ne.nIter0 ) THEN       &                               and.myIter.ne.nIter0 ) THEN
# Line 59  C     Dump files and restart average com Line 59  C     Dump files and restart average com
59  C      Normalize by integrated time  C      Normalize by integrated time
60         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
61          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
62           CALL TIMEAVE_NORMALIZ(SURave,dic_timeave, 1 ,           CALL TIMEAVE_NORMALIZE(SURave,dic_timeave, 1 ,
63       &                 bi,bj,myThid)       &                 bi,bj,myThid)
64           CALL TIMEAVE_NORMALIZ(SURCave,dic_timeave, 1 ,           CALL TIMEAVE_NORMALIZE(SURCave,dic_timeave, 1 ,
65       &                 bi,bj,myThid)       &                 bi,bj,myThid)
66           CALL TIMEAVE_NORMALIZ(SUROave,dic_timeave, 1 ,           CALL TIMEAVE_NORMALIZE(SUROave,dic_timeave, 1 ,
67       &                 bi,bj,myThid)       &                 bi,bj,myThid)
68  #ifdef pH_3D  #ifdef pH_3D
69           CALL TIMEAVE_NORMALIZ(pCO2ave,dic_timeave, Nr,           CALL TIMEAVE_NORMALIZE(pCO2ave,dic_timeave, Nr,
70       &                 bi,bj,myThid)       &                 bi,bj,myThid)
71           CALL TIMEAVE_NORMALIZ(pHave,dic_timeave, Nr,           CALL TIMEAVE_NORMALIZE(pHave,dic_timeave, Nr,
72       &                 bi,bj,myThid)       &                 bi,bj,myThid)
73  #else  #else
74           CALL TIMEAVE_NORMALIZ(pCO2ave,dic_timeave, 1 ,           CALL TIMEAVE_NORMALIZE(pCO2ave,dic_timeave, 1 ,
75       &                 bi,bj,myThid)       &                 bi,bj,myThid)
76           CALL TIMEAVE_NORMALIZ(pHave,dic_timeave, 1 ,           CALL TIMEAVE_NORMALIZE(pHave,dic_timeave, 1 ,
77       &                 bi,bj,myThid)       &                 bi,bj,myThid)
78  #endif  #endif
79          ENDDO          ENDDO
# Line 144  C      Reset averages to zero Line 144  C      Reset averages to zero
144           CALL TIMEAVE_RESET(pCO2ave,1,bi,bj,myThid)           CALL TIMEAVE_RESET(pCO2ave,1,bi,bj,myThid)
145           CALL TIMEAVE_RESET(pHave,1,bi,bj,myThid)           CALL TIMEAVE_RESET(pHave,1,bi,bj,myThid)
146  #endif  #endif
147           do k=1,Nr           DIC_TimeAve(bi,bj)=0.d0
            DIC_TimeAve(bi,bj,k)=0.d0  
          enddo  
148          ENDDO          ENDDO
149         ENDDO         ENDDO
150    
151        ENDIF        ENDIF
152    
153  #endif /*TIMEAVE*/  #endif /* ALLOW_TIMEAVE */
154    #endif /* ALLOW_CARBON */
155    #endif /* ALLOW_DARWIN */
156    #endif /* ALLOW_PTRACERS */
157    
158        RETURN        RETURN
159        END        END
 #endif  /*ALLOW_CARBON*/  
160    
 #endif  /*DARWIN*/  
 #endif  /*ALLOW_PTRACERS*/  
161  c ==================================================================  c ==================================================================
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22