anchor.diff:          allocate(gf_stress%array(gf_stress%nrecords,gf_stress%nfields))
       endif
   endif
  +if (inversion_mode.eq.'anneal-psc') then
  +    gf_stress%nfields = 2*fault%nrecords
  +    gf_stress%nrecords = 2*fault%nrecords
  +    if (allocated(gf_stress%array)) then
  +        deallocate(gf_stress%array)
  +    endif
  +    allocate(gf_stress%array(gf_stress%nrecords,gf_stress%nfields))
  +endif
   
   !----
   ! Smoothing
  
positive.code:         else
              write(ounit,5001) fault_slip(i,1), fault_slip(i,2)
          endif
      else
          call print_usage('!! Error: frankly, I do not know how you got this far using an '//&
                           'inversion mode that does not seem to exist...')
  
negative.code:     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
  