anchor.diff:  
   
           call init_lindblad_fit()
  +        call init_nakajima_zwanzig_shared(Ham)
   
           !write(*,*) indices_to_superindex(1,2,2,1,351)
           !call superindex_to_indices(indices_to_superindex(1,2,2,1,351), r,s,k,l,i)
           !write(*,*) r,s,k,l,i
           !write(*,*) 1,2,2,1,351
  +
  +        call open_files('E')
  +        call read_evops('E', Uee)
  +        call close_files('E')
   
           call create_hustomatrix(A,B)
   
  
positive.code:       do Utnemele=1,N1_from_type(type)
        do Utnemele2=1,N2_from_type(type)
  
  
          write(ind(Uelement,Uelement2,Utnemele,Utnemele2,.false.) &
          ,*) timeStep*i,' ',real(actual_U(Uelement,Uelement2,Utnemele,Utnemele2,i)),' ',&
                               aimag(actual_U(Uelement,Uelement2,Utnemele,Utnemele2,i))
  
        end do
        end do
        end do
        end do
        end do
      end subroutine write_evolution_operators
  
      subroutine open_files(type)
        character, intent(in) :: type
  
negative.code:           A(super1,super2) = A(super1,super2) + calc(super1,i,j,k,l)*conjg(calc(super2,i,j,k,l))
  
            if(super1 == 1) then
              B(super2) = B(super2) + Uee(i,j,k,l,tind)*conjg(calc(super2,i,j,k,l))
            end if
          end do
          end do
  