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:     if (ierr.ne.0) then
          write(stderr,*) 'read_inputs: problem reading slip constraint file'
          write(stderr,*) 'Format:'
          call usage(     '    slip_ss slip_ds')
      endif
  
      ! There can be 1 value for all faults or 1 value for each fault
      if (slip_constraint%nrows.ne.1.and.slip_constraint%nrows.ne.fault%nrows) then
          call usage('read_inputs: number of slip constraints must be 1 or number of faults')
      endif
  
      ! If there is only 1 value, set the full constraint array to that value
  
negative.code:         endif
          allocate(gf_los%array(gf_los%nrows,gf_los%ncols),stat=ierr)
          if (ierr.ne.0) then
              call usage('read_inputs: error allocating memory to gf_los%array (usage:none)')
          endif
      endif
  endif
  