anchor.diff:  ! Add equations to model matrix to minimize the curvature (second derivative)
   ! of the solution to Ax = b with weighting factor smooth.
   !----
  -use command_line, only : rake_file
  +use command_line, only : verbosity, rake_file
   use arrays, only : nfaults, nsmooth, smoothing_pointers, smoothing_neighbors
  +use io, only : stderr
   implicit none
   ! I/O variables
   integer :: nrows, ncols, ptr_smooth
   double precision :: A(nrows,ncols), b(nrows,1), smoothing_constant
   ! Local variables
   integer :: i, j, ifault, nneighbor, ineighbor, neighbor
  +
  +if (verbosity.ge.2) then
  +    write(stderr,'(A)') "Starting lsqr_load_smooth()"
  +endif
   
   ! Read smoothing element linking file and load into A
   do i = 1,nsmooth
  
positive.code: enddo
  
  if (verbosity.ge.2) then
      write(stderr,'("Finished damping_linear_ls()")')
  endif
  
  return
  
negative.code: integer :: i, j
  
  if (verbosity.ge.2) then
      write(stderr,'(A)') "Starting lsqr_load_damp()"
  endif
  
  do i = 0,m-1
  