anchor.diff:            write(*,*)
         endif
         write(*,*)
  -     1 '-shear SHEARFILE        Shear stress (',
  +     1 '-shear[max] SHEARFILE   Shear stress (',
        2                               'requires target fault kinematics)'
         if (str.eq.'long') then
             write(*,*)
  
positive.code:               knor = 1
                i = i + 1
                call getarg(i,norf)
            elseif (tag(1:7).eq.'-shearu') then
                if (kstn.eq.0) kstn = 2
                if (ksts.eq.0) ksts = 2
                kshr = 2
  
negative.code:           write(*,2003)
            write(*,2004) 'HORIZONTAL'
            write(*,2005) adist
            write(*,2008) 'W: ',W
            write(*,2008) 'E: ',E
            write(*,2008) 'S: ',S
            write(*,2008) 'N: ',N
            write(*,2009) incr
            write(*,2010) 'NX: ',int((E-W)/incr)+1
            write(*,2010) 'NY: ',int((N-S)/incr)+1
            call hgrid(staf,W,E,S,N,incr,adist)
        elseif (auto.eq.2.or.auto.eq.3) then
            call autovlims(lmin,lmax,zmin,zmax,x0,y0,z0,str0,
  