/[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.5 by jmc, Tue Feb 5 15:31:19 2008 UTC revision 1.6 by jmc, Mon Jul 14 23:26:03 2014 UTC
# Line 65  C---+----1----+----2----+----3----+----4 Line 65  C---+----1----+----2----+----3----+----4
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 |',  c        WRITE(nUnit,'(6A)')     ' k |',
69       &             ' --  Average --  |', ' --  Std.Dev  -- |',  c    &             ' --  Average  -- |', ' --  Std.Dev  -- |',
70       &             ' --     min  --  |', ' --     max   -- |',' Vol'  c    &             ' --    min    -- |', ' --    max    -- |',' --  Vol'
71             WRITE(nUnit,'(6A)') ' k |',
72         &     ' --   Average    -- |', ' --   Std.Dev    -- |',
73         &     ' --     min      -- |', ' --     max      -- |',' --    Vol'
74           DO k=0,klev           DO k=0,klev
75  C     full precision, do not write the volume:  C     full precision, do not write the volume:
76  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)
77  C     reduced precision + write the volume (usefull for testing):  C     reduced precision + write the volume (usefull for testing):
78            WRITE(nUnit,'(I3,1P5E18.10)')  c         WRITE(nUnit,'(I3,1P5E18.10)')
79    c    &       k, (statGlob(i,k,j),i=1,nStats), statGlob(0,k,j)
80    C     full precision + write the volume:
81              WRITE(nUnit,'(I3,1P5E21.13)')
82       &       k, (statGlob(i,k,j),i=1,nStats), statGlob(0,k,j)       &       k, (statGlob(i,k,j),i=1,nStats), statGlob(0,k,j)
83           ENDDO           ENDDO
84          ENDIF          ENDIF

Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6

  ViewVC Help
Powered by ViewVC 1.1.22