anchor.diff:      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
  
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:     elseif (trim(tag).eq.'-gf:stress') then
          i = i + 1
          call get_command_argument(i,gf_stress%file)
      elseif (trim(tag).eq.'-xy'.or.trim(tag).eq.'-cartesian') then
          coord_type = 'cartesian'
      elseif (trim(tag).eq.'-geo'.or.trim(tag).eq.'-geographic') then
          coord_type = 'geographic'
  
      ! Half-space options
  