/[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.1 by jmc, Fri May 20 07:28:50 2005 UTC revision 1.5 by jmc, Tue Feb 5 15:31:19 2008 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
       CHARACTER*(MAX_LEN_MBUF) msgBuf  
44    
45  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
46    
# Line 62  C---+----1----+----2----+----3----+----4 Line 61  C---+----1----+----2----+----3----+----4
61         DO j=0,nRegions         DO j=0,nRegions
62          IF ( diagSt_region(j,listId).GE.1 ) THEN          IF ( diagSt_region(j,listId).GE.1 ) THEN
63           WRITE(nUnit,'(3A,I10,2(A,I4))') ' field : ', cdiag(ndId),           WRITE(nUnit,'(3A,I10,2(A,I4))') ' field : ', cdiag(ndId),
64       &     ' ; Iter =', myIter,' ; region #', j, ' ; nb.Lev =', klev         &    ' ; Iter =',myIter,' ; region #',j, ' ; nb.Lev =',kdiag(ndId)
65  c        WRITE(nUnit,'(5A)')       ' k |',  c        WRITE(nUnit,'(5A)')       ' k |',
66  c    &             ' --   Average  --  |', ' --   Std.Dev   -- |',  c    &             ' --   Average  --  |', ' --   Std.Dev   -- |',
67  c    &             ' --      min   --  |', ' --      max    -- |'  c    &             ' --      min   --  |', ' --      max    -- |'
68           WRITE(nUnit,'(6A)')       ' k |',           WRITE(nUnit,'(6A)')       ' k |',
69       &             ' --  Average --  |', ' --  Std.Dev  -- |',       &             ' --  Average --  |', ' --  Std.Dev  -- |',
70       &             ' --     min  --  |', ' --     max   -- |',' Vol'       &             ' --     min  --  |', ' --     max   -- |',' Vol'
71           DO k=0,klev           DO k=0,klev
72  C     full precision, do not write the volume:  C     full precision, do not write the volume:
73  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)
74  C     reduced precision + write the volume (usefull for testing):  C     reduced precision + write the volume (usefull for testing):
75            WRITE(nUnit,'(I3,1P4E18.10,1X,1PE12.5)')            WRITE(nUnit,'(I3,1P5E18.10)')
76       &       k, (statGlob(i,k,j),i=1,nStats), statGlob(0,k,j)       &       k, (statGlob(i,k,j),i=1,nStats), statGlob(0,k,j)
77           ENDDO           ENDDO
78          ENDIF          ENDIF
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22