integer i, k, md

      integer ki, i, md

            do md = 1, mdiff

               dkmp5         =      casea(i)  * diffus(i,ki,md) +

     1                         (1.- casea(i)) * blmc  (i,ki,md)

               dstar         = (1.- delta)**2 * dkm1(i,md)

               blmc(i,ki,md) = (1.- delta)*diffus(i,ki,md)

      do md = 1, mdiff

             if(k.ge.kmtj(i))  diffus(i,k,md) = 0.0

         mda(k2) = (md(k2)*drF(k2)+md(k2+1)*drF(k2+1))/

                mda(k2) = (md(k2)*dz(k2)+md(k2+1)*dz(k2+1))/

      integer md, dd, sd, ld, wd

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

      print *,'md,dd,sd       ', md,dd,sd

      integer md, dd, sd, ld, wd

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

      print *,'md,dd,sd       ', md,dd,sd

      integer md, dd, sd, ld, wd

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

      print *,'md,dd,sd       ', md,dd,sd

     I                        listId, md, ndId, ip, im, lm,

     &      ' :',listId, md, ndId, ' ', cdiag(ndId)

      INTEGER listId, md, ndId, ip, im, lm

         DO md = 1,nfields(ld)

          nd = ABS(jdiag(md,ld))

      INTEGER ld,md,nd

     &      'fails to interpolate diag.(#', ndId,'): ',flds(md,listId)

     I                        listId, md, ndId, ip, im, lm,

      INTEGER listId, md, ndId, ip, im, lm

     I                       ld, md,

      INTEGER ld, md, nbLists

      m1 = md

       DO md = 1,nfields(listId)

        ndId = ABS(jdiag(md,listId))

           IF ( mdiag(md,listId).NE.0 ) mDbl = hdiag(nn)

        IF ( idiag(md,listId).NE.0 .AND. gcode(5:5).NE.'D' ) THEN

          ip = ABS(idiag(md,listId)) + kdiag(ndId)*(lm-1)

          im = mdiag(md,listId)

     &       '- WARNING -   diag.#',ndId, ' : ',flds(md,listId),

     &       ' (#',md,' ) in outp.Stream: ',fnames(listId)

     I                         listId, md, ndId, ip, im, lm,

     I                         listId, md, ndId, ip, im, lm,

             IF ( flds(md,listId).EQ.'PhiVEL  '

     &       .OR. flds(md,listId).EQ.'PsiVEL  '

     I                         listId, md, ndId, ip, im, lm,

           nRec = md + (lm-1)*nfields(listId)

      INTEGER md, ndId, nn, ip, im

         DO md = 1,nfields(ld),10

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

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

       DO md=1,numperList

        idiag(md,ld) = 0

        jdiag(md,ld) = 0

        mdiag(md,ld) = 0

       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'

       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

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

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

          DO md=1,nfields(ld)

            nd = ABS(jdiag(md,ld))

          DO md=1,nfields(ld)

            nd = ABS(jdiag(md,ld))

      INTEGER md,ld,nd

       DO md=1,nfields(ld)

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

         IF ( j.GE.1 .AND. j.NE.md-1 ) THEN

           IF ( j.LT.md-1 ) THEN

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

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

             DO i=j,md-2

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

           ELSEIF ( j.GT.md ) THEN

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

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

             DO i=j,md+1,-1

             flds(md,ld) = cdiag(mate)

     I                        listId, md, ndId, ip, im, lm,

      INTEGER listId, md, ndId, ip, im, lm

        wFac = jdiag(md,listId).LT.0

         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))

            im = mdiag(md,ld)

         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

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

      INTEGER md,ld,nd

       DO md=1,numperlist

        iSdiag(md,ld) = 0

        jSdiag(md,ld) = 0

        mSdiag(md,ld) = 0

       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'

      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