anchor.diff:        do Utnemele2=1,Nl2
   
   
  -        write(ind(Uelement,Uelement2,Utnemele,Utnemele2,'w'),*, IOSTAT=file_ios) timeStep*(i-1),              &
  +        write(ind(Uelement,Uelement2,Utnemele,Utnemele2,code),*, IOSTAT=file_ios) timeStep*(i-1),              &
                           real(Evops(Uelement,Uelement2,Utnemele,Utnemele2,i)),                                 &
                           aimag(Evops(Uelement,Uelement2,Utnemele,Utnemele2,i))
   
  
positive.code:         call open_files('D')
          call write_devops()
          call close_files('D')
  
      end subroutine do_lindblad_fit_work
  
  
negative.code:       end do
      end subroutine read_evops
  
      subroutine write_evops(code)
        integer (i4b)          :: i, file_ios
        integer(i4b)           :: Uelement, Uelement2,Utnemele,Utnemele2
        character, intent(in)  :: code
  
        do i=1,size(Evops,5)
        do Uelement=1,Nl1
  