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

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

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

revision 1.3 by jmc, Tue Feb 21 01:23:57 2006 UTC revision 1.4 by jmc, Tue Feb 5 15:31:19 2008 UTC
# Line 59  C-      get a free unit number as the I/ Line 59  C-      get a free unit number as the I/
59       &          diagSt_Fname(n)(1:iL), '.', nIter0, '.txt'       &          diagSt_Fname(n)(1:iL), '.', nIter0, '.txt'
60            OPEN( nUnit, FILE=dataFName, STATUS='unknown' )            OPEN( nUnit, FILE=dataFName, STATUS='unknown' )
61    
62            WRITE(msgBuf,'(4A,I6)') 'DIAGSTATS_INI_IO: ',            WRITE(msgBuf,'(4A,I6)') 'DIAGSTATS_INI_IO: ',
63       &         'open file: ',dataFName(1:iL+15), ' , unit=', nUnit       &         'open file: ',dataFName(1:iL+15), ' , unit=', nUnit
64            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
65       &                        SQUEEZE_RIGHT , myThid)       &                        SQUEEZE_RIGHT , myThid)
# Line 74  C-      write a kind of header: Line 74  C-      write a kind of header:
74            DO j=0,nRegions            DO j=0,nRegions
75              IF (diagSt_region(j,n).GE.1 .AND.iL+3.LE.MAX_LEN_MBUF) THEN              IF (diagSt_region(j,n).GE.1 .AND.iL+3.LE.MAX_LEN_MBUF) THEN
76                tmpBuf(1:iL) = msgBuf(1:iL)                tmpBuf(1:iL) = msgBuf(1:iL)
77                WRITE(msgBuf,'(A,I3)') tmpBuf(1:iL),j                            WRITE(msgBuf,'(A,I3)') tmpBuf(1:iL),j
78                iL = iL+3                iL = iL+3
79              ENDIF              ENDIF
80            ENDDO            ENDDO
81            WRITE(nUnit,'(A)') msgBuf(1:iL)            WRITE(nUnit,'(A)') msgBuf(1:iL)
82            DO j=1,diagSt_nbFlds(n),10            DO j=1,diagSt_nbFlds(n),10
83              WRITE(nUnit,'(A,20A)')      '# Fields       :',              WRITE(nUnit,'(A,20A)')      '# Fields       :',
84       &        (' ', diagSt_Flds(m,n), m=j,MIN(diagSt_nbFlds(n),j+9) )       &        (' ', diagSt_Flds(m,n), m=j,MIN(diagSt_nbFlds(n),j+9) )
85            ENDDO            ENDDO
86            DO j=1,diagSt_nbFlds(n),50            DO j=1,diagSt_nbFlds(n),50
87              WRITE(nUnit,'(A,50I4)')   '# Nb of levels : ',              WRITE(nUnit,'(A,50I4)')   '# Nb of levels : ',
88       &         ( kdiag(jSdiag(m,n)), m=j,MIN(diagSt_nbFlds(n),j+49) )       &         ( kdiag(jSdiag(m,n)), m=j,MIN(diagSt_nbFlds(n),j+49) )
89            ENDDO            ENDDO
90            WRITE(nUnit,'(2A)') '# end of header ----------------------',            WRITE(nUnit,'(2A)') '# end of header ----------------------',
# Line 92  C-      write a kind of header: Line 92  C-      write a kind of header:
92            WRITE(nUnit,'(A)') ' '            WRITE(nUnit,'(A)') ' '
93    
94          ENDDO          ENDDO
95          
96        ENDIF        ENDIF
97    
98        _END_MASTER( myThid )        _END_MASTER( myThid )
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22