anchor.diff:            write(*,2003)
             write(*,2004) 'HORIZONTAL'
             write(*,2005) adist
  -          write(*,2008) 'WEST: ',W,'EAST: ',E
  -          write(*,2008) 'SOUT: ',S,'NORT: ',N
  +          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,
  -     1                    'NY: ',int((N-S)/incr)+1
  +          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,
  
positive.code:           write(*,*)
       1 '    onto the rake vector of the target fault.'
            write(*,*)
       1 '    To have output be absolute value of shear traction, ',
       1       'use -shearu.'
            write(*,*)
            write(*,9999)
            write(*,*)
  
negative.code:  2005 format('  DEPTH:    ',F10.2)
   2006 format('  GRID_AZ:  ',F10.2)
   2007 format('  HOR_SHFT: ',F10.2)
   2008 format(A5,F10.2)
   2009 format('  GRID_INCR:',F10.2)
   2010 format(A6,I8)
        RETURN
        END
  
  