anchor.diff:  c
   c.... close varts file for probes
   c
  -      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
  +      call finalizeTimeSeries()
   
         if (numpe > 1) call MPI_BARRIER(MPI_COMM_WORLD, ierr)
         if(myrank.eq.0)  then
  
positive.code:       subroutine dumpTimeSeries()
        use timedata   !allows collection of time series
        include "common.h"
         character*60    fvarts
         character*10    cname2
     
  
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)
  
  