/[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.20 by jmc, Mon Aug 27 23:04:44 2007 UTC revision 1.21 by jmc, Tue Feb 5 15:31:19 2008 UTC
# Line 227  C-     Only lists with non-empty file na Line 227  C-     Only lists with non-empty file na
227             averagePhase(n) = averagingPhase(l)             averagePhase(n) = averagingPhase(l)
228             averageCycle(n) = repeatCycle(l)             averageCycle(n) = repeatCycle(l)
229           ELSEIF (averagingFreq(l).NE.0. .OR. repeatCycle(l).NE.0) THEN           ELSEIF (averagingFreq(l).NE.0. .OR. repeatCycle(l).NE.0) THEN
230             WRITE(msgBuf,'(2A,F18.6,I3)') 'DIAGNOSTICS_READPARMS: ',             WRITE(msgBuf,'(2A,F18.6,I4)') 'DIAGNOSTICS_READPARMS: ',
231       &       'unvalid Average-Freq & Cycle:',       &       'unvalid Average-Freq & Cycle:',
232       &       averagingFreq(l), repeatCycle(l)       &       averagingFreq(l), repeatCycle(l)
233             CALL PRINT_ERROR( msgBuf , myThid )             CALL PRINT_ERROR( msgBuf , myThid )
234             WRITE(msgBuf,'(2A,I3,2A)') 'DIAGNOSTICS_READPARMS: ',             WRITE(msgBuf,'(2A,I4,2A)') 'DIAGNOSTICS_READPARMS: ',
235       &         ' for list l=', l, ', filename: ', filename(l)       &         ' for list l=', l, ', filename: ', filename(l)
236             CALL PRINT_ERROR( msgBuf , myThid )             CALL PRINT_ERROR( msgBuf , myThid )
237             STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'             STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'
# Line 252  C-     Only lists with non-empty file na Line 252  C-     Only lists with non-empty file na
252                 nlevels(n) = nlevels(n) + 1                 nlevels(n) = nlevels(n) + 1
253                 levs(nlevels(n),n) = levels(k,l)                 levs(nlevels(n),n) = levels(k,l)
254               ELSEIF ( levels(k,l).NE.UNSET_RL ) THEN               ELSEIF ( levels(k,l).NE.UNSET_RL ) THEN
255                WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',                WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
256       &         'Exceed Max.Num. of Levels numLevels=', numLevels       &         'Exceed Max.Num. of Levels numLevels=', numLevels
257                CALL PRINT_ERROR( msgBuf , myThid )                CALL PRINT_ERROR( msgBuf , myThid )
258                WRITE(msgBuf,'(2A,I3,A,F8.0)') 'DIAGNOSTICS_READPARMS: ',                WRITE(msgBuf,'(2A,I4,A,F8.0)') 'DIAGNOSTICS_READPARMS: ',
259       &         'when trying to add level(k=', k, ' )=', levels(k,l)       &         'when trying to add level(k=', k, ' )=', levels(k,l)
260                CALL PRINT_ERROR( msgBuf , myThid )                CALL PRINT_ERROR( msgBuf , myThid )
261                WRITE(msgBuf,'(2A,I3,2A)') 'DIAGNOSTICS_READPARMS: ',                WRITE(msgBuf,'(2A,I4,2A)') 'DIAGNOSTICS_READPARMS: ',
262       &         ' for list l=', l, ', filename: ', filename(l)       &         ' for list l=', l, ', filename: ', filename(l)
263                CALL PRINT_ERROR( msgBuf , myThid )                CALL PRINT_ERROR( msgBuf , myThid )
264                STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'                STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'
# Line 275  C-       will set levels later, once the Line 275  C-       will set levels later, once the
275               nfields(n) = nfields(n) + 1               nfields(n) = nfields(n) + 1
276               flds(nfields(n),n) = fields(m,l)               flds(nfields(n),n) = fields(m,l)
277             ELSEIF ( fields(m,l).NE.blk8c ) THEN             ELSEIF ( fields(m,l).NE.blk8c ) THEN
278               WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
279       &        'Exceed Max.Num. of Fields/list numperlist=', numperlist       &        'Exceed Max.Num. of Fields/list numperlist=', numperlist
280               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
281               WRITE(msgBuf,'(2A,I3,3A,I3,2A)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4,3A,I4,2A)') 'DIAGNOSTICS_READPARMS: ',
282       &        'when trying to add field (m=', m, ' ): ',fields(m,l)       &        'when trying to add field (m=', m, ' ): ',fields(m,l)
283               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
284               WRITE(msgBuf,'(2A,I3,2A)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4,2A)') 'DIAGNOSTICS_READPARMS: ',
285       &        ' in list l=', l, ', filename: ', filename(l)       &        ' in list l=', l, ', filename: ', filename(l)
286               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
287               STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'               STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'
# Line 290  C-       will set levels later, once the Line 290  C-       will set levels later, once the
290           nlists = nlists + 1           nlists = nlists + 1
291  c        write(6,*) 'list summary:',n,nfields(n),nlevels(n)  c        write(6,*) 'list summary:',n,nfields(n),nlevels(n)
292         ELSEIF ( iLen.GE.1 ) THEN         ELSEIF ( iLen.GE.1 ) THEN
293           WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
294       &            'Exceed Max.Num. of list numlists=', numlists       &            'Exceed Max.Num. of list numlists=', numlists
295           CALL PRINT_ERROR( msgBuf , myThid )           CALL PRINT_ERROR( msgBuf , myThid )
296           WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
297       &    'when trying to add list l=', l       &    'when trying to add list l=', l
298           CALL PRINT_ERROR( msgBuf , myThid )           CALL PRINT_ERROR( msgBuf , myThid )
299           WRITE(msgBuf,'(2A,F18.6,2A)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,F18.6,2A)') 'DIAGNOSTICS_READPARMS: ',
# Line 347  C     for now, simpler just to read it f Line 347  C     for now, simpler just to read it f
347         ENDIF         ENDIF
348        ENDDO        ENDDO
349        IF ( j.GT.nRegions ) THEN        IF ( j.GT.nRegions ) THEN
350           WRITE(msgBuf,'(2A,I3,A)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,I4,A)') 'DIAGNOSTICS_READPARMS: ',
351       &   'set_regMask & val_regMask lists assume at least',j,' regions'       &   'set_regMask & val_regMask lists assume at least',j,' regions'
352           CALL PRINT_ERROR( msgBuf , myThid )           CALL PRINT_ERROR( msgBuf , myThid )
353           WRITE(msgBuf,'(2A)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A)') 'DIAGNOSTICS_READPARMS: ',
# Line 376  C-     Only lists with non-empty file na Line 376  C-     Only lists with non-empty file na
376               diagSt_region(j,n) = 1               diagSt_region(j,n) = 1
377               regionCount = regionCount + 1               regionCount = regionCount + 1
378              ELSE              ELSE
379               WRITE(msgBuf,'(2A,I3,2A)')               WRITE(msgBuf,'(2A,I4,2A)')
380       &        'DIAGNOSTICS_READPARMS:',       &        'DIAGNOSTICS_READPARMS:',
381       &        ' in list l=', l, ', stat_fname: ', stat_fname(l)       &        ' in list l=', l, ', stat_fname: ', stat_fname(l)
382               CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,               CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
383       &                           SQUEEZE_RIGHT , myThid )       &                           SQUEEZE_RIGHT , myThid )
384               WRITE(msgBuf,'(A,I3,A)')               WRITE(msgBuf,'(A,I4,A)')
385       &        'DIAGNOSTICS_READPARMS: region=',j,       &        'DIAGNOSTICS_READPARMS: region=',j,
386       &        ' can only be selected once => ignore 2nd selection'       &        ' can only be selected once => ignore 2nd selection'
387               CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,               CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
388       &                           SQUEEZE_RIGHT , myThid )       &                           SQUEEZE_RIGHT , myThid )
389              ENDIF              ENDIF
390             ELSEIF ( j.NE.UNSET_I ) THEN             ELSEIF ( j.NE.UNSET_I ) THEN
391               WRITE(msgBuf,'(A,I3,A,I3,2A)')               WRITE(msgBuf,'(A,I4,A,I4,2A)')
392       &       'DIAGNOSTICS_READPARMS: region=',j,       &       'DIAGNOSTICS_READPARMS: region=',j,
393       &         ' in list l=', l, ', stat_fname: ', stat_fname(l)       &         ' in list l=', l, ', stat_fname: ', stat_fname(l)
394               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
395               WRITE(msgBuf,'(2A,I3,A,I3,2A)')               WRITE(msgBuf,'(2A,I4,A,I4,2A)')
396       &       'DIAGNOSTICS_READPARMS: ==> exceed Max.Nb of regions',       &       'DIAGNOSTICS_READPARMS: ==> exceed Max.Nb of regions',
397       &       '(=',nRegions,' )'       &       '(=',nRegions,' )'
398               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
# Line 410  C-       no region selected => default i Line 410  C-       no region selected => default i
410               diagSt_nbFlds(n) = diagSt_nbFlds(n) + 1               diagSt_nbFlds(n) = diagSt_nbFlds(n) + 1
411               diagSt_Flds(diagSt_nbFlds(n),n) = stat_fields(m,l)               diagSt_Flds(diagSt_nbFlds(n),n) = stat_fields(m,l)
412             ELSEIF ( stat_fields(m,l).NE.blk8c ) THEN             ELSEIF ( stat_fields(m,l).NE.blk8c ) THEN
413               WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
414       &        'Exceed Max.Num. of Fields/list numperlist=', numperlist       &        'Exceed Max.Num. of Fields/list numperlist=', numperlist
415               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
416               WRITE(msgBuf,'(2A,I3,3A,I3,2A)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4,3A,I4,2A)') 'DIAGNOSTICS_READPARMS: ',
417       &        'when trying to add stat_field (m=', m,       &        'when trying to add stat_field (m=', m,
418       &        ' ): ',stat_fields(m,l)       &        ' ): ',stat_fields(m,l)
419               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
420               WRITE(msgBuf,'(2A,I3,2A)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4,2A)') 'DIAGNOSTICS_READPARMS: ',
421       &        ' in list l=', l, ', stat_fname: ', stat_fname(l)       &        ' in list l=', l, ', stat_fname: ', stat_fname(l)
422               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
423               STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'               STOP 'ABNORMAL END: S/R DIAGNOSTICS_READPARMS'
# Line 426  C-       no region selected => default i Line 426  C-       no region selected => default i
426           diagSt_nbLists = diagSt_nbLists + 1           diagSt_nbLists = diagSt_nbLists + 1
427  c        write(6,*) 'stat-list summary:',n,diagSt_nbFlds(n),regionCount  c        write(6,*) 'stat-list summary:',n,diagSt_nbFlds(n),regionCount
428         ELSEIF ( iLen.GE.1 ) THEN         ELSEIF ( iLen.GE.1 ) THEN
429           WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
430       &            'Exceed Max.Num. of list numlists=', numlists       &            'Exceed Max.Num. of list numlists=', numlists
431           CALL PRINT_ERROR( msgBuf , myThid )           CALL PRINT_ERROR( msgBuf , myThid )
432           WRITE(msgBuf,'(2A,I3)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
433       &    'when trying to add stat_list l=', l       &    'when trying to add stat_list l=', l
434           CALL PRINT_ERROR( msgBuf , myThid )           CALL PRINT_ERROR( msgBuf , myThid )
435           WRITE(msgBuf,'(2A,F18.6,2A)') 'DIAGNOSTICS_READPARMS: ',           WRITE(msgBuf,'(2A,F18.6,2A)') 'DIAGNOSTICS_READPARMS: ',
# Line 457  C     Echo History List Data Structure Line 457  C     Echo History List Data Structure
457          WRITE(msgBuf,'(2(A,F18.6))') 'Output Frequency:', freq(n),          WRITE(msgBuf,'(2(A,F18.6))') 'Output Frequency:', freq(n),
458       &                               ' ; Phase: ', phase(n)       &                               ' ; Phase: ', phase(n)
459          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
460          WRITE(msgBuf,'(2(A,F18.6),A,I3)')          WRITE(msgBuf,'(2(A,F18.6),A,I4)')
461       &    ' Averaging Freq.:', averageFreq(n),       &    ' Averaging Freq.:', averageFreq(n),
462       &    ' , Phase: ', averagePhase(n), ' , Cycle:', averageCycle(n)       &    ' , Phase: ', averagePhase(n), ' , Cycle:', averageCycle(n)
463          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)          CALL PRINT_MESSAGE( msgBuf, stdUnit,SQUEEZE_RIGHT, myThid)
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22