integer          ido, m, n, nev, ncv, lworkl, info, ierr,

      n = admtlmrec

      if ( n .gt. maxn ) then

         CALL DSAUPD ( ido, bmat, n, which, nev, tol, resid, 

            do l = 1, n

            do l = 1, n

     &        bmat, n, which, nev, tol, resid, ncv, v, ldv, 

               do l = 1, n

         print *, ' Size of the matrix is ', n

      INTEGER n

      do i=1,4*n,4

      SUBROUTINE MDS_BYTESWAPR8( n, arr )

      INTEGER n

      do i=1,8*n,8

      SUBROUTINE MDS_BYTESWAPR4( n, arr )

        n=1

     &                          n,ijROatm(n),ijROocn(n),arROmap(n)

        n=nROmap

     &                          n,ijROatm(n),ijROocn(n),arROmap(n)

      INTEGER n, ijo, ija

        DO n=1,ROsize

          ijROocn(n)=0

          ijROatm(n)=0

          arROmap(n)=0.

        DO n=1,LEN( runoffmapFile )

         IF ( runoffmapFile(n:n).NE.' ' ) lengthName=n

        DO n=1,nROmap

         iRec = n

         tmpfld(1,n)=r8seg(1)

         tmpfld(2,n)=r8seg(2)

         tmpfld(3,n)=r8seg(3)

        DO n=1,nROmap

          ija=NINT(tmpfld(1,n))

          ijo=NINT(tmpfld(2,n))

          ijROatm(n)=ija

            arROmap(n)=tmpfld(3,n)/rAc(ijo)

          ijROocn(n)=ijo

      INTEGER n, ijo, ija

        DO n=1,nROmap

          ijo=ijROocn(n)

          ija=ijROatm(n)

     &                -fld_atm(ija)*arROmap(n)

      INTEGER  n, iL, nUnit

        DO n=1,diagSt_nbLists

          nUnit = diagSt_ioUnit(n)

          iL = ILNBLNK(diagSt_Fname(n))

     &          diagSt_Fname(n)(1:iL), '.', nIter0, '.txt'

          DO n = 1,nlists

            DO m = 1,nfields(n)

              WRITE(dUnit,'(I10)') ndiag(jdiag(m,n))

      INTEGER dUnit, n, m

      INTEGER n,m

      DO n = 2,nlists

       DO m = 1,n-1

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

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

      DO n = 2,diagSt_nbLists

       DO m = 1,n-1

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

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

        DO n=1,ndiagt

         IF ( MOD(n,100).EQ.0 ) THEN

         l = ILNBLNK(tdiag(n))

           WRITE(dUnit,'(I4,3A,I3,6A)') n,' |',cdiag(n),'|',

     &           kdiag(n),' |',gdiag(n),'|',udiag(n),'|',tdiag(n)(1:l)

           WRITE(dUnit,'(I4,3A,I3,6A)') n,' |',cdiag(n),'|',

     &           kdiag(n),' |',gdiag(n),'|',udiag(n),'|'

      DO n = 2,ndiagt

       IF ( cdiag(n).NE.blk8c ) THEN

        DO l = 1,n-1

         IF ( cdiag(l).EQ.cdiag(n) ) THEN

     &       'diag.Name: ',cdiag(n),' registered 2 times :'

     &       ' 2nd (n=', n, ' ), title= ',tdiag(n)

      DO n=1,nlists

        IF ( nlevels(n).GT.nlevs ) THEN

     &       ' in list n=', n, ', filename: ', fnames(n)

     &       ' number of lev= ', nlevels(n), ' exceeds Max=',nlevs

      INTEGER l, n, nlevs, nGroundLev

      DO n = 1,ndiagt

        IF     (gdiag(n)(10:10) .EQ. '0') THEN

          kdiag(n) = 0

        ELSEIF (gdiag(n)(10:10) .EQ. '1') THEN

          kdiag(n) = 1

        ELSEIF (gdiag(n)(10:10) .EQ. 'R') THEN

          kdiag(n) = Nr

        ELSEIF (gdiag(n)(10:10) .EQ. 'L') THEN

          kdiag(n) = nlevs

        ELSEIF (gdiag(n)(10:10) .EQ. 'M') THEN

          kdiag(n) = nlevs - 1

        ELSEIF (gdiag(n)(10:10) .EQ. 'G') THEN

          kdiag(n) = nGroundLev

        ELSEIF (gdiag(n)(10:10) .NE. ' ') THEN

          kdiag(n) = 1

          DO n = 1,nlists

            DO m = 1,nfields(n)

              READ(dUnit,'(I10)') ndiag(jdiag(m,n))

      INTEGER dUnit, n, m

      DO n=0,nRegions

       IF (region2fill(n).NE.0) THEN

         kRegMsk = diagSt_kRegMsk(n)

     I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),

     I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),

     I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),

     I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),

     I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),

     I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),

     I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),

     I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),

            IF ( statFld(0,n).LE.0. ) THEN

              statFld(im,n) = statLoc(im)

              statFld(ix,n) = statLoc(ix)

              statFld(im,n) = MIN( statFld(im,n), statLoc(im) )

              statFld(ix,n) = MAX( statFld(ix,n), statLoc(ix) )

              statFld(i,n) = statFld(i,n) + statLoc(i)

      INTEGER i, n, km, kFr, kRegMsk

      integer n

      do n=1,ndiagMax

        kdiag(n) = 0

        cdiag(n) = '        '

        gdiag(n) = '                '

        udiag(n) = '                '

        tdiag(n) = ' '

      INTEGER i, j, n

      DO n=0,nStats

        statArr(n) = 0.

           ndId = jdiag(m,n)

           ipt = ipt + pdiag(n,bi,bj)*kdiag(ndId)

           IF ( mdiag(m,n).EQ.0 ) mate = 0

      DO n=1,diagSt_nbLists

       DO m=1,diagSt_nbActv(n)

        IF ( chardiag.EQ.diagSt_Flds(m,n) .AND. iSdiag(m,n).GT.0 ) THEN

         iSp = iSdiag(m,n)

           ndId = jSdiag(m,n)

           IF ( mSdiag(m,n).EQ.0 ) mate = 0

            region2fill(j) = diagSt_region(j,n)

      INTEGER m, n, j, k, l, bi, bj

      DO n=1,nlists

       DO m=1,nActive(n)

        IF ( chardiag.EQ.flds(m,n) .AND. idiag(m,n).GT.0 ) THEN

         ipt = idiag(m,n)

      INTEGER m, n

      DO n=1,nlists

       DO m=1,nActive(n)

        IF ( chardiag.EQ.flds(m,n) .AND. idiag(m,n).GT.0 ) THEN

         ipt = idiag(m,n)

      INTEGER i,j,n

      n = LEN(diagUnitsInPieces)

      DO i=1,n

      INTEGER i,j,n

      n = LEN(diagTitleInPieces)

      DO i=1,n

      INTEGER n,m

      DO n=1,nlists

       DO m=1,nActive(n)

        IF ( diagName.EQ.flds(m,n) .AND. idiag(m,n).NE.0 ) THEN

          ipoint = abs(idiag(m,n))

          jpoint = jdiag(m,n)

      DO n=1,diagSt_nbLists

       DO m=1,diagSt_nbActv(n)

        IF ( chardiag.EQ.diagSt_Flds(m,n) .AND. iSdiag(m,n).GT.0 ) THEN

         iSp = iSdiag(m,n)

           ndId = jSdiag(m,n)

            region2fill(j) = diagSt_region(j,n)

      INTEGER m, n, j, k, l, bi, bj

      DO n=1,nlists

       DO m=1,nActive(n)

        IF ( chardiag.EQ.flds(m,n) .AND. idiag(m,n).GT.0 ) THEN

         ipt = idiag(m,n)

           ndId = jdiag(m,n)

           ipt = ipt + pdiag(n,bi,bj)*kdiag(ndId)

      DO n=1,diagSt_nbLists

       DO m=1,diagSt_nbActv(n)

        IF ( chardiag.EQ.diagSt_Flds(m,n) .AND. iSdiag(m,n).GT.0 ) THEN

         iSp = iSdiag(m,n)

           ndId = jSdiag(m,n)

            region2fill(j) = diagSt_region(j,n)

      INTEGER m, n, j, k, l, bi, bj

      DO n=1,nlists

       DO m=1,nActive(n)

        IF ( chardiag.EQ.flds(m,n) .AND. idiag(m,n).GT.0 ) THEN

         ipt = idiag(m,n)

           ndId = jdiag(m,n)

           ipt = ipt + pdiag(n,bi,bj)*kdiag(ndId)

      integer n, ndiagnum, ipointer

      do n = 1,ndiagt

        if(chardiag.eq.cdiag(n)) then

          ndiagnum = n

      DO n = 1,numlists

        freq(n) = 0.

        phase(n) = 0.

        averageFreq(n)  = 0.

        averagePhase(n) = 0.

        averageCycle(n) = 1

        nlevels(n) = 0

        nfields(n) = 0

        fnames(n) = blkFilName

          levs(k,n) = 0

          flds(m,n) = blk8c

        fflags(n)   = blk8c

         n = nlists + 1

         freq(n)    = frequency(l)

           phase(n) = timePhase(l)

           phase(n) = -0.5 _d 0 * frequency(l)

           averageFreq(n)  = averagingFreq(l)

           averagePhase(n) = averagingPhase(l)

           averageCycle(n) = repeatCycle(l)

           averageFreq(n)  = nTimeSteps*deltaTClock

           averagePhase(n) = phase(n)

           averageFreq(n)  = frequency(l)

           averagePhase(n) = phase(n)

         fnames(n)  = filename (l)

         fflags(n)  = fileflags(l)

         nlevels(n) = 0

     &            nlevels(n).LT.numLevels ) THEN

               nlevels(n) = nlevels(n) + 1

               levs(nlevels(n),n) = levels(k,l)

           nlevels(n) = -1

         nfields(n) = 0

     &          nfields(n).LT.numperlist ) THEN

             nfields(n) = nfields(n) + 1

             flds(nfields(n),n) = fields(m,l)

      DO n = 1,numlists

        diagSt_freq(n) = 0.

        diagSt_phase(n) = 0.

        diagSt_nbFlds(n) = 0

        diagSt_ioUnit(n) = 0

        diagSt_Fname(n) = blkFilName

          diagSt_region(j,n) = 0

          diagSt_Flds(m,n) = blk8c

         n = diagSt_nbLists + 1

         diagSt_freq(n) = stat_freq(l)

           diagSt_phase(n) = stat_phase(l)

           diagSt_phase(n) = -0.5 _d 0 * stat_freq(l)

         diagSt_Fname(n)  = stat_fname(l)

            IF ( diagSt_region(j,n).EQ.0 ) THEN

             diagSt_region(j,n) = 1

           diagSt_region(0,n) = 1

         diagSt_nbFlds(n) = 0

     &          diagSt_nbFlds(n).LT.numperlist ) THEN

             diagSt_nbFlds(n) = diagSt_nbFlds(n) + 1

             diagSt_Flds(diagSt_nbFlds(n),n) = stat_fields(m,l)

      DO n = 1,nlists

        WRITE(msgBuf,'(2a)') 'Creating Output Stream: ', fnames(n)

        WRITE(msgBuf,'(2(A,F17.6))') 'Output Frequency:', freq(n),

     &                               ' ; Phase: ', phase(n)

     &    ' Averaging Freq.:', averageFreq(n),

     &    ' , Phase: ', averagePhase(n), ' , Cycle:', averageCycle(n)

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

         DO l=1,nlevels(n),20

          m = MIN(nlevels(n),l+19)

          WRITE(msgBuf,'(A,20F5.0)') ' Levels:    ', (levs(k,n),k=l,m)

        DO nf = 1,nfields(n),10

        m = MIN(nfields(n),nf+9)

        WRITE(msgBuf,*) 'Fields:   ',(' ',flds(l,n),l=nf,m)

      DO n = 1,diagSt_nbLists

     &                       diagSt_Fname(n)

     &               diagSt_freq(n), ' ; Phase: ', diagSt_phase(n)

         IF ( diagSt_region(j,n).GE.1 ) THEN

     &                 (' ',diagSt_Flds(l,n),l=1,diagSt_nbFlds(n))

      INTEGER j,k,l,n,m,nf

       IF ( averageCycle(n).GT.1 ) THEN

        xInterval = myTime + deltaTclock*0.5 _d 0 - averagePhase(n)

        xInterval = xInterval / averageFreq(n)

          nInterval = 1 + INT( -xInterval/FLOAT(averageCycle(n)) )

          nInterval = nInterval*averageCycle(n) + INT(xInterval) - 1

        nInterval = MOD(nInterval,averageCycle(n))

        IF (dBugFlag.AND.pdiag(n,1,1).NE.nInterval)

     &      myIter, ', list:', n, ' switch',

     &      pdiag(n,1,1),' ->', nInterval, xInterval

        IF ( nInterval.LT.0 .OR. nInterval.GE.averageCycle(n) ) THEN

     &       ' error setting pdiag(n=',n,') to:', nInterval

     &       ' cycle=', averageCycle(n), ', xInt=', xInterval

           pdiag(n,bi,bj) = nInterval

      DO n = 1,diagSt_nbLists

       IF ( diagSt_freq(n).LT.0. ) THEN