anchor.diff:      endif
   
   enddo
  -write(stdout,*) 'anneal_post: finished reading log file'
  +
  +
  +if (verbosity.ge.1) then
  +    write(stdout,*) 'read_models: finished'
  +endif
  +
  +
   1002 if (ios.ne.0) then
       write(stderr,*) 'anneal_post: reached end of anneal log file before finished reading faults'
       write(stderr,*) 'iteration: ',it
  
positive.code: 
  subroutine read_header(nit,nflt,npoles,log_format)
  
  use io, only: stderr, stdout
  
  use anneal_post, only: lu_log_file, &
                         burnin_iterations
  
negative.code:         allocate(poles(nit,npoles,3))
      endif
  
      ! Initialize the random number seed
      if (iseed.eq.0) then
          iseed = timeseed()
      endif
  
      ! Resample the ensemble
      if (verbosity.ge.2) then
          write(stdout,*) 'anneal_post: resampling models'
      endif
      allocate(samp(nit))
      allocate(obj_resamp(nit))
      max_obj = maxval(obj(1:nit))
  