anchor.diff:    end subroutine
   
     subroutine finalize_work
  -    integer :: istat=0
  +    !integer :: istat=0
       call finalize(psi0)
       call finalize(hamilt)
       call finalize(c_ops)
  
positive.code:     logical :: states
      type(operat) :: rho_sparse
      integer :: istat=0,istat2=0,traj,progress
      integer :: i,j,l,m,n
      real(wp) :: mu,nu,S
      real(wp), allocatable :: p(:)
      ! ITASK  = An index specifying the task to be performed.
  
negative.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
  