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

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

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

revision 1.1 by jmc, Fri May 20 07:28:52 2005 UTC revision 1.2 by jmc, Sun Jun 26 16:51:49 2005 UTC
# Line 11  C     !INTERFACE: Line 11  C     !INTERFACE:
11        SUBROUTINE DIAGSTATS_GLOBAL(        SUBROUTINE DIAGSTATS_GLOBAL(
12       O                    qtmp1, qtmp2,       O                    qtmp1, qtmp2,
13       I                    undef, nLev, jReg,       I                    undef, nLev, jReg,
14       I                    ndId, mate, myThid )       I                    ndId, mate, iSp, iSm, myThid )
15    
16  C     !DESCRIPTION:  C     !DESCRIPTION:
17  C     Retrieve averaged model diagnostic  C     Retrieve averaged model diagnostic
# Line 24  C     !USES: Line 24  C     !USES:
24  #include "DIAGNOSTICS.h"  #include "DIAGNOSTICS.h"
25    
26  C     !INPUT PARAMETERS:  C     !INPUT PARAMETERS:
27  C     undef  ..... Undefined value  C     undef     :: Undefined value
28  C     nLev   ....  2nd Dimension (max Nb of levels) of qtmp1,2 arrays  C     nLev      :: 2nd Dimension (max Nb of levels) of qtmp1,2 arrays
29  C     jReg   ..... region Index to be process.  C     jReg      :: region Index to be process.
30  C     ndId     ... diagnostic Id number (in diagnostics long list)  C     ndId      :: diagnostic Id number (in available diagnostics list)
31  C     mate   ..... counter diagnostic Id number if any ; 0 otherwise  C     mate      :: counter mate Id number if any ; 0 otherwise
32  C     myThid ..... my thread Id number  C     iSp       :: diagnostics  pointer to storage array
33    C     iSm       :: counter-mate pointer to storage array
34    C     myThid    :: my thread Id number
35        _RL undef        _RL undef
36        INTEGER nLev, jReg, ndId, mate        INTEGER nLev, jReg, ndId, mate, iSp, iSm
37        INTEGER myThid        INTEGER myThid
38    
39  C     !OUTPUT PARAMETERS:  C     !OUTPUT PARAMETERS:
# Line 68  C--     Retrieve tile statistics first Line 70  C--     Retrieve tile statistics first
70           DO bi=myBxLo(myThid),myBxHi(myThid)           DO bi=myBxLo(myThid),myBxHi(myThid)
71    
72            DO k=1,klev            DO k=1,klev
73              kd = iSdiag(ndId) + k - 1              kd = iSp + k - 1
74              IF ( qSdiag(0,jReg,kd,bi,bj).GT.0. ) THEN              IF ( qSdiag(0,jReg,kd,bi,bj).GT.0. ) THEN
75               IF ( qtmp1(0,k).LE.0. ) THEN               IF ( qtmp1(0,k).LE.0. ) THEN
76                DO i=0,nStats                DO i=0,nStats
# Line 85  C--     Retrieve tile statistics first Line 87  C--     Retrieve tile statistics first
87            ENDDO            ENDDO
88            IF ( mate.GT.0 ) THEN            IF ( mate.GT.0 ) THEN
89             DO k=1,kMlev             DO k=1,kMlev
90              kd = iSdiag(mate) + k - 1              kd = iSm + k - 1
91              IF ( qSdiag(0,jReg,kd,bi,bj).GT.0. ) THEN              IF ( qSdiag(0,jReg,kd,bi,bj).GT.0. ) THEN
92               IF ( qtmp1(0,k).LE.0. ) THEN               IF ( qtmp1(0,k).LE.0. ) THEN
93                DO i=0,nStats                DO i=0,nStats
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22