anchor.diff:      elseif (trim(tag).eq.'-gf:stress') then
           i = i + 1
           call get_command_argument(i,gf_stress%file)
  -    elseif (trim(tag).eq.'-xy') then
  +    elseif (trim(tag).eq.'-xy'.or.trim(tag).eq.'-cartesian') then
           coord_type = 'cartesian'
  -    elseif (trim(tag).eq.'-geo') then
  +    elseif (trim(tag).eq.'-geo'.or.trim(tag).eq.'-geographic') then
           coord_type = 'geographic'
   
       ! Half-space options
  
positive.code: call init_fltinv_data(gf_disp)
  call init_fltinv_data(gf_los)
  call init_fltinv_data(gf_stress)
  coord_type = 'cartesian'
  
  halfspace_file = 'none'
  
  
negative.code:     write(stdout,*) 'read_inputs: starting'
  endif
  
  
  !----
  ! Check that coordinate type is specified
  !----
  if (coord_type.ne.'cartesian'.and.coord_type.ne.'geographic') then
      call usage('read_inputs: no coordinate type specified: use -xy or -geo option')
  endif
  
  !----
  ! Observations are required to constrain the fault slip
  