integer md, dd, sd, ld, wd

      call CAL_CONVDATE( daydate,yday,md,dd,sd,ld,wd,mythid )

      print *,'ld,wd          ', ld,wd

      integer md, dd, sd, ld, wd

      call CAL_CONVDATE( daydate,yday,md,dd,sd,ld,wd,mythid )

      print *,'ld,wd          ', ld,wd

      integer md, dd, sd, ld, wd

      call CAL_CONVDATE( daydate,yday,md,dd,sd,ld,wd,mythid )

      print *,'ld,wd          ', ld,wd

      DO ld = 2,nlists

       DO m = 1,ld-1

        IF ( fnames(ld).EQ.fnames(m) ) THEN

     &    ' 2nd (n=', ld, ' ): ', fnames(ld)

      DO ld = 2,diagSt_nbLists

       DO m = 1,ld-1

        IF ( diagSt_Fname(ld).EQ.diagSt_Fname(m) ) THEN

     &    ' 2nd (n=', ld, ' ): ', diagSt_Fname(ld)

      DO ld = 1,nlists

       IF ( fflags(ld)(2:2).EQ.'P' ) THEN

          DO k=1,nlevels(ld)

           IF ( levs(k,ld)-MAX(rF(1),rF(Nr+1)).GT.margin

     &     .OR. levs(k,ld)-MIN(rF(1),rF(Nr+1)).LT.-margin

     &     .OR. levs(k,ld).LE.0. ) THEN

     &       'Vertical Interp. for list l=', ld,

     &       ', filename: ', fnames(ld)

     &       ' lev(k=', k, ') p=', levs(k,ld),

     &       ' for list l=', ld, ', filename: ', fnames(ld)

         DO md = 1,nfields(ld)

          nd = ABS(jdiag(md,ld))

          CALL DIAGNOSTICS_GET_POINTERS( 'RSURF   ', ld,

           CALL DIAGNOSTICS_GET_POINTERS('FIZPRES ', ld,

           CALL DIAGNOSTICS_GET_POINTERS('RCENTER ', ld,

     &      'to interpolate diags from output list:', ld

      INTEGER ld,md,nd

     I                       ld, md,

      INTEGER ld, md, nbLists

      DO n = ld,nbLists

     &             'Updated list in Outp.Stream #', ld, ' :'

         DO md = 1,nfields(ld),10

           j = MIN(nfields(ld),md+9)

           WRITE(msgBuf,'(21A)') ' Fields:   ',(' ',flds(i,ld),i=md,j)

      DO ld=1,numLists

        idiag(md,ld) = 0

        jdiag(md,ld) = 0

        mdiag(md,ld) = 0

      DO ld=1,nlists

       nActive(ld) = nfields(ld)

       DO md=1,nfields(ld)

          IF ( flds(md,ld).EQ.cdiag(nd) ) THEN

            CALL DIAGNOSTICS_SETDIAG(mate,ndiagcount,md,ld,nd,myThid)

     I                      ld, md, nlists, nfields, flds, myThid )

     &                      flds(md,ld),' is not a Diagnostic'

       nActiveMax = MAX(nActive(ld),nActiveMax)

      DO ld=1,nlists

       DO md=1,nActive(ld)

        IF (mdiag(md,ld).EQ.0 ) THEN

         k = SIGN(1,jdiag(md,ld))

         nd = ABS(jdiag(md,ld))

            IF ( mdiag(md,ld).EQ.0 .AND. (k*jdiag(i,j)).EQ.mate ) THEN

             IF ( freq(j).EQ.freq(ld) .AND. phase(j).EQ.phase(ld)

     &           .AND. averageFreq(j) .EQ.averageFreq(ld)

     &           .AND. averagePhase(j).EQ.averagePhase(ld)

     &           .AND. averageCycle(j).EQ.averageCycle(ld) )

     &          mdiag(md,ld) = ABS(idiag(i,j))

         IF ( mdiag(md,ld).NE.0 ) THEN

      DO ld=1,nlists

        IF ( nlevels(ld).EQ.-1 ) THEN

          DO md=1,nfields(ld)

            nd = ABS(jdiag(md,ld))

     &      ' setting levs of list l=',ld,', fnames=', fnames(ld)

     &      ' setting levs of list l=',ld,', fnames=', fnames(ld)

          nlevels(ld) = kLev

           levs(k,ld) = k

     &      'Set levels for Outp.Stream: ',fnames(ld)

          IF ( fflags(ld)(2:2).EQ.'I' ) suffix = ' Sum Levels:'

          DO k1=1,nlevels(ld),20

            k2 = MIN(nlevels(ld),k1+19)

            WRITE(msgBuf,'(A,20F5.0)') suffix, (levs(k,ld),k=k1,k2)

        ELSEIF ( fflags(ld)(2:2).NE.'P' ) THEN

          DO k=1,nlevels(ld)

            kLev = MAX(NINT(levs(k,ld)),kLev)

          DO md=1,nfields(ld)

            nd = ABS(jdiag(md,ld))

     &         ' in list l=', ld, ', filename: ', fnames(ld)

      INTEGER md,ld,nd

      DO ld=1,nlists

       DO md=1,nfields(ld)

          IF ( nd.EQ.0 .AND. flds(md,ld).EQ.cdiag(i) ) nd = i

               DO i=1,nfields(ld)

                 IF ( j.EQ.0 .AND. flds(i,ld).EQ.cdiag(mate) ) j = i

     &             'Re-Order Diags in Outp.Stream: ',fnames(ld)

     &         ' move ',flds(j,ld),' from ',j,' down to',md-1,

     &         ' just before ',flds(md,ld)

               flds(i,ld) = flds(i+1,ld)

             flds(md-1,ld) = cdiag(mate)

     &         ' move ',flds(j,ld),' from ',j,'  up to ',md,

     &         ' just before ',flds(md,ld)

               flds(i,ld) = flds(i-1,ld)

             flds(md,ld) = cdiag(mate)

         IF ( fflags(ld)(2:2).EQ.'P' ) THEN

          DO k1=1,nlevels(ld),10

           k2 = MIN(nlevels(ld),k1+9)

           WRITE(dUnit,'(A,1P10E10.3)')' interp:', (levs(k,ld),k=k1,k2)

          DO k1=1,nlevels(ld),25

           k2 = MIN(nlevels(ld),k1+24)

           WRITE(dUnit,'(A,25I4)')' levels:',(NINT(levs(k,ld)),k=k1,k2)

         DO md=1,nActive(ld)

           ndId = ABS(jdiag(md,ld))

     &       jdiag(md,ld),' |', flds(md,ld),'|',idiag(md,ld),' |',

     &       mdiag(md,ld),' |', kdiag(ndId),' |'

           ipt = ABS(idiag(md,ld))

           IF (ipt.NE.0 .AND. averageCycle(ld).GT.1) THEN

             DO l=1,averageCycle(ld)

               IF (l.EQ.averageCycle(ld)) THEN

            im = mdiag(md,ld)

        DO ld=1,diagSt_nbLists

         iLen = ILNBLNK(diagSt_Fname(ld))

         WRITE(dUnit,'(A,I4,2A)') 'listId=', ld,

     &                     ' ; file name: ',diagSt_Fname(ld)(1:iLen)

     &    diagSt_nbFlds(ld), '  |',diagSt_nbActv(ld), '  |',

     &    diagSt_freq(ld), diagSt_phase(ld), ' |'

           IF ( diagSt_region(j,ld).GE.1

         DO md=1,diagSt_nbActv(ld)

     &       jSdiag(md,ld),' |', diagSt_Flds(md,ld),'|',iSdiag(md,ld),

     &       ' |', mSdiag(md,ld),' |'

           ipt = ABS(iSdiag(md,ld))

           im = mSdiag(md,ld)

      INTEGER md, ld, ndId, ipt, im

        DO ld=1,nlists

         iLen = ILNBLNK(fnames(ld))

         WRITE(dUnit,'(A,I5,2A)') 'listId=', ld,

     &                            ' ; file name: ',fnames(ld)(1:iLen)

     &    nfields(ld), '  |',nActive(ld), '  |',

     &    freq(ld), phase(ld), ' |', nlevels(ld)

            mm = diagSt_nbActv(ld) + 1

             iSdiag(mm,ld) = SIGN(mSdiag(md,ld),mate)

             jSdiag(mm,ld) = mate

             diagSt_Flds(mm,ld) = cdiag(mate)

            diagSt_nbActv(ld) = mm

       nActiveMax = MAX(diagSt_nbActv(ld),nActiveMax)

      INTEGER md,ld,nd

      DO ld=1,numlists

        iSdiag(md,ld) = 0

        jSdiag(md,ld) = 0

        mSdiag(md,ld) = 0

      DO ld=1,diagSt_nbLists

       diagSt_nbActv(ld) = diagSt_nbFlds(ld)

       DO md=1,diagSt_nbFlds(ld)

          IF ( diagSt_Flds(md,ld).EQ.cdiag(nd) ) THEN

            CALL DIAGSTATS_SETDIAG(mate,ndiagcount,md,ld,nd,myThid)

            jSdiag(md,ld) = nd

     I                      ld,md, diagSt_nbLists,

     &                     diagSt_Flds(md,ld),' is not a Diagnostic'

           DO l=1,ld-1

     &                 .OR. (diagSt_region(j,l).LT.diagSt_region(j,ld))

      integer md, dd, sd, ld, wd

           call CAL_CONVDATE( daydate,yday,md,dd,sd,ld,wd,mythid )

             call CAL_CONVDATE( daydate,yday,md,dd,sd,ld,wd,mythid )

             call CAL_CONVDATE( daydate,yday,md,dd,sd,ld,wd,mythid )

      integer md, dd, sd, ld, wd

      call CAL_CONVDATE( daydate,yday,md,dd,sd,ld,wd,mythid )

      integer md, dd, sd, ld, wd

      integer md, dd, sd, ld, wd