anchor.diff:            trgmatchsta = .true.
         else
             write (*,8999)
  -          call usage()
  +          call usage(' ')
         endif
    8997 format ('NL in targetparams.txt = 1')
    8998 format ('NL in targetparams.txt = NL in stations.txt')
  
positive.code: 
  C----------------------------------------------------------------------C
  
        SUBROUTINE usage()
        IMPLICIT none
        write (*,*)
       1 'Usage: flt2coulomb -V -h/-?'
        write (*,*)
       1 '  -V (default false) turn on verbose operation'
        write (*,*)
       1 '  -h/-?              help'
        write (*,*) ''
        write (*,*)
       1 '  flt2coulomb calculates strain, normal, shear, and ',
       2   'coulomb stresses from faults at user-defined locations.'
        write (*,*)
       1 '  Output files: strain.out, norml.out (pos. = dilation), ',
       2    'shear.out, coul.out'
        write (*,*)
       1 '  Required input files:'
        write (*,*) ''
        write (*,*)
       1 '    faults.txt: list of dislocation sources'
        write (*,*)
       1 '      evlo evla evdp str dip rak slip width length fault_type'
        write (*,*)
       1 '                (km)             (m)  (km)   (km)  (0=pt/1=fin)'
        write (*,*) ''
        write (*,*)
       1 '    stations.txt: list of station locations and depths'
        write (*,*)
       1 '      stlo stla stdp'
        write (*,*)
       1 '                (km)'
        write (*,*) ''
        write (*,*)
       1 '    structure.txt: vp vs dens (only vp/vs ratio matters)'
        write (*,*) ''
        write (*,*)
       1 '    targetparams.txt: str dip rak coeff_frict'
        write (*,*) ''
        stop
        END
  
negative.code:           write (*,8887)
            read *,haffile
            if (haffile.eq.'quit') stop
            if (haffile.eq.'help') call usage(' ')
            goto 13
        else
             if (verbose) write(*,9999)
  