anchor.diff:        endif
         if (ref(5).gt.8.5d1) then
             call usage('!! Error: fault dip greater than 85')
  +      endif
  + 9001 if (datain.eq.1.and.clip.eq.'none') then
  +          call usage('!! Error: can only check if input is inside '//
  +     1               'polygon along with -clip')
         endif
   C If output file unspecified, print result to standard output
         if (ofile.eq.'none'.and.p.eq.0) p = 1
  
positive.code:       xsec = 0
        xz = 0
        clip = 'none'
        do 100 i = 1,5
            ref(i) = -1.0d0
    100 continue
  
negative.code:       REAL*8 pi,d2r
        PARAMETER (pi=4.0d0*datan(1.0d0),d2r=pi/180.0d0)
        CHARACTER*40 ofile
        CHARACTER*200 clip,line
        REAL*8 x1,x2,y1,y2,dx,dy,dz,dd,ref(5),dist,az,lo,la
        REAL*8 x,y,z,xx(1000),yy(1000),zsave
        INTEGER i,j,nx,ny,p,xsec,xz,nn,io,datain
  
  C----
  C Get parameters from the command line
  