anchor.diff:                      BB = 1.0/300,    &
                       CC = 1.0/100
   
  -        allocate_SVD = .false.
           STEPS = 500
           call init_lindblad_fit()
   
  
positive.code:         allocate(prhodx1(Nl1,Nl2))
  
          allocate(Evops(Nl1,Nl2,Nl1,Nl2,STEPS) )
  
          rho1    = 0.0_dp
          prhox1  = 0.0_dp
          prhodx1 = 0.0_dp
  
          Evops = 0.0_dp
  
      end subroutine init_lindblad_fit
  
  
negative.code:     ! declarations
      real(dp), dimension(:,:), allocatable  :: Ham
      complex(dpc), allocatable:: rho1(:,:), prhox1(:,:), prhodx1(:,:)
      complex(dpc), dimension(:,:,:,:,:), allocatable      :: Evops, Devops
      complex(dpc), dimension(:,:), allocatable            :: Evops2, Devops2, AAA, VT, U
      real(dp), dimension(:), allocatable                  :: eigval
      integer(i4b)         :: Nl1, Nl2, Nl
      character, parameter :: type = 'E'
      real(dp)             :: timeStep = 0
  