anchor.diff:          call read_evops()
           call close_files('r')
   
  +        if(to_exciton_at_output) then
  +        do i=1, size(Evops,5)
  +          call superops_to_exc(Evops(:,:,:,:,i), type)
  +        end do
  +        end if
  +
           call evops_derivative()
           call calculate_coeff()
   
  
positive.code: 
            AAA = matmul(dEvops2,AAA)
  
            call superops_2indexed_to_4indexed(AAA,Evops(:,:,:,:,tind),type)
  
  
          end do
  
negative.code:         allocate(Evops(Nl1,Nl2,Nl1,Nl2,STEPS) )
          allocate(DEvops(Nl1,Nl2,Nl1,Nl2,STEPS) )
  
          rho1    = 0.0_dp
          prhox1  = 0.0_dp
          prhodx1 = 0.0_dp
  
          Evops = 0.0_dp
          DEvops = 0.0_dp
  
      end subroutine init_lindblad_fit
  
  