anchor.diff:                                                   'subset of faults'
   write(stderr,*) '-obj OBJ_FILE                   Iteration and objective function'
   write(stderr,*) '-uncertainty UNCERT_FILE        Uncertainty distribution'
  -write(stderr,*) '-resample [SEED]                Resample models proportional to probability (objective)'
  +write(stderr,*) '-resample NRESAMPLE             Resample models proportional to probability (objective)'
  +write(stderr,*) '-seed SEED                      Set random number seed'
   write(stderr,*)
   
   call error_exit(1)
  
positive.code:             i = i - 1
          else
              call get_command_argument(i,arg)
              read(arg,*,iostat=ios) iseed
              if (ios.ne.0) then
                  iseed = 0
                  i = i - 1
              endif
          endif
  
      else
  
negative.code: logical, allocatable :: isFltInSubset(:)
  integer :: it, nit
  integer :: i, ios, iflt, nflt, ipole, npoles, nsubset
  integer, allocatable :: samp_tmp(:)
  integer, allocatable :: locked(:,:), plocked(:)
  double precision, allocatable :: temp(:)
  double precision, allocatable :: obj(:), obj_resamp(:)
  