anchor.diff:      ! declarations
       real(dp), dimension(:,:), allocatable  :: Ham
       complex(dpc), allocatable:: rho1(:,:), prhox1(:,:), prhodx1(:,:)
  -    !complex(dpc), allocatable:: rho2(:,:), prhox2(:,:), prhodx2(:,:)
       complex(dpc), dimension(:,:,:,:,:), allocatable      :: Ueg, Uee, Ugg
       complex(dpc), dimension(:,:), private, allocatable   :: A, U, VT
       complex(dpc), dimension(:), private, allocatable     :: B
  
positive.code:     private::read_config_file
      private::Lmult1
      private::propagate1
      !private::Lmult2
      !private::propagate2
  
      private::read_Nsys
      private::ind
  
negative.code:         end do
        end if
  
        result = result * arbitrary_basis_element(t)
  
      end subroutine Lmult1
  
  