anchor.diff:                         nmarg1d, &
                          marg1d_flt, &
                          slip_marg_output_file, &
  +                       mean_pole_output_file, &
  +                       nmarg_pole, &
  +                       marg_pole, &
  +                       pole_marg_output_file, &
                          nflocked, &
                          flocked_subset_file, &
                          flocked_output_file
  
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:     close(23)
      deallocate(plocked)
  endif
  
  
  ! Mean slip (AND STANDARD DEVIATION???)
  if (mean_slip_output_file.ne.'') then
  