anchor.diff:      ! Do all the simulations within the module
       !
       subroutine do_lindblad_fit_work()
  -        integer(i4b) :: r,s,k,l,i
  +        integer(i4b) :: i
   
           call init_lindblad_fit()
   
  
positive.code:         call read_evops()
          call close_files('r')
  
          call evops_derivative()
          call calculate_coeff()
  
  
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
  
  