anchor.diff:        subroutine dumpTimeSeries()
         use timedata   !allows collection of time series
         include "common.h"
  +      include "mpif.h"
          character*60    fvarts
          character*10    cname2
      
  
positive.code:        character*60    fvarts
         character*10    cname2
  
  
          inquire(file='xyzts.dat',exist=exts)
          if(exts) then
             
             open(unit=626,file='xyzts.dat',status='old')
             read(626,*) ntspts, freq, tolpt, iterat, varcod
             call sTD             ! sets data structures
  
negative.code: c
        return
        end subroutine
  
        subroutine finalizeTimeSeries()
        use timedata   !allows collection of time series
        include "common.h"
        if(exts) then
          do jj=1,ntspts
            if (myrank == iv_rank(jj)) then
              close(1000+jj)
            endif
          enddo
          call dTD   ! deallocates time series arrays
        endif
        return
        end subroutine
  
  
  
         subroutine initEQS(iBC, rowp, colm)
  
  