anchor.diff:  !----
   ! Best fitting model
   if (best_output_file.ne.'') then
  +    write(*,*) 'anneal_post: working on best-fitting output file'
       ! Find minimum chi-squared model
       min_chi2 = 1e10
       it = 1
  
positive.code:     deallocate(mean_ds)
  endif
  
  ! 1-D marginal distributions for slip components
  if (nmarg1d.ge.1) then
      do i = 1,nmarg1d
          open(unit=25,file=slip_marg_output_file(i),status='unknown')
          do it = nburnin+1,nit
              write(25,*) slip(it,marg1d_flt(i),1:2)
  
negative.code: character(len=512) :: input_line
  logical :: fileExists
  logical, allocatable :: isFltInSubset(:)
  integer :: i, ios, iflt, nflt, ipole, npoles, it, nit, nburnin, nsubset, itmp
  integer, allocatable :: locked(:,:), plocked(:)
  double precision :: temp0, obj0, min_chi2, fraction_locked
  double precision, allocatable :: obj(:), slip(:,:,:), mean_ss(:), mean_ds(:)
  double precision :: px, py, pz, plon, plat, pmag
  double precision, allocatable :: poles(:,:,:), mean_px(:), mean_py(:), mean_pz(:)
  
  
  ios = 0
  