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

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

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

revision 1.23 by jmc, Fri May 30 01:58:51 2008 UTC revision 1.24 by jmc, Mon Jan 11 19:44:07 2010 UTC
# Line 68  C     stat_fname  :: List of statistics Line 68  C     stat_fname  :: List of statistics
68        CHARACTER*8 fileflags(ldimLoc)        CHARACTER*8 fileflags(ldimLoc)
69        CHARACTER*8 blk8c        CHARACTER*8 blk8c
70        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
71          CHARACTER*12 suffix
72        INTEGER stat_region(rdimLoc,ldimLoc)        INTEGER stat_region(rdimLoc,ldimLoc)
73        INTEGER set_regMask(rdimLoc)        INTEGER set_regMask(rdimLoc)
74        _RS     val_regMask(rdimLoc)        _RS     val_regMask(rdimLoc)
# Line 111  C     Initialize and Read Diagnostics Na Line 112  C     Initialize and Read Diagnostics Na
112          filename(l)   = blkFilName          filename(l)   = blkFilName
113          missing_value(l)     = UNSET_RL          missing_value(l)     = UNSET_RL
114          missing_value_int(l) = UNSET_I          missing_value_int(l) = UNSET_I
 C       eight spaces:        12345678  
 c       fileflags(l)(1:8) = '        '  
115          fileflags(l)  = blk8c          fileflags(l)  = blk8c
116          DO k = 1,kdimLoc          DO k = 1,kdimLoc
117            levels(k,l) = UNSET_RL            levels(k,l) = UNSET_RL
# Line 251  C-     Only lists with non-empty file na Line 250  C-     Only lists with non-empty file na
250             averageFreq(n)  = frequency(l)             averageFreq(n)  = frequency(l)
251             averagePhase(n) = phase(n)             averagePhase(n) = phase(n)
252           ENDIF           ENDIF
253           IF ( missing_value(l) .NE. UNSET_RL )           IF ( missing_value(l) .NE. UNSET_RL )
254       &        misvalFlt(n) = missing_value(l)       &        misvalFlt(n) = missing_value(l)
255           IF ( missing_value_int(l) .NE. UNSET_I )           IF ( missing_value_int(l) .NE. UNSET_I )
256       &        misvalInt(n) = missing_value_int(l)       &        misvalInt(n) = missing_value_int(l)
257           fnames(n)  = filename (l)           fnames(n)  = filename (l)
258           fflags(n)  = fileflags(l)           fflags(n)  = fileflags(l)
# Line 474  C     Echo History List Data Structure Line 473  C     Echo History List Data Structure
473       &    ' Averaging Freq.:', averageFreq(n),       &    ' Averaging Freq.:', averageFreq(n),
474       &    ' , Phase: ', averagePhase(n), ' , Cycle:', averageCycle(n)       &    ' , Phase: ', averagePhase(n), ' , Cycle:', averageCycle(n)
475          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
476          WRITE(msgBuf,'(A,1PE20.12,A,I12)')          IF ( fflags(n).EQ.blk8c ) THEN
477              WRITE(msgBuf,'(A,1PE20.12,A,I12,3A)')
478       &       ' missing value:',  misvalFlt(n),       &       ' missing value:',  misvalFlt(n),
479       &       ' ; for integers:', misvalInt(n)       &       ' ; for integers:', misvalInt(n)
480            ELSE
481              WRITE(msgBuf,'(A,1PE20.12,A,I12,3A)')
482         &       ' missing value:',  misvalFlt(n),
483         &       ' ; for integers:', misvalInt(n),
484         &       ' ; F-Flags="', fflags(n),'"'
485            ENDIF
486          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
487          IF ( nlevels(n).EQ.-1 ) THEN          IF ( nlevels(n).EQ.-1 .AND. fflags(n)(2:2).EQ.'I' ) THEN
488              WRITE(msgBuf,'(A)') ' Cumulate all Levels (to be set later)'
489              CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
490            ELSEIF ( nlevels(n).EQ.-1 ) THEN
491            WRITE(msgBuf,'(A,A)') ' Levels:    ','will be set later'            WRITE(msgBuf,'(A,A)') ' Levels:    ','will be set later'
492            CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)            CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
493          ELSEIF ( fflags(n)(2:2).EQ.'P' ) THEN          ELSEIF ( fflags(n)(2:2).EQ.'P' ) THEN
# Line 488  C     Echo History List Data Structure Line 497  C     Echo History List Data Structure
497            CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)            CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
498           ENDDO           ENDDO
499          ELSE          ELSE
500             suffix = ' Levels:    '
501             IF ( fflags(n)(2:2).EQ.'I' ) suffix = ' Sum Levels:'
502           DO l=1,nlevels(n),20           DO l=1,nlevels(n),20
503            m = MIN(nlevels(n),l+19)            m = MIN(nlevels(n),l+19)
504            WRITE(msgBuf,'(A,20F5.0)') ' Levels:    ', (levs(k,n),k=l,m)            WRITE(msgBuf,'(A,20F5.0)') suffix, (levs(k,n),k=l,m)
505            CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)            CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
506           ENDDO           ENDDO
507          ENDIF          ENDIF
Legend:
Removed from v.1.23  
changed lines
  Added in v.1.24
  ViewVC Help
Powered by ViewVC 1.1.22