anchor.diff:            elseif (tag(1:2).eq.'-o') then
                 i = i + 1
                 call getarg(i,ofile)
  +          elseif (tag(1:2).eq.'-d') then
  +              write(*,*) 'Running with default file names'
  +          elseif (tag(1:2).eq.'-u') then
  +              user = 1
             elseif (tag(1:2).eq.'-h'.or.tag(1:2).eq.'-?') then
                 call usage()
             endif
  
positive.code:       IMPLICIT none
  
        write(*,*)
       1 'Usage: sphfinrot -f [IFILE] -o [OFILE] -h/-?'
        write(*,*)
       1 '  -f [IFILE] (Default none) name of input file'
        write(*,*)
       1 '                            if no file name given, fault ',
       2                            'parameters prompted for manual input'
        write(*,*)
       1 '  -o [OFILE] (Default sphfinrot.out) name of output file'
        write (*,*)
       1 '  -h/-?        help'
        write (*,*) ''
        write (*,*)
       1 '  sphfinrot makes a counter-clockwise finite rotation of a',
       2     ' point around a pole'
        write (*,*) ''
        write (*,*)
       1 '    input file'
        write (*,*)
       1 '        lon_init lat_init pole_lon pole_lat angle'
        write (*,*)
       1 '         :  :'
        write (*,*)
       1 '    output file'
        write (*,*)
       1 '        lon_final lat_final'
        write (*,*)
  
negative.code: 
        IMPLICIT none
        REAL*8 pi,d2r
        PARAMETER (pi=4.0d0*datan(1.0d0),d2r=pi/1.8d2)
        REAL*8 lon,lat,x,y,z
  
        x = dcos(lon*d2r)*dcos(lat*d2r)
        y = dsin(lon*d2r)*dcos(lat*d2r)
        z = dsin(lat*d2r)
  
        RETURN
        END
  
  C----------------------------------------------------------------------C
  
        SUBROUTINE gcmdln(ifile,ofile,user)
  
        IMPLICIT NONE
        CHARACTER*30 ifile,ofile,tag
        INTEGER i,narg,user
  
        user = 0
        ifile = 'sphfinrot.in'
        ofile = 'sphfinrot.out'
  
        narg = iargc()
        if (narg.eq.0) call usage()
  
        i = 0
   9998 i = i + 1
        if (i.gt.narg) goto 9999
  