anchor.diff:            A(super2,super1) = calc(super1,i,j,k,l)
   
             if(super1 == 1) then
  -            B(super2) = Uee(i,j,k,l,tind)
  +            B(super2) = Evops(i,j,k,l,tind)
             end if
           end do
           end do
  
positive.code:         allocate(rho1(Nl1,Nl2))
          allocate(prhox1(Nl1,Nl2))
          allocate(prhodx1(Nl1,Nl2))
  
          allocate(A(Nl1*Nl2*Nl1*Nl2*STEPS,Nl1*Nl2*Nl1*Nl2*Nbasis))
          allocate(U(Nl1*Nl2*Nl1*Nl2*STEPS,Nl1*Nl2*Nl1*Nl2*STEPS))
  
negative.code:     ! declarations
      real(dp), dimension(:,:), allocatable  :: Ham
      complex(dpc), allocatable:: rho1(:,:), prhox1(:,:), prhodx1(:,:)
      complex(dpc), dimension(:,:,:,:,:), allocatable      :: Evops
      complex(dpc), dimension(:,:), private, allocatable   :: A, U, VT
      complex(dpc), dimension(:), private, allocatable     :: B
      real(dp), dimension(:), private, allocatable         :: EIGVAL
  