anchor.diff:              ! 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)*slip_constraint%array(i,2)
  +                                     - 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
  
positive.code:     ! Apply constraints
      if (slip_constraint%file.ne.'none') then
          call load_slip_constraints()
      endif
  
      if (verbosity.ge.2) then
  
negative.code: 
      subroutine load_prestresses()
      use io_module, only: stderr, verbosity
      use variable_module, only : fault, prestress, rake_constraint, gf_stress, stress_weight
      implicit none
      ! Local variables
      integer :: i, j, nflt
  