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

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

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

revision 1.6 by jmc, Mon Jul 14 23:26:03 2014 UTC revision 1.7 by jmc, Tue Dec 8 19:18:31 2015 UTC
# Line 40  CEOP Line 40  CEOP
40  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
41        INTEGER im, ix, iv        INTEGER im, ix, iv
42        PARAMETER ( iv = nStats - 2 , im = nStats - 1 , ix = nStats )        PARAMETER ( iv = nStats - 2 , im = nStats - 1 , ix = nStats )
43        INTEGER i, j, k, klev, nUnit        INTEGER i, j, k, kLev, nUnit
44    
45  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
46    
# Line 53  C---+----1----+----2----+----3----+----4 Line 53  C---+----1----+----2----+----3----+----4
53  #endif  #endif
54    
55         nUnit = diagSt_ioUnit(listId)         nUnit = diagSt_ioUnit(listId)
56         klev = kdiag(ndId)         kLev = kdiag(ndId)
57  C-     single level field: Vertical Integral (k=0) & 1rst level are identical  C-     single level field: Vertical Integral (k=0) & 1rst level are identical
58  C      => write only 1 of the 2:  C      => write only 1 of the 2:
59         IF ( klev.EQ.1 ) kLev = 0         IF ( kLev.EQ.1 ) kLev = 0
60           kLev = MIN( kLev, nLev )
61  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
62         DO j=0,nRegions         DO j=0,nRegions
63          IF ( diagSt_region(j,listId).GE.1 ) THEN          IF ( diagSt_region(j,listId).GE.1 ) THEN
# Line 71  c    &             ' --    min    -- |', Line 72  c    &             ' --    min    -- |',
72           WRITE(nUnit,'(6A)') ' k |',           WRITE(nUnit,'(6A)') ' k |',
73       &     ' --   Average    -- |', ' --   Std.Dev    -- |',       &     ' --   Average    -- |', ' --   Std.Dev    -- |',
74       &     ' --     min      -- |', ' --     max      -- |',' --    Vol'       &     ' --     min      -- |', ' --     max      -- |',' --    Vol'
75           DO k=0,klev           DO k=0,kLev
76  C     full precision, do not write the volume:  C     full precision, do not write the volume:
77  c         WRITE(nUnit,'(I3,1P4E20.12)') k,(statGlob(i,k,j),i=1,nStats)  c         WRITE(nUnit,'(I3,1P4E20.12)') k,(statGlob(i,k,j),i=1,nStats)
78  C     reduced precision + write the volume (usefull for testing):  C     reduced precision + write the volume (usefull for testing):
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22