/[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.18 by jmc, Mon Jun 5 18:15:53 2006 UTC revision 1.24 by jmc, Mon Jan 11 19:44:07 2010 UTC
# Line 35  C     timePhase :: phase (in s) within t Line 35  C     timePhase :: phase (in s) within t
35  C     averagingFreq  :: frequency (in s) for periodic averaging interval  C     averagingFreq  :: frequency (in s) for periodic averaging interval
36  C     averagingPhase :: phase     (in s) for periodic averaging interval  C     averagingPhase :: phase     (in s) for periodic averaging interval
37  C     repeatCycle    :: number of averaging intervals in 1 cycle  C     repeatCycle    :: number of averaging intervals in 1 cycle
38    C     mising_value     :: missing value for floats   in output
39    C     mising_value_int :: missing value for integers in output
40  C     levels    :: List Output Levels  C     levels    :: List Output Levels
41  C     fields    :: List Output Fields  C     fields    :: List Output Fields
42  C     filename  :: List Output Filename  C     filename  :: List Output Filename
# Line 55  C     stat_fname  :: List of statistics Line 57  C     stat_fname  :: List of statistics
57        _RL         frequency(ldimLoc), timePhase(ldimLoc)        _RL         frequency(ldimLoc), timePhase(ldimLoc)
58        _RL         averagingFreq(ldimLoc), averagingPhase(ldimLoc)        _RL         averagingFreq(ldimLoc), averagingPhase(ldimLoc)
59        INTEGER     repeatCycle(ldimLoc)        INTEGER     repeatCycle(ldimLoc)
60          _RL         missing_value(ldimLoc)
61          INTEGER     missing_value_int(ldimLoc)
62        _RL         levels(kdimLoc,ldimLoc)        _RL         levels(kdimLoc,ldimLoc)
63        _RL         stat_freq(ldimLoc), stat_phase(ldimLoc)        _RL         stat_freq(ldimLoc), stat_phase(ldimLoc)
64        CHARACTER*8 fields(fdimLoc,ldimLoc)        CHARACTER*8 fields(fdimLoc,ldimLoc)
# Line 64  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 77  C--   full level output: Line 82  C--   full level output:
82        NAMELIST / DIAGNOSTICS_LIST /        NAMELIST / DIAGNOSTICS_LIST /
83       &     frequency, timePhase,       &     frequency, timePhase,
84       &     averagingFreq, averagingPhase, repeatCycle,       &     averagingFreq, averagingPhase, repeatCycle,
85         &     missing_value, missing_value_int,
86       &     levels, fields, filename, fileflags,       &     levels, fields, filename, fileflags,
87       &     dumpAtLast, diag_mnc,       &     dumpAtLast, diag_mnc,
88       &     diag_pickup_read,     diag_pickup_write,       &     diag_pickup_read,     diag_pickup_write,
# Line 104  C     Initialize and Read Diagnostics Na Line 110  C     Initialize and Read Diagnostics Na
110          averagingPhase(l)= 0.          averagingPhase(l)= 0.
111          repeatCycle(l)   = 0          repeatCycle(l)   = 0
112          filename(l)   = blkFilName          filename(l)   = blkFilName
113  C       eight spaces:        12345678          missing_value(l)     = UNSET_RL
114  c       fileflags(l)(1:8) = '        '          missing_value_int(l) = UNSET_I
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 183  C     Initialise DIAG_SELECT common bloc Line 189  C     Initialise DIAG_SELECT common bloc
189          nlevels(n) = 0          nlevels(n) = 0
190          nfields(n) = 0          nfields(n) = 0
191          fnames(n) = blkFilName          fnames(n) = blkFilName
192            misvalFlt(n) = UNSET_RL
193            misvalInt(n) = UNSET_I
194          DO k = 1,numLevels          DO k = 1,numLevels
195            levs(k,n) = 0            levs(k,n) = 0
196          ENDDO          ENDDO
# Line 227  C-     Only lists with non-empty file na Line 235  C-     Only lists with non-empty file na
235             averagePhase(n) = averagingPhase(l)             averagePhase(n) = averagingPhase(l)
236             averageCycle(n) = repeatCycle(l)             averageCycle(n) = repeatCycle(l)
237           ELSEIF (averagingFreq(l).NE.0. .OR. repeatCycle(l).NE.0) THEN           ELSEIF (averagingFreq(l).NE.0. .OR. repeatCycle(l).NE.0) THEN
238             WRITE(msgBuf,'(2A,F17.6,I3)') 'DIAGNOSTICS_READPARMS: ',             WRITE(msgBuf,'(2A,F18.6,I4)') 'DIAGNOSTICS_READPARMS: ',
239       &       'unvalid Average-Freq & Cycle:',       &       'unvalid Average-Freq & Cycle:',
240       &       averagingFreq(l), repeatCycle(l)       &       averagingFreq(l), repeatCycle(l)
241             CALL PRINT_ERROR( msgBuf , myThid )             CALL PRINT_ERROR( msgBuf , myThid )
242             WRITE(msgBuf,'(2A,I3,2A)') 'DIAGNOSTICS_READPARMS: ',             WRITE(msgBuf,'(2A,I4,2A)') 'DIAGNOSTICS_READPARMS: ',
243       &         ' for list l=', l, ', filename: ', filename(l)       &         ' for list l=', l, ', filename: ', filename(l)
244             CALL PRINT_ERROR( msgBuf , myThid )             CALL PRINT_ERROR( msgBuf , myThid )
245             STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'             STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'
# Line 242  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 )
254         &        misvalFlt(n) = missing_value(l)
255             IF ( missing_value_int(l) .NE. UNSET_I )
256         &        misvalInt(n) = missing_value_int(l)
257           fnames(n)  = filename (l)           fnames(n)  = filename (l)
258           fflags(n)  = fileflags(l)           fflags(n)  = fileflags(l)
259           nlevels(n) = 0           nlevels(n) = 0
# Line 252  C-     Only lists with non-empty file na Line 264  C-     Only lists with non-empty file na
264                 nlevels(n) = nlevels(n) + 1                 nlevels(n) = nlevels(n) + 1
265                 levs(nlevels(n),n) = levels(k,l)                 levs(nlevels(n),n) = levels(k,l)
266               ELSEIF ( levels(k,l).NE.UNSET_RL ) THEN               ELSEIF ( levels(k,l).NE.UNSET_RL ) THEN
267                WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',                WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
268       &         'Exceed Max.Num. of Levels numLevels=', numLevels       &         'Exceed Max.Num. of Levels numLevels=', numLevels
269                CALL PRINT_ERROR( msgBuf , myThid )                CALL PRINT_ERROR( msgBuf , myThid )
270                WRITE(msgBuf,'(2A,I3,A,F3.0)') 'DIAGNOSTICS_READPARMS: ',                WRITE(msgBuf,'(2A,I4,A,F8.0)') 'DIAGNOSTICS_READPARMS: ',
271       &         'when trying to add level(k=', k, ' )=', levels(k,l)       &         'when trying to add level(k=', k, ' )=', levels(k,l)
272                CALL PRINT_ERROR( msgBuf , myThid )                CALL PRINT_ERROR( msgBuf , myThid )
273                WRITE(msgBuf,'(2A,I3,2A)') 'DIAGNOSTICS_READPARMS: ',                WRITE(msgBuf,'(2A,I4,2A)') 'DIAGNOSTICS_READPARMS: ',
274       &         ' for list l=', l, ', filename: ', filename(l)       &         ' for list l=', l, ', filename: ', filename(l)
275                CALL PRINT_ERROR( msgBuf , myThid )                CALL PRINT_ERROR( msgBuf , myThid )
276                STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'                STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'
# Line 275  C-       will set levels later, once the Line 287  C-       will set levels later, once the
287               nfields(n) = nfields(n) + 1               nfields(n) = nfields(n) + 1
288               flds(nfields(n),n) = fields(m,l)               flds(nfields(n),n) = fields(m,l)
289             ELSEIF ( fields(m,l).NE.blk8c ) THEN             ELSEIF ( fields(m,l).NE.blk8c ) THEN
290               WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
291       &        'Exceed Max.Num. of Fields/list numperlist=', numperlist       &        'Exceed Max.Num. of Fields/list numperlist=', numperlist
292               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
293               WRITE(msgBuf,'(2A,I3,3A,I3,2A)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4,3A,I4,2A)') 'DIAGNOSTICS_READPARMS: ',
294       &        'when trying to add field (m=', m, ' ): ',fields(m,l)       &        'when trying to add field (m=', m, ' ): ',fields(m,l)
295               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
296               WRITE(msgBuf,'(2A,I3,2A)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4,2A)') 'DIAGNOSTICS_READPARMS: ',
297       &        ' in list l=', l, ', filename: ', filename(l)       &        ' in list l=', l, ', filename: ', filename(l)
298               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
299               STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'               STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'
# Line 290  C-       will set levels later, once the Line 302  C-       will set levels later, once the
302           nlists = nlists + 1           nlists = nlists + 1
303  c        write(6,*) 'list summary:',n,nfields(n),nlevels(n)  c        write(6,*) 'list summary:',n,nfields(n),nlevels(n)
304         ELSEIF ( iLen.GE.1 ) THEN         ELSEIF ( iLen.GE.1 ) THEN
305           WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
306       &            'Exceed Max.Num. of list numlists=', numlists       &            'Exceed Max.Num. of list numlists=', numlists
307           CALL PRINT_ERROR( msgBuf , myThid )           CALL PRINT_ERROR( msgBuf , myThid )
308           WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
309       &    'when trying to add list l=', l       &    'when trying to add list l=', l
310           CALL PRINT_ERROR( msgBuf , myThid )           CALL PRINT_ERROR( msgBuf , myThid )
311           WRITE(msgBuf,'(2A,F17.6,2A)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,F18.6,2A)') 'DIAGNOSTICS_READPARMS: ',
312       &    ' Frq=', frequency(l), ', filename: ', filename(l)       &    ' Frq=', frequency(l), ', filename: ', filename(l)
313           CALL PRINT_ERROR( msgBuf , myThid )           CALL PRINT_ERROR( msgBuf , myThid )
314           STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'           STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'
# Line 347  C     for now, simpler just to read it f Line 359  C     for now, simpler just to read it f
359         ENDIF         ENDIF
360        ENDDO        ENDDO
361        IF ( j.GT.nRegions ) THEN        IF ( j.GT.nRegions ) THEN
362           WRITE(msgBuf,'(2A,I3,A)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,I4,A)') 'DIAGNOSTICS_READPARMS: ',
363       &   'set_regMask & val_regMask lists assume at least',j,' regions'       &   'set_regMask & val_regMask lists assume at least',j,' regions'
364           CALL PRINT_ERROR( msgBuf , myThid )           CALL PRINT_ERROR( msgBuf , myThid )
365           WRITE(msgBuf,'(2A)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A)') 'DIAGNOSTICS_READPARMS: ',
# Line 376  C-     Only lists with non-empty file na Line 388  C-     Only lists with non-empty file na
388               diagSt_region(j,n) = 1               diagSt_region(j,n) = 1
389               regionCount = regionCount + 1               regionCount = regionCount + 1
390              ELSE              ELSE
391               WRITE(msgBuf,'(2A,I3,2A)')               WRITE(msgBuf,'(2A,I4,2A)')
392       &        'DIAGNOSTICS_READPARMS:',       &        'DIAGNOSTICS_READPARMS:',
393       &        ' in list l=', l, ', stat_fname: ', stat_fname(l)       &        ' in list l=', l, ', stat_fname: ', stat_fname(l)
394               CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,               CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
395       &                           SQUEEZE_RIGHT , myThid )       &                           SQUEEZE_RIGHT , myThid )
396               WRITE(msgBuf,'(A,I3,A)')               WRITE(msgBuf,'(A,I4,A)')
397       &        'DIAGNOSTICS_READPARMS: region=',j,       &        'DIAGNOSTICS_READPARMS: region=',j,
398       &        ' can only be selected once => ignore 2nd selection'       &        ' can only be selected once => ignore 2nd selection'
399               CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,               CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
400       &                           SQUEEZE_RIGHT , myThid )       &                           SQUEEZE_RIGHT , myThid )
401              ENDIF              ENDIF
402             ELSEIF ( j.NE.UNSET_I ) THEN             ELSEIF ( j.NE.UNSET_I ) THEN
403               WRITE(msgBuf,'(A,I3,A,I3,2A)')               WRITE(msgBuf,'(A,I4,A,I4,2A)')
404       &       'DIAGNOSTICS_READPARMS: region=',j,       &       'DIAGNOSTICS_READPARMS: region=',j,
405       &         ' in list l=', l, ', stat_fname: ', stat_fname(l)       &         ' in list l=', l, ', stat_fname: ', stat_fname(l)
406               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
407               WRITE(msgBuf,'(2A,I3,A,I3,2A)')               WRITE(msgBuf,'(2A,I4,A,I4,2A)')
408       &       'DIAGNOSTICS_READPARMS: ==> exceed Max.Nb of regions',       &       'DIAGNOSTICS_READPARMS: ==> exceed Max.Nb of regions',
409       &       '(=',nRegions,' )'       &       '(=',nRegions,' )'
410               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
# Line 410  C-       no region selected => default i Line 422  C-       no region selected => default i
422               diagSt_nbFlds(n) = diagSt_nbFlds(n) + 1               diagSt_nbFlds(n) = diagSt_nbFlds(n) + 1
423               diagSt_Flds(diagSt_nbFlds(n),n) = stat_fields(m,l)               diagSt_Flds(diagSt_nbFlds(n),n) = stat_fields(m,l)
424             ELSEIF ( stat_fields(m,l).NE.blk8c ) THEN             ELSEIF ( stat_fields(m,l).NE.blk8c ) THEN
425               WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
426       &        'Exceed Max.Num. of Fields/list numperlist=', numperlist       &        'Exceed Max.Num. of Fields/list numperlist=', numperlist
427               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
428               WRITE(msgBuf,'(2A,I3,3A,I3,2A)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4,3A,I4,2A)') 'DIAGNOSTICS_READPARMS: ',
429       &        'when trying to add stat_field (m=', m,       &        'when trying to add stat_field (m=', m,
430       &        ' ): ',stat_fields(m,l)       &        ' ): ',stat_fields(m,l)
431               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
432               WRITE(msgBuf,'(2A,I3,2A)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4,2A)') 'DIAGNOSTICS_READPARMS: ',
433       &        ' in list l=', l, ', stat_fname: ', stat_fname(l)       &        ' in list l=', l, ', stat_fname: ', stat_fname(l)
434               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
435               STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'               STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'
# Line 426  C-       no region selected => default i Line 438  C-       no region selected => default i
438           diagSt_nbLists = diagSt_nbLists + 1           diagSt_nbLists = diagSt_nbLists + 1
439  c        write(6,*) 'stat-list summary:',n,diagSt_nbFlds(n),regionCount  c        write(6,*) 'stat-list summary:',n,diagSt_nbFlds(n),regionCount
440         ELSEIF ( iLen.GE.1 ) THEN         ELSEIF ( iLen.GE.1 ) THEN
441           WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
442       &            'Exceed Max.Num. of list numlists=', numlists       &            'Exceed Max.Num. of list numlists=', numlists
443           CALL PRINT_ERROR( msgBuf , myThid )           CALL PRINT_ERROR( msgBuf , myThid )
444           WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
445       &    'when trying to add stat_list l=', l       &    'when trying to add stat_list l=', l
446           CALL PRINT_ERROR( msgBuf , myThid )           CALL PRINT_ERROR( msgBuf , myThid )
447           WRITE(msgBuf,'(2A,F17.6,2A)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,F18.6,2A)') 'DIAGNOSTICS_READPARMS: ',
448       &    ' Frq=', stat_freq(l), ', stat_fname: ', stat_fname(l)       &    ' Frq=', stat_freq(l), ', stat_fname: ', stat_fname(l)
449           CALL PRINT_ERROR( msgBuf , myThid )           CALL PRINT_ERROR( msgBuf , myThid )
450           STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'           STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'
# Line 452  C     Echo History List Data Structure Line 464  C     Echo History List Data Structure
464       & '-----------------------------------------------------'       & '-----------------------------------------------------'
465        CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)        CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
466        DO n = 1,nlists        DO n = 1,nlists
467          WRITE(msgBuf,'(2a)') 'Creating Output Stream: ', fnames(n)          WRITE(msgBuf,'(2A)') 'Creating Output Stream: ', fnames(n)
468          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
469          WRITE(msgBuf,'(2(A,F17.6))') 'Output Frequency:', freq(n),          WRITE(msgBuf,'(2(A,F18.6))') 'Output Frequency:', freq(n),
470       &                               ' ; Phase: ', phase(n)       &                               ' ; Phase: ', phase(n)
471          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
472          WRITE(msgBuf,'(2(A,F17.6),A,I3)')          WRITE(msgBuf,'(2(A,F18.6),A,I4)')
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          IF ( nlevels(n).EQ.-1 ) THEN          IF ( fflags(n).EQ.blk8c ) THEN
477              WRITE(msgBuf,'(A,1PE20.12,A,I12,3A)')
478         &       ' missing value:',  misvalFlt(n),
479         &       ' ; 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)
487            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
494             DO l=1,nlevels(n),10
495              m = MIN(nlevels(n),l+9)
496              WRITE(msgBuf,'(A,1P10E10.3)')' interp:  ', (levs(k,n),k=l,m)
497              CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
498             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
508          DO nf = 1,nfields(n),10          DO nf = 1,nfields(n),10
509          m = MIN(nfields(n),nf+9)            m = MIN(nfields(n),nf+9)
510          WRITE(msgBuf,*) 'Fields:   ',(' ',flds(l,n),l=nf,m)            WRITE(msgBuf,'(21A)') ' Fields:   ',(' ',flds(l,n),l=nf,m)
511          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)            CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
512          ENDDO          ENDDO
513        ENDDO        ENDDO
514        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
# Line 484  C     Echo History List Data Structure Line 518  C     Echo History List Data Structure
518       &     ' DIAGNOSTICS_READPARMS: statistics diags. summary:'       &     ' DIAGNOSTICS_READPARMS: statistics diags. summary:'
519        CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)        CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
520        DO n = 1,diagSt_nbLists        DO n = 1,diagSt_nbLists
521          WRITE(msgBuf,'(2a)') 'Creating Stats. Output Stream: ',          WRITE(msgBuf,'(2A)') 'Creating Stats. Output Stream: ',
522       &                       diagSt_Fname(n)       &                       diagSt_Fname(n)
523          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
524          WRITE(msgBuf,'(2(A,F17.6))') 'Output Frequency:',          WRITE(msgBuf,'(2(A,F18.6))') 'Output Frequency:',
525       &               diagSt_freq(n), ' ; Phase: ', diagSt_phase(n)       &               diagSt_freq(n), ' ; Phase: ', diagSt_phase(n)
526          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
527          WRITE(msgBuf,'(A)') ' Regions : '          WRITE(msgBuf,'(A)') ' Regions: '
528          l = 12          l = 10
529          DO j=0,nRegions          DO j=0,nRegions
530           IF ( diagSt_region(j,n).GE.1 ) THEN           IF ( diagSt_region(j,n).GE.1 ) THEN
           IF (l+3.LE.MAX_LEN_MBUF) WRITE(msgBuf,'(A,I3)') msgBuf(1:l),j  
531            l = l+3            l = l+3
532              IF (l.LE.MAX_LEN_MBUF) WRITE(msgBuf(l-2:l),'(I3)') j
533           ENDIF           ENDIF
534          ENDDO          ENDDO
535          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
536          WRITE(msgBuf,*) 'Fields:   ',          DO nf = 1,diagSt_nbFlds(n),10
537       &                 (' ',diagSt_Flds(l,n),l=1,diagSt_nbFlds(n))            m = MIN(diagSt_nbFlds(n),nf+9)
538          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)            WRITE(msgBuf,'(21A)') ' Fields:   ',
539         &                 (' ',diagSt_Flds(l,n),l=nf,m)
540              CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
541            ENDDO
542        ENDDO        ENDDO
543        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
544       & '-----------------------------------------------------'       & '-----------------------------------------------------'
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.24
  ViewVC Help
Powered by ViewVC 1.1.22