anchor.diff:  burnin_iterations = 0
   best_output_file = ''
   plocked_output_file = ''
  +
   mean_slip_output_file = ''
   nmarg1d = 0
   marg1d_flt = 0
   slip_marg_output_file = ''
  +
   nflocked = 0
   flocked_subset_file = ''
   flocked_output_file = ''
  +
  +mean_pole_output_file = ''
  +nmarg_pole = 0
  +marg_pole = 0
  +pole_marg_output_file = ''
   
   ! Number of arguments
   narg = command_argument_count()
  
positive.code:     endif
      do iflt = 1,nflt
          read(21,'(A)',end=1003,iostat=ios) input_line
          if (it.gt.0) then
              read(input_line,*,end=1004,err=1004,iostat=ios) locked(it,iflt),slip(it,iflt,1:2)
          endif
      enddo
      it = it + 1
  enddo
  ios = 0
  write(*,*) 'anneal_post progress: ',100*it/nit,'%'
  
  1001 if (it.lt.nit) then
      write(stderr,*) 'anneal_post: you requested more iterations (',nit,') than are in log file'
  
negative.code: integer :: marg1d_flt(1000)
  character(len=64) :: slip_marg_output_file(1000)
  
  character(len=512) :: mean_pole_output_file
  integer :: nmarg_pole
  integer :: marg_pole(20)
  character(len=64) :: pole_marg_output_file(20)
  
  integer :: nflocked
  character(len=64) :: flocked_subset_file(10)
  character(len=64) :: flocked_output_file(10)
  