anchor.diff:      implicit none
       ! Local variables
       integer :: i, j, n, dn, ndsp, nflt, nsmo
  +    integer :: cols2delete(ncols), rows2delete(nrows)
  +    character(len=256) :: fmt
   
       if (verbosity.ge.2) then
           write(stderr,'(A)') "load_slip_constraints says: starting"
       endif
  +
  +    if (verbosity.ge.3) then
  +        write(stderr,'(A)') 'load_arrays says: arrays before modifying topology'
  +        write(stderr,'("A =")')
  +        write(fmt,9001) ncols
  +        do i = 1,nrows
  +            if (i.eq.ptr_disp) write(stderr,'("displacements")')
  +            if (i.eq.ptr_stress) write(stderr,'("pre-stresses")')
  +            if (i.eq.ptr_damp) write(stderr,'("damping")')
  +            if (i.eq.ptr_smooth) write(stderr,'("smoothing")')
  +            write(0,fmt) (A(i,j),j=1,ncols)
  +        enddo
  +        write(stderr,*)
  +        write(stderr,'("b =")')
  +        do i = 1,nrows
  +            if (i.eq.ptr_disp) write(stderr,'("displacements")')
  +            if (i.eq.ptr_stress) write(stderr,'("pre-stresses")')
  +            if (i.eq.ptr_damp) write(stderr,'("damping")')
  +            if (i.eq.ptr_smooth) write(stderr,'("smoothing")')
  +            write(stderr,'(1PE12.4)') b(i,1)
  +        enddo
  +        write(stderr,*)
  +    endif
  +    9001 format('(1P',I3,'E12.4)')
   
       nflt = fault%nrecords
       ndsp = displacement%nrecords
       nsmo = smoothing%nrecords
  -
  -    ! Move contribution from constrained fault segments to observation side of equation (do not
  -    ! renumber yet.....will do that later in subroutine)
  +    rows2delete = 0
  +    cols2delete = 0
  +
  +    ! Move contribution from constrained fault segments to observation side of equation
  +    ! (do not renumber yet.....will do that later in subroutine)
       do i = 1,nflt
           ! Is the strike-slip (or constant rake) component of this fault constrained?
           if (dabs(slip_constraint%array(i,1)).lt.99998.0d0) then
  
positive.code:     ! Apply constraints
      if (slip_constraint%file.ne.'none') then
          call load_slip_constraints()
      endif
  
      if (verbosity.ge.2) then
  
negative.code:             ! Move stresses corresponding to this slip component to b vector
              if (prestress%file.ne.'none') then
                  do j = 1,nflt
                      b(ptr_stress+j-1,1) = b(ptr_stress+j-1,1) &
                                              - A(ptr_stress+j-1,i+nflt)*slip_constraint%array(i,2)
                      b(ptr_stress+nflt+j-1,1) = b(ptr_stress+nflt+j-1,1) &
                                       - A(ptr_stress+nflt+j-1,i+nflt)*slip_constraint%array(i,2)
                  enddo
                  rows2delete(ptr_stress+i+nflt-1) = 1
              endif
  
              ! Move smoothing corresponding to this slip component to b vector
  