anchor.diff:          i = i + 1
           call get_command_argument(i,tag)
           read(tag,*) disp_components
  +    elseif (trim(tag).eq.'-disp:misfit') then
  +        i = i + 1
  +        call get_command_argument(i,disp_misfit_file)
       elseif (trim(tag).eq.'-prests') then
           i = i + 1
           call get_command_argument(i,prestress%file)
  +    elseif (trim(tag).eq.'-prests:weight') then
  +        i = i + 1
  +        call get_command_argument(i,tag)
  +        read(tag,*) stress_weight
       elseif (trim(tag).eq.'-flt') then
           i = i + 1
           call get_command_argument(i,fault%file)
  
positive.code:     elseif (trim(tag).eq.'-lsqr:mode') then
          i = i + 1
          call get_command_argument(i,lsqr_mode)
      elseif (trim(tag).eq.'-lsqr:stress_weight') then
          i = i + 1
          call get_command_argument(i,tag)
          read(tag,*) stress_weight
  
      ! Simulated annealing options
      elseif (trim(tag).eq.'-anneal:init_mode') then
  
negative.code:         call get_command_argument(i,gf_stress%file)
  
      ! Inversion options
      elseif (trim(tag).eq.'-damp') then
          i = i + 1
          call get_command_argument(i,tag)
  