/[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.8 by jmc, Fri May 20 07:28:50 2005 UTC revision 1.11 by heimbach, Sat Nov 4 14:28:09 2006 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
# Line 13  C     !INTERFACE: Line 13  C     !INTERFACE:
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
19  #include "SIZE.h"  #include "SIZE.h"
# Line 27  C     !INPUT PARAMETERS: Line 27  C     !INPUT PARAMETERS:
27  CEOP  CEOP
28    
29  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
30        INTEGER i,j,k,n,bi,bj        INTEGER i,j,k,bi,bj
31    
32  C     Zero out the qdiag array which accumulates during integration  C     Zero out the qdiag array which accumulates during integration
33        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
# Line 38  C     Zero out the qdiag array which acc Line 38  C     Zero out the qdiag array which acc
38                  qdiag(i,j,k,bi,bj) = 0. _d 0                  qdiag(i,j,k,bi,bj) = 0. _d 0
39                ENDDO                ENDDO
40              ENDDO              ENDDO
41    C     Zero out the counters for the qdiag array
42                ndiag(k,bi,bj) = 0
43              ENDDO
44              DO k = 1,numlists
45    C     Zero out the index array for periodic averaging diagnostic
46                pdiag(k,bi,bj) = 0
47            ENDDO            ENDDO
48          ENDDO          ENDDO
49        ENDDO        ENDDO
50    
 C     Zero out the counters for the qdiag array  
       do n=1,ndiagMax  
         ndiag(n) = 0  
       enddo  
   
51  C     Zero out the qSdiag array (statistics) which accumulates during integration  C     Zero out the qSdiag array (statistics) which accumulates during integration
52        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
53          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
54            DO k = 1,diagSt_size            DO k = 1,diagSt_size
55              DO j = 0,nRegions              DO j = 0,nRegions
56                DO i = 0,nStats                DO i = 0,nStats
57                  qSdiag(i,j,k,bi,bj) = 0. _d 0                  qSdiag(i,j,k,bi,bj) = 0. _d 0
# Line 62  C     Zero out the qSdiag array (statist Line 63  C     Zero out the qSdiag array (statist
63    
64        CALL DIAGNOSTICS_READ_PICKUP( myThid )        CALL DIAGNOSTICS_READ_PICKUP( myThid )
65    
66        CALL DIAGSTATS_INI_IO( myThid )        CALL DIAGNOSTICS_SUMMARY( startTime, nIter0, myThid )
67    
68        RETURN        RETURN
69        END        END

Legend:
Removed from v.1.8  
changed lines
  Added in v.1.11

  ViewVC Help
Powered by ViewVC 1.1.22