anchor.diff:            write(*,*)
        1 '    onto the rake vector of the target fault.'
             write(*,*)
  -     1 '    To have output be absolute value of shear traction, ',
  -     1       'use -shearu.'
  +     1 '    To have output be maximum value of shear traction, ',
  +     1       'use -shearmax.'
             write(*,*)
             write(*,9999)
             write(*,*)
  
positive.code:               write(*,2006) str0
            endif
            write(*,2007) adist
            write(*,2008) 'LMIN: ',lmin,'LMAX: ',lmax
            write(*,2008) 'ZMIN: ',zmin,'ZMAX: ',zmax
            write(*,2009) incr
            write(*,2010) 'NL: ',int((lmax-lmin)/incr)+1,
       1               'NZ: ',int((zmax-zmin)/incr)+1
            call vgrid(staf,lmin,lmax,zmin,zmax,incr,x0,y0,str0,xy,auto,
       1               adist)
        else
  
negative.code:      1 '               -disp DSPFILE [-az] -strain STNFILE -stress ',
       2                             'STSFILE -estress ESTSFILE'
        write(*,*)
       1 '               -normal NORFILE -shear[max] SHRFILE -coul ',
       2                        'COULFILE'
        write(*,*)
       1 '               [-auto h|d|s DIST INCR] [-long] [-prog] ',
       2                   '[-gmt GMTFILE]'
  