/[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.33 by jmc, Wed Feb 6 21:25:26 2013 UTC revision 1.34 by jmc, Wed Aug 14 00:57:33 2013 UTC
# Line 104  C--   per level statistics output: Line 104  C--   per level statistics output:
104       &     set_regMask, val_regMask,       &     set_regMask, val_regMask,
105       &     diagSt_regMaskFile, nSetRegMskFile       &     diagSt_regMaskFile, nSetRegMskFile
106    
107  C     Initialize and Read Diagnostics Namelist  C-    PKG diagnostics is not used
108          IF ( .NOT.useDiagnostics ) THEN
109            _BEGIN_MASTER(myThid)
110    C-    Track diagnostics pkg activation status:
111             diag_pkgStatus = -1
112    C     print a (weak) warning if data.diagnostics is found
113            _END_MASTER(myThid)
114            RETURN
115          ENDIF
116    
117    C-    Initialize and Read Diagnostics Namelist
118        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
119    
120        blk8c  = '        '        blk8c  = '        '
# Line 130  c       missing_value(l)     = UNSET_RL Line 140  c       missing_value(l)     = UNSET_RL
140            levels(k,l) = UNSET_RL            levels(k,l) = UNSET_RL
141          ENDDO          ENDDO
142          DO m = 1,fdimLoc          DO m = 1,fdimLoc
143            fields(m,l) = blk8c            fields(m,l) = blkName
144          ENDDO          ENDDO
145        ENDDO        ENDDO
146        diagLoc_ioUnit = 0        diagLoc_ioUnit = 0
       settingDiags = .FALSE.  
147        dumpAtLast   = .FALSE.        dumpAtLast   = .FALSE.
148        diag_mnc     = useMNC        diag_mnc     = useMNC
149        useMissingValue = .FALSE.        useMissingValue = .FALSE.
# Line 164  c       missing_value(l)     = UNSET_RL Line 173  c       missing_value(l)     = UNSET_RL
173            stat_region(k,l) = UNSET_I            stat_region(k,l) = UNSET_I
174          ENDDO          ENDDO
175          DO m = 1,fdimLoc          DO m = 1,fdimLoc
176            stat_fields(m,l) = blk8c            stat_fields(m,l) = blkName
177          ENDDO          ENDDO
178        ENDDO        ENDDO
179    C-    Track diagnostics pkg activation status:
180          diag_pkgStatus = 1
181    
182        WRITE(msgBuf,'(2A)')        WRITE(msgBuf,'(2A)')
183       &     ' DIAGNOSTICS_READPARMS: opening data.diagnostics'       &     ' DIAGNOSTICS_READPARMS: opening data.diagnostics'
# Line 222  c       misValInt(n) = UNSET_I Line 233  c       misValInt(n) = UNSET_I
233            levs(k,n) = 0            levs(k,n) = 0
234          ENDDO          ENDDO
235          DO m = 1,numperList          DO m = 1,numperList
236            flds(m,n) = blk8c            flds(m,n) = blkName
237          ENDDO          ENDDO
238          fflags(n)   = blk8c          fflags(n)   = blk8c
239        ENDDO        ENDDO
# Line 311  C-       will set levels later, once the Line 322  C-       will set levels later, once the
322           ENDIF           ENDIF
323           nfields(n) = 0           nfields(n) = 0
324           DO m=1,fdimLoc           DO m=1,fdimLoc
325             IF ( fields(m,l).NE.blk8c .AND.             IF ( fields(m,l).NE.blkName .AND.
326       &          nfields(n).LT.numperList ) THEN       &          nfields(n).LT.numperList ) THEN
327               nfields(n) = nfields(n) + 1               nfields(n) = nfields(n) + 1
328               flds(nfields(n),n) = fields(m,l)               flds(nfields(n),n) = fields(m,l)
329             ELSEIF ( fields(m,l).NE.blk8c ) THEN             ELSEIF ( fields(m,l).NE.blkName ) THEN
330               WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
331       &        'Exceed Max.Num. of Fields/list numperList=', numperList       &        'Exceed Max.Num. of Fields/list numperList=', numperList
332               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
# Line 367  C-    Initialise DIAG_STATIS common bloc Line 378  C-    Initialise DIAG_STATIS common bloc
378            diagSt_region(j,n) = 0            diagSt_region(j,n) = 0
379          ENDDO          ENDDO
380          DO m = 1,numperList          DO m = 1,numperList
381            diagSt_Flds(m,n) = blk8c            diagSt_Flds(m,n) = blkName
382          ENDDO          ENDDO
383        ENDDO        ENDDO
384    
# Line 446  C-       no region selected => default i Line 457  C-       no region selected => default i
457           ENDIF           ENDIF
458           diagSt_nbFlds(n) = 0           diagSt_nbFlds(n) = 0
459           DO m=1,fdimLoc           DO m=1,fdimLoc
460             IF ( stat_fields(m,l).NE.blk8c .AND.             IF ( stat_fields(m,l).NE.blkName .AND.
461       &          diagSt_nbFlds(n).LT.numperList ) THEN       &          diagSt_nbFlds(n).LT.numperList ) THEN
462               diagSt_nbFlds(n) = diagSt_nbFlds(n) + 1               diagSt_nbFlds(n) = diagSt_nbFlds(n) + 1
463               diagSt_Flds(diagSt_nbFlds(n),n) = stat_fields(m,l)               diagSt_Flds(diagSt_nbFlds(n),n) = stat_fields(m,l)
464             ELSEIF ( stat_fields(m,l).NE.blk8c ) THEN             ELSEIF ( stat_fields(m,l).NE.blkName ) THEN
465               WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',               WRITE(msgBuf,'(2A,I4)') 'DIAGNOSTICS_READPARMS: ',
466       &        'Exceed Max.Num. of Fields/list numperList=', numperList       &        'Exceed Max.Num. of Fields/list numperList=', numperList
467               CALL PRINT_ERROR( msgBuf , myThid )               CALL PRINT_ERROR( msgBuf , myThid )
Legend:
Removed from v.1.33  
changed lines
  Added in v.1.34
  ViewVC Help
Powered by ViewVC 1.1.22