anchor.diff:  integer :: i, j
   
   if (verbosity.ge.2) then
  -    write(stderr,'("Starting damping_linear_ls()")')
  +    write(stderr,'(A)') "Starting lsqr_load_damp()"
   endif
   
   do i = 0,m-1
  
positive.code: ! Local variables
  integer :: i, j, n
  
  if (verbosity.ge.4) then
      write(stderr,'(A)') 'Starting lsqr_load_disp()'
  endif
  
  do i = 1,ndisplacements
      do j = 1,nfaults
          ! For selected components of displacements, load strike-slip source GFs into model matrix
          if (disp_comp.eq.'a') then
              A(i,j)                    = disp_gfs(i,j,1) ! ssx
              A(i+ndisplacements,j)     = disp_gfs(i,j,2) ! ssy
              A(i+2*ndisplacements,j)   = disp_gfs(i,j,3) ! ssz
          elseif (disp_comp.eq.'h') then
              A(i,j)                    = disp_gfs(i,j,1) ! ssx
              A(i+ndisplacements,j)     = disp_gfs(i,j,2) ! ssy
          elseif (disp_comp.eq.'v') then
              A(i,j)                    = disp_gfs(i,j,3) ! ssz
          endif
  
          ! Choose columns to load selected fault slip components
          if (rake_file.eq.'none') then
              n = nfaults
          else
              exit
          endif
  
          ! For selected components of displacements, load dip-slip source GFs into model matrix
          if (disp_comp.eq.'a') then
              A(i,j+n)                  = disp_gfs(i,j,4) ! dsx
              A(i+ndisplacements,j+n)   = disp_gfs(i,j,5) ! dsy
              A(i+2*ndisplacements,j+n) = disp_gfs(i,j,6) ! dsz
          elseif (disp_comp.eq.'h') then
              A(i,j+n)                  = disp_gfs(i,j,4) ! dsx
              A(i+ndisplacements,j+n)   = disp_gfs(i,j,5) ! dsy
          elseif (disp_comp.eq.'v') then
              A(i,j+n)                  = disp_gfs(i,j,6) ! dsz
          endif
      enddo
  
      ! Load selected components of observed displacements into constraint vector
      if (disp_comp.eq.'a') then
          b(i,1)                  = displacements(i,4) ! x
          b(i+ndisplacements,1)   = displacements(i,5) ! y
          b(i+2*ndisplacements,1) = displacements(i,6) ! z
      elseif (disp_comp.eq.'h') then
          b(i,1)                  = displacements(i,4) ! x
          b(i+ndisplacements,1)   = displacements(i,5) ! y
      elseif (disp_comp.eq.'v') then
          b(i,1)                  = displacements(i,6) ! z
      endif
  enddo
  
  if (verbosity.ge.4) then
      write(stderr,'(A)') 'Finished lsqr_load_disp()'
      write(stderr,*)
  endif
  
negative.code: ! Local variables
  integer :: i, j
  
  if (verbosity.ge.2) then
      write(stderr,'(A)') 'Finished lsqr_load_stress()'
  endif
  
  