--- MITgcm/pkg/diagnostics/diagstats_ascii_out.F 2006/02/21 01:23:57 1.4 +++ MITgcm/pkg/diagnostics/diagstats_ascii_out.F 2008/02/05 15:31:19 1.5 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/diagnostics/diagstats_ascii_out.F,v 1.4 2006/02/21 01:23:57 jmc Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/diagnostics/diagstats_ascii_out.F,v 1.5 2008/02/05 15:31:19 jmc Exp $ C $Name: $ #include "DIAG_OPTIONS.h" @@ -40,7 +40,7 @@ C !LOCAL VARIABLES: INTEGER im, ix, iv PARAMETER ( iv = nStats - 2 , im = nStats - 1 , ix = nStats ) - INTEGER i, j, k, klev, nUnit + INTEGER i, j, k, klev, nUnit C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| @@ -64,15 +64,15 @@ & ' ; Iter =',myIter,' ; region #',j, ' ; nb.Lev =',kdiag(ndId) c WRITE(nUnit,'(5A)') ' k |', c & ' -- Average -- |', ' -- Std.Dev -- |', -c & ' -- min -- |', ' -- max -- |' +c & ' -- min -- |', ' -- max -- |' WRITE(nUnit,'(6A)') ' k |', & ' -- Average -- |', ' -- Std.Dev -- |', - & ' -- min -- |', ' -- max -- |',' Vol' + & ' -- min -- |', ' -- max -- |',' Vol' DO k=0,klev C full precision, do not write the volume: c WRITE(nUnit,'(I3,1P4E20.12)') k,(statGlob(i,k,j),i=1,nStats) C reduced precision + write the volume (usefull for testing): - WRITE(nUnit,'(I3,1P5E18.10)') + WRITE(nUnit,'(I3,1P5E18.10)') & k, (statGlob(i,k,j),i=1,nStats), statGlob(0,k,j) ENDDO ENDIF