anchor.diff:        REAL*8 uLOS,stlo,stla,stdp,uN,uE,uZ
         REAL*8 az,incl,cosaz,cosinc,sinaz,sininc
         CHARACTER*40 ifile,ofile
  +      CHARACTER*256 azincfile
         LOGICAL ex
         INTEGER uin,uout
   
   C----
   C Parse command line
   C----
  -      call gcmdln(az,incl,ifile,ofile)
  -      if (az.lt.-1000.0d0.or.incl.lt.-1000.0d0) then
  +      call gcmdln(az,incl,ifile,ofile,azincfile)
  +      if ((az.lt.-1000.0d0.or.incl.lt.-1000.0d0)
  +     1                                .and.azincfile.eq.'none') then
             write(*,*) '!! Error: azimuth or inclination are undefined'
  -          call usage('!! Use -a AZ or -i INCL')
  +          call usage('!! Use -a AZ, -i INCL, or -azinc FILE')
         endif
         if (ifile.eq.'none') then
             write(*,*) '!! Error: Input displacement vector file ',
  
positive.code:           write(*,*) '!! Error: No output file specified'
            call usage('!! Use -o OFILE to specify LOS output file')
        endif
        cosaz = dcos(az*d2r)
        sinaz = dsin(az*d2r)
        cosinc = dcos(incl*d2r)
        sininc = dsin(incl*d2r)
        !print *,'AZ',az
        !print *,'INCL',incl
        !print *,'IFILE ',ifile 
        !print *,'OFILE ',ofile
  
  C----
  C Control files
  
negative.code:           write(*,*)
        endif
        write(*,*)
       1 'Usage: vec2los -a AZ -i INC | -azinc file [-f IFILE] ',
       2                 '[-o OFILE] [-h]'
        write(*,*)
        write(*,*)
       1 '-a AZ     Look azimuth (CW from N)'
        write(*,*)
       1 '-i INC    Look inclination (from horizontal=0 to vertical=90)'
        write(*,*)
       1 '-azinc FILE Azimuth and inclination in a file'
        write(*,*)
       1 '-f IFILE  Input vector displacement file (default: stdin)',
       2            ' (stlo stla stdp uE uN uZ)'
  