anchor.diff:      write(stderr,'("    gf_disp%file:           ",A)') trim(gf_disp%file)
       write(stderr,'("    gf_stress%file:         ",A)') trim(gf_stress%file)
       write(stderr,*)
  -    write(stderr,'("    inversion_mode:         ",A)') trim(inversion_mode)
       write(stderr,'("    damping_constant:       ",1PE14.6)') damping_constant
       write(stderr,'("    smoothing_constant:     ",1PE14.6)') smoothing_constant
       write(stderr,'("    smoothing_file:         ",A)') trim(smoothing%file)
  
positive.code: 
      call get_command_argument(i,tag)
  
      ! Output options
      if (trim(tag).eq.'-o') then
          i = i + 1
          call get_command_argument(i,output_file)
  
  
negative.code: gf_type = 'okada_rect'
  coord_type = 'cartesian'
  disp_components = '123'
  disp_misfit_file = 'none'
  
  ! Initialize regularization variables
  damping_constant = -1.0d0
  