/[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.2 by jmc, Sun Jun 26 16:51:49 2005 UTC revision 1.3 by jmc, Sun Jul 10 00:52:12 2005 UTC
# Line 89  C--     Retrieve tile statistics first Line 89  C--     Retrieve tile statistics first
89             DO k=1,kMlev             DO k=1,kMlev
90              kd = iSm + 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 ( qtmp2(0,k).LE.0. ) THEN
93                DO i=0,nStats                DO i=0,nStats
94                  qtmp2(i,k) = qSdiag(i,jReg,kd,bi,bj)                  qtmp2(i,k) = qSdiag(i,jReg,kd,bi,bj)
95                ENDDO                ENDDO
# Line 178  C            Standard deviation : Line 178  C            Standard deviation :
178               IF (qtmp1(iv,k).GT.0.) qtmp1(iv,k) = SQRT(qtmp1(iv,k))               IF (qtmp1(iv,k).GT.0.) qtmp1(iv,k) = SQRT(qtmp1(iv,k))
179             ENDIF             ENDIF
180            ENDDO            ENDDO
181    C         return global (& vertically integrated) volume in qtmp2(0,0):
182              qtmp2(0,0) = qtmp1(0,0)
183          ELSE          ELSE
184  C       With counter diagnostics => average = Sum / Sum(counter) :  C       With counter diagnostics => average = Sum / Sum(counter) :
185            DO k=0,klev            DO k=0,klev
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22