anchor.diff:      call invert_lsqr()
   elseif (inversion_mode.eq.'anneal') then
       call invert_anneal()
  +elseif (inversion_mode.eq.'anneal-psc') then
  +    call invert_anneal_pseudocoupling()
   else
       call print_usage('!! Error: no inversion mode named '//trim(inversion_mode))
   endif
  
positive.code: 
      ! Read the fault data
      call read_program_data_file(fault)
  endif
  
  ! If pre-stresses are defined, calculate the shear stresses on the faults
  
negative.code: 
  subroutine calc_greens_functions()
  use io_module, only: stderr, verbosity
  use variable_module, only: inversion_mode, displacement, prestress, gf_type, gf_disp, gf_stress
  use okada_module, only: calc_gf_disp_okada_rect, calc_gf_stress_okada_rect
  implicit none
  
  