pkg/diagnostics/diagnostics_init_varia.F

09ceb40cd6 Jean*0001 #include "DIAG_OPTIONS.h"
3ae5f90260 Jean*0002
09ceb40cd6 Jean*0003 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0004 CBOP 0
                0005 C     !ROUTINE: DIAGNOSTICS_INIT_VARIA
                0006
                0007 C     !INTERFACE:
41c4545f8f Jean*0008       SUBROUTINE DIAGNOSTICS_INIT_VARIA( myThid )
09ceb40cd6 Jean*0009
                0010 C     !DESCRIPTION:
                0011 C     Initialize the qdiag array which accumulates during integration
3ae5f90260 Jean*0012
09ceb40cd6 Jean*0013 C     !USES:
                0014       IMPLICIT NONE
                0015 #include "SIZE.h"
ee2e7fad64 Ed H*0016 #include "EEPARAMS.h"
                0017 #include "PARAMS.h"
09ceb40cd6 Jean*0018 #include "DIAGNOSTICS_SIZE.h"
                0019 #include "DIAGNOSTICS.h"
                0020
                0021 C     !INPUT PARAMETERS:
                0022       INTEGER myThid
                0023 CEOP
                0024
                0025 C     !LOCAL VARIABLES:
3ae5f90260 Jean*0026       INTEGER i,j,k,bi,bj
3871a20d99 Ed H*0027
ab43bc12c4 Jean*0028 C--   Track diagnostics pkg activation status:
e01144d3ad Jean*0029       _BARRIER
ab43bc12c4 Jean*0030       _BEGIN_MASTER( myThid )
                0031 c      IF ( diag_pkgStatus.NE.3 ) STOP
                0032        diag_pkgStatus = 10
                0033       _END_MASTER( myThid )
e01144d3ad Jean*0034       _BARRIER
ab43bc12c4 Jean*0035
                0036 C--   Zero out the qdiag array which accumulates during integration
09ceb40cd6 Jean*0037       DO bj = myByLo(myThid), myByHi(myThid)
                0038         DO bi = myBxLo(myThid), myBxHi(myThid)
e129400813 Jean*0039           DO k = 1,numDiags
ab43bc12c4 Jean*0040             DO j = 1-OLy,sNy+OLy
                0041               DO i = 1-OLx,sNx+OLx
3e5de6a370 Jean*0042                 qdiag(i,j,k,bi,bj) = 0. _d 0
09ceb40cd6 Jean*0043               ENDDO
                0044             ENDDO
ab43bc12c4 Jean*0045 C-    Zero out the counters for the qdiag array
3ae5f90260 Jean*0046             ndiag(k,bi,bj) = 0
09ceb40cd6 Jean*0047           ENDDO
ab43bc12c4 Jean*0048           DO k = 1,numLists
                0049 C-    Zero out the index array for periodic averaging diagnostic
666b944083 Jean*0050             pdiag(k,bi,bj) = 0
                0051           ENDDO
09ceb40cd6 Jean*0052         ENDDO
                0053       ENDDO
                0054
ab43bc12c4 Jean*0055 C--   Zero out the qSdiag array (statistics) which accumulates during integration
3e5de6a370 Jean*0056       DO bj = myByLo(myThid), myByHi(myThid)
                0057         DO bi = myBxLo(myThid), myBxHi(myThid)
3ae5f90260 Jean*0058           DO k = 1,diagSt_size
3e5de6a370 Jean*0059             DO j = 0,nRegions
                0060               DO i = 0,nStats
                0061                 qSdiag(i,j,k,bi,bj) = 0. _d 0
                0062               ENDDO
                0063             ENDDO
                0064           ENDDO
                0065         ENDDO
                0066       ENDDO
                0067
987bbc7d38 Ed H*0068       CALL DIAGNOSTICS_READ_PICKUP( myThid )
e797dd72d4 Ed H*0069
41c4545f8f Jean*0070       CALL DIAGNOSTICS_SUMMARY( 0, startTime, nIter0, myThid )
3ae5f90260 Jean*0071
                0072       RETURN
09ceb40cd6 Jean*0073       END