anchor.diff:  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)
  
positive.code: 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 = ''
  
  ! Number of arguments
  narg = command_argument_count()
  
negative.code: 
  write(stderr,*) 'Usage: anneal_post -f LOG_FILE -nit NITER [-nburn NBURN]'
  write(stderr,*) '       [-best BEST_FILE] [-plocked PROB_FILE] [-slip:mean SLIP_FILE]'
  write(stderr,*) '       [-slip:marg IFLT FILE] [-pole:mean POLE_FILE] [-pole:marg IPOL FILE]'
  write(stderr,*) '       [-flocked[:subset FLT_FILE] FILE]'
  write(stderr,*)
  write(stderr,*) '-f LOG_FILE                     Annealing log file from fltinv'
  write(stderr,*) '-nit NITER                      Number of iterations in search'
  write(stderr,*) '-nburn NBURN                    Number of burn-in iterations in search'
  write(stderr,*) '-best BEST_FILE                 Best fitting model (locked ss ds): also pole_BEST_FILE'
  write(stderr,*) '-plocked PROB_FILE              Probability locked (plocked)'
  write(stderr,*) '-slip:mean SLIP_FILE            Mean slip in each fault'
  write(stderr,*) '-slip:marg IFLT FILE            1-D marginal for IFLT slip (repeat for other faults)'
  write(stderr,*) '-pole:mean POLE_FILE            Mean Euler poles'
  write(stderr,*) '-pole:marg IPOL FILE            Marginals for IPOL rotation pole'
  write(stderr,*) '-flocked FILE                   Fraction of locked faults in each iteration'
  write(stderr,*) '-flocked:subset FLT_FILE FILE   Fraction of locked faults in each iteration from a ',&
                                                   'subset of faults'
  