anchor.diff:      !     smoothing_neighbors(nneighbors): single column array with neighbors only
   
       if (.not.fileExists(smoothing_file)) then
  -        call usage('read_inputs: no smoothing file found named '//trim(smoothing_file))
  +        call usage('read_inputs: no smoothing file found named '//trim(smoothing_file)// &
  +                   '(usage:lsqr)')
       endif
       nsmooth = line_count(smoothing_file)
       if (nsmooth.gt.fault%nrows) then
  -        call usage('read_inputs: number of faults to smooth is larger than total number of faults')
  +        call usage('read_inputs: number of faults to smooth is larger than total number of '// &
  +                   'faults (usage:none)')
       endif
   
       ! Read smoothing file and calculate pointers to neighbor array
       allocate(smoothing_pointers(nsmooth,3),stat=ierr)
       if (ierr.ne.0) then
  -        call usage('read_inputs: error allocating memory to smoothing_pointers')
  +        call usage('read_inputs: error allocating memory to smoothing_pointers (usage:none)')
       endif
       open(unit=34,file=smoothing_file,status='old')
       do i = 1,nsmooth
  
positive.code:         endif
          allocate(gf_stress%array(gf_stress%nrows,gf_stress%ncols),stat=ierr)
          if (ierr.ne.0) then
              call usage('read_inputs: error allocating memory to gf_stress%array')
          endif
      endif
  endif
  
negative.code:     allocate(smoothing_neighbors(smoothing_pointers(nsmooth,3)+smoothing_pointers(nsmooth,2)),&
               stat=ierr)
      if (ierr.ne.0) then
          call usage('read_inputs: error allocating memory to smoothing_neighbors (usage:none)')
      endif
      rewind(34)
      do i = 1,nsmooth
  