anchor.diff:      logical :: states
       type(operat) :: rho_sparse
       integer :: istat=0,istat2=0,traj,progress
  -    integer :: i,j,l,m,n
  +    integer :: i,j,l
  +    !integer :: m,n
       real(wp) :: mu,nu,S
       real(wp), allocatable :: p(:)
       ! ITASK  = An index specifying the task to be performed.
  
positive.code:     double precision, optional :: atol(1),rtol(1)
      real(wp), optional :: norm_tol
      integer, intent(in), optional :: ml,mu
      integer :: istat
  
      ode%neq = neq
      if (lzw.ne.0) then
  
negative.code:   end subroutine
  
    subroutine finalize_work
      !integer :: istat=0
      call finalize(psi0)
      call finalize(hamilt)
      call finalize(c_ops)
  