anchor.diff:  
       ! Read the fault data
       call read_program_data_file(fault)
  +endif
  +
  +! Make coordinates Cartesian
  +if (coord_type.eq.'geographic') then
  +    reflon = fault%array(1,1)
  +    reflat = fault%array(1,2)
  +    do i = 1,fault%nrecords
  +        call ddistaz(dist,az,reflon,reflat,fault%array(i,1),fault%array(i,2))
  +        dist = dist*6.371d6
  +        fault%array(i,1) = dist*dsin(az)
  +        fault%array(i,2) = dist*dcos(az)
  +    enddo
  +    if (displacement%file.ne.'none') then
  +        do i = 1,displacement%nrecords
  +            call ddistaz(dist,az,reflon,reflat,displacement%array(i,1),displacement%array(i,2))
  +            dist = dist*6.371d6
  +            displacement%array(i,1) = dist*dsin(az)
  +            displacement%array(i,2) = dist*dcos(az)
  +        enddo
  +    endif
   endif
   
   ! If pre-stresses are defined, calculate the shear stresses on the faults
  
positive.code: use io_module, only: stderr, verbosity
  use variable_module, only: inversion_mode
  use lsqr_module, only: invert_lsqr
  use anneal_module, only: invert_anneal
  implicit none
  
  if (verbosity.ge.1) then
  
negative.code:                            displacement, prestress, fault, &
                             gf_type, gf_disp, gf_stress, &
                             smoothing, rake_constraint, slip_constraint, &
                             halfspace, coord_type
  use elast_module, only: calc_plane_unit_vectors, calc_traction, calc_traction_components
  implicit none
  ! Local variables
  integer :: i
  double precision :: stress(3,3), nor(3), str(3), upd(3), traction(3), traction_comp(3)
  double precision :: reflon, reflat, dist, az
  
  if (verbosity.ge.1) then
      write(stderr,'(A)') 'read_fltinv_inputs says: starting'
  