anchor.diff:        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
  +        read(ind(Uelement,Uelement2,Utnemele,Utnemele2,.false.),*, IOSTAT=file_ios) time, a, b
  +        actual_U(Uelement,Uelement2,Utnemele,Utnemele2,i) = a + b * cmplx(0,1)
  +
  +        if(i == 1 .and. Uelement == 1 .and. Uelement2 == 1 .and. Utnemele == 1 .and. Utnemele2 == 1) then
  +          timeStep = time
  +        elseif(i == 2 .and. Uelement == 1 .and. Uelement2 == 1 .and. Utnemele == 1 .and. Utnemele2 == 1) then
  +          timeStep = time - timeStep
  +          write(*,*) 'timeStep', timeStep
  +        end if
  +
  +      end do
  +      end do
  +      end do
  +      end do
  +      end do
  +    end subroutine read_evops
   
       subroutine open_files(type)
         character, intent(in) :: type
  
positive.code: 
      private::init_lindblad_fit
      private::clean_lindblad_fit
      private::write_evolution_operators
      private::open_files
      private::close_files
      private::random_test
  
negative.code: 
  
          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)
  
  