anchor.diff:  integer, allocatable :: locked(:,:), plocked(:)
   double precision, allocatable :: temp(:)
   double precision, allocatable :: obj(:), obj_resamp(:)
  +double precision :: max_obj
   double precision, allocatable :: uncert(:)
   double precision :: min_obj, fraction_locked
   double precision, allocatable :: slip(:,:,:), mean_ss(:), mean_ds(:)
   double precision :: px, py, pz, plon, plat, pmag
   double precision, allocatable :: poles(:,:,:), mean_px(:), mean_py(:), mean_pz(:)
   double precision :: ddum
  +integer :: idum
   logical :: ldum
   
   ios = 0
  
positive.code:     ! Sort the sampled models for more efficient reading
      call heapsort(samp,nresamp)
  
      ! Reset file to read again
      rewind(lu_log_file)
      call read_header(nit,nflt,npoles,log_format)
  
  	! Read search results, saving only sampled models
      call read_resampled_models(nit,nflt,npoles,log_format,locked,obj,slip,poles,temp,uncert,isAccepted)
  
  else
  
  
negative.code:         call usage('')
      endif
      1001 if (ios.ne.0) then
          write(stderr,*) 'anneal_post: more iterations in header (',nit,') than found in file (',it-1,')'
          write(stderr,*) 'Continuing with iterations that have been read'
          ios = 0
          nit = it-1
          exit
      endif
  
      ! Read fault information defined by log file format
  