anchor.diff:  
           ! Slip magnitudes
           do iflt = 1,nflt
  -            read(21,'(A)',end=1002,iostat=ios) input_line
  +            read(lu_log_file,'(A)',end=1002,iostat=ios) input_line
               if (it.gt.0) then
                   read(input_line,*,end=1003,err=1003,iostat=ios) slip(it,iflt,1)
               endif
  
positive.code: 
      elseif (log_format.eq.'anneal_psc') then
          do iflt = 1,nflt
              read(21,'(A)',end=1002,iostat=ios) input_line
              if (it.gt.0) then
                  read(input_line,*,end=1003,err=1003,iostat=ios) locked(it,iflt), slip(it,iflt,1:2)
              endif
  
negative.code: character(len=32) :: log_format
  character(len=512) :: input_line
  logical :: fileExists
  logical, allocatable :: isAccepted(:)
  logical, allocatable :: isFltInSubset(:)
  integer :: it, it0, nit
  