/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_init_varia.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_init_varia.F

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

revision 1.1 by jmc, Mon Dec 13 21:43:54 2004 UTC revision 1.13 by jmc, Wed Aug 14 00:57:33 2013 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "DIAG_OPTIONS.h"  #include "DIAG_OPTIONS.h"
5          
6  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7  CBOP 0  CBOP 0
8  C     !ROUTINE: DIAGNOSTICS_INIT_VARIA  C     !ROUTINE: DIAGNOSTICS_INIT_VARIA
9    
10  C     !INTERFACE:  C     !INTERFACE:
11        SUBROUTINE DIAGNOSTICS_INIT_VARIA(myThid)        SUBROUTINE DIAGNOSTICS_INIT_VARIA(
12         I     myThid )
13    
14  C     !DESCRIPTION:  C     !DESCRIPTION:
15  C     Initialize the qdiag array which accumulates during integration  C     Initialize the qdiag array which accumulates during integration
16          
17  C     !USES:  C     !USES:
18        IMPLICIT NONE        IMPLICIT NONE
 #include "EEPARAMS.h"  
19  #include "SIZE.h"  #include "SIZE.h"
20    #include "EEPARAMS.h"
21    #include "PARAMS.h"
22  #include "DIAGNOSTICS_SIZE.h"  #include "DIAGNOSTICS_SIZE.h"
23  #include "DIAGNOSTICS.h"  #include "DIAGNOSTICS.h"
24    
# Line 25  C     !INPUT PARAMETERS: Line 27  C     !INPUT PARAMETERS:
27  CEOP  CEOP
28    
29  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
30        INTEGER i,j,n,bi,bj        INTEGER i,j,k,bi,bj
31    
32    C--   Track diagnostics pkg activation status:
33          _BEGIN_MASTER( myThid )
34    c      IF ( diag_pkgStatus.NE.3 ) STOP
35           diag_pkgStatus = 10
36          _END_MASTER( myThid )
37    
38    C--   Zero out the qdiag array which accumulates during integration
39          DO bj = myByLo(myThid), myByHi(myThid)
40            DO bi = myBxLo(myThid), myBxHi(myThid)
41              DO k = 1,numDiags
42                DO j = 1-OLy,sNy+OLy
43                  DO i = 1-OLx,sNx+OLx
44                    qdiag(i,j,k,bi,bj) = 0. _d 0
45                  ENDDO
46                ENDDO
47    C-    Zero out the counters for the qdiag array
48                ndiag(k,bi,bj) = 0
49              ENDDO
50              DO k = 1,numLists
51    C-    Zero out the index array for periodic averaging diagnostic
52                pdiag(k,bi,bj) = 0
53              ENDDO
54            ENDDO
55          ENDDO
56    
57  C     Zero out the qdiag array which accumulates during integration  C--   Zero out the qSdiag array (statistics) which accumulates during integration
58        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
59          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
60            DO n = 1,numdiags            DO k = 1,diagSt_size
61              DO j = 1-Oly,sNy+Oly              DO j = 0,nRegions
62                DO i = 1-Olx,sNx+Olx                DO i = 0,nStats
63                  qdiag(i,j,n,bi,bj) = 0. _d 0                  qSdiag(i,j,k,bi,bj) = 0. _d 0
64                ENDDO                ENDDO
65              ENDDO              ENDDO
66            ENDDO            ENDDO
67          ENDDO          ENDDO
68        ENDDO        ENDDO
69    
70        RETURN        CALL DIAGNOSTICS_READ_PICKUP( myThid )
71    
72          CALL DIAGNOSTICS_SUMMARY( startTime, nIter0, myThid )
73    
74          RETURN
75        END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22