pkg/diagnostics/diagnostics_init_varia.F

C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_init_varia.F,v 1.7 2005/05/06 20:34:16 molod Exp $
C $Name:  $

#include "DIAG_OPTIONS.h"
      
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
CBOP 0
C     !ROUTINE: DIAGNOSTICS_INIT_VARIA

C     !INTERFACE:
      SUBROUTINE DIAGNOSTICS_INIT_VARIA(
     I     myThid )

C     !DESCRIPTION:
C     Initialize the qdiag array which accumulates during integration
      
C     !USES:
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DIAGNOSTICS_SIZE.h"
#include "DIAGNOSTICS.h"

C     !INPUT PARAMETERS:
      INTEGER myThid
CEOP

C     !LOCAL VARIABLES:
      INTEGER i,j,k,n,bi,bj

C     Zero out the qdiag array which accumulates during integration
      DO bj = myByLo(myThid), myByHi(myThid)
        DO bi = myBxLo(myThid), myBxHi(myThid)
          DO k = 1,numdiags
            DO j = 1-Oly,sNy+Oly
              DO i = 1-Olx,sNx+Olx
                qdiag(i,j,k,bi,bj) = 0. _d 0
              ENDDO
            ENDDO
          ENDDO
        ENDDO
      ENDDO

C     Zero out the counters for the qdiag array
      do n=1,ndiagMax
        ndiag(n) = 0
      enddo

C     Zero out the qSdiag array (statistics) which accumulates during integration
      DO bj = myByLo(myThid), myByHi(myThid)
        DO bi = myBxLo(myThid), myBxHi(myThid)
          DO k = 1,diagSt_size 
            DO j = 0,nRegions
              DO i = 0,nStats
                qSdiag(i,j,k,bi,bj) = 0. _d 0
              ENDDO
            ENDDO
          ENDDO
        ENDDO
      ENDDO

      CALL DIAGNOSTICS_READ_PICKUP( myThid )

      CALL DIAGSTATS_INI_IO( myThid )

      RETURN 
      END
1	jmc	1.8	C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_init_varia.F,v 1.7 2005/05/06 20:34:16 molod Exp $
2	jmc	1.1	C $Name: $
3
4			#include "DIAG_OPTIONS.h"
5
6			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
7			CBOP 0
8			C !ROUTINE: DIAGNOSTICS_INIT_VARIA
9
10			C !INTERFACE:
11	edhill	1.2	SUBROUTINE DIAGNOSTICS_INIT_VARIA(
12			I myThid )
13	jmc	1.1
14			C !DESCRIPTION:
15			C Initialize the qdiag array which accumulates during integration
16
17			C !USES:
18			IMPLICIT NONE
19	edhill	1.2	#include "SIZE.h"
20	jmc	1.1	#include "EEPARAMS.h"
21	edhill	1.2	#include "PARAMS.h"
22	jmc	1.1	#include "DIAGNOSTICS_SIZE.h"
23			#include "DIAGNOSTICS.h"
24
25			C !INPUT PARAMETERS:
26			INTEGER myThid
27			CEOP
28
29			C !LOCAL VARIABLES:
30	jmc	1.8	INTEGER i,j,k,n,bi,bj
31	edhill	1.4
32	jmc	1.1	C Zero out the qdiag array which accumulates during integration
33			DO bj = myByLo(myThid), myByHi(myThid)
34			DO bi = myBxLo(myThid), myBxHi(myThid)
35	jmc	1.8	DO k = 1,numdiags
36	jmc	1.1	DO j = 1-Oly,sNy+Oly
37			DO i = 1-Olx,sNx+Olx
38	jmc	1.8	qdiag(i,j,k,bi,bj) = 0. _d 0
39	jmc	1.1	ENDDO
40			ENDDO
41			ENDDO
42			ENDDO
43			ENDDO
44
45	molod	1.7	C Zero out the counters for the qdiag array
46			do n=1,ndiagMax
47			ndiag(n) = 0
48			enddo
49
50	jmc	1.8	C Zero out the qSdiag array (statistics) which accumulates during integration
51			DO bj = myByLo(myThid), myByHi(myThid)
52			DO bi = myBxLo(myThid), myBxHi(myThid)
53			DO k = 1,diagSt_size
54			DO j = 0,nRegions
55			DO i = 0,nStats
56			qSdiag(i,j,k,bi,bj) = 0. _d 0
57			ENDDO
58			ENDDO
59			ENDDO
60			ENDDO
61			ENDDO
62
63	edhill	1.6	CALL DIAGNOSTICS_READ_PICKUP( myThid )
64	edhill	1.5
65	jmc	1.8	CALL DIAGSTATS_INI_IO( myThid )
66
67	jmc	1.1	RETURN
68			END