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-|--+----|
acacc28f7f Jean*0004 CBOP
09ceb40cd6 Jean*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
acacc28f7f Jean*0036 C--   Set I/O unit and open diag-stats ASCII output file:
                0037       CALL DIAGSTATS_INI_IO( myThid )
                0038
ab43bc12c4 Jean*0039 C--   Zero out the qdiag array which accumulates during integration
09ceb40cd6 Jean*0040       DO bj = myByLo(myThid), myByHi(myThid)
                0041         DO bi = myBxLo(myThid), myBxHi(myThid)
e129400813 Jean*0042           DO k = 1,numDiags
ab43bc12c4 Jean*0043             DO j = 1-OLy,sNy+OLy
                0044               DO i = 1-OLx,sNx+OLx
3e5de6a370 Jean*0045                 qdiag(i,j,k,bi,bj) = 0. _d 0
09ceb40cd6 Jean*0046               ENDDO
                0047             ENDDO
ab43bc12c4 Jean*0048 C-    Zero out the counters for the qdiag array
3ae5f90260 Jean*0049             ndiag(k,bi,bj) = 0
09ceb40cd6 Jean*0050           ENDDO
ab43bc12c4 Jean*0051           DO k = 1,numLists
                0052 C-    Zero out the index array for periodic averaging diagnostic
666b944083 Jean*0053             pdiag(k,bi,bj) = 0
                0054           ENDDO
09ceb40cd6 Jean*0055         ENDDO
                0056       ENDDO
                0057
ab43bc12c4 Jean*0058 C--   Zero out the qSdiag array (statistics) which accumulates during integration
3e5de6a370 Jean*0059       DO bj = myByLo(myThid), myByHi(myThid)
                0060         DO bi = myBxLo(myThid), myBxHi(myThid)
3ae5f90260 Jean*0061           DO k = 1,diagSt_size
3e5de6a370 Jean*0062             DO j = 0,nRegions
                0063               DO i = 0,nStats
                0064                 qSdiag(i,j,k,bi,bj) = 0. _d 0
                0065               ENDDO
                0066             ENDDO
                0067           ENDDO
                0068         ENDDO
                0069       ENDDO
                0070
987bbc7d38 Ed H*0071       CALL DIAGNOSTICS_READ_PICKUP( myThid )
e797dd72d4 Ed H*0072
41c4545f8f Jean*0073       CALL DIAGNOSTICS_SUMMARY( 0, startTime, nIter0, myThid )
3ae5f90260 Jean*0074
                0075       RETURN
09ceb40cd6 Jean*0076       END