anchor.diff:  
       subroutine load_prestresses()
       use io_module, only: stderr, verbosity
  -    use variable_module, only : fault, prestress, rake_constraint, gf_stress
  +    use variable_module, only : fault, prestress, rake_constraint, gf_stress, stress_weight
       implicit none
       ! Local variables
       integer :: i, j, nflt
  
positive.code:     integer :: nrhs
      integer :: lwork, info
      character(len=1) :: trans
      double precision :: work(100000)
      double precision :: btmp(nrows,1)
  
      if (verbosity.ge.2) then
          write(stderr,'(A)') 'lsqr_solve_dgels says: starting'
      endif
  
      nrhs = 1
      trans = 'N'      ! A has form (nrows x ncols), i.e. not transposed
  
negative.code:                 enddo
              endif
  
              ! Indicate this column has been adjusted by changing cols2delete
              cols2delete(i) = 1
          endif
  
          ! Is the dip-slip component of this fault constrained?
  