anchor.diff:  
       elseif (log_format.eq.'anneal_psc') then
           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) locked(it,iflt), slip(it,iflt,1:2)
               endif
  
positive.code: 
      elseif (log_format.eq.'anneal_psc_euler') 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
          enddo
          do ipole = 1,npoles
              do i = 1,3
                  read(21,'(A)',end=1004,iostat=ios) input_line
                  if (it.gt.0) then
                      read(input_line,*,end=1005,err=1005,iostat=ios) poles(it,ipole,i)
                  endif
  
negative.code: 
          ! Rake angles
          do iflt = 1,nflt
              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) rake
                  ! Convert to strike- and dip-slip
  